直杆碰撞刚性壁弹塑性动力后屈曲有限元分析

利用显式动力学有限元方法对直杆弹塑性动力后屈曲进行了分析,模拟了直杆轴向碰撞动力屈曲的变形及发展过程。分析中在直杆碰撞端局部临界屈曲长度范围内引入半正弦波形式的初始缺陷,计算结果与文献中的实验数据获得了很好的一致。分析结果表明,随着碰撞过程中所产生的应力波逐渐向前传播,后屈曲变形过程中所呈现的多个半波形式的高阶屈曲模态由初始具有单个半波形式的简单屈曲模态迅速演变而成。分析结果同时也揭示了直杆动力屈曲变形发展的机理,以及轴向应力波和屈曲变形的相互作用规律。     

电磁内爆套筒动力屈曲的实验和数值模拟

用实验及数值模拟方法研究了电磁内爆套筒的屈曲响应规律。用电容器组脉冲发生器装置作为驱动源,设计4种不同材料、不同尺寸的金属套筒,通过调节电容器组的充电电压得到不同的加载电流,研究材料及几何参数、加载脉冲特性对套筒屈曲的影响;采用瞬态非线性有限元方法对实验结果进行了数值模拟。研究表明,在同样加载电流作用下相同材料套筒形成的屈曲波数随半径/厚度比增大而增大,不同材料套筒的屈曲波数随屈服强度/塑性强化模量比值增大而增大,而相同材料的屈曲波数不随加载电流的大小而改变;模拟计算结果与实验结果基本吻合。     

压电板屈曲和后屈曲的有限元分析

采用板的一附剪切理论,并考虑结构的几何非线性,基于Ｔｏｔａｌ Ｌａｇｒａｎｇｅ方法,建立了在弱非线性耦合假设下压电智能的有限元控制方法,分析了没边界条件下压电板的屈曲和后屈曲,计算结果表明,在单向压力作用下,材料的压电效应和外加电压对板的屈曲载荷及后屈曲影响很小;而电场对位移加载简支板的屈曲影响较显著。该文为压电材料的工程应用提供了理论指导,同时提供了一种有效的有限元分析方法。     

圆柱约束下细长压杆屈曲及后屈曲行为

研究了轴向力作用下受圆柱横向约束的弹性杆的屈曲和后屈曲行为。通过有限元模拟,分析了细长压杆发生正弦和螺旋屈曲的轴向临界力,提出了正弦和螺旋屈曲临界点判定方法,且与文献结果比较验证了分析的正确性。同时考察了长细比和边界条件等因素对临界力的影响。结果表明正弦屈曲临界力和螺旋屈曲临界力随杆的长细比减小而增大。杆在一定长度范围内,端部约束条件对临界屈曲载荷及模态的影响不可忽略;当杆的长度足够长时,可以忽略边界条件对临界载荷的影响。     

海底管道的弹塑性稳定性和屈曲传播

本文根据我国南海海底输油气铺管的要求系统分析了海底管道在弯曲与外水压力共同作用下的弹塑性稳定性和屈曲传播.研究了管道的极值型屈曲和分枝型屈曲.在考虑管道的初始非圆度和材料的物理非线性的情况下提出了临界屈曲载荷的计算方法.综述和评论了屈曲传播现象的本质和各种计算方法.介绍了我们所进行的全尺寸管道实验.在分析理论结果时与现行的有关设计规程进行了比较和评论.     

薄膜褶皱的非线性屈曲有限元分析

有限元模拟已经成为目前薄膜褶皱预测的重要工具。本文采用ANSYS SHELL63单元对薄膜受剪情况下的褶皱形变进行了非线性屈曲分析。通过本征屈曲分析得到的模态模拟了薄膜的初始缺陷。利用本文的非线性有限元模型分析得到了薄膜褶皱的波长和幅度,并与理论分析结果进行了比较,结果相近。     

组合圆柱壳静态屈曲的几个影响因素分析

本文的研究是以“九五”国家科技攻关专题“快堆主钠池堆芯抗震性能的安全评价方法研究、中的屈曲研究子课题为背景展开的,所用模型简化自中国实验快堆钠池主容,是由薄壁圆柱壳和多个加劲肋结构而成的组合圆柱壳体,本文分别采用大型有限元程度ANSYS 5．4和LGOR FEAS（SUPER SAP 93）,对该壳体进行了常温时水平,轴向荷载共同作用静态屈曲的计算,同时考虑了诸如塑性,边界条件及初始缺陷等因素的影响,并进行了相关实验研究,最后将有限元计算结果与实验所获得的静屈曲荷载进行比较,结果吻合较好。     

