含局部减薄弯头塑性极限载荷的数值分析

局部减薄是一种压力管道表面常见的体积型缺陷，它不仅会降低管道的承载能力，而且还可能引起管道破坏，导致严重的事故．采用数值分析方法，对内压和面内弯矩作用下含局部减薄弯头的极限载荷进行了研究，分析了载荷组合、缺陷形状、位置、尺寸对弯头极限承载能力的影响，讨论了导致弯头破坏的典型失效模式．提出了含局部减薄缺陷弯头的塑性极限载荷工程估算式．计算结果为含缺陷弯头的安全评定提供了理论依据．     

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

为研究承压热冲击下含缺陷的核反应堆压力容器(RPV)的极限承载力,借助ABAQUS软件建立含表面裂纹的RPV的3D有限元模型．通过热-力耦合分析,分别得到筒体区和喷嘴区的瞬态温度场和应力场．考虑到实际RPV断裂试验工作具有相当大的危险性,使用扩展有限元法(XFEM)对裂纹的扩展情况进行数值模拟．通过对含有不同位置和尺寸裂纹、具有不同壁厚容器的极限承载力分析,得到满足强度条件的各临界裂纹尺寸．结果表明,当表面裂纹距喷嘴口的距离达到一定程度,其长度对RPV结构的极限承载力不构成影响．在喷嘴区,随裂纹距喷嘴口距离的减小,临界裂纹尺寸迅速减小,喷嘴口上仅允许深长比较小的浅表面裂纹出现．     

含裂纹薄壁圆柱壳弹塑性失稳极限荷载的解析解

为了得到含环向裂纹压弯薄壁圆柱壳失稳的极限荷载,基于Je?ek解析法,先利用裂纹面的平衡方程、变形几何关系和物理方程,得到在压弯荷载作用下圆柱壳轴力与裂纹面挠度的关系式;再利用极值条件,推导出在轴压和弯矩共同作用下含裂纹薄壁圆柱壳失稳时极限荷载的解析表达式。利用数值计算分析裂纹长度、长细比以及端部施加弯矩对圆柱壳极限荷载、弹性区高度和圆柱壳裂纹截面处挠度的影响,通过有限元数值结果验证极限荷载解析解的准确性。计算结果表明:与无裂纹完善构件相比,弯矩对含裂纹圆柱壳的极限荷载Pu的影响最为显著,其降低幅值可达3.78%～23.18%;弯矩的增加也会导致构件弹性区高度逐渐减小,裂纹面处的挠度呈非线性增长。本文研究弥补了目前该类问题缺乏理论解析解的不足。     

混凝土中爆炸模拟的数值方法比较

应用AUTODYN商业软件对装药在混凝土靶介质中的爆炸过程进行数值模拟,对比了混凝土中爆炸模拟的不同方法。将各数值方法模拟得到的混凝土毁伤特征尺寸与相关实验数据进行比较,研究了各数值方法的准确性及描述混凝土毁伤特征的可靠性参数。结果发现：Euler算法计算的压力最大,SPH和Lagrange耦合算法计算的压力最小;Euler算法和SPH算法计算的压力衰减过程会发生比较显著的扰动。以不同的平均网格尺寸对数值模拟结果的影响分析了各种数值方法的稳定性。比较靶体上同一测点的压力时间历程表明网格变化对测点处压力影响不大。通过应用各种数值方法模拟再现混凝土中的爆炸过程,揭示了各种方法的优缺点,明确了现有商业软件模拟混凝土中爆炸的适用性。     

共线裂纹体尺寸对应力强度因子的影响

本文在文献[1]的基础之上,用文献[2]中提出的奇异准谐调元,进一步对二维裂纹体几何尺寸对一对共线裂纹应力强度因子的影响进行了大量的数值计算和研究,给出了计算表格、拟合曲线以及经验公式,定量地分析了裂纹体几何尺寸对共线裂纹应力强度因子的影响。为工程中的裂纹体的脆性断裂分析与评定,提供了必要的计算依据。     

S型功能梯度扁球壳非线性屈曲的渐近分析

