CN110096772B - 一种面向航空航天薄壳结构的形位偏差特征库建立方法 - Google Patents

Abstract

本发明属于航空航天结构主承力构件设计技术领域，涉及一种面向航空航天薄壳结构的形位偏差特征库建立方法。包括以下步骤：1)给出适合现场工人测量的壳体质量检测卡，采点测量并填写壳体质量检测卡，获取壳体关键结构参数和形位偏差特征数据；2)对第一步得到的形位偏差数据，进行形位偏差的特征梳理和数学描述，进行缺陷成分分析；3)使用第一步和第二步对于多个典型航空航天薄壁结构的形位偏差数据进行收集和分析，建立航空航天薄壳结构形位偏差特征库。本发明专利有助于建立拥有过我国自主知识产权的航空航天薄壳结构形位偏差特征库，能够有效服务于我国航空航天装备的研制，缩短设计周期，为主要承力薄壳结构提供指导和设计规范。

Description

一种面向航空航天薄壳结构的形位偏差特征库建立方法

技术领域

本发明属于航空航天结构主承力构件设计技术领域，涉及一种面向航空航天薄壳结构的形位偏差特征库建立方法。

背景技术

薄壳结构具有高比刚度和高比强度，因此被广泛应用于航空航天结构的主承力构件中。航空航天中主承力薄壳结构主要承受轴压载荷，破坏的主要模式为屈曲失稳，因此薄壳结构在轴压工况下的极限承载能力是航空航天结构设计过程中的关键考核指标。然而通过大量轴压圆柱壳试验发现，试验结果和经典理论存在巨大的差异，试验得到的屈曲载荷往往只有经典理论解得20％-60％，并且实验结果存在巨大的分散和差异性。经过数十年的研究发现，初始几何缺陷的存在是试验结果和理论结果存在巨大差异的主要原因之一。而几何缺陷是加工制造、运输装配、服役等过程中不可避免的，因此如何精确地考虑不同种类的形位偏差(加工制造几何缺陷)是准确预估薄壳结构极限承载力并进行精细化设计的有效途径之一。欧盟采用欧洲钢壳规范的方法，根据初始缺陷的情况对已加工薄壳结构进行评估分级，该方法被认为是圆柱壳稳定设计方面的权威性规范。美国NASA使用实验设备测量圆柱壳蒙皮的完整初始几何缺陷并建立初始几何缺陷数据库，为航天薄壳结构设计起到指导作用。我国运载火箭主承力结构采用网格加筋壳结构形式，建立能够充分考虑加筋薄壳结构中不同种类形位偏差的数值模型，是准确预测薄壁筒壳结构实际承载力的有效途径，也是我国运载火箭简体筒壳结构轻量化设计的基础，建立拥有我国自主知识产权的航空航天薄壳结构形位偏差特征库，更是我国重型运载火箭和大型飞机制造的坚实基础。

航天航空主承力薄壳结构一般包括筒壳结构，锥壳结构，球壳结构和复杂异性曲面薄壳等结构形式。随着有限元等数值分析方法的快速发展，充分考虑形位偏差，建立能够体现不同种类和位置制造几何缺陷的有限元数值模型，进而精确分析计算筒壳结构极限承载能力是航天主承力结构精细化分析与设计的必要途径之一。针对网格加筋薄壳结构，主要的形位偏差包括壳体外表面母线偏移值、端面平整度、蒙皮径向凹陷、焊缝宽度和高度、筋条厚度不均匀度等制造几何缺陷形式。工程中一般采用较为保守的折减因子方法预估薄壳结构的极限承载力，这样限制了航空航天主承力结构的精细化和轻量化设计。因此，亟需建立一种能充分考虑不同种类形位偏差，实现精细化模型的建立方法，并开展一种面向航空航天薄壳结构的形位偏差特征库建立方法，为我国未来重型运载火箭和大型飞机制造提供指导和规范。

发明内容

本发明主要解决现有航空航天薄壳结构数值分析模型难以充分考虑不同类型形位偏差，难以精细化建模并针对形位偏差进行有效收集和分析的问题。提出一种面向航空航天薄壳结构的形位偏差特征库建立方法，通过填写壳体质量检测卡完成制造特征形位偏差的信息汇总和梳理，在完美有限元网格模型的基础上，采用插值和概率估计的方法对有限元模型中节点坐标进行修调，进而获得考虑不同种类形位偏差的精细化有限元模型，达到精细化建模并有效梳理制造特征，指导优化设计的目的。这种方法方便使用，可实现多种形位偏差的特征梳理并在数值模型中有效体现，并且可以完成航空航天薄壳结构形位偏差特征库的建立，为后续多种型号设计提供规范和设计基础，缩短设计周期。

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