非线性地基上弹性梁屈曲传播抑制的动力有限元分析

基于非线性地基弹性梁屈曲传播的动态控制方程，考虑梁质量惯性的影响，用有限元和时间积分相结合的方法，计算了分析了止屈器对梁屈曲传播抑制的作用，止屈器是通过局部加大地基刚度加以模拟。研究表明：考虑动态得到的止屈效果和前人以准静态为前提的结果相比有明显不同。     

无拉力弹性地基上矩形板屈曲/后屈曲问题的辛求解方法

无拉力弹性地基上矩形薄板的屈曲/后屈曲问题是板壳力学中一类重要课题,在工程中有着大量应用.因涉及接触非线性,目前主要采用数值方法对该类问题进行求解,发展具有重要基准价值的解析方法是当前面临的一项挑战.针对上述问题,本文将板划分为若干包含强制边界条件的板,形成子问题,在辛空间下利用分离变量与辛本征展开对子问题进行解析求解,通过子问题边界处的连续条件确定板与地基的接触状态;通过迭代求解上述过程,获得子问题划分的收敛结果,并得到最终屈曲载荷及模态.结果表明,无拉力弹性地基与Winkler地基上板的屈曲行为存在显著差异,且无拉力弹性地基的刚度对板的屈曲载荷与屈曲模态均有重要影响.在此基础上,结合Koiter摄动法与辛方法,对无拉力弹性地基上矩形板的后屈曲问题进行求解,获得板的后屈曲平衡路径.所得到的屈曲与后屈曲分析结果均与有限元计算结果吻合良好,确认了本文结果的正确性.由于本文方法数学推导严格,求解效率高,因此不仅为研究无拉力弹性地基上矩形薄板的屈曲/后屈曲行为提供了一种有价值的理论工具,更有望拓展至更多复杂板壳力学问题的求解.     

基于弧长法的加筋板后屈曲特性分析及试验

作为飞机上重要的承载部件,加筋壁板在发生初始屈曲后仍具有较强的后屈曲承载能力,因此研究其后屈曲特性对于确定破坏载荷具有重要意义。传统的特征值屈曲分析是以小位移小应变的线弹性理论为基础的,且不考虑结构在受载过程中结构构形的变化,因此误差较大。本文采用Riks弧长法,结合材料弹塑性理论对铝合金整体加筋壁板轴压加载后的屈曲破坏过程、传载机制、极限载荷进行了研究,并进行了轴压加载的试验验证,得到了加载过程中的应力、应变曲线以及极限载荷,还对后屈曲破坏形式进行了分析。数值模拟结果表明:本文研究的整体加筋板初始屈曲发生在蒙皮,后屈曲过程筋条是主要的承载部位,与试验中观察到的现象一致;试验中加筋板最终破坏部位发生在筋与蒙皮连接处,有限元模拟结果与试验中加筋板的最终破坏部位一致;数值模拟得到的极限载荷与试验的相对误差在5%以内。这表明基于弧长法的后屈曲计算能够准确跟踪整体加筋板的后屈曲平衡路径和预测极限载荷。     

热过屈曲功能梯度壁板的气动弹性颤振

功能梯度材料的宏观物理性能随空间位置连续变化,能充分减少不同组份材料结合部位界面性能的不匹配因素.功能梯度壁板用作高速飞行器的热防护结构,能有效消除气动加热带来的壁板内部热应力集中.本文考虑热过屈曲变形引入的结构几何非线性,分析功能梯度壁板的气动弹性颤振边界.基于幂函数材料分布假设,采用混合定律计算功能梯度材料的等效力学性能.根据一阶剪切变形板理论、冯·卡门应变-位移关系和一阶活塞理论,基于虚功原理建立超声速气流中受热功能梯度壁板的非线性气动弹性有限元方程.采用牛顿-拉弗森迭代法数值求解壁板的热屈曲变形,分析超声速气流对热屈曲变形的影响机理.在壁板热过屈曲的静力平衡位置分析动态稳定性,确定了壁板的颤振边界.研究表明,当陶瓷-金属功能梯度壁板的组份材料沿厚度方向梯度分布时,会破坏结构的对称性导致壁板在面内热应力作用下发生指向金属侧的热屈曲变形.超声速气流中壁板热屈曲变形最大的位置随气流速压增大向下游推移,并伴随屈曲变形量的减小.热过屈曲壁板的几何非线性效应会提高壁板的颤振边界,这种影响在高温、低无量纲速压且壁板发生大挠度热屈曲变形时表现显著.较高无量纲气流速压下由于壁板的热屈曲变形被气动力限定在小挠度范围,几何非线性效应不明显.      

