热超弹性圆筒的不稳定性

应用有限变形弹性理论分析了受内压和轴向拉伸作用的不可压热超弹性圆筒发生非均 匀变形的不稳定性问题. 受内压和轴向拉力作用的薄壁圆筒，当内压较小时，圆筒发生稳定 的均匀膨胀变形；当内压大于某一临界值时，圆筒产生复杂的非均匀变形，其一部分膨胀变 形很大，形如``灯泡'状，而另一部分仅仅是轻微膨胀，且此时的变形是不稳定的. 但对厚 壁圆筒而言，不论压力如何，总是发生稳定的均匀膨胀变形. 根据圆筒的变形曲线，给出了 圆筒可以发生不稳定变形的临界厚度. 同时，讨论了轴向拉伸和温度场对圆筒变形的影响.     

复合材料层合板的二次屈曲和二次分枝点分析

为了研究复合材料层合板的二次分叉特性 ,利用能量变分原理和非线性几何方程建立了具有弹性约束的复合材料层合板在面内载荷作用下的非线性稳定性控制方程组。控制方程组用广义傅立叶级数法进行求解 ,并得到载荷 -挠度曲线。基于分叉理论中的 Lerray-Schaulder定理 ,采用小挠动法 ,直接导出了复合材料层合板的二次失稳方程。研究结果表明 ,非对称层板也可能存在分叉 ,弹性转动支持系数和铺层等因素对二次分叉有很重要的影响。随着弹性系数的增大 ,二次失稳载荷值与初次失稳载荷值之比下降     

注浆加固和可压缩层联合作用下隧道力学行为研究

根据高地应力深埋隧道的变形特征,以围岩加固和变形能量释放为原则,提出了注浆加固与可压缩层联合支护的方法,建立了基于注浆加固与可压缩层联合作用下的隧道力学模型;推导了在联合作用下隧道应力场和位移场的解析解,并与数值计算结果进行了对比验证。结果表明:围岩的注浆加固和可压缩层的设置均能有效改善衬砌的受力情况。在其他条件保持不变时,衬砌的压力和位移,以及围岩的位移都随着注浆加固层厚度的增大而减小,但压力和位移的减小趋势随着其厚度的增加而趋于平缓。同样地,可压缩层厚度增大可导致围岩变形能量更充分地耗散,衬砌的压力及位移均随着可压缩层厚度的增大而不断减小,该减小趋势随可压缩层厚度的变化趋于平缓。但是,随着可压缩层厚度的增加,围岩的位移却不断增加。另外,隧道变形的改善效果与可压缩层填充材料的力学参数密切相关,当其弹性模量由26MPa增至160MPa时,衬砌位移迅速增加。因此,该填充材料应具备足够的变形能力才能达到吸收围岩变形、改善衬砌受力的目的。     

?????????????????е???????

研究了埋置于弹性地基内充液压力管道中非线性波的传播. 假设管壁是线弹 性的,地基反力采用Winkler线性地基模型,管中流体为不可压缩理想流体. 假定系统初始 处于内压为$P_0$的静力平衡状态,动态的位移场及内压和流速的变化是叠加在静 力平衡状态上的扰动. 基于质量守恒和动量定理,建立了管壁和流体耦合作用的非 线性运动方程组; 进而用约化摄动法, 在长波近似情况下得到了KdV方程,表征 着系统有孤立波解.     

闭孔泡沫铝泊松比及三轴压缩变形模式