一种面向航空航天薄壳结构的形位偏差特征库建立方法，包括以下步骤：

第一步，为搭建航空航天薄壳结构形位偏差特征库，为航空航天主承力薄壳结构的精细化建模提供数据基础，建立航空航天薄壳结构的设计规范，设计给出适合现场工人测量的壳体质量检测卡，采点测量并填写壳体质量检测卡，获取壳体关键结构参数和形位偏差特征数据。具体为：

针对网格加筋筒壳结构，所述的壳体关键结构参数包括壳体直径、高度、蒙皮厚度，筋条高度及厚度、加紧类型、焊接类型等结构参数和描述。所述的形位偏差特征数据包括壳体外表面母线偏移值、端面平整度、蒙皮径向凹陷、焊缝宽度和高度、筋条厚度不均匀度等形式，各种形位偏差数据测量方法如下：(1)壳体外表面母线偏移值：沿环向每隔一定角度在壳体外表面取一条测量线，测量线尽量避开纵向焊缝，沿测量线长度方向等分形成若干个测量点，利用靠尺测量并记录其径向偏移值，其中内凹为正，外凸为负。(2)端面平整度：沿环向每隔一定角度，利用塞尺测量下端面与地面间的距离Z_下，翻转壳体，按照同样环向位置，利用塞尺测量上端面与地面间的距离Z_上。(3)蒙皮径向凹陷：通过目测，将蒙皮凹陷由大到小依次排序，测量并记录每个凹陷的位置(轴向坐标和环向角度)、幅值(凹陷深度)。其中内凹为正，外凸为负。(4)焊缝宽度和高度：对每个焊缝，沿轴向或环向等分若干个点，作为测量点，测量并记录每个测量点处的焊缝宽度和高度。(5)筋条厚度不均匀度：沿环向每隔一定角度取一条纵向筋条(或有纵向分量)，沿纵向筋条的长度方向等分形成若干个测量点，利用游标卡尺测量并记录其厚度，并统计其不均匀度。需要说明，这些形位偏差数据类型和测量方法同样适用于锥壳、球壳、异性薄壳等其他薄壁结构形式。

通过对单一构件形位偏差数据的统计分析，可以用于实际结构的精细化数值分析，而数据的积累可实现航空航天薄壳结构形位偏差特征库的搭建，服务于我国火箭主承力构件精细化分析和设计，为火箭的进一步减重起到重要作用。所提出的壳体质量检测卡模板起到标准示范作用。

第二步，对第一步得到的形位偏差数据，进行形位偏差的特征梳理和数学描述，进行缺陷成分分析。具体为：

采用数据拟合、概率分布、数据插值等数学方法对壳体质量检测卡中收集的形位偏差特征(制造几何特征)数据进行分析，针对不同类型的形位偏差数据分别建立合适的数学模型，完成制造几何特征的特征梳理和描述。具体为：对于壳体外表面母线偏移值和端面平整度采用数据插值方法或数据拟合方法进行数学建模，所述的数据插值和数据拟合方法包括双线性、最临近、三次样条、多项式、傅里叶等不同形式的插值和函数拟合方法。对于蒙皮径向凹陷，采用凹陷形式的数学描述函数。对于焊缝宽度和高度、筋条厚度不均匀度采用数学概率统计的方法对尺寸的离散程度进行数学描述。

进一步的，对于壳体外表面母线偏移值和端面平整度，最有效的方法是使用傅里叶级数进行数据拟合，将形位偏差在空间频域范围内展开半波余弦或半波正弦的双重傅里叶级数是一种高效准确的三维形貌表征方法。所述的傅里叶半波余弦和半波正弦的双重傅里叶级数分别表示如下：

其中，

为无量纲离面位移，x为形位偏差环向坐标，y为形位偏差轴向坐标；L和R为柱壳的名义高度和名义半径；A_kl、B_kl为半波余弦傅里叶级数的系数；C_kl、D_kl为半波正弦傅里叶级数的系数；k和l为傅里叶级数系数的编号，n1和n2为傅里叶系数的个数。将式(1)转化为广义多变元线性拟合问题，然后使用最小二乘法进行系数的求解。

第三步，采用第一步和第二步的方法对于多个典型航空航天薄壁结构的形位偏差数据进行收集和分析，建立航空航天薄壳结构形位偏差特征库。具体为：

