一种离散的飞机气动载荷参数化的方法技术

本发明专利技术涉及飞机结构件的强度计算领域，尤其涉及一种将离散的飞机气动载荷参数化并用于有限元结构分析的一种方法。本发明专利技术采用专业的数值工具和有限元软件对离散气动载荷压力系数进行参数化处理，使其以节点坐标对应的函数表示出来。通过该函数计算出来的气动载荷压力系数和原始的气动载荷压力系数之间的相对误差不超过10％，可以作为一种方便快捷的工程方法应用在各种飞机结构件的力学计算中。本发明专利技术中的载荷回归函数是坐标的函数，不会因为待分析结构网格节点的改变而发生变化。当有限元模重新划分网格模型时，同种载荷工况下的气动载荷压力系数回归函数可以重复使用，大大减少了有限元模型调试时的工作量,极大的提高了工作效率。

【技术实现步骤摘要】
一种离散的飞机气动载荷参数化的方法
本专利技术涉及飞机结构件的强度计算领域，尤其涉及一种将离散的飞机气动载荷参数化并用于有限元结构分析的一种方法。
技术介绍
气动载荷是飞机机头、机身、机翼等航空器部件所承受的主要载荷。在进行结构件或者组件静强度分析时，由于其结构较为复杂，通常会采用相应的有限元软件进行分析计算，例如ABAQUS、ANSYS、Nastran等。气动载荷一般会以节点压力的形式给出，例如飞机结构表面某点的压力系数。对于某个结构件或者组件而言，就是一系列带节点坐标的离散载荷形式。以往在利用有限元软件分析过程中，因为要对结构件重新建立三维模型，这时的结构网格节点很难和离散气动载荷的网格节点一一对应，一般的处理方法是选用一个最严苛的压力施加到结构件上，此计算结构得到的结果偏于保守；尤其对于类似直升机有多种飞行姿态对应多种气动载荷的情况下，并不能通过这种计算方法得到最为严苛的飞行工况。因此，需要引入一种新的更为方便、快捷的气动载荷施加方法得到较为精确的强度计算结果。
技术实现思路
本专利技术的目的在于提供一种飞机载荷的参数化方法来克服或至少减轻现有技术的至少一个上述缺陷。本专利技术的技术解决方案是，所述飞机离散气动载荷的处理方法包括如下步骤：步骤1：根据已有的飞机载荷信息文件，提取待分析结构的节点坐标范围内的离散气动载荷压力系数以及各离散气动载荷压力系数对应的节点坐标值，形成离散气动载荷压力系数信息文件；步骤2：对离散气动载荷压力系数信息文件中节点坐标范围内的数据进行分析，找出离散气动载荷压力系数和各节点坐标值之间的线性和非线性分布规律，设定自变量和因变量，其中节点坐标值为自变量，节点处的气动载荷压力系数为因变量，形成一个多变量函数，赋予待定系数初始值，初始值大于零或者小于零，通过数值分析工具对离散气动载荷压力系数进行数值回归，得到气动载荷压力系数关于坐标值的回归函数，通过节点坐标值得到气动载荷压力系数的函数值；步骤3：计算气动载荷压力系数的函数值和离散气动载荷信息文件中的气动载荷压力系数的相对误差，使其在10％的范围内，若相对误差大于10％，对函数进行修正，重新进行数值回归，确定针对待分析结构的气动载荷压力系数回归函数；使离散气动载荷参数化，将参数化后的离散气动载荷压力系数施加在有限元模型上，进行力学性能计算。所述步骤1采用有限元软件提取待分析结构的节点坐标以及对应的气动载荷压力系数，形成气动载荷压力系数信息文件。所述步骤2以节点坐标值为自变量、气动载荷压力系数为因变量的函数为多项式，每项前有一个待定系数。所述数值分析工具为MATLAB，所述有限元软件为ABAQUS。本专利技术的有益效果是：1、本专利技术中的离散飞机气动载荷参数化方法采用专业的数值工具和有限元软件对离散气动载荷压力系数进行参数化处理，使其以节点坐标对应的函数表示出来。通过该函数计算出来的气动载荷压力系数和原始的气动载荷压力系数之间的相对误差不超过10％，可以作为一种方便快捷的工程方法应用在各种飞机结构件的力学计算中。2、本专利技术中的载荷回归函数是坐标的函数，不会因为待分析结构网格节点的改变而发生变化。当有限元模重新划分网格模型时，同种载荷工况下的气动载荷压力系数回归函数可以重复使用，大大减少了有限元模型调试时的工作量。3、本专利技术中的主要工作步骤由MATLAB程序实现，不同的气动载荷压力系数的参数化工作用同样的程序文件实现。极大的提高了工作效率。附图说明图1是根据本专利技术一实施例的飞机气动载荷处理方法的流程示意图。图2为直升机俯仰角α和侧滑角β的示意图。图3是直升机风挡蒙皮三维网格示意图。图4是直升机风挡原始气动载荷压力系数的分布图。图5是直升机风挡回归气动载荷压力系数在的分布图。具体实施方式下面结合附图和具体的实施实例对本专利技术技术方案作进一步的详细描述。所描述的实施例是本专利技术的一部分实施例，而不是全部的实施例。下面通过参考附图描述的实施例是示例性的，旨在用于解释本专利技术，而不能理解为对本专利技术的限制。计算待分析结构的力学性能前，需要得到用于表征气动载荷压力系数的函数。图1所示的飞机气动载荷处理方法旨在获得这种函数，其包括如下步骤：步骤1：已有的载荷信息文件包含多种飞机结构表面蒙皮的网格节点和节点对应的气动载荷压力系数。首先要从原始载荷信息文件中提取待分析结构的气动载荷信息和网格节点坐标。采用有限元软件输出蒙皮节点坐标以及气动载荷压力系数信息，以文件的形式输出，气动载荷信息包括各个节点坐标以及各个节点所对应的气动载荷压力系数。步骤2：对离散气动载荷压力系数和节点坐标值之间的对应关系进行数值分析，找出离散气动载荷压力系数和坐标值之间的线性和非线性分布规律，假设一个多变量函数，该函数以三个坐标(X、Y、Z)为自变量，坐标值对应的气动载荷压力系数为因变量，形成一个多变量函数，一般设置为多项式，每项前有一个待定系数。赋予待定系数一系列不等于零的初始值，利用已有的坐标值与气动载荷压力系数值对设定的未知函数进行回归，目的是确定待定系数。从而得到回归函数。函数回归工作在matlab等数值分析软件中进行。如图4为原始气动载荷压力系数的分布规律图，图5为回归函数计算出来的气动载荷压力系数的分布规律图。步骤3：计算气动载荷压力系数回归值和气动载荷压力系数原始值的相对误差，使平均相对误差不大于10％。当使其平均误差值在允许范围内，确定该回归函数。当相对误差超过允许范围，可对多项式函数进行修正，重复步骤2。直至得到满足误差允许条件的回归函数。确定了待定系数的回归函数，通过任意坐标值可以计算出该点的气动载荷压力系数值。有限元软件中以ABAQUS为例，可以施加以坐标为变量的载荷函数。前面得到的气动载荷压力系数函数，可以施加在有限元模型表面。对于本领域技术人员而言，本专利技术不仅限于上述示范性实例的细节，而且在不背离本专利技术的精神和基本特征的情况下，能够以其他的具体形式实现本专利技术中的
技术实现思路
，都应包含在本专利技术所要求保护的范围内。实施例：根据直升机的俯仰角α和侧滑角β不同，可以将直升机分为若干种飞行工况。图2所示为直升机俯仰角α和侧滑角β的示意图。下面选取其中一种飞行工况的气动载荷压力载荷系数参数化作为实施例进行说明。图3是某种直升机风挡蒙皮的3D网格图；图4则是离散的气动载荷压力系数在风挡蒙皮节点上的分布图。(1)提取风挡蒙皮的气动载荷压力系数信息文件，文件中包括蒙皮节点坐标以及各个节点所对应的气动载荷压力系数。(2)假设一个以x、y、z三个坐标为自变量的函数：设定待定系数初始值：a0＝1，a1＝1，a2＝1，a3＝1，a4＝1，a5＝1，a6＝1，a7＝1，a8＝1，a9＝1利用Matlab里的nlinfit回归函数确定系数a0～a9的值如下：a0＝0.2721，a1＝0.3118，a2＝-3.1148，a3＝14.8490，a4＝6.3261，a5＝-417.6138，a6＝0.0043，a7＝-0.0109，a8＝0.0223，a9＝133.2567。将上述a0～a9的值带入公式(1)中得到通过该公式(2)计算得到的气动载荷压力系数分布图如图5所示。可以看出其与原始气动载荷分布规律(图4)相同，但有一定偏差。(3)计算得出的气动载荷压力系数相对于原始气动载荷压力系数的相对误差平均本文档来自技高网...