Analysis model to simulate the cracked pipe buried in soil

This paper describes the use of the finite element method to simulate crack propagation in a pipeline in the presence of soil backfill. Since typical values for soil modulus are known, a spring element can be used to simulate the effects of backfilled soil on crack propagation. This consists of the relationship between the soil property and the stiffness of the spring element, and the effect of backfill depth on the crack propagation in the pipes. By equating crack driving force,G, in the presence of soil, to the dynamic fracture toughness of the pipeline material, a theoretical predictive capability is obtained for the first time for the gas transmission and distribution pipelines. Numerical calculations have also been carried out for the cracked pipe with backfill and without backfill.     

承压热冲击下RPV承载力影响因素的数值分析

为研究承压热冲击下含缺陷的核反应堆压力容器(RPV)的极限承载力,借助ABAQUS软件建立含表面裂纹的RPV的3D有限元模型.通过热-力耦合分析,分别得到简体区和喷嘴区的瞬态温度场和应力场.考虑到实际RPV断裂试验工作具有相当大的危险性,使用扩展有限元法(XFEM)对裂纹的扩展情况进行数值模拟.通过对含有不同位置和尺寸...     

Chater梁屈曲传播的动态分析

本文从非线性有限元的一般方法出发,结合应用加权残数技术和ALE边界元法,对Chater梁的屈曲传播进行了完全的动态分析,抛弃了前人一直采用的准静态和动稳态假设;另外,基于势流理论和非线性液面条件,还考虑了结构充液后由于流固耦合作用对屈曲传播所造成的影响,大量算例揭示了结构的菲线性、作用载荷、屈曲传播速度以及液体介质之间的内在关系,从而对屈曲传播有了更深入的认识。     

Material combinations and parametric study of thermal and mechanical performance of pyramidal core sandwich panels used for hypersonic aircrafts

A novel kind of lightweight integrated thermal protection system, named pyramidal core sandwich panel, is proposed to be a good safeguard for hypersonic aircrafts in the current study. Such system is considered as not only an insulation structure but also a load-bearing structure. In the context of design for hypersonic aircrafts, an efficient optimization should be paid enough attention. This paper concerns with the homogenization of the proposed pyramidal sandwich core panel using two-dimensional model in subsequent research for material selection. According to the required insulation performance and thermal–mechanical properties, several suitable material combinations are chosen as candidates for the pyramidal core sandwich panel by adopting finite element analysis and approximate response surface. To obtain lightweight structure with an excellent capability of heat insulation and load-bearing, an investigation on some specific design variables, which are significant for thermal–mechanical properties of the structure, is performed. Finally, a good balance between the insulation performance, the capability of load-bearing and the lightweight has attained.     

SBFEM与非局部宏微观损伤模型相结合的准脆性材料开裂模拟

混凝土是一种被广泛应用于土木和水利工程中的准脆性材料, 在各种内外部因素的作用下, 开裂是混凝土结构最为普遍的破坏形式, 准确模拟结构的开裂过程, 对于结构的安全评估至关重要. 将比例边界有限元与非局部宏微观损伤模型相结合提出一种准脆性材料开裂模拟新方法. 以比例边界有限元子域的比例中心作为物质点, 通过两比例中心(物质点对)之间的物质键的正伸长率来定义微细观损伤, 将某点影响域内物质键的微细观损伤加权平均得到该点的宏观拓扑损伤. 再引入能量退化函数, 将宏观拓扑损伤嵌入到比例边界有限元的基本框架中. 充分利用比例边界有限元网格允许存在悬挂节点的优势, 采用四叉树网格离散技术进行快速、高质量的多级网格划分与过渡. 通过一个I型开裂与一个混合型开裂的两个典型算例, 验证了该方法可捕获结构裂纹扩展路径与荷载变形曲线. 与现有的方法相比, 本文的损伤模型可得到更准确的局部开裂损伤带, 结果更为合理, 且具有更高的计算精度和计算效率. 当损伤过程区网格尺寸小于影响域半径的1/5时, 计算结果不存在网格敏感性问题.      

Initial thermo-mechanical post-buckling of beams with temperature-dependent physical properties

