基于拓扑优化技术的集中力扩散结构设计

对于运载飞行器一类由多个部段装配而成的复杂结构, 外载荷往往被转化为多点集中载荷传递到结构某一部段, 为此, 结构上需要采取使集中力扩散的措施. 基于连续体拓扑优化, 以结构最小柔顺性为目标函数, 同时考虑设计域内材料用量约束和考查区域(承载结构)内力均匀性约束, 提出了集中力扩散结构优化设计的理论模型, 并给出了平面和三维两种不同的优化算例, 均得到了合理的优化设计结果. 最后, 针对运载火箭燃料贮箱短壳, 考虑工程实际的受力及约束情况, 给出了可行的概念性设计方案.      

薄壁结构多层级并发加筋拓扑优化研究

新一代航天装备的主承力薄壁舱段在追求极致轻量化的同时,还具有更高的刚度和抗屈曲等设计指标.传统结构形式和设计方法难以满足轻质高承载的设计要求.为此,本文提出了一种薄壁结构多层级并发加筋拓扑优化方法,通过构建主层级稀疏加筋和次层级密集点阵增强结构整体和局部力学性能,扩展结构设计空间,有效提升材料利用率.其中,主层级加筋布局通过变密度拓扑优化方法获得,次层级点阵构型通过基于改进的渐进均匀化方法提出的两种设计方法获得,并基于材料插值模型,建立了多层级并发加筋拓扑优化框架,实现在一次拓扑优化求解中同时获得主层级加筋布局和次层级单胞拓扑构型.基于上述方法,本文分别给出了考虑结构刚度和稳定性设计需求的优化算例,并与传统单层级加筋拓扑优化进行了对比.结果 表明,多层级并发加筋方法可以根据承载边界和设计目标寻找优化的结构形式,且在相同质量下,其优化构型相比传统单层级拓扑优化结果表现出更高的承载性能,证明了本方法在薄壁结构设计上的优势.     

薄壁加筋肋圆柱壳稳定性分析的参数化研究

针对在轴向载荷作用下的正置、正交网格形式的薄壁加筋肋圆柱壳结构，利用有限元程序，对薄壁加筋肋圆柱壳稳定性分析进行了参数化研究，得到了进行结构优化设计的准则，对于给定的设计载荷，当结构参数位于某一个局部失稳与整体失稳的临界区域时，结构的重量最轻。提出了基于有限元分析进行结构优化设计的策略，利用优化策略，获得了一薄壁加筋肋圆柱壳结构的优化设计结果，同时给出了粘合刚度简化模型与有限元计算结果的比较。     

多室弯箱力学参数的非线性动态Markov优化识别研究

对于多室弯箱结构,基于误差理论建立了多室弯箱力学参数的非线性动态Markov误差函数和方差计算公式,并根据正交函数系理论推导了固体壳单元计算列式;利用直接优化理论研究了多室弯箱力学参数的动态识别方法,其中最优步长采用一维自动寻优方案确定。同时给出了多室弯箱力学参数的非线性动态Markov优化识别的具体步骤,并研制了相应的计算程序。研究结果表明:基于正交函数系理论构造的固体壳单元是正确的;动态误差函数能同时考虑多位移考察点和多量测次数的系统响应信息;直接理论用于多室弯箱力学参数的非线性Markov优化识别时表现出良好的数值稳定性和收敛性,具有较高的计算效率。     

基于缺陷敏感性分析的加筋圆柱壳结构设计

为了对加筋圆柱壳进行面向低缺陷敏感度的结构设计,针对一个直径为3m 的正置正交加筋圆柱壳,基于显式后屈曲数值分析,研究了预制不同幅度的内凹、外凸双曲母线形状设计的加筋圆柱壳轴压临界荷载对初始模态缺陷的敏感性。结果表明：外凸双曲母线形状可在幅度为20mm时提高含缺陷结构的轴压临界荷载,可达6.6%,表现为对缺陷的低敏感性;内凹双曲母线则可在幅度为-10mm时小幅提高完善结构的轴压临界荷载,可达1.7%;当加筋圆柱壳应用于运载火箭的燃料贮箱时,燃料加注会升高贮箱内压,因此可通过内压控制形成外凸形状来降低结构的缺陷敏感性。     

