小样本数据下圆柱薄壳初始缺陷不确定性量化的极大熵方法

具有不确定性特征的初始缺陷被认为是导致薄壳结构实际临界载荷值与理论解不相符并呈现分散特征的主要原因.对实际薄壳结构初始缺陷的建模至少需要考虑两个方面的不确定性量化,一是对缺陷分布形式和幅值等固有随机性的量化,二是对小样本量和不准确测量所导致缺陷统计量的不确定性的量化.本文在利用随机场的Karhunen-Loeve展开法对薄壳初始几何缺陷建模的基础上,提出一种基于极大熵原理的缺陷建模方法.首先,采用极大熵分布来估计Karhunen-Loeve随机变量的概率密度函数,以适应不能使用高斯随机场进行缺陷随机场建模的情况.随后,通过将经典的等式约束极大熵模型扩展为区间约束极大熵模型,实现对实际工程中仅能获得少量薄壳结构几何缺陷样本数据所导致的认知不确定性的量化.最后,将所提方法用于对国际缺陷数据库的A-Shell进行缺陷建模和临界载荷预测.研究表明,所提基于区间约束极大熵原理的随机场建模方法在能够有效表征实测数据高阶矩信息的同时,还具备量化小样本数据导致的认知不确定性的能力,并且高斯随机场模型和基于等式约束极大熵原理的随机场模型是本文所提建模方法的两种特殊情况.     

考虑周长约束的圆柱薄壳初始几何缺陷随机场建模方法

利用随机场对圆柱薄壳结构的初始几何缺陷进行建模,并据此建立了一种用于含初始几何缺陷轴压圆柱薄壳屈曲分析的随机分析方法。首先,指出已有将圆柱薄壳初始几何缺陷表征为二维高斯随机场的方法会导致与实际不相符的初始几何缺陷,如圆柱周长显著增大或缩小的几何缺陷。其次,提出一种考虑周长不变约束的随机场建模方法,以剔除与实际不相符的随机几何缺陷。最后,基于所建立的初始几何缺陷随机场模型,利用非干涉多项式混沌展开法进行圆柱薄壳的随机屈曲分析,给出临界屈曲载荷的概率分布。数值试验结果表明,基于随机场理论的初始几何缺陷建模方法可有效刻画几何缺陷对结构承载能力的影响,而提出的约束随机场建模方法又能有效减小结果的分散性。     

关于圆筒壳稳定性中的初始缺陷

本文综述了圆筒壳稳定性中初始几何缺陷的几个重要问题。讨论了初始几何缺陷对圆筒壳临界载荷的基本效应,它的随机性、测量以及相应的数据库的建立。最后介绍了失稳破坏占主要地位的壳体结构的交互设计的基本思想和方法,并对进一步的研究工作提出了看法。      

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

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

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

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

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

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

含轴对称初始缺陷的具有正交各向异性表层的夹层圆柱壳的非线性屈曲

夹层圆柱壳具有很高的结构效能。在许多工程结构中被广泛采用。本文研究分析了含有轴对称初始缺陷的夹层圆柱壳在轴压下的非线性屈曲问题。该夹层壳具有正交各向异性表层和各向同性可承剪的夹心．利用Stein的前屈曲一致理论得出了前屈曲挠度随轴向载荷及缺陷参数的变化情况，运用Galerkin法导出了屈曲控制方程，并进行了数值计算，得到了屈曲载荷、缺陷幅值、缺陷波数、夹心模量等参量之间的关系．结果表明与壳体实际屈曲模态相同的初始缺陷具有很大的危害性，可以通过增加壳体表层的轴向弹性模量或优化夹心的有关参数等途径来提高屈曲载荷，改善壳体屈曲性能。     

基于弧长法的受压球壳稳定性分析

对受压球壳进行特征值屈曲分析,得到了前6阶屈曲模态及线性屈曲临界载荷;采用弧长法进行非线性有限元分析,对理想球壳施加初始扰动,通过2次扰动值折半的方法求得引起结构屈曲的最小扰动值,追踪到了屈曲分支点和全过程载荷-位移路线。基于前6阶屈曲模态位移,在受压球壳中分别引入2.5mm和1mm两种缺陷值,分析缺陷对球壳屈曲特性的影响。结果表明:取壳厚的0.5%即0.05mm时,得最小扰动值,近似模型与完善结构极值载荷的差值为0.93%;球壳是缺陷敏感性结构,缺陷的幅值和分布都对其极限载荷有影响,缺陷幅值与厚度比为0.1时,缺陷球壳承载力相对理想结构下降了约11%,缺陷幅值与厚度比为0.25时,承载力相对下降了约30%,说明提高球壳稳定性需要提高球壳加工精度。     