【技术保护点】
1.一种飞机离散气动载荷的参数化方法，其特征在于，所述飞机离散气动载荷的处理方法包括如下步骤：步骤1：根据已有的飞机气动载荷信息文件，提取待分析结构的节点坐标范围内的离散气动载荷压力系数以及各离散气动载荷压力系数对应的节点坐标值，形成离散气动载荷压力系数信息文件；步骤2：对离散气动载荷压力系数信息文件中节点坐标范围内的数据进行分析，找出离散气动载荷压力系数和各节点坐标值之间的线性和非线性分布规律，设定自变量和因变量，其中坐标值为自变量，气动载荷压力系数为因变量，形成一个多变量函数，赋予待定系数初始值，初始值大于零或者小于零，通过数值分析工具对气动载荷压力系数进行数值回归，得到气动载荷压力系数关于坐标值的回归函数，通过节点坐标值得到气动载荷压力系数的函数值；步骤3：计算气动载荷压力系数的函数值和离散气动载荷信息文件中的气动载荷压力系数的相对误差，使其在10％的范围内，若相对误差大于10％，对函数进行修正，重新进行数值回归，确定针对待分析结构的气动载荷压力系数回归函数；使离散气动载荷参数化，将参数化后的离散气动载荷压力系数施加在有限元模型上，进行力学性能计算。

【技术特征摘要】
1.一种飞机离散气动载荷的参数化方法，其特征在于，所述飞机离散气动载荷的处理方法包括如下步骤：步骤1：根据已有的飞机气动载荷信息文件，提取待分析结构的节点坐标范围内的离散气动载荷压力系数以及各离散气动载荷压力系数对应的节点坐标值，形成离散气动载荷压力系数信息文件；步骤2：对离散气动载荷压力系数信息文件中节点坐标范围内的数据进行分析，找出离散气动载荷压力系数和各节点坐标值之间的线性和非线性分布规律，设定自变量和因变量，其中坐标值为自变量，气动载荷压力系数为因变量，形成一个多变量函数，赋予待定系数初始值，初始值大于零或者小于零，通过数值分析工具对气动载荷压力系数进行数值回归，得到气动载荷压力系数关于坐标值的回归函数，通过节点坐标值得到气动载荷压力系数的函数值；步骤3：计算气动载荷压力系数的函数值和离散气...

【专利技术属性】
技术研发人员：赵景云，颜悦，霍钟祺，姜良宝，
申请(专利权)人：中国航发北京航空材料研究院，
类型：发明
国别省市：北京,11