本文利用渐近迭代法获得了边界弹性支撑的功能梯度扁球壳的非线性屈曲问题的理论解．假设材料组分体积分数沿壳体厚度方向呈sigmoid幂函数变化,边界上考虑一般的弹性支撑约束．基于经典的薄壳理论和几何非线性关系,导出了S型功能梯度扁球壳的非线性屈曲问题的控制方程．采用渐近迭代法通过两次迭代得到了无量纲挠度和均布荷载之间的非线性特征关系．通过与已有有限元方法和其他数值方法的结果对比,验证了本文解的有效性和高精度．同时,通过计算阐述了缺陷因子、材料参数、边界约束系数及特征几何参数对壳体临界屈曲荷载的影响．     

含有初始缺陷双曲扁壳的固有振动特性分析

以工程中常用的双曲壳结构如圆柱壳、球壳、双曲抛物壳为研究对象,利用薄扁壳理论,基于瑞利-里兹法和切比雪夫多项式求得了几种边界条件下的双曲扁壳的自由振动固有频率,并与ANSYS分析结果进行了对比,验证了该方法的适用性。详细研究了在不同边界条件下的双曲扁壳的几何参数、初始几何缺陷尺寸、初始几何缺陷密集程度对频率大小、频率转向、振型变化的影响。结果表明:随着壳体结构厚度的增加及曲率半径的减小,壳体的固有频率会增加;几何缺陷半波数及缺陷尺寸对频率影响情况较为复杂,并且会使系统发生频率转向问题,这些结果对于工程实际具有重要的理论指导意义。     

高层框筒和筒中筒结构的整体稳定分析

本文将有限元法用于高层筒 体结构的整体稳定分析。导出了截面含双对称轴的矩形筒的筒单元刚度矩阵和单元几何刚度矩阵。可用于框筒和筒中筒结构的整体稳定计算。计算结果验证了本文方法的精度。     

基于模态缺陷的受压球壳屈曲特性研究

针对受压球壳非线性屈曲过程，对含初始缺陷受压球壳的稳定性进行研究。根据EN1993-1-6（2007）规范，给出不同制造等级壳体的等效缺陷值计算方法；基于线性特征值分析的模态构型给出初始缺陷的分布。利用非线性有限元弧长法对球壳受压失稳过程进行数值模拟，得到屈曲前后球壳变形情况及全过程载荷-位移曲线。计算一致缺陷模态法和N阶缺陷模态法对应的球壳屈曲临界载荷，结果表明，受压球壳对缺陷较敏感，承载能力随缺陷值增大而降低；一致缺陷模态法计算便捷，在工程应用上具有合理性，N阶缺陷模态法考虑高阶模态缺陷构型，结果更加全面，可以为工程中缺陷结构稳定性设计提供参考。     

均布压力下正交异性圆锥扁壳的后屈曲分析

基于Ｋａｒｍａｎ型大挠度方程,用修正迭代法分析了均布压力下夹支正交异性圆锥扁壳的几何非线性的后屈曲行为,给出二阶近似的荷载挠度特征关系式及临界荷载,给出了三种正交异性参数对应数值结果,分析了正交异性参数对壳体变形和屈曲荷载的影响。     

NUMERICAL INVESTIGATION OF THE LIMIT LOADS FOR PRESSURE VESSELS WITH PART-THROUGH SLOTS

A numerical investigation of the limit loads is carried out for pressurevessels with part-through slots using a general computational method for the limitanalysis of 3-D structures.The limit pressures are given for a comprehensive range ofgeometric parameters.Some of the calculated results are compared with the results of3-D elastic-plastic finite element analysis and existing numerical solutions.The effectsof various shapes and sizes of part-through slots on the load carrying capacity ofcylindrical shells are investigated and evaluated.Two kinds of typical failure modescorresponding to different dimensions of slots are studied.Based on the numericalresults,a geometric parameter G which combines the slot dimensions and the cylindergeometry is presented.It reasonably reflects the overall effect of slots on the limit loadsof cylinders.An empirical formula for estimating the limit pressures of cylindricalshells with part-through slots is obtained.     

易拉罐在轴-侧-扭-内压联合作用下的屈曲地貌