不确定初始几何缺陷壳的动态屈曲分析

研究了动载作用下含有不确定初始几何缺陷壳的动态屈曲问题.给出了不确定初始几何缺陷在椭球描述和区间描述下,壳的动态响应和安全因子的界限估计.从数学证明和数值算例两方面,讨论了基于不确定性初始几何缺陷两种描述的凸模型方法和区间分析方法所得壳动态响应及安全因子的界限的关系,为判断壳的动态屈曲失效提供了依据.     

含初缺陷裂纹损伤梁的冲击动力屈曲

由Hamilton原理导出考虑初始缺陷及横向剪切变形时裂纹梁的动力屈曲控制方程;应用断裂力学中常用的线弹簧模型将裂纹引入到屈曲控制方程中;基于B-R动力屈曲判断准则,采用数值方法求解了受轴向冲击载荷作用时裂纹梁的动力屈曲;对比讨论了不同冲击速度、初始几何缺陷大小以及分布形式等因素对梁冲击动力屈曲的影响。     

Effect of multiple localized geometric imperfections on stability of thin axisymmetric cylindrical shells under axial compression

Stability of imperfect elastic cylindrical shells which are subjected to uniform axial compression is analyzed by using the finite element method. Multiple interacting localized axisymmetric initial geometric imperfections, having either triangular or wavelet shapes, were considered. The effect of a single localized geometric imperfection was analyzed in order to assess the most adverse configuration in terms of shell aspect ratios. Then two or three geometric imperfections of a given shape and which were uniformly distributed along the shell length were introduced to quantify their global effect on the shell buckling strength. It was shown that with two or three interacting geometric imperfections further reduction of the buckling load is obtained. In the ranges of parameters that were investigated, the imperfection wavelength was found to be the major factor influencing shell stability; it is followed by the imperfection amplitude, then by the interval distance separating the localized imperfections. In a wide range of parameters this last factor was recognized to have almost no effect on buckling stresses.     

Modified Monte Carlo method for buckling analysis of nonlinear imperfect structures

In this paper, we propose a modified Monte Carlo method for analysis of buckling of an imperfect beams on softening nonlinear elastic foundation. Such structures exhibit considerable imperfection sensitivity, i.e. reduction in the maximum load that the structure is able to support in contrast to classical buckling load of the perfect structure. The initial imperfections are treated as random functions of axial coordinate. In order to reduce the needed number of simulations, the Monte Carlo method is coupled with maximum likelihood methodology and the Kolmogorov–Smirnov test.     

Non-linear stability analysis of imperfect thin-walled composite beams

The static non-linear behavior of thin-walled composite beams is analyzed considering the effect of initial imperfections. A simple approach is used for determining the influence of imperfection on the buckling, prebuckling and postbuckling behavior of thin-walled composite beams. The fundamental and secondary equilibrium paths of perfect and imperfect systems corresponding to a major imperfection are analyzed for the case where the perfect system has a stable symmetric bifurcation point. A geometrically non-linear theory is formulated in the context of large displacements and rotations, through the adoption of a shear deformable displacement field. An initial displacement, either in vertical or horizontal plane, is considered in presence of initial geometric imperfection. Ritz's method is applied in order to discretize the non-linear differential system and the resultant algebraic equations are solved by means of an incremental Newton-Rapshon method. The numerical results are presented for a simply supported beam subjected to axial or lateral load. It is shown in the examples that a major imperfection reduces the load-carrying capacity of thin-walled beams. The influence of this effect is analyzed for different fiber orientation angle of a symmetric balanced lamination. In addition, the postbuckling response obtained with the present beam model is compared with the results obtained with a shell finite element model (Abaqus).     

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.     

初始几何缺陷对网壳截面优化结果影响研究

单层网壳结构是缺陷敏感结构,初始几何缺陷不同对"网壳截面优化结果"有显著影响,该文研究了网壳杆件截面取值与缺陷分布、结构极限承载力之间的关系.结果表明:经截面优化设计的网壳结构最不利应力最早并且一直出现在初始几何缺陷较大处;优化寻优结果与初始几何缺陷分布有直接关系,增大缺陷较大处杆件截面能显著提高结构承载力;由于缺陷随机分布,经截面优化设计的网壳结构须校核其在不同缺陷下的稳定承载力以确保结构安全.     

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

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