一种快速求解含缺陷薄壳结构屈曲载荷的数值振动方法与流程

本发明属于航空航天中主承力薄壁构件快速设计领域，提供一种快速求解含缺陷薄壳结构屈曲载荷的数值振动方法。

背景技术：

薄壳结构广泛用于航空航天承力部件，如运载火箭燃料贮箱、飞机机身、航空发动机加力筒体等。典型的薄壳结构包括柱壳、锥壳、球壳、组合壳等。典型的服役工况包括轴压、外压、扭转、轴内压、轴外压、轴扭转等。在上述工况下，薄壳结构的主要失效模式为屈曲失稳。在薄壳结构的加工、装配、运输等过程中，不可避免地产生几何缺陷，其会对薄壳结构的屈曲载荷值造成较大折减，影响薄壳结构的服役性能。

实验技术虽然可以用来准确的评估薄壳结构屈曲载荷，但实验周期长、成本高，无法用于结构设计阶段，更适合用于后期验证。近年来，随着计算力学的快速发展，数值分析技术在薄壳结构屈曲载荷预测方面发挥了越来越重要的作用。代表性的方法包括特征值屈曲分析方法、弧长法、牛顿-拉夫森方法、显式动力学方法等，上述方法各有优劣。例如，特征值屈曲分析方法计算效率高，但只能分析薄壳结构的线性屈曲载荷、无法获得后屈曲极限承载载荷、且无法考虑缺陷的影响；弧长法及牛顿-拉夫森方法在求解含缺陷薄壳结构后屈曲极限承载载荷时存在收敛问题，常常导致分析失败；显式动力学方法可以稳定、准确地计算含缺陷薄壳结构后屈曲极限承载，与实验结果误差较小，但计算效率过低，无法满足结构设计大量迭代计算需求。因此，发明适用于含缺陷薄壳结构的准确高效的屈曲载荷预测技术，具有重要研究意义和实用价值。

技术实现要素：

本发明主要解决含缺陷薄壳结构屈曲载荷预测精度差、分析耗时长的问题，提出一种基于数值振动方法的薄壳结构屈曲载荷快速求解方法。计算不同外载荷水平下含缺陷薄壳结构的一阶负载固有频率值，并与对应的外载荷幅值进行拟合，将拟合值与完美薄壳结构线性屈曲载荷值相乘，即可预测出含缺陷薄壳结构的屈曲载荷值。

为了达到上述目的，本发明的技术方案为：

一种快速求解含缺陷薄壳结构屈曲载荷的数值振动方法，包括以下几个步骤：

第一步，计算完美薄壳结构的线性屈曲载荷值pcr，并将其作为后续计算的参考值。其中，可采用特征值屈曲分析、瑞利-里兹法等方法计算薄壳结构线性屈曲载荷值。所述薄壳结构的构型包括柱壳、锥壳、球壳、组合壳等；材料包括金属、复合材料等。

第二步，以一定比例α％的线性屈曲载荷值pcr作为外载荷上限值，以0作为外载荷下限值，对含缺陷薄壳结构逐步施加不同幅值的外载荷pi(共分n步进行加载，i＝1,2,…,n)，并分别计算含缺陷薄壳结构的一阶负载固有频率代表外载荷幅值为0时含缺陷薄壳结构的一阶负载固有频率。

其中，外载荷比例值α％的取值范围建议为50％-80％，逐步加载可以为等比例加载(例如，逐步加载0％*pcr、10％*pcr、20％*pcr、30％*pcr、40％*pcr、50％*pcr、60％*pcr)和非等比例加载(例如，逐步加载0％*pcr、35％*pcr、50％*pcr、65％*pcr)。一阶负载固有频率的计算可以采用精细模型的频率分析，也可以采用本征正交分解法(properorthogonaldecomposition，pod)、组合逼近法(combinedapproximation，ca)、并行计算等方法提高频率分析的效率。

所述含缺陷薄壳结构的缺陷形式包括实测缺陷、单点凹陷、多点凹陷、一阶模态缺陷等。所述外载荷工况包括轴压、外压、扭转、轴内压、轴外压、轴扭转等。

第三步，根据施加的外载荷和计算获得的一阶负载固有频率结果，对(1-pi/pcr)²与间的二次关系进行拟合，即确定的待定系数a、b、c，进而，求得的极小值

第四步，根据求得的pcr和ξ，可以得到数值振动方法预测的含缺陷薄壳结构的屈曲载荷pvct＝(1-ξ)pcr。