柱壳结构广泛应用于各个领域, 但由于其对初始缺陷较为敏感, 容易发生灾难性的屈曲失稳. 本文利用非线性有限元分析程序ABAQUS研究了柱壳屈曲问题, 并应用到了易拉罐的屈曲分析. 首先采用数值模拟的方法验证了Virot等学者的易拉罐屈曲试验结果, 然后为了获得屈曲的一些普适结果, 进一步考察了柱壳的屈曲表现. 对柱壳结构在不同载荷组合、不同几何参数作用下进行了细致分析. 为了讨论的直观, 本文绘制了柱壳结构在受到侧压-轴压载荷作用下外力-屈曲载荷-位移三维屈曲地貌图(称为landscape). 结果表明: 在侧压-轴压-扭转载荷作用下, 试件力-位移曲线出现了"cliff"(断崖)现象; 扭转载荷的施加不利于试件整体稳定性, 并造成了试件对初始缺陷的敏感性; 对于受到轴压-扭转载荷作用的试件, 本文定义承载力为零的平面为"sea level"(海平面)来区分试件破坏模式; 通过对不同边界条件的试件进行分析, 发现试件两端固定可以有效地增加结构的承载能力, 提高稳定性. 对柱壳结构内部充气可以大幅度提升结构的承载能力和稳定性, 减小对初始缺陷的敏感度.      

含缺陷轴对称体的安定与极限分析

利用静力安全定理得到了计算轴对称体安定与极限载荷的统一格式，采用温度参数法构造安定分析所需的残余应力场，为了克服对工程实际问题进行安定分析时解题规模与计算精度的矛盾，针对轴对称体的特点，采用两种线性化方案对屈服面进行线性化处理，即直接内接法和在降维应力偏量空间中对屈服面的线性化处理，使安定分析转化为一线性规划问题，在简化过程中合理选择线性化方案以便使应力校核点接近精确的屈服面；为了减小计算量，在求     

Axial load capacity of cylindrical shells with local geometric defects

The effect of local geometric defects on the buckling load of axially compressed thin circular cylindrical shells is investigated experimentally. Defects, in the form of diamond-shaped local dimples, similar to the buckles of the Yoshimura pattern, were introduced in otherwise near perfect isotropic epoxy shells by locally heating the shell wall. The behavior of the defects under load was monitored optically using a special whole-field grid-reflection technique. The effects of variations in shell geometry and defect size are also investigated. In general, the results indicate that the effect of local diamond-shaped defects on the stability of the axially loaded cylinder is not as detrimental as that of the global initial imperfections hitherto investigated.     

考虑材料损伤累积单层柱面网壳在强震下的失效研究

采用考虑钢材损伤累积的本构模型,对跨度15米单层柱面网壳在强震下的失效进行了系统的研究。在大量算例的基础上,考察了屋面质量、矢跨比、初始几何缺陷和长宽比等结构几何参数的变化对柱面网壳失效特征和失效极限荷载的影响,探讨了结构的强震失效模式及失效机理。对大量强度破坏算例在失效时刻的特征响应进行统计分析,建立了考虑材料损伤累积的单层柱面网壳动力损伤模型,可对结构在不同荷载强度下的损伤程度进行评估;提出了单层柱面网壳强震失效判别准则,用于判别单层柱面网壳在强震下的失效极限荷载。     

考虑形位公差的薄壁筒壳折减因子预测方法研究

缺陷敏感性是薄壁筒壳结构设计所面临的主要问题之一,通常利用折减因子来量化筒壳结构的缺陷敏感性程度。然而现有缺陷敏感性分析方法大多以预测筒壳折减因子下限为目的,未考虑不同形位公差水平对筒壳折减因子的影响。针对此问题,本文提出了一种考虑形位公差的薄壁筒壳折减因子预测方法。该方法基于多点最不利扰动载荷法进行最不利缺陷搜索,获得筒壳不同形位公差下的折减因子下限值,从而确定考虑形位公差的薄壁筒壳折减因子参考值,并利用不完全折减刚度法对计算过程进行加速。算例结果表明本文提出的筒壳折减因子预测方法可在保证安全可靠的前提下,有效提高折减因子预测精度,消除筒壳结构设计过程中不必要的安全裕度,对结构减重有积极意义。未来可基于本方法展开我国新一代航天结构薄壁筒壳折减因子设计规范的研究工作,进一步提高我国航天筒壳结构设计的精细化和轻量化水平。     

Buckling analysis of cylindrical shells with random geometric imperfections

