夹层板的边界单元法

本文提出了Ｒｅｉｓｓｎｅｒ型和Ｈｏｆｆ型夹层板的直接边界单元法。从夹层弯曲问题的基本方程出发，利用偏微分算子导出夹层板偏微分方程的基本解，收此建立起问题的边界积分方程组，采用边界元法求解，所得结果具有较高精度。     

Hoff型夹层板的边界积分方程法

本文利用Hormander算子法构造了Hoff型夹层板的基本解,建立了Hoff型夹层板的边界积分方程并分析了一些算例。数值结果表明本文方法具有较高的精度。     

夹层板的有限单元法

本文根据Hoff型夹层板弯曲理论,采用有限单元法导出了夹层板单元的刚度矩阵、质量矩阵、一致的载荷矢量和温度载荷矢量;用以计算一般夹层板的变形、应力或自由振动。其结果与理论解相比,具有较满意的精度。     

阻尼夹层筋板结构有限元动力分析

在小变形线弹性理论基础上，导出了粘弹性阻尼夹层筋板壳结构的有限元动力方程，构造了壳单元、阻尼夹层壳单元和开口薄壁梁单元。采用这些单元可以很好地模拟具有不连续阻尼夹层处理的筋板壳结构。数值分析采用［３］中提出的两次渐近法。通过一块阻尼加筋板的特征分析和振动测量，证实了所建立的有限元两次渐近法的良好精度和广泛的工程适用性。     

考虑横法向热应变的C~0型Reddy板理论和三角形板单元

考虑横法向热变形,建议了C0型Reddy理论,并用于分析复合材料层合/夹层板热膨胀问题。虽然考虑了横法向热应变,但不增加额外的位移变量。此理论位移场不含有横向位移一阶导数,构造有限元时仅需C0插值函数。基于这一模型,运用虚位移原理推导了复合材料板平衡方程以及构造了6节点三角形板单元,并分析了简支复合材料层合/夹层板的热膨胀问题。数值结果表明,建立的模型能准确分析复合材料层合/夹层板热膨胀问题,而忽略横法向热应变的理论分析热膨胀问题误差较大。     

用于多种夹层板等效的有限元分析法

针对夹层板力学性能解析法难于计算复杂结构的夹层板且通用性差的问题,本文采用有限元分析法研究了夹层板性能的等效方法。对夹层板的代表体单元模型施加位移约束,模拟弯曲变形时线性独立的应变分量和弯曲内力;根据夹层板内力与应变的本构关系,求出刚度矩阵;最后由刚度矩阵得出宏观等效弹性常数,从而把夹层板等效成连续材料的单层板单元。将该方法与解析法计算结果进行比较得到的夹层板单元四个主要弹性常数误差在0.2%以内,验证了该方法的有效性;另外采用该方法等效三种典型结构夹层板,比较实际模型和等效模型的弯曲响应,得到的误差均在1.4%以内,表明该方法在不考虑复杂多变的夹芯结构时具有通用性。     

利用Ferguson曲面构造广义协调板单元

根据广义协调原理，首先利用Ｆｅｒｇｕｓｏｎ曲面构造出薄板弯曲单元，将中厚度板视为双向深梁，由Ｔｉｍｏｓｈｅｎｋｏ理论拟合单元边界，利用Ｆｅｒｇｕｓｏｎ曲面的张量积性质，将薄板单元推广到中厚度板。数值结果表明此单元精度高，适应性强，且不出现剪切闭锁现象。     

复合材料面层夹层板中转动一致有效理论

采用分层位移模式，引入夹层板层间和表面剪应力协调条件，导出位移场的修正形式。根据能量误差一致原则，应用最小位能原理建立了复合材料层夹层板的基本方程和边界条件。     

扁球壳的边界元几何非线性分析

本文从壳体位移的三个微分方程出发，采用付立叶积分变换的基本解，利用加权残值法推导了几何非线性边界积分方程。这种基本解的壳体边界元法类似于板的非线性边界元法，各种变量物理意义明确，能方便地处理各种复杂边界条件及有开口情况。文末算例说明本文方法的可行性、收敛性和精确性，并与二变量边界单元法或有限元结果相比较，吻合较好。     

矩形夹层板稳定性分析

本文在Reissner 理论基础上，应用功的互等定理推导了夹层矩形板稳定问题的基本解，在已推导出的夹层板基本解的基础上，利用功的互等法求解了两对边固定一边简支一边自由、两邻边简支另两邻边自由且角点支承、两邻边简支另两邻边自由且角点悬空三种不同边界条件下夹层板的稳定问题，给出了挠曲面方程及其对应的执行方程；进行了数值计算，并与有限元结果进行对照分析．结果表明：本文方法求解过程更简单，提供了一种求解夹层板稳定问题的新方法，计算结果对解决工程实际问题具有一定的参考价值．     

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

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

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