张拉结构非线性分析两节点曲线单元有限元法

提出适合张拉结构几何非线性分析的两节点曲线单元有限元方法．假定索元的初始形状呈二次抛物线,根据单索的平衡条件、几何和物理关系建立了索元的位移函数;由拉格朗日应变的定义建立了可以考虑任意次高阶位移影响的索元轴向应变的精确表达式,并基于拉格明日描述方法和虚功原理得到了索元的非线性平衡方程与切线刚度矩阵．采用荷载增量法与Newton－Raphson法相结合的混合法进行了实例计算,结果表明：本文方法的精度明显优于两节点直线索单元,适合于大跨度索阿、索穹顶等张拉结构的几何非线性分析．     

一种连续索滑移的处理方法

目的是建立一种快速简便的拉索滑移问题的有限元计算方法。采用基于弹性悬链线精确解的两节点悬链线单元模拟索结构中的拉索，利用求解非线性方程组的延拓法，再配以常用的牛顿法，建立了快速的给定索原长情况下的拉索刚度计算方法。此法对于一般的给定初值，能很快逼近真实解，从而扩大了解的收敛范围，更重要的是与牛顿法相比更加高效。进而采用二分法建立了索结构中拉索滑移问题处理的算法，编制了相应的程序段，并通过算例证明了其正确性。本文的拉索滑移处理方法可用于索结构的设计与施工。     

基于欧拉描述的两节点索单元非线性有限元法

本文针对柔性悬索结构几何非线性分析的特点,提出了一种用欧拉描述来表示的两节点索单元非线性有限元模型,在索元变形后的位置上由虚功能建立了非线有限元基本方程及切线刚度矩阵。这样建立的非线性有限元分析方法可充分考虑拉索的几何非线性特性的影响并给悬索结构的初始平衡分析带来方便,算例结果表明,本文方法是精确有效的。     

非线性梁索单元理论在轻型桅杆设计中的应用

导出了轻型桅杆设计计算的非线性梁与索单元理论和方法，求出的桅杆截面内力考察了索单元的非线性影响，计入了桅杆杆身在轴压下弯曲失稳的效应。在计算桅杆时采用了增量Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ迭代法，如果在设计荷载作用下的迭代过程收敛，则桅杆不会发生绕两个主轴的整体弯曲失稳定。本文方法解决了轻型桅杆设计中的强度、整体和层间稳定及刚度计算的难题。     

非线性梁索单元理论在轻型桅杆设计中的应用

导出了轻型桅杆设计计算的非线性梁与索单元理论和方法，求出的桅杆截面内力考虑了索单元的非线性影响，计入了桅杆杆身在轴压下弯曲失稳的效应。在计算桅杆时采用了增量Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ迭代法，如果在设计荷载作用下的迭代过程收敛，则桅杆不会发生绕两个主轴的整体弯曲失稳。本文方法解决了轻型桅杆设计中的强度、整体和层间稳定及刚度计算的难题。     

桅杆风致破坏机理研究

桅杆风致破坏的情况屡见不鲜。桅杆的非线性是相当复杂的,它是由拉索的几何非线性和拉索的张力产生的,拉索的张力又依赖于拉索的几何位置。研究指出由于在经常性风荷载作用下,桅杆产生空间混沌振动,并由此产生疲劳破坏的结果,但这些研究仅限于确定性情况。本文中将研究扩展到随机问题,数值计算结果表明该系统存在随机混沌现象,其响应平均值与平均应力指出混沌性质。为了进一步研究桅杆破坏机理,因此在本文中,既考虑突发性故障问题,又考虑随机疲劳寿命问题,并给出时间与应力极值概率。     

基于初弯曲单元的某弦支穹顶非线性稳定承载力分析

在单元的随动坐标系下,建立了初始弯曲对杆单元轴向刚度影响的计算公式,导出了考虑初始弯曲的非线性杆单元刚度矩阵;从经典的梁-柱理论出发,给出了考虑初始弯曲的非线性空间梁单元切线刚度矩阵推导过程。将建立的初始弯曲单元应用于弦支穹顶结构算例的非线性稳定承载力分析,研究了杆件初始弯曲对结构整体非线性稳定性能的影响,结果表明,杆件初始弯曲的存在会降低结构的整体刚度和极限承载力,而在不同的初始预应力状态下,极限荷载随初始弯曲增加都基本呈线性下降,并且其关系曲线基本平行。本文建立的基于初弯曲单元的非线性分析方法,可推广应用于各种类型空间网格结构的非线性分析,为结构设计方案的选择提供更为科学的依据。     

张拉结构非线性分析的五节点等参单元

本文针对张拉结构的特点，提出了一种五节点等参数单元有限元模型，采用四次多项式作为位移插值函数及单元初始形状函数，并假定索是理想柔性的且满足虎克定律，基于修正的Ｌａｇｒａｎｇｉａｎ坐标描述法，建立了非线性有限元基本方程和切线刚度矩阵，利用Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ法进行了实例计算。结果表明：本文方法精度极高，可供大跨度索网，索穹顶，拉线塔等张拉结构分析，设计时采用。     

一种考虑初始垂度影响的非线性索单元

本文首先运用Mathematic的符号运算功能，通过解偏微分方程给出了不等高单索的显式解析解，然后把该解析解应用于索微元的应变计算中，在此基础上用虚功原理推导了考虑初始垂度影响的两节点非线性索单元，并给出了索单元的单刚矩阵的具体形式，通过与两节点直线索单元的单铡矩阵的形式的比较，明确了该曲线非线性索单元的修正项，并给出了垂度影响因子的变化曲线。比较直观的给出了垂度对索单元刚度各项的影响程度。该非线性索单元既有多节点索单元精度高的特点，又有节点少，刚度元素较易求解以及有限元列式简洁等特点。本文通过两个数值算例表明，本文的非线性索单元是正确的，也表明了所编制的非线性计算程序是正确和可靠的。     

基于正交基无单元Galerkin法和非线性规划的安定分析方法

基于安定分析的下限定理,用正交基无单元Galerkin法建立了交交载荷作用下理想弹塑性结构安定分析的下限计算格式.在给定载荷域的载荷角点所对应的载荷作用下,采用正交基无单元Galerkin法计算相应的虚拟弹性应力场.并且利用结构在正交基无单元Galerkin法弹塑性增量分析中平衡迭代结果计算得到自平衡应力场基矢量,然后由这些基矢量的线性组合模拟自平街应力场.安定分析问题最终被归结为一系列未知变量较少的非线性数学规划子问题,通过复合形法求解.算例表明该方法的计算结果是令人满意的,并且对初始复合形顶点和用于构造自平衡应力场基矢量的载荷增量是非常不敏感的.     

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

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

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


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).     

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

正~~     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).