In this paper the effect of random geometric imperfections on the limit loads of isotropic, thin-walled, cylindrical shells under deterministic axial compression is presented. Therefore, a concept for the numerical prediction of the large scatter in the limit load observed in experiments using direct Monte Carlo simulation technique in context with the Finite Element method is introduced. Geometric imperfections are modeled as a two dimensional, Gaussian stochastic process with prescribed second moment characteristics based on a data bank of measured imperfections. (The initial imperfection data bank at the Delft University of Technology, Part 1. Technical Report LR-290, Department of Aerospace Engineering, Delft University of Technology). In order to generate realizations of geometric imperfections, the estimated covariance kernel is decomposed into an orthogonal series in terms of eigenfunctions with corresponding uncorrelated Gaussian random variables, known as the Karhunen-Loéve expansion. For the determination of the limit load a geometrically non-linear static analysis is carried out using the general purpose code STAGS (STructural Analysis of General Shells, user manual, LMSC P032594, version 3.0, Lockheed Martin Missiles and Space Co., Inc., Palo Alto, CA, USA). As a result of the direct Monte Carlo simulation, second moment characteristics of the limit load are presented. The numerically predicted statistics of the limit load coincide reasonably well with the actual observations, particularly in view of the limited data available, which is reflected in the statistical estimators.     

Slot/slots air jet impinging cooling of a cylinder for different jets–cylinder configurations

Heat transfer characteristics of a slot/slots jet air impinging on a cylinder have been experimentally investigated for two different orientations of slot/slots jet plan with respect to cylinder axis. The experiments were carried out to study the effects of orientation of slot/slots jet plan with respect to cylinder axis on the rate of heat transfer from the cylinder. Two different jet–cylinder configurations were studied: (1) single slot jet aligned with cylinder axis (slot length = cylinder length), and (2) multiple slot jets equally spaces distributed orthogonal to cylinder axis (each slot length = cylinder diameter and sum of slots lengths = cylinder length). For each configuration, parametric effects of Reynolds numbers (Re) ranging from 1,000 to 10,000, dimensionless slot widths (W/d) ranging from 0.125 to 0.5, and dimensionless slot orifice-to-cylinder spacing (Z/W) ranging from 1 to 12 on local and average Nusselt numbers around cylinder surface have been investigated. The results showed that: (1) cooling the cylinder by multiple slots jets situated orthogonal to cylinder axis gave more uniform surface temperature distributions and higher heat transfer rate than the case of cooling the cylinder by single slot jet aligned with cylinder axis, (2) for both configurations the Nusselt number around the cylinder increased with increasing Re and W/d, and (3) for both configurations there was a certain Z/W in the range 4<Z/W<6 at which the stagnation and mean Nusselt number were maximum. Correlations for the mean Nusselt numbers around cylinder surface have been presented for both configurations. Comparisons between the correlations predictions and the present and other previous experimental data have been conducted.     

Repeated buckling and its influence on the geometrical imperfections of stiffened cylindrical shells under combined loading

The present experimental study aims at providing better inputs for improvement of the buckling load predictions of stiffened cylindrical shells subjected to combined loading. The work focuses on two main factors which considerably affect the combined buckling load of stiffened shells, namely geometric imperfections and boundary conditions. Six shells with nominal simple supports were tested under various combinations of axial compression and external pressure. The vibration correlation technique is employed to define the real boundary conditions. The geometric imperfections of the integrally stiffened shells are measured in the present experiments in situ and are used as inputs to a multimode analysis which yields the corresponding “knockdown” factor for various combinations of loading. Thus, when employing the repeated buckling procedure for obtaining interaction curves, each point on the curve is adjusted (using the multimode analysis) for the measured “new” surface of the shell and this results in more realistic interaction curves. The geometrical imperfections of the preloaded shells can also serve as an input to the International Imperfection Data Bank for future studies on the correlation between the manufacturing method of the shell and their geometric imperfections.     

A boundary layer theory for the buckling of thin cylindrical shells under external pressure

Based on the boundary layer theory for the buckling of thin elastic shells suggested in ref. [14]. the buckling and postbuckling behavior of clamped circular cylindrical shells under lateral or hydrostatic pressure is studied applying singular perturbation method by taking deflection as perturbation parameter. The effects of initial geometric imperfection are also considered. Some numerical results for perfect and imperfect cylindrical shells are given. The analytical results obtained are compared with some experimental data in detail, which shows that both are rather coincident.