在闭孔泡沫铝的唯象本构模型中, 泊松比是一个非常关键的参数, 为了探究闭孔泡沫铝泊松比变化规律研究结果存在分歧的原因, 认识闭孔泡沫铝泊松比变化规律中特征点的物理意义, 采用数值模拟方法, 建立了闭孔泡沫铝的3D-Voronoi模型及2D-Voronoi模型, 对模型进行侧面位移耦合单轴压缩边界条件下的仿真分析; 基于闭孔泡沫铝本构模型的唯象特性, 对闭孔泡沫铝变形模式的研究同样十分重要, 为明确其三轴压缩下的变形模式, 对闭孔泡沫铝的3D-Voronoi模型进行侧面位移受限轴向压缩边界条件下的仿真分析. 研究结果表明, 常规壳单元接触中的厚度减薄特性是闭孔泡沫铝泊松比变化规律的研究结论存在分歧的原因, 但厚度减薄不影响泡沫铝模型致密前胞孔结构的变形模式; 闭孔泡沫铝泊松比的准确变化规律为“增高?降低?再增高”的“S”型曲线, 并且, 曲线极大值对应闭孔泡沫铝吸能效率的增速下降点; 等比压缩应力状态下, 闭孔泡沫铝存在四种侧面变形模式, 分别为“(短期)压缩变形→膨胀变形”、“压缩变形→膨胀变形→压缩变形→膨胀变形”、“压缩变形→(短期)膨胀变形”及“压缩变形”.      

平面应变不可压缩橡胶圆柱的大变形

针对一类新的橡胶材料应变能函数，推导了受内压作用下不可压缩橡胶圆柱的大变形公式，给出了位移、应力的解析表达式（用积分形式表示）。建立了适用于分析非线性不可压缩橡胶材料的罚有限元列式，算例表明：位移与应力能很好地与理论解吻合，并提出了控制计算稳定的方法，特别是详细地讨论罚因子的选取及对计算结果的影响。     

生长对超弹性球壳变形和稳定性的影响

应用连续介质力学有限变形理论建立受内压作用不可压超弹性球壳大变形问题的力学模型, 且运用基于变形梯度张量极分解的弹性体积生长理论分析生长对不可压超弹性球壳变形和稳定性的影响. 通过对球壳变形与内压关系式的数值计算得到不同生长条件下球壳的变形曲线和应力分布曲线及由生长引起的残余变形和残余应力分布. 计算结果表明生长对球壳变形特性有明显的影响, 生长影响球壳可产生不稳定变形的临界壁厚和临界内压, 且在某些情况下可改变球壳的稳定性. 生长在球壳中可产生一定的残余变形和残余应力, 对球壳中的应力分布有一定的影响. 另外当生长的程度足够大时, 即便没有外力作用,球壳仅在生长引起的残余应力作用下也可产生不稳定变形.      

薄壁细胞受微吸管与探针共同作用接触模型及数值模拟

建立了以典型的薄壁球型植物细胞为原型的细胞、微吸管及探针接触模型.模型的细胞壁采用封闭球形薄膜,其本构关系为体积不可压超弹性,膜球内充满有压流体以模拟细胞质.应用轴对称几何非线性方法得出了基本微分方程组,并应用龙格-库塔法进行了求解;同时,应用流固耦合有限元法进行了数值模拟以资比较.两种方法得出了较为一致的变形和应力分...     

含孔边裂纹球壳的断裂分析

在依据Reissner理论得出的球壳裂纹尖端应力应变场展开式基础上,采用局部—整体分析法和权函数方法分别计算承受内压的含孔边裂纹球壳的应力强度因子.在有限元的模式中考虑剪切变形的影响,并对奇异元模式的应力应变场展开式的项数选择、奇异元最佳尺寸的选取进行了分析.本文计算和分析了在不同几何尺寸,不同开孔大小以及不同剪切参量条件下承受内压的含孔边裂纹球壳的应力强度因子及其变化规律.     

反平面弹性中边缘内分叉裂纹的奇异积分方程解法

利用复变函数和奇异积分方程方法,求解反平面弹性中半平面边缘内分叉裂纹问题。提出了满足半平面边界自由的由分布位错密度表示的半平面中单裂纹的基本解,此基本解由主要部分和辅助部分组成。将半平面边缘内分叉裂纹问题看作是许多单裂纹问题的叠加,建立了以分布位错密度为未知函数的Cauchy型奇异积分方程组。然后,利用半开型积分法则求解奇异积分方程,得到了裂纹端处的应力强度因子。文中给出两个数值算例的计算结果。     

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

人造地震动合成中的位移误差分析

三角函数级数法是合成人工地震动常用算法之一,但是通过对加速度积分求取位移时,却存在与零线漂移相类似的位移漂移现象。