辐射换热下瞬态热-结构分析的一种空间薄壁杆单元

发展了一种用于辐射换热条件下瞬态热-结构分析的空间薄壁杆单元,其截面形式可以是任意形状的闭口截面和单支开口截面。该单元温度场分解为平均温度和多谐摄动温度,沿杆轴方向采用两结点线性插值,沿杆截面周向用三角函数展开,每结点含多个解耦的自由度,其中结点平均温度方程同传统一维温度有限元方程为非线性,各谐摄动温度方程为线性,然后利用Wilson-θ法求解结构的瞬态温度场。本文选择了两节点Ber-noulli直梁单元得到准静态热弹性有限元方程并求解,针对非对称开口截面考虑了翘曲变形及弯扭耦合的影响。温度场引起的等效热载荷不仅包括常规的热轴力,还包括热弯矩以及热双力矩。本文针对不同截面形式的梁单元给出了瞬态温度场以及热变形的验证算例,并通过与商业程序中二维壳元计算结果的比较说明了本文所提出方法的正确性和高效性。     

大型空间结构的热-动力学耦合问题及其有限元分析

论文对辐射换热条件下闭口薄壁杆件与单枝开口薄壁杆件的瞬态温度场问题,提出了一种一维傅立叶温度有限元,克服了传统一维温度单元只能计算薄壁杆截面平均温度的缺点,通过增加结点摄动温度自由度的方法,该一维单元能计算杆截面的温度分布.在此一维温度单元与梁位移单元相协调的基础上,进一步发展了大型空间结构热诱发振动稳定性判据与热颤振响应有限元计算方法.对于柔性空间结构发展了考虑几何非线性的热-结构动力学耦合有限元计算方法,成功地对这类结构的热动力屈曲问题进行了数值模拟.     

空间薄壁管结构瞬态温度场、热变形有限元分析

对于辐射换热的薄壁杆件航天结构，为分析其在太空中不同时刻的姿态下的温度场和热变形，构造了一种相对自由度矩形管单元。其基本思想是从2维非线性瞬态热传导方程出发，假设沿矩形管横截面上每边的温度为线性分布，用4个角点的平均温度和温差来表示矩形管横截面上的温度分布，构造了一种1维2结点温度杆单元，该单元每个结点包含平均温度，上下面温差，左右面温差3个自由度；在计算热变形时，此三个广义温度参数分别对应热轴力和2个热弯矩载荷。经与三维有限元计算结果的比较，证明用该单元计算的矩形管温度场和位移场是可靠的。利用这种新型的薄壁矩形管单元和本文作者在其他文章中提出的薄壁圆管单元，可以对非线性换热条件下的复杂空间结构进行比较准确的温度场和热变形有限元分析，最后本文计算了考虑遮挡的太阳能帆板的瞬态温度场和热变形以说明其应用价值。     

大型空间结构的热-结构动力学分析

空间结构在辐射换热条件下的热诱发振动是导致空间结构失效的一种典型模式。弄清热诱发振动的机理是理解热诱发振动失效的基础。本文针对常见的空间薄壁杆件结构，提出了一种能够对复杂结构及加热条件进行比较准确的温度场和热诱发振动分析的有限元方法。首先利用一种Founer-有限元方法，同时考虑杆截面内平均温度和温差，求解了包含辐射非线性的瞬态热传导问题，并推荐了一种有效降低求解规模的减缩近似方法-Lanczos方法。在此基础上，用有限元法求解了杆件结构的热诱发振动问题，并就杆截面内平均温度和温差对结构振动的影响以及最大动静态响应的比值分别进行了讨论．合理地解释了一类常见的热诱发振动现象，本文的数值算例说明了这点。     

TEC结构的三维非线性瞬态温度场分析