复合材料燃料贮箱渐进损伤分析与剩余强度预报

基于复合材料渐进失效准则和刚度降阶模型,建立了复合材料燃料贮箱性能演化的 表征模型和分析模型,借助有限元软件ABAQUS对内压作用下复合材料贮箱结构的渐进损伤过 程进行了分析,实现了复合材料贮箱结构剩余强度的预报,通过与文献数据的比较,验证了 该分析方法的有效性.     

带频率约束的浸水圆柱壳结构优化设计

与干结构相比，流固耦合系统属于非经典问题，其动态设计要复杂得多。本文基于流固耦合分析，主要探讨了浸水圆柱壳在频率，重量等约束条件下的优化设计问题。针对圆柱壳的结构特点，对壳体，流场分别采用有限条，边界元法进行分析以减小计算量。充分利用系统的灵敏数度信息，采用多目标优化方法对浸水圆柱壳进行了优化设计。通过具体算例证明了本文方法的可行性。     

充液矩形贮箱弹性箱底板应力与变形分析

基于相容拉格朗日-欧拉法,通过对流场与弹性固体间流固耦合作用的分析,得到了矩形贮箱弹性底板流固耦合系统的自由振动微分方程。将伯努利方程与外加激励条件、速度势函数耦合到自由振动方程中,采用迦辽金积分法,给出了矩形贮箱在流体作用下的应力与变形的解析解。讨论了弹性底板的抗弯刚度、结构尺寸、底板材料参数及流体深度等因素对底板应力与变形的影响。研究结果表明:在液体晃动非线性激励作用下,贮箱底板的应力和变形随着液体深度、板长的增大而增大,随着板厚的减小而增大,且成非线性变化关系;底板的变形及应力与底板的材料常数相关,其中板厚的变化对其变形和应力影响要比板长及液体深度的影响显著得多。本文结果可为工程实际中矩形贮箱的设计提供参考。     

箱型结构内部爆炸破坏研究进展

结构内部爆炸破坏机理和规律是常规武器毁伤效能预测与评估、建筑物和舰船抗爆防护设计的重要支撑。基于结构内爆炸载荷、内部爆炸作用下结构塑性响应、内部爆炸作用下箱壁结构破坏模式、内部爆炸作用下多箱型结构破坏模式和分布四个方面详细论述了箱型结构内部爆炸破坏的研究现状及存在的问题,并对内部爆炸后续研究给出了建议。建议研究并建立更加复杂的结构内部爆炸载荷和破坏效应描述模型、内部爆炸作用下箱壁的动力响应机理、多箱型结构与内部爆炸波产生的耦合效应、内部爆炸作用下结构的破坏模式和破坏范围的快速准确预测方法等。     

双层环肋圆柱壳受多个物体撞击下的结构动响应

针对双层环肋圆柱壳受到多个物体的撞击问题,采用MSC.Dytran软件对受撞过程中的结构损伤变形、撞击力变化和能量转换进行数值模拟,并与模型试验相对比后发现:双层环肋圆柱壳结构同时受多物体撞击是一个瞬态动响应过程,在巨大瞬时冲击载荷作用下,受撞区壳板会迅速超越弹性变形而产生塑性变形;多物体撞击会造成外壳板一定区域的损伤变形,撞击力会相互干扰,导致其非线性特征更明显。结果表明,双层圆柱壳的外壳能对内壳起到较好的防护作用,在外壳没被撞穿的情况下,其结构变形会吸收绝大部分的撞击动能,可以通过优化外壳的吸能效率来达到双层壳体结构物内壳防撞的目的。     

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

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

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

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

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

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).