本发明的有益效果为：本发明提供一种快速求解含缺陷薄壳结构屈曲载荷的数值振动方法，区别于传统的特征值屈曲分析、显式动力学等屈曲分析方法，突破了含缺陷薄壳结构屈曲载荷预测传统方法精度差、分析耗时长的缺点。本发明提出的技术对于多种类型的薄壳结构构型(柱壳、锥壳、球壳、组合壳等)、不同的载荷工况(轴压、外压、扭转、轴内压、轴外压、轴扭转等)、不同类型的材料属性(金属、复合材料等)和不同类型的缺陷形式(实测缺陷、单点凹陷、多点凹陷、一阶模态缺陷等)，表现出较高的适用性，取得了较高的预测精度和效率。该技术可有效缩短航空航天薄壳结构研发周期。

附图说明

图1为一种快速求解含缺陷薄壳结构屈曲载荷的数值振动方法的实现流程图；

图2为外载荷与一阶负载固有频率结果二次拟合示意图。

具体实施方式

为使本发明解决的方法问题、采用的方法方案和达到的方法效果更加清楚，下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部内容。

图1为本发明实施例提供的一种快速求解含缺陷薄壳结构屈曲载荷的数值振动方法的实现流程图。如图1所示，本发明实施例提供的一种快速求解含缺陷薄壳结构屈曲载荷的数值振动方法包括：

第一步，建立含实测缺陷圆柱壳有限元模型。壳体模型直径d＝1000mm，壳体长度l＝600mm，蒙皮厚度ts＝1.5mm。为了便于加载和固定，在圆柱壳的两端设计t形环。t形环的厚度tring＝6.0mm，宽度wring＝30.0mm。圆柱壳采用各向同性的铝合金材料。两端t形环翼板采用固支边界条件。有限元模型的网格单元类型采用具有4个节点、每个节点具有5个自由度和缩减积分的四边形壳单元，简称为s4r5单元。经网格收敛性分析后确定网格大小的全局尺寸为8.0mm。所使用的计算工作站性能为intelxeone5-2687w@3.10ghz的cpu和64g的内存。

基于瑞利-里兹法计算完美圆柱壳结构的线性屈曲载荷值pcr＝596kn，并将其作为后续计算的参考值。

第二步，以一定比例70％的线性屈曲载荷值pcr作为外载荷上限值，以0作为外载荷下限值，对含缺陷圆柱壳结构逐步施加不同幅值的外载荷pi(共分8步进行加载，i＝1,2,…,8)，并分别计算含缺陷圆柱壳结构的一阶负载固有频率代表外载荷幅值为0时含缺陷圆柱壳结构的一阶负载固有频率。

其中，逐步加载采用等比例加载，即逐步加载0％*pcr、10％*pcr、20％*pcr、30％*pcr、40％*pcr、50％*pcr、60％*pcr、70％*pcr。一阶负载固有频率的计算采用基于本征正交分解法(properorthogonaldecomposition，pod)的降阶频率分析，可以提高频率分析的效率。首先，针对外载荷为0％*pcr和50％*pcr的含缺陷圆柱壳结构，开展基于精细模型的频率分析，获得其一阶负载固有频率及特征向量。然后，通过本征正交分解法将上述特征向量进行正交分解，获得降阶基，构建降阶频率分析方程。最后，基于降阶频率分析方程快速计算出其他幅度外载荷下含缺陷圆柱壳结构的一阶负载固有频率，从而获得全部外载荷下含缺陷圆柱壳结构的一阶负载固有频率

第三步，根据施加的外载荷和计算获得的一阶负载固有频率结果，对(1-pi/pcr)²与间的二次关系进行拟合，如图2所示。确定二次方程中待定系数，a＝2.7449、b＝-3.2529、c＝0.9971。进而，求得上式的极小值

第四步，根据求得的pcr和ξ，可以得到数值振动方法预测的含缺陷圆柱壳结构的屈曲载荷pvct＝(1-ξ)pcr＝487.3kn。含缺陷圆柱壳结构的屈曲实验结果为482.2kn。对比实验结果可知，本发明的预测误差仅为1％，验证了本技术的高预测精度。基于显式动力学方法计算的含缺陷圆柱壳结构的屈曲载荷为479.3kn，计算耗时为40分钟。对比显式动力学分析耗时，本发明的计算耗时为1.6分钟，节省了96％的计算耗时，验证了本技术的高预测效率。

最后应说明的是：以上各实施例仅用以说明本发明的方法方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通方法人员应当理解：其对前述各实施例所记载的方法方案进行修改，或者对其中部分或者全部方法特征进行等同替换，并不使相应方法方案的本质脱离本发明各实施例方法方案的范围。