The objective of this work is to analyze the elastic buckling and initial post-buckling behavior of slender beams subjected to uniform heating. The beams are assumed to be double-hinged with fixed ends, preventing thermal expansion. Consequently, destabilizing compressive forces arise that may lead to beam buckling. When the temperature is further increased, the beam experiences finite displacements, with the result that the analysis is geometrically non-linear. The modulus of elasticity and the thermal induced strain, key material properties for this problem, are temperature-dependent. Thus, the coefficients of the governing equations are not constant. This suggests the physical non-linearity of the mathematical model. Hence, the analysis is geometrically and physically non-linear. The analysis is sensitive to the beam initial temperature, as the thermal strain is a function of the initial and final temperatures. The material is considered to be linear elastic, and consequently viscoelastic and plastic effects are not taken into account. Furthermore, the beam cross-section properties are assumed to be constant, which is consistent with the small strain formulation. A perturbation method is applied to the governing non-linear differential equations so that the initial post-buckling behavior may be analytically determined when temperatures above the critical temperature are applied to the beam. To illustrate the application of the formulation we present a case study for the aluminum 7075-T6 alloy, a material commonly used in aerospace and naval industries. Nonetheless, it is expected similar behavior for other metallic materials. The curves that define the variation of the modulus of elasticity, the thermal strain and the yield stress with temperature are considered in our analysis. The change in length, reaction forces at the supports and geometric configurations are obtained as a function of temperature and the beam slenderness ratio. The critical buckling loads and temperatures and the initial post-buckling analysis are also calculated in the context of the temperature-independent physical properties. Our results emphasize the importance of modeling the material's non-linearity if accuracy is required. However, from a practical application point of view results are acceptable if temperature-independent physical properties are employed, especially for large slenderness ratios.     

Asymptotic solutions for buckling delamination induced crack propagation in the thin film-compliant substrate system

In a thin film-substrate system in-plane compressive stress is commonly generated in the film due to thermal mismatch in operation or fabrication process. If the stress exceeds a critical value, part of the film may buckle out of plane along the defective interface. After buckling delamination, the interface crack at the ends may propagate. In the whole process, the compliance of the substrate compared with the film plays an important role. In this work, we study a circular film subject to compressive stress on an infinitely thick substrate. We study the effects of compliance of the substrate by modeling the system as a plate on an elastic foundation. The critical buckling condition is formulated. The asymptotic solutions of post-buckling deformation and the corresponding energy release rate of the interface crack are obtained with perturbation methods. The results show that the more compliant the substrate is, the easier for the film to buckle and easier for the interface crack to propagate after buckling.     

Finite element analysis of post-buckling dynamics in plates—Part I: An asymptotic approach

Various static and dynamic aspects of post-buckled thin plates, including the transition of buckled patterns, post-buckling dynamics, secondary bifurcation, and dynamic snapping (mode jumping phenomenon), are investigated systematically using asymptotical and non-stationary finite element methods. In part I, the secondary dynamic instability and the local post-secondary buckling behavior of thin rectangular plates under generalized (mechanical and thermal) loading is investigated using an asymptotic numerical method which combines Koiter’s nonlinear instability theory with the finite element technique. A dynamic multi-mode reduction method—similar to its static single-mode counterpart: Liapunov–Schmidt reduction—is developed in this perturbation approach. Post-secondary buckling equilibrium branches are obtained by solving the reduced low-dimensional parametric equations and their stability properties are determined directly by checking the eigenvalues of the resulting Jacobian matrix. Typical post-secondary buckling forms—transcritical, supercritical and subcritical bifurcations are observed according to different combinations of boundary conditions and load types. Geometric imperfection analysis shows that not only the secondary bifurcation load but also changes in the fundamental post-secondary buckling behavior are affected. The post-buckling dynamics and the global analysis of mode jumping of the plates are addressed in part II.     

Incubation time fracture criterion for FEM simulations

The paper is discussing problems connected with embedment of the incubation time criterion for brittle fracture into finite element computational schemes.Incubation time fracture criterion is reviewed;practical questions of its numerical implementation are extensively discussed.Several examples of how the incubation time fracture criterion can be used as fracture condition in finite element computations are given.The examples include simulations of dynamic crack propagation and arrest,impact crater formation(i.e.fracture in initially intact media),spall fracture in plates,propagation of cracks in pipelines.Applicability of the approach to model initiation,development and arrest of dynamic fracture is claimed.