机翼快速变外形参数建模方法与流程

本发明涉及一种机翼快速变外形参数建模方法。
背景技术：
：对于一个初始机翼模型，其外形参数等还存在着一定的可变化空间，但伴随外形参数的修改，结构有限元模型及气动模型也需作出相应的更新，这样便面临着需要重新建模的问题。尤其在优化过程中，机翼外形参数需要进行多次反复设计修改，大量的设计时间就会消耗在修改模型上，对于较复杂的机翼三维板杆有限元模型，难度尤甚。有限元模型的建立通常包括几何模型建立、有限元网格划分，定义材料、约束、载荷等必要的过程。几何模型的建立可以在cad(如catia，pro-e等)软件中实现，也可以在cae(如msc.patran、ansys)软件中实现。虽然cad软件在几何建模方面在几何建模方面具有较高的建模效率，也提供了较为方便的参数化建模方法，但若进行有限元分析还需要将几何模型导入cae软件，并进行网格划分、定义材料、约束、载荷等工作。目前较为普遍的参数化建模方法都是基于cad软件实现模型变参重建，再通过cae软件前处理模块识别才能进入到后续的结构有限元分析。技术实现要素：根据本发明的一个方面，提供了一种机翼快速变外形参数建模方法，其特征在于包括：a)机翼外形参数作为决定机翼外形的参数，在变参前给定，变参后得到相应的机翼模型；初始机翼的关键点参数作为实现变参所需的辅助参数，在变参前根据初始模型给定；其中外形参数包括：内翼展弦比λ1，外翼展弦比λ2，内翼根梢比η1，外翼根梢比η2，内翼上反角τ1，外翼上反角τ2，后掠角χ，前梁在翼根位置的弦向占比m1，前梁在机翼转折处的弦向占比m2，前梁在翼尖位置的弦向占比m3，后梁在翼根位置的弦向占比n1，后梁在机翼转折处的弦向占比n2，后梁在翼尖位置的弦向占比n3，机翼转折处即内翼和外翼的相交处，此处亦为机翼的后缘连线发生转折的位置，内翼即机翼模型的平面投影在a1a2a3a4梯形区域内的部分，外翼即机翼模型的平面投影在a3a4a5a6梯形区域的部分，初始关键点参数包括：翼根前缘点a1，翼根后缘点a2，机翼转折处前缘点a3，机翼转折处后缘点a4，翼稍前缘点a5，翼稍后缘点a6，翼根处的前梁点位置b1，机翼转折处的前梁点位置b3，翼尖处的前梁点位置b5，翼根处的后梁点的位置b2，机翼转折处的后梁点的位置b4，翼尖处的的后梁点的位置b6，b1)计算关键点位置，包括：把机翼三维板杆模型的平面形状视为两个分别由a1a2a3a4和a3a4a5a6围成的梯形区域，使翼根处的前、后缘点a1、a2的位置固定不变，在改变机翼的十三个外形参数λ1、λ2、η1、η2、τ1、τ2、χ、m1、m2、m3、n1、n2、n3后，计算出其他关键点a3、a4、a5、a6的新位置，b2)计算初始模型任意一点g0与点g0在横截面的前后缘点p0,q0的坐标值，对初始模型所有点依次执行这一操作,b3)确定初始模型各点的比例系数，包括：通过对模型的nastran文件中的grid卡片(包含模型所有点的位置信息)的读取，读入初始模型内翼处任何一点g0的x,y,z轴的坐标值，确定g0点对应的展向比例系数k1，计算弦向比例系数k2和厚度方向比例系数k3,b4)计算变参后新模型任意一点g0′与点g0′在横截面的前后缘点p0′,q0′的坐标值，对变参后新模型所有点依次执行这一操作,b5)进行比例系数的逆向计算，确定g0点变参后的点g0′的位置，包括：模型任一点在变参前后的展向、弦向和厚度方向的比例系数不变，即g0点在初始机翼中的比例系数和g0′点在新机翼中的比例系数一致，从已知g0点展向、弦向和厚度方向的比例系数以及g0,p0,q0三点坐标得到g0′,p0′,q0′三点坐标后，对比例系数k1,k2.k3进行逆向计算，得出点g0变参后的g0′点坐标。附图说明图1是机翼有限元模型坐标系示意图。图2是机翼有限元模型简化后的平面投影与外形参数说明的示意图。图3是说明机翼有限元模型内部翼梁占比参数的示意图。图4是阐明机翼关键点及比例系数计算方法的示意图。图5是阐明翼盒内翼部分新角点计算方法的示意图。图6是根据本发明的一个实施例的机翼的变参方法的流程图。图7是图6中“单个梯形区域的变参”步骤的子流程图。图8是验证变参方法所使用的初始机翼。图9是变后掠角后模型2的示意图。图10是变内翼展弦比后模型3的示意图。图11是变外翼上反角后模型4的示意图。具体实施方式表1：对变参过程涉及的总体外形参数的说明在上述现有技术的参数化建模方法中，由于整个过程需要采用两种软件，会给分析带来不便，因此，本发明人提出了本发明的建模方法；本发明的该方法即基于这样的情况——在实际优化的过程中往往已经具备了初始模型。根据本发明的一个具体实施例，基于msc.patran软件平台，使用fortran语言编程修改初始模型点的坐标，建立了变化外形参数后的新模型。根据本发明的其他实施例适用于其他软件平台的机翼各类结构模型及气动模型，在改变初始模型内翼展弦比、外翼展弦比、内翼根梢比、外翼根梢比、内翼上反角、外翼上反角、后掠角、翼根位置翼尖位置和机翼转折位置这三个位置的前梁和后梁弦向占比(计六个参数)，总计十三个参数后，通过修改初始模型点的坐标，可快速新建对应的模型。机翼外形参数改变后，虽然模型发生了改变，但机翼的基本结构形式保持不变。所以变参只需改变点的位置即可，在msc.patran&nastran的软件中便体现为对于所有grid卡片坐标值的修改。所以，变参后重建模型的核心在于两点：一是建立初始模型关键点的位置与外形参数的函数关系；二是找到初始模型的所有点和新模型点之间的位置函数关系，生成新的grid卡片。这两点即为本方法的核心，前者在cad软件参数化建模中被广泛应用，后者则是本发明中首先提出和明确的。在介绍根据本发明的建模方法前，首先对其机翼变参过程中涉及到的参数进行说明。参数包括两类：机翼外形参数和初始关键点参数。外形参数有：内翼展弦比；外翼展弦比；内翼根梢比；外翼根梢比；内翼上反角；外翼上反角；后掠角；前梁和后梁分别在翼根位置、翼尖位置、机翼转折位置这三个位置的弦向占比(计六个参数)；总计十三个参数。各外形参数的介绍以初始模型为例进行说明，见表1。这十三项外形参数在变参前由使用者决定输入，或在优化过程中随机产生。初始关键点参数有：翼根前缘点a1；翼根后缘点a2；机翼转折位置前缘点a3；机翼转折位置后缘点a4；翼稍前缘点a5；翼稍后缘点a6；翼根处、翼尖处和机翼转折位置处的前梁点的位置(依次为b1、b3、b5)；翼根处、翼尖处和机翼转折位置处的后梁点的位置(依次为b2、b4、b6)；总计12个点的x,y,z轴位置参数。这些参数是初始模型的参数，a1-a6的位置如图2所示，b1-b6的位置如图3所示。在整个变参过程中a1、a2不发生改变。初始关键点位置如图2、图3所示。根据本发明的一个实施例的模型坐标系的选取如图1所示，机翼有限元模型外部轮廓在x_y的简化平面投影如图2所示。介绍完各参数后，以下详细描述根据本发明的机翼快速变外形参数建模方法，其包括：第一步，计算新的关键点位置翼根处前后缘点a1、a2固定，机翼外形参数改变后，其余的关键点位置随之改变，需计算出其他关键点的新位置。如图2所示，机翼可以从机翼转折处分为内翼和外翼。这样，机翼的平面外形可以看做两个连接着的梯形，分别对应机翼的内翼和外翼。例如，机翼内翼可以等效为如图4所示的梯形。就一个梯形区域而言，如图4所示，所给的输入参数相当于已知了这个梯形四个角点a1、a2、a3、a4的原始位置坐标，a1、a2两点不变，新角点a′3、a′4(分别对应a3、a4)位置及新展长l′1计算公式如下(带上标与不带上标分别对应变参后机翼和初始机翼)：l′1＝b0·λ1′其中，带上标的符号对应变参后机翼参数，不带上标的符号对应初始机翼参数，所有其余点坐标值的表示方式以此类推。l1′,λ1′,η1′,τ1′,χ′分别表示变参后机翼内翼展长、内翼展弦比、内翼根梢比、内翼上反角和后掠角；l1,λ1,η1,τ1,χ分别表示初始的机翼内翼展长、内翼展弦比、内翼根梢比、内翼上反角和后掠角；b0为初始机翼根部弦长，变参后不发生改变，分别为点a3′的x,y,z三个方向坐标值，分别点a4′的x,y,z三个方向坐标值，分别点a1′的x,y,z三个方向坐标值，即点a2′的z轴坐标值。第二步，如图4所示，计算初始模型任意一点g0与点g0在横截面的前后缘点p0,q0的坐标值。根据三角形的相似性，计算公式如下：其中，其中，分别为点p0,q0,g0,a1,a2的y轴坐标值，分别为点a1,a2,a3,a4,p0,q0的x轴坐标值，分别为点a1,a2,p0,q0的z轴坐标值，l1为初始机翼内翼展长，τ1为初始机翼内翼上反角。第三步，计算比例系数在经由第一步计算得到关键点位置和第二步计算得到p0,q0,g0三点坐标后，需要建立模型中任意一点与关键点位置之间的函数关系。如图4所示，方框中为机翼中间任意位置的横截面简图。就一个梯形区域而言，读入初始模型内翼处任何一点g0的x,y,z轴的坐标值，计算其展向比例系数为弦向比例系数和厚度方向比例系数其中，分别为点g0,p0,q0的x轴坐标，为点a1的y轴坐标，分别为点g0,p0的z轴坐标，l1为初始机翼内翼展长。第四步，计算得出变参后点p0′、q0′点(对应初始机翼的p0、q0点)的坐标，方法和第二步一致，计算公式如下：其中，其中，分别为点p0′,q0′,g0′,a1′,a2′,a1的y轴坐标值，分别为点a1,a2,a3′,a4′,p0′,q0′的x轴坐标值，分别为点a1,a2,p0′,q0′的z轴坐标值，l1′为变参后新机翼的内翼展长，k1为展向比例系数，τ1′为变参后新机翼的内翼上反角。第五步，通过比例系数的逆向计算，输出变参后点g0′的全部位置信息，计算公式如下：其中，分别为点g0′,p0′,q0′的x轴坐标，为点a1的y轴坐标，分别为点g0′,p0′的z轴坐标，l1′为变参后新机翼内翼展长。以上即根据本发明的一个实施例的具体变参建模方法。对于如图1整个机翼包括两段梯形区域的情况而言，先计算内翼，内翼计算完毕后，内翼稍处的a′3、a′4点作为外翼根处不变的点，同理计算出外翼稍新的关键点a′5、a′6位置。外翼部分计算方法完全相同。若机翼带小翼，则为三段梯形区域；同样的，若机翼为飞翼构型，可能会增加机翼转折处的位置，相应地增加梯形区域即可。以下为前后梁在翼根位置、翼尖位置、机翼转折位置的弦向占比这六个参数变参后点的位置计算方法。根据本发明的一个具体实施例，此步骤被作为变参的第一步，即先改变前后梁占比，得到整个模型新的点坐标文件，在此文件基础上再改变其他外形参数，得到最终的点坐标文件。如图3所示，b1b3b4b2、b3b5b6b4分别标出了两个梯形区域，这两部分组成机翼有限元模型的翼盒部分。给定六个占比参数后，投影在翼盒内的所有点的z,y轴坐标不变，相当于仅对翼盒部分作了x方向的拉伸收缩变换，各点x轴新坐标值的计算方法和前文相同。需要注意的是，对机翼外形要求较精细的气动模型一般不存在此部分步骤，改变前后梁占比一般仅针对结构有限元模型，外形要求不需要非常精细。另外，在实际应用中各占比参数的变化范围一般较小，例如前后梁新位置与原位置在弦向的距离，不超过此处弦长的5％(此比例根据使用的精度要求界定)，故仅对点的x轴坐标进行改变，z，y轴坐标不变。(如果各占比参数的变化范围较大或者有限元模型外形要求较精细，则不适用此方法进行前后梁弦向占比的变参。)以翼盒内翼部分平面投影后的梯形区域为例，与前文所述方法相比，区别在于第一步，梯形四个角点b1b3b4b2的坐标，均由新的占比参数m1′、m2′、n1′、n2′(定义见表1，加上标与初始机翼区分)给定。如图5所示，新的四个角点b1′b3′b4′b2′的x轴坐标值计算公式如下(在初始机翼的基础上改变前后梁占比)：接下来计算此区域内任一点新坐标，如图5，方框中为翼盒任意位置的横截面简图，r0、j0为上下两短线中点。点r0、j0、t0分别对应前文的p0、q0、g0点，相同的方法计算弦向比例系数k2，再逆向计算出变占比参数后的点t0的x轴坐标即可。整个机翼(以常规布局的机翼进行说明，即仅有内翼和外翼两段梯形区域的机翼)的变参流程如图6所示，其中“单个梯形区域的变参”步骤包含“单个梯形区域的变参”子流程。图6中“单个梯形区域的变参”这一步骤的子流程图如图7所示。根据本发明的变参建模方法的实例验证：图8为初始机翼有限元模型，为方便对比，模型1-3(如图9-11所示)仅改变一个外形参数，具体如表2所示：表2：模型1-3的变参详情模型1：机翼后掠角初始53°，变为65°。模型2：机翼内翼展弦比初始1.61，变为1.4。模型3：机翼外翼上反角初始1.1°，变为5.0°。图8-11中，所有的虚线部分勾勒出了初始机翼的外轮廓。图8为初始机翼，改变后掠角后如图9所示，从图中可以看出整个机翼后掠角(后掠角定义如图2中的χ角所示)按照给定的新参数发生变化。图10为给定新的内翼展弦比后的模型，虚线轮廓为初始机翼，新机翼与初始机翼相比内翼展向明显变短，整个机翼也较初始机翼展向变短。图11为给定新的外翼上反角后的模型，左侧的虚线轮廓为初始机翼，右侧的虚线显示了初始机翼外翼上反程度，从图中可以看出整个机翼外翼上反角按照给定的新参数发生了变化。本发明确立了在具有初始模型的前提下，无需接入cad软件，仅对初始模型文件的点的坐标进行更改的情况下，改变初始模型内翼展弦比、外翼展弦比、内翼根梢比、外翼根梢比、内翼上反角、外翼上反角、后掠角、翼根位置翼尖位置和机翼转折位置这三个位置的前梁和后梁弦向占比(计六个参数)，总计十三个参数后，得到新模型的一套快捷有效的方法。本发明消除了更改参数后重新建模困难或接入cad软件进行跨平台变参的繁琐操作，通过对整个机翼模型所有点的坐标进行伸缩变换，得到了与初始机翼构型相同的新机翼。本发明不仅可以用作迅速生成新机翼(新机翼与原始机翼构型需一致)，而且可以接入机翼优化的过程中，对初始机翼的外形参数进行优化，得到机翼的最佳外形。当前第1页12