一种薄壁复合材料舱体负刚度稳定性分析方法技术

随着对临近空间飞行器性能需求的提升，结构系统轻质化逐渐成为关键技术。为了实现飞行器结构系统轻质化，薄壁舱体结构被广泛应用，在满足强度需求的前提下，实现结构轻质化设计。然而薄壁舱段结构在大弯矩载荷工况及横法向振动条件下，蒙皮容易发生大变形，形成负刚度形式的全新几何包络。薄壁结构的负刚度变形导致舱体稳定特性发生大幅改变，传统基于线性化的稳定性计算方法不再适用。本发明专利技术一种薄壁复合材料舱体负刚度稳定性分析方法，建立基于几何非线性假设的有限元修正模型，采用弧长法对产生负刚度变形的薄壁结构进行稳定性预示，避免薄壁结构发生失稳带来的气动外形变化及强度破坏。强度破坏。强度破坏。

【技术实现步骤摘要】
一种薄壁复合材料舱体负刚度稳定性分析方法


[0001]本专利技术是一种通过识别薄壁舱体在振动条件下的负刚度变形分布，建立基于几何非线性假设的有限元修正模型，采用弧长法对产生负刚度变形的薄壁结构进行稳定性求解。

技术介绍

[0002]临近空间飞行器由于其具有重大的战略意义，当前已经成为国内外研究的热点之一。随着对飞行器性能及突防能力的需求，结构系统轻质化逐渐成为关键技术。
[0003]为了实现飞行器结构系统轻质化，薄壁舱体结构被广泛应用，尤其随着复合材料的日渐成熟与发展，薄壁结构的厚度尺寸可以做到更小，异形舱体结构的蒙皮壁厚甚至可以达到1mm，在满足强度需求的前提下，实现结构轻质化设计。
[0004]然而薄壁复合材料舱段结构在大弯矩载荷工况及横法向振动条件下，蒙皮容易发生大幅变形，存在形成负刚度形式的全新几何包络。薄壁结构的负刚度变形导致舱体稳定特性发生大幅改变，传统基于线性化的稳定性计算方法不再适用。

技术实现思路