热电制冷器(TEC)以其体积小、作用速度快及无噪音等机械制冷无法替代的优点在航空航天和电子工业等领域得到了越来越广泛的应用。本文根据TEC的导热特点,推导了TEC结构稳态温度场的解析解,建立了其瞬态非线性温度场分析的微分方程。利用伽辽金法导出TEC结构热分析的有限元方程,对非线性热分析的有限元方程进行了求解,得到了TEC的稳态温度场和瞬态响应温度场。算例结果表明,本文提出的TEC结构热分析有限元模型具有较高的精度,能够有效地分析TEC的非线性瞬态温度场。     

热流密度/对流换热边界下双向梯度板的瞬态温度场分析

采用有限单元-有限差分法研究了热流密度/对流换热边界条件下双向梯度板的瞬态热传导问题。采用细观力学方法结合混合律准则描述了材料的热物理属性,通过推导一种8节点高阶双向梯度单元建立了结构的连续梯度有限元模型。计算给出了在考虑组份属性的温度效应下,温度场的时间响应历程以及不同时刻温度场的空间分布形式,并与材料属性温度无关时的计算结果进行了比较,最后讨论了相关参数对瞬态温度场的影响规律。结果表明：温度较低时,组份属性的温度效应对瞬态温度场影响很小;在 y 方向热流密度载荷的作用下,温度场沿 x、y 方向均存在明显的梯度;x 方向组份体积分布系数的增大,延长了温度场达到稳态需要的时间,绝对温度梯度沿 x、y 方向均增大,稳态温度场升高;增大 y 方向组份体积分布系数的值,情况相反。     

四边形等参单元温度场有限元分析中的混合边界条件

在温度场有限元分析中，边界条件的合理确定是一个非常重要的问题。本文以四边形等参单元为基本单元，采用基于微分方程等效积分原理的Galerkin加权余量法，建立了热传导问题的有限元方程，推导出混合边界条件的有限元计算公式。最后，根据求得的温度场可得截面温度应力分布。     

一期水管冷却效应的数值模拟新方法

采用有限单元离散空间域，视冷却水管所在单元为子结构，其它单元为常规单元，然后，把混凝土结构瞬态温度场的有限元控制方程表示成状态方程，最后，利用精细积分法求解．结果表明，子结构、状态方程和精细积分技术的引用降低了对计算机贮存量的要求和提高了冷却水管附近温度梯度的精度，显示出极大的优越性．     

Timoshenko-Euler楔形梁有限元

本文首先建立楔形梁包含轴力和剪切变形效应的平衡微分方程，由于该方法是二阶变系数微分方程，其争析解很难得到，本文通过将该方程中的变系数和方程的解用Chebyshev多项式逼近得到了Timoshenko-Euler楔形梁的单元刚度方程，最后通过算例检验了所得单元刚度方程的对称性，以及验证了计算悬臂梁挠度和悬臂柱弹性临界力的正确性及其收敛性，本文提出的方法可适用于任意变截面Timoshenko-Euler梁单元刚度方程的求解，运用此方法，除可以考虑轴力和剪切变形的影响外，还可以减少结构分析中的单元数和自由度，提高包含楔形构件的结构分析的精度和速度。     

复合多层混合结构三维热传导分析

复合多层混合结构由不同热导率材料混杂组成。对于分析复合多层混合结构内部热传导机理以及非均匀温度场的复合混合结构热力耦合计算,常需要采用三维模型以得到较为精确的结构温度场分布结果。针对包含复合材料的多层混合结构的热传导问题,本文通过基于散度表示的一般瞬态三维热传导方程,利用Galikin法建立了复合多层混合结构瞬态三维热传导的有限元方程。对于复合多层混合结构单层材料过薄,引起的沿厚度方向划分单元的困难,采用等效三维单元使一个单元内包含一个或数个材料铺层,并在沿厚度方向热导率变化剧烈的情况对单元进行细化。数值仿真结果表明该方法网格划分灵活,计算结果较为理想,同时由于其拓扑结构包含的信息较为充分,能够实现多物理场条件下的耦合分析。     

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

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

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

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

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

正~~     

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

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