热结构稳态响应的耦合灵敏度分析方法

研究结构稳态热变形和热应力的灵敏度分析方法，给出了直接法和伴随法两种算法。考虑了温度场的耦合作用，在直接法中需要计算温度场对设计变量的导数，在伴随法中需要计算热载荷对温度场的导数。对尺寸和形状两类设计变量的灵敏度分析算例，验证了本文方法的精度。伴随法在应用程序中的实现，为大型结构优化提供了高效率的灵敏度计算方法。     

航空发动机叶片-机匣碰摩热-结构耦合仿真

针对转静子碰摩过程中的摩擦热效应现象,建立了叶片-机匣的碰摩热-结构耦合模型。采用热-结构耦合单元,对航空发动机叶片-机匣进行了碰摩热效应特性研究。考虑摩擦因数及旋转过程中的离心力作用,对瞬态-热结构耦合场进行了有限元分析,研究了结构在温度场作用下的热-结构耦合应力和温度分布。与纯机械载荷下的应力分布对比,发现了碰摩热效应在叶片与机匣上的扩散规律。研究结果表明,考虑摩擦热效应时由于热应力的作用,导致结构的总应力水平升高进而产生热变形,从而使故障进一步恶化。由此可见,碰摩热效应的影响在实际航空发动机振动分析中不能忽视。     

热结构瞬态响应的耦合灵敏度分析方法与优化设计

研究结构瞬态热变形和热应力的灵敏度分析方法及其优化设计,灵敏度计算给出了直接法和伴随法两种算法.考虑了温度场的耦合作用,在直接法中需要计算温度场对设计变量的导数,在伴随法中需要计算热载荷对温度场的导数.数值算例验证了该方法的精度.伴随法在应用程序中的实现,为大型结构优化提供了高效率的灵敏度计算方法.     

采用包络-准则法优化考虑瞬态传热的薄板结构

刚度和强度是薄板结构的两个主要性能。在瞬态传热过程中,考虑热-力耦合,随时间和空间变化的非均匀温度场在结构中会引起热变形和热应力,温度场随时间变化的规律和空间分布依赖于板的厚度变化,进而影响板的刚度和强度。因此,考虑瞬态传热的薄板优化问题具有更强的非线性,更加难以求解。本文给出一种包络-准则方法处理这类结构优化问题。首先,针对外力荷载,进行一个结构柔顺性的优化设计;以这一设计为基础,通过瞬态热-力耦合分析及优化准则,计算多个时刻的优化设计变量并取其包络,对上述优化结果进行迭代修正,以消除瞬态温度场作用下较高的局部应力。优化算例表明,该方法对于考虑瞬态传热薄板优化问题有效。     

中心舱体与薄壁梁刚-柔耦合系统的热弹性-结构动力学分析

空间柔性结构受太阳热流冲击而诱发的振动是导致航天器失效的典型模式之一,准确预测结构热致振动的响应及稳定性是卫星设计的基础。针对常见的中心舱体与附属薄壁杆件组成的空间结构,提出了考虑刚-柔耦合、耦合热弹性和耦合热-结构三重耦合效应的热致振动分析理论模型。其中,刚-柔耦合是指舱体姿态角、顶端集中质量转动与柔性附件运动的耦合;耦合热弹性是指应变率与温度场的耦合;耦合热-结构是指舱体转动及结构变形与薄壁杆件吸收太阳热流的耦合。基于热弹性理论和Lagrange方程,推导了传热和运动的耦合方程;采用Laplace变换方法并使用Routh-Hurwitz稳定判据推导了稳定性边界方程。结果表明,该模型能够更为准确的给出热致振动响应及稳定性预测。     

考虑热流固耦合作用的多孔介质孔隙尺度两相流动模拟