采用第一步获得的壳体质量检测卡对多组典型薄壳结构进行形位偏差数据收集，采用第二步的方法进行形位偏差的特征梳理和数学描述，进行缺陷成分分析。基于薄壳结构的关键几何参数和形位偏差特征进行数据分组归类、无量纲处理、数据存储，完成航空航天薄壳结构形位偏差特征库的建立。基于航空航天薄壳结构形位偏差特征库，可进行薄壳结构的质量评级和缺陷敏感性分析，从而建立不同形式薄壳结构的设计规范，为我国未来重型运载火箭和大飞机制造打下坚实的基础。

第四步，基于第三步建立的航空航天薄壳结构形位偏差特征库，建立充分考虑形位偏差的有限元数值模型。具体为：

首先，采用有限元分析软件ANSYS、ABAQUS或MSC.Patran，针对特定的薄壳结构建立完美的三维有限元模型，该模型为不含有制造几何缺陷的完美三维有限元模型。然后，在航空航天薄壳结构形位偏差特征库中搜索与指定模型具有相同结构形式和相近几何尺寸的形位偏差特征数据；基于该形位偏差的数学特征描述，对形位偏差缺陷进行合理的放缩以适应指定的结构尺寸。最后，采用插值方法对完美结构三维有限元模型整体或局部的节点坐标进行修调或者基于几何参数的概率分布函数为焊缝和筋条选取截面几何尺寸，达到将不同真实制造几何缺陷引入完美结构三维有限元模型的目的，进而获得充分考虑制造几何缺陷的精细化有限元模型。通过对缺陷幅值放大的方法可实现形位偏差数据的可视化，为设计人员提供参考。

另外，第一步中所述的适合现场工人测量的壳体质量检测卡，该模板同样适用于锥壳、球壳等结构形式，该模板可服务于搭建不同结构形式的形位偏差特征库，为建立精细化、真实化数值模型和建立行业规范起到重要作用。

本发明的有益效果为：本发明针对现有的航空航天薄壳结构数值分析模型难以考虑加工制造过程中产生的不同类型形位偏差问题，为搭建航空航天薄壳结构的形位偏差特征库，建立精细化薄壳结构数值模型并进行精细化数值分析，给出了适合现场工人测量使用的壳体质量检测卡，搭建了航空航天薄壳结构的形位偏差特征库，并给出了建立考虑多种形位偏差的精细化有限元模型的方法。可以通过逐点测量的方式对筒壳主要的制造几何缺陷特征信息进行采集，包括：壳体外表面母线偏移值、端面平整度、蒙皮径向凹陷、焊缝宽度和高度、筋条厚度不均匀度。通过对壳体质量检测卡所测量制造几何缺陷的特征梳理和数学描述，建立航空航天薄壳结构的形位偏差特征库，可以通过移动有限元模型节点的方式对完美有限元模型进行修调，从而获得考虑多种制造几何缺陷的精细化有限元数值模型。本发明所建立的航空航天薄壳结构的形位偏差特征库能够充分利用实验构件所体现的真实形位偏差信息，建立考虑不同种类制造几何缺陷的精细化数值分析模型，增强了有限元预测模型的真实性，进一步提高了有限元预测的准确性。建立拥有过我国自主知识产权的航空航天薄壳结构形位偏差特征库能够有效服务于我国航空航天装备的研制，缩短设计周期，为主要承力薄壳结构提供指导和设计规范。

附图说明

图1为本发明提供的面向航空航天薄壳结构的形位偏差特征库建立方法的实现流程图；

图2为典型形位偏差示意图；其中，图2(a)表示壳体外表面母线偏移，图2(b)表示端面平整度，图2(c)表示蒙皮径向凹陷，图2(d)表示焊缝宽度和厚度；

图3(a)为壳体外表面母线偏移形位偏差特征梳理和数学描述效果示意图；

图3(b)为端面平整度形位偏差特征梳理和数学描述效果示意图；

图3(c)为焊缝厚度形位偏差特征梳理和数学描述效果示意图；

图3(d)为筋条厚度形位偏差特征梳理和数学描述效果示意图；

图4为轴压筒壳结构的形位偏差放大可视化示意图。

具体实施方式

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

图1为本发明实例提供的一种面向航空航天薄壳结构的形位偏差特征库建立方法的实现流程图。如图1所示，本发明实施例提供的一种面向航空航天薄壳结构的形位偏差特征库建立方法包括：

第一步，给出适合现场工人测量的壳体质量检测卡，采点测量并填写壳体质量检测卡，获取壳体关键结构参数和形位偏差特征数据。采用逐点测量的方式，对直径2米左右的29个加筋筒壳的主要形位偏差数据进行测量统计，并填写壳体质量检测卡，获得网格加筋筒壳结构试验件的关键几何参数和形位偏差数据。测量的形位偏差主要包括：壳体外表面母线偏移值、端面平整度、蒙皮径向凹陷、焊缝宽度和高度、筋条厚度不均匀度。图2为测量的典型形位偏差示意图。这里采用的壳体质量检测卡模板如下表1(a)、表1(b)所示。