[0005]本专利技术解决的技术问题是：克服现有技术的不足，提供一种薄壁复合材料舱体负刚度稳定性分析方法，避免薄壁结构发生失稳带来的气动外形变化及强度破坏。
[0006]本专利技术的技术方案是：一种薄壁复合材料舱体负刚度稳定性分析方法，步骤如下：
[0007]1)利用正弦扫频分析，获得结构变形分布U；
[0008]2)建立等效静力分析模型，根据步骤1)中获得的结构变形分布场U，获得等效静力过载；
[0009]3)根据步骤2)获得的等效静力过载，结合非线性静力分析模型，判断结构是否出现负刚度；若出现则继续步骤4)，若未出现则转至步骤6)；
[0010]4)提取非线性静力分析模型变形结果，修正结构稳定性分析模型；
[0011]5)对稳定性分析模型进行求解，获得薄壁结构发生失稳时的载荷比值和后屈曲的破坏状态，本方法结束；
[0012]6)建立小变形假设下的有限元分析模型，并对薄壁结构稳定性进行求解。
[0013]所述步骤1)获得结构变形分布的过程为：建立薄壁舱体有限元分析模型，对薄壁舱体有限元模型进行基于频域方法的正弦扫频计算分析，获得正弦扫频分析结果，并重点获取结构整体最大变形分布场U。
[0014]采用模态叠加法的稳态动力学分析方法，对薄壁舱体有限元模型进行基于频域方法的正弦扫频计算分析。
[0015]所述步骤2)获得等效静力过载的过程为：基于小变形假设，建立薄壁舱段稳态静力学有限元分析模型，对模型施加等效过载，使薄壁结构变形分布满足差值条件，以此确认等效过载量值。
[0016]所述差值条件为：minσ|U
g
‑
U|，即等效过载下的结构各有限元节点变形与正弦扫频分析获得的有限元节点变形相减后的方差值最小，以此确认等效过载量值；其中U
g
为静力学分析模型中的结构位移场，σ为方差。
[0017]所述步骤3)结合非线性静力分析模型，判断结构是否出现负刚度的过程为：基于几何非线性假设，即考虑应变二阶微分项，重新建立薄壁舱段稳态静力学有限元分析模型，加载确定后的等效过载，提取稳态静力结果中的薄壁结构节点位移场，通过变形分布识别是否出现负刚度变形。
[0018]所述步骤3)中负刚度变形判定方法如下：
[0019]31)以舱体经线为基准，按指定间距逐条提取经线上节点变形信息，获取经线上变形后的位置结果；
[0020]32)利用舱体旋成体构型特点，建立经线局部坐标系；
[0021]33)在经线局部坐标系下，以提取的经线两端节点变形后位置构造线性辅助线，并形成函数：f x,y＝ax+by+c；式中，a、b、c为线性方程系数，x、y为经线局部坐标系下坐标；
[0022]34)通过经线上的各节点位置信息到辅助线的矢量距离，判断结构是否发生负刚度变形。
[0023]所述步骤34)中，通过经线上的各节点位置信息x0,y0到辅助线的矢量距离判断结构是否发生负刚度变形，若始终d＞0，则结构未发生负刚度变形，若存在d＜0，则结构发生负刚度变形；式中x0、y0为结构壁面在经线坐标系下位置。
[0024]所述步骤4)修正结构稳定性分析模型的过程为：利用几何非线性假设，将薄壁结构节点位移场的影响施加至结构刚度矩阵中；基于更改后的结构刚度矩阵，建立稳定性分析模型。
[0025]所述步骤6)中，利用模态叠加法对薄壁结构稳定性进行求解。
[0026]本专利技术与现有技术相比的优点在于：
[0027]1、薄壁结构的负刚度变形导致舱体稳定特性发生大幅改变，传统基于线性化的稳定性计算方法无法捕获负刚度变形，因此不再适用，本专利技术基于几何非线性假设对薄壁结构负刚度引起的稳定性问题进行准确预示。
[0028]2、本专利技术通过建立等效静力过载的方法，将结构动力学问题等效为静力学问题，准确识别负刚度现象的产生并大幅降低问题求解难度。
[0029]3、本专利技术将线性分析模型与非线性分析模型相结合，通过预分析的方式，快速建立结构动力学模型与静力学模型之间的关系，从而提高分析效率。
附图说明
[0030]图1为本专利技术方法流程图。
[0031]图2为负刚度识别方法示意图。
[0032]图3为弧长与加载载荷比例关系图。
具体实施方式
[0033]如图1所示，本专利技术一种薄壁复合材料舱体负刚度稳定性分析方法，步骤如下：
[0034]1、建立薄壁舱体有限元分析模型，采用模态叠加法的稳态动力学分析方法，对薄壁舱体有限元模型进行基于频域方法的正弦扫频计算分析，获得正弦扫频分析结果，重点获取结构整体最大变形分布场U；
[0035]2、基于小变形假设(即)，建立薄壁舱段稳态静力学有限元分析模型，对模型施加等效过载，使薄壁结构变形分布满足如下条件：minσ|U
g
‑
U|，即等效过载下的结构各有限元节点变形与正弦扫频分析获得的有限元节点变形相减后的方差值最小，以此确认等效过载量值；
[0036]3、基于几何非线性假设(即考虑应变二阶微分项)，重新建立薄壁舱段稳态静力学有限元分析模型，加载确定后的等效过载，提取稳态静力结果中考虑应变二阶项的薄壁结构节点更新位移场，通过变形分布识别是否出现负刚度变形，若出现则继续以下步骤，若未出现则转至步骤6。
[0037]如图2所示，负刚度变形判定方法如下：
[0038]31)以舱体经线为基准，按指定间距逐条提取经线上节点变形信息，获取经线上变形后的位置结果；
[0039]32)利用舱体旋成体构型特点，建立经线局部坐标系；
[0040]33)在局部坐标系下，以提取的经线两端节点变形后位置构造线性辅助线，并形成函数：f x,y＝ax+by+c；
[0041]34)通过经线上的各节点位置信息x0,y0到辅助线的矢量距离判断结构是否发生负刚度变形，若始终d＞0，则结构未发生负刚度变形，若存在d＜0，则结构发生负刚度变本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种薄壁复合材料舱体负刚度稳定性分析方法，其特征在于包括：1)利用正弦扫频分析，获得结构变形分布U；2)建立等效静力分析模型，根据步骤1)中获得的结构变形分布场U，获得等效静力过载；3)根据步骤2)获得的等效静力过载，结合非线性静力分析模型，判断结构是否出现负刚度；若出现则继续步骤4)，若未出现则转至步骤6)；4)提取非线性静力分析模型变形结果，修正结构稳定性分析模型；5)对修正后的稳定性分析模型进行求解，获得薄壁结构发生失稳时的载荷比值和后屈曲的破坏状态，本方法结束；6)建立小变形假设下的有限元分析模型，并对薄壁结构稳定性进行求解。2.根据权利要求1所述的一种薄壁复合材料舱体负刚度稳定性分析方法，其特征在于：所述步骤1)获得结构变形分布的过程为：建立薄壁舱体有限元分析模型，对薄壁舱体有限元模型进行基于频域方法的正弦扫频计算分析，获得正弦扫频分析结果，并重点获取结构整体最大变形分布场U。3.根据权利要求2所述的一种薄壁复合材料舱体负刚度稳定性分析方法，其特征在于：采用模态叠加法的稳态动力学分析方法，对薄壁舱体有限元模型进行基于频域方法的正弦扫频计算分析。4.根据权利要求1所述的一种薄壁复合材料舱体负刚度稳定性分析方法，其特征在于：所述步骤2)获得等效静力过载的过程为：基于小变形假设，建立薄壁舱段稳态静力学有限元分析模型，对模型施加等效过载，使薄壁结构变形分布满足差值条件，以此确认等效过载量值。5.根据权利要求4所述的一种薄壁复合材料舱体负刚度稳定性分析方法，其特征在于：所述差值条件为：minσ|U
g
‑
U|，即等效过载下的结构各有限元节点变形与正弦扫频分析获得的有限元节点变形相减后的方差值最小，以此确认等效过载量值；其中U
g
为静力学分析模型中的结构位移...

【专利技术属性】
技术研发人员：王一凡，孙晓玲，肖振，张宏宇，崔慧永，张皓，王捷冰，周启超，高鹏，涂积奇，李萌萌，
申请(专利权)人：北京临近空间飞行器系统工程研究所，
类型：发明
国别省市：