为了研究深层油气资源在岩石多孔介质内的运移过程, 使用一种基于Darcy-Brinkman-Biot的流固耦合数值方法, 结合传热模型, 完成了Duhamel-Neumann热弹性应力的计算, 实现了在孔隙模拟多孔介质内的考虑热流固耦合作用的两相流动过程. 模型通过求解Navier-Stokes方程完成对孔隙空间内多相流体的计算, 通过求解Darcy方程完成流体在岩石固体颗粒内的计算, 二者通过以动能方式耦合的形式, 计算出岩石固体颗粒质点的位移, 从而实现了流固耦合计算. 在此基础上, 加入传热模型考虑温度场对两相渗流过程的影响. 温度场通过以产生热弹性应力的形式作用于岩石固体颗粒, 总体上实现热流固耦合过程. 基于数值模型, 模拟油水两相流体在二维多孔介质模型内受热流固耦合作用的流动过程. 研究结果表明: 热应力与流固耦合作用产生的应力方向相反, 使得总应力比单独考虑流固耦合作用下的应力小; 温度的增加使得模型孔隙度增加, 但当注入温差达到150 K后, 孔隙度不再有明显增加; 温度的增加使得水相的相对渗流能力增加, 等渗点左移.      

热沥青混合料摊铺后基层温度应力瞬态分析

根据传热学原理,利用有限元方法对沥青混合料摊铺后的路面结构瞬态温度场进行数值模拟.在瞬态温度场分析过程中,不仅考虑了热沥青混合料面层对下承层结构的热传导作用,而且也考虑了面层表面与外界环境的热交换,包括:太阳辐射、空气对流换热和路表辐射换热等.将各不同时刻的路面结构温度场作为荷载施加到结构上,计算出半刚性基层结构的瞬态温度应力场.计算结果表明:热沥青混合料摊铺引起的基层最大温度应力明显小于基层的极限强度,不会造成基层的破坏.     

金属互连结构的电迁移失效分析新算法

综合考虑了电子风、温度梯度、应力梯度和原子密度梯度四种电迁移驱动机制, 基于ANSYS软件平台和FORTRAN程序提出一种新的电迁移失效分析算法．通过ANSYS电－热－结构耦合分析获得模型的电流密度分布、温度分布和应力分布, 基于FORTRAN编写的原子密度重分布算法获得不同时刻的原子密度, 依据空洞生成和扩展失效准则进行电迁移动态空洞演化模拟并得到失效寿命．最后, SWEAT结构与CSP结构的应用算例验证了算法的精度．     

含损伤演化的TM耦合数值模型及其应用研究

从岩石材料的细观结构层次出发，应用损伤力学和热弹性理论，对热力耦合作用下岩石破裂 过程中热－应力相互作用关系进行了分析. 初步建立了细观岩石热－力(TM)耦合数值模型， 探讨了TM耦合作用下岩石材料的细观结构损伤及其诱发的材料力学性能演化机制，把岩石 热固耦合问题的研究从应力状态分析深入到损伤、破坏过程分析之中. 运用该数值模型对某 硬岩实验室开展的原位尺度实验中的废料处理井间柱稳定性进行了模拟分析，其应力场、岩 石剥离破坏形态及破坏诱发的AE特性等均与实验监测结果表现出了较好的一致性，证 明了该数值模型的合理性和有效性.     

热力耦合结构的弹性支撑分析与拓扑优化设计

弹性支撑是保证热力耦合载荷作用下结构有效承载的一种支撑设计方案, 可以在满足结构刚度设计要求的同时有效降低结构内的应力集中, 保证结构的热变形协调. 从理论分析与数值计算两个方面, 研究了热力耦合载荷作用下结构的弹性支撑优化设计. 首先以弹性支撑梁模型为例, 通过给出热力耦合载荷下应力计算的解析表达式, 从理论上说明弹性支撑对结构应力峰值的影响, 阐述了弹性支撑的承载与热变形协调的双重效果. 在此基础上, 提出了解决该类问题的弹性支撑通用优化数学模型和拓扑构型优化方法, 通过数值算例与优化结果展示了方法的有效性.     

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

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

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