表1：壳体质量检测卡模板(a)

表1：壳体质量检测卡模板(b)

第二步，对第一步得到的形位偏差数据，进行形位偏差的特征梳理和数学描述，进行缺陷成分分析。使用数据拟合，概率分布，数据插值等数学方法对壳体质量检测卡中收集的形位偏差特征(制造几何特征)数据进行分析，针对不同类型的形位偏差数据分别建立合适的数学模型，完成制造几何特征的特征梳理和描述。对于本发明实例，壳体外表面母线偏移值使用样条曲线插值；端面平整度使用傅里叶函数拟合；焊缝宽度和高度和筋条厚度不均匀度使用数学概率统计的方法对尺寸的离散程度进行数学描述。图3为缺陷特征梳理和数学描述的效果示意图。

第三步，使用第一步和第二步对于多个典型航空航天薄壁结构的形位偏差数据进行收集和分析，建立航空航天薄壳结构形位偏差特征库。

使用第一步获得的壳体质量检测卡对多组典型薄壳结构进行形位偏差数据收集，使用第二步进行形位偏差的特征梳理和数学描述，进行缺陷成分分析。基于薄壳结构的关键几何参数和形位偏差特征进行数据分组归类、无量纲处理、数据存储，完成航空航天薄壳结构形位偏差特征库的建立。

第四步，基于航空航天薄壳结构形位偏差特征库，建立充分考虑形位偏差的有限元数值模型。首先，采用有限元分析软件ABAQUS针对特定的薄壳结构建立完美的三维有限元模型，该模型为不含有制造几何缺陷的完美三维有限元模型。然后，在航空航天薄壳结构形位偏差特征库中搜索与指定模型具有相同结构形式和相近几何尺寸的形位偏差特征数据；基于该形位偏差的数学特征描述，对形位偏差缺陷进行合理的放缩以适应指定的结构尺寸。最后，采用插值方法对完美结构三维有限元模型整体或局部的节点坐标进行修调或者基于几何参数的概率分布函数为焊缝和筋条选取截面几何尺寸，达到将不同真实制造几何缺陷引入完美结构三维有限元模型的目的，进而获得充分考虑制造几何缺陷的精细化有限元模型。对于端面平整度，蒙皮径向凹陷可通过准确的数学函数计算节点坐标的偏移量。基于壳体外表面母线偏移值数据插值获得的空间散点数据很难与有限元网格逐一对应，这里采用基于网格节点逐点插值的方法逐一修调蒙皮节点坐标。通过对缺陷幅值放大的方法可实现，形位偏差数据的可视化(如图4，形位偏差缺陷数据方法200倍)，为设计人员提供参考。

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

Claims (5)

1.一种面向航空航天薄壳结构的形位偏差特征库建立方法，其特征在于，所述的形位偏差特征库建立方法包括以下步骤：

第一步，设计给出适合现场工人测量的壳体质量检测卡，采点测量并填写壳体质量检测卡，获取壳体关键结构参数和形位偏差特征数据，形位偏差特征数据包括壳体外表面母线偏移值、端面平整度、蒙皮径向凹陷、焊缝宽度和高度、筋条厚度不均匀度；

第二步，对第一步得到的形位偏差数据，进行形位偏差的特征梳理和数学描述，进行缺陷成分分析；

采用数学方法对壳体质量检测卡中收集的形位偏差特征数据进行分析，针对不同类型的形位偏差数据分别建立合适的数学模型，完成形位偏差特征的特征梳理和描述：对于壳体外表面母线偏移值和端面平整度采用数据插值方法或数据拟合方法进行数学建模；对于蒙皮径向凹陷，采用凹陷形式的数学描述函数；对于焊缝宽度和高度、筋条厚度不均匀度采用数学概率统计的方法对尺寸的离散程度进行数学描述；

第三步，采用第一步和第二步的方法对于多个典型航空航天薄壁结构的形位偏差数据进行收集和分析，建立航空航天薄壳结构形位偏差特征库；

采用第一步获得的壳体质量检测卡对多组典型薄壳结构进行形位偏差数据收集，采用第二步的方法进行形位偏差的特征梳理和数学描述，进行缺陷成分分析；基于薄壳结构的关键几何参数和形位偏差特征进行数据分组归类、无量纲处理、数据存储，完成航空航天薄壳结构形位偏差特征库的建立；基于航空航天薄壳结构形位偏差特征库，进行薄壳结构的质量评级和缺陷敏感性分析，建立不同形式薄壳结构的设计规范；

第四步，基于第三步建立的航空航天薄壳结构形位偏差特征库，建立充分考虑形位偏差的有限元数值模型；

首先，针对特定的薄壳结构建立完美的三维有限元模型，该模型为不含有制造几何缺陷的完美三维有限元模型；然后，在航空航天薄壳结构形位偏差特征库中搜索与指定模型具有相同结构形式和相近几何尺寸的形位偏差特征数据；基于该形位偏差的数学特征描述，对形位偏差缺陷进行合理的放缩以适应指定的结构尺寸；最后，采用插值方法对完美结构三维有限元模型整体或局部的节点坐标进行修调或者基于几何参数的概率分布函数为焊缝和筋条选取截面几何尺寸，达到将不同真实制造几何缺陷引入完美结构三维有限元模型的目的，进而获得充分考虑制造几何缺陷的精细化有限元模型；通过对缺陷幅值放大的方法可实现形位偏差数据的可视化，为设计人员提供参考。

2.根据权利要求1所述的一种面向航空航天薄壳结构的形位偏差特征库建立方法，其特征在于，第一步中所述的适合现场工人测量的壳体质量检测卡，形位偏差数据类型和测量方法同样适用于锥壳、球壳、异性薄壳或其他薄壁结构形式，用于搭建不同结构形式的形位偏差特征库，为建立精细化、真实化数值模型和建立行业规范起到重要作用。

3.根据权利要求1或2所述的一种面向航空航天薄壳结构的形位偏差特征库建立方法，其特征在于，第一步中形位偏差数据测量方法如下：

(1)壳体外表面母线偏移值：沿环向每隔一定角度在壳体外表面取一条测量线，测量线避开纵向焊缝，沿测量线长度方向等分形成若干个测量点，利用靠尺测量并记录其径向偏移值，其中内凹为正，外凸为负；

(2)端面平整度：沿环向每隔一定角度，利用塞尺测量下端面与地面间的距离Z_下，翻转壳体，按照同样环向位置，利用塞尺测量上端面与地面间的距离Z_上；

(3)蒙皮径向凹陷：通过目测，将蒙皮凹陷由大到小依次排序，测量并记录每个凹陷的位置、幅值；其中内凹为正，外凸为负；

(4)焊缝宽度和高度：对每个焊缝，沿轴向或环向等分若干个点，作为测量点，测量并记录每个测量点处的焊缝宽度和高度；

(5)筋条厚度不均匀度：沿环向每隔一定角度取一条纵向筋条，沿纵向筋条的长度方向等分形成若干个测量点，利用游标卡尺测量并记录其厚度，并统计其不均匀度。

4.根据权利要求1或2所述的一种面向航空航天薄壳结构的形位偏差特征库建立方法，其特征在于，第二步中，对于壳体外表面母线偏移值和端面平整度，优选采用傅里叶级数进行数据拟合，将形位偏差在空间频域范围内展开半波余弦或半波正弦的双重傅里叶级数是一种高效准确的三维形貌表征方法；所述的傅里叶半波余弦和半波正弦的双重傅里叶级数分别表示如下：

其中，

为无量纲离面位移，x为形位偏差环向坐标，y为形位偏差轴向坐标；L和R为柱壳的名义高度和名义半径；A_kl、B_kl为半波余弦傅里叶级数的系数；C_kl、D_kl为半波正弦傅里叶级数的系数；k和l为傅里叶级数系数的编号，n1和n2为傅里叶系数的个数；

将公式(1)转化为广义多变元线性拟合问题，采用最小二乘法求解系数。

5.根据权利要求3所述的一种面向航空航天薄壳结构的形位偏差特征库建立方法，其特征在于，第二步中，对于壳体外表面母线偏移值和端面平整度，优选采用傅里叶级数进行数据拟合，将形位偏差在空间频域范围内展开半波余弦或半波正弦的双重傅里叶级数是一种高效准确的三维形貌表征方法；所述的傅里叶半波余弦和半波正弦的双重傅里叶级数分别表示如下：

其中，

为无量纲离面位移，x为形位偏差环向坐标，y为形位偏差轴向坐标；L和R为柱壳的名义高度和名义半径；A_kl、B_kl为半波余弦傅里叶级数的系数；C_kl、D_kl为半波正弦傅里叶级数的系数；k和l为傅里叶级数系数的编号，n1和n2为傅里叶系数的个数；

将公式(1)转化为广义多变元线性拟合问题，采用最小二乘法求解系数。

Images