CN112434457B - 一种气动载荷转换为结构载荷的方法 - Google Patents

一种气动载荷转换为结构载荷的方法 Download PDF

Info

Publication number
CN112434457B
CN112434457B CN202011416150.2A CN202011416150A CN112434457B CN 112434457 B CN112434457 B CN 112434457B CN 202011416150 A CN202011416150 A CN 202011416150A CN 112434457 B CN112434457 B CN 112434457B
Authority
CN
China
Prior art keywords
load
structural
data
pneumatic
node
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011416150.2A
Other languages
English (en)
Other versions
CN112434457A (zh
Inventor
易俊杰
洪清泉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Boke Industrial Co ltd
Original Assignee
Shanghai Boke Industrial Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Boke Industrial Co ltd filed Critical Shanghai Boke Industrial Co ltd
Priority to CN202011416150.2A priority Critical patent/CN112434457B/zh
Publication of CN112434457A publication Critical patent/CN112434457A/zh
Application granted granted Critical
Publication of CN112434457B publication Critical patent/CN112434457B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Fluid Mechanics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明涉及计算机仿真技术领域，尤其涉及一种气动载荷转换为结构载荷的方法包括：气动面载荷预处理，形成气动面Cp数据；结构数据预处理，形成结构节点空间位置数据；根据形成的气动面Cp数据和结构节点空间位置数据，插值计算得到结构节点Cp数据；获得结构节点Cp数据后，在单元面上积分形成单元载荷数据，对单元载荷数据在结构节点上进行分载计算，得到结构节点力载荷本发明通过读取气动载荷数据，继而形成压力场，通过读取结构有限元数据，并在形成的压力场内进行插值，形成单元面压力载荷，并最终转化为节点力载荷。本发明提高了气动载荷处理的效率和准确性。

Description

一种气动载荷转换为结构载荷的方法
技术领域
本发明涉及计算机仿真技术领域，尤其涉及一种气动载荷转换为结构载荷的方法。
背景技术
飞行器在实际使用过程中承受着各种各样的载荷，大体可以分为两大类：一类为空气动力载荷，另一类为惯性力载荷。高速运行飞行器受到气动载荷的影响比较大，在飞机结构有限元分析中，空气动力载荷最终在结构设计中往往需要转化为结构上的受力，也即通常需要将气动载荷转换为结构有限元的节点载荷，基于该种需求，提供一种将气动载荷转换为结构节点载荷的方法具有非常重要的工程意义。
发明内容
针对上述现有技术的缺点，本发明的目的旨在提供一种高效、快速、准确的气动载荷转换为结构载荷的方法。
本发明实施例提供的一种气动载荷转换为结构载荷的方法，该方法包括：
气动面载荷预处理，形成气动面Cp数据；结构数据预处理，形成结构节点空间位置数据；
根据形成的气动面Cp数据和结构节点空间位置数据，插值计算得到结构节点Cp数据；
获得结构节点Cp数据后，在单元面上积分形成单元载荷数据，对单元载荷数据在结构节点上进行分载计算，得到结构节点力载荷。
优选地，上述气动载荷转换为结构载荷的方法中，所述气动面类型包括扁平面和曲面；
所述气动面为扁平面时，气动面载荷预处理流程包括：
读取输入的气动面轮廓数据，形成分载网格面；
读取气动分布载荷Cp，根据双线性插值方法，形成分载网格面Cp数据；
所述气动面为曲面时，气动面载荷预处理流程包括：
读取预存的气动面节点Cp数据。
优选地，上述气动载荷转换为结构载荷的方法中，所述根据双线性插值方法，形成分载网格面Cp数据的流程如下：
假设分载网格面的四个角点的坐标和Cp数据分别为：(x1，y1，cp1)、(x2，y2，cp2)、(x3，y3，cp3)和(x4，y4，cp4)；
由四个角点插值生成双线性曲面，则该面内任一点的参数坐标和Cp数据为：
x(u，w)＝x1+u(x2-x1)+[x3+u(x4-x3)-x1-u(x2-x1)]w
y(u，w)＝y1+u(y2-y1)+[y3+u(y4-y3)-y1-u(y2-y1)]w
cp(u，w)＝cp1(1-u)(1-w)+cp3(1-u)w+cp2u(1-w)+cp4uw
式中：u∈[0，1]；w∈[0，1]，u、w均为参数坐标。
优选地，上述气动载荷转换为结构载荷的方法中，所述插值计算得到结构节点Cp数据的具体流程包括：
4.1)根据结构节点的空间位置数据，确定在设定范围内与所有结构节点关联的气动Cp节点；
4.2)根据步骤4.1)确定的气动Cp节点的空间位置数据和Cp数据，通过样条基函数法或径向基函数法，计算得到局部气动Cp数据分布场的插值系数；
4.3)根据结构节点的空间位置数据和步骤4.2)中得到的局部气动Cp数据分布场的插值系数，计算得到结构节点Cp数据。
优选地，上述气动载荷转换为结构载荷的方法中，所述对单元载荷数据在结构节点上进行分载计算的具体流程包括：
5.1)获取结构单元的节点编号及其空间坐标；
5.2)根据结构单元的类型，计算结构单元法向量；
5.3)设置结构单元积分点数目；
5.4)根据积分点参数坐标和结构单元节点的物理坐标，计算得到结构单元高斯点的物理坐标；
5.5)根据高斯点物理坐标和插值系数计算得到高斯点的Cp数据；
5.6)根据结构单元法向量计算任一空间平面旋转至XY平面的旋转矩阵，并将空间单元旋转到XY平面；
5.7)通过高斯积分计算得到空间单元旋转后的压心和压力；
5.8)在旋转后的单元上进行分载；
5.9)将分载得到的基于XY平面单元的节点力旋转回原空间；
5.10)对属于不同结构单元的同一节点的节点载荷进行叠加，得到结构节点力载荷；
5.11)判断所有结构单元是否均完成处理，如没完成处理，则继续重复步骤5.1)-5.10)。
优选地，上述气动载荷转换为结构载荷的方法中，所述结构单元的类型包括三角形和四边形，所述结构单元的类型为三角形时，通过平面单元两个相连的单元边向量叉乘获得垂直于平面的结构单元法向量；所述结构单元的类型为四边形时，通过对角线向量叉乘获得结构单元法向量。
优选地，上述气动载荷转换为结构载荷的方法中，所述在旋转后的单元上进行分载计算具体采用最小优化法，或非线性规划法，或多点排直接分配法进行计算。
优选地，上述气动载荷转换为结构载荷的方法中，所述通过高斯积分计算得到空间单元旋转后的压心和压力的公式如下：
Figure BDA0002818723570000041
其中，F为计算得到的压力，xF，yF，zF是压力作用点，即压心的坐标。
优选地，上述气动载荷转换为结构载荷的方法中，还包括在得到结构节点力载荷后，判断载荷转换前后总压和压心是否一致，如一致，则输出节点载荷文件，如不一致，则对节点力载荷进行修正。
优选地，上述气动载荷转换为结构载荷的方法中，所述对节点力载荷进行修正的具体流程包括：
10.1)读取气动面提供的总压和压心数据；
10.2)基于步骤10.1)的总压和压心数据，计算各个方向上的合力矩信息；
10.3)针对单个方向建立线性方程组，并求解系数；
10.4)根据步骤10.3)获得的系数，结合结构节点坐标计算得到修正后的单方向分量载荷；
10.5)重复步骤10.1)-10.4)，直到X、Y、Z三个方向分量载荷均处理完毕；
10.6)输出修正后的节点力载荷。
与现有技术相比，本发明气动载荷转换为结构载荷的方法包括：气动面载荷预处理，形成气动面Cp数据；结构数据预处理，形成结构节点空间位置数据；根据形成的气动面Cp数据和结构节点空间位置数据，插值计算得到结构节点Cp数据；获得结构节点Cp数据后，在单元面上积分形成单元载荷数据，对单元载荷数据在结构节点上进行分载计算，得到结构节点力载荷本发明通过读取气动载荷数据，继而形成压力场，通过读取结构有限元数据，并在形成的压力场内进行插值，形成单元面压力载荷，并最终转化为节点力载荷。本发明通过读取气动载荷数据，继而形成压力场，通过读取结构有限元数据，并在形成的压力场内进行插值，形成单元面压力载荷，并最终转化为节点力载荷。本发明能够快速有效地转化气动载荷为结构载荷，并保证转换前后的总载荷大小和压心相同。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简要介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明提供的一种气动载荷转换为结构载荷的方法流程图；
图2为本发明提供的插值计算得到结构节点Cp数据的方法流程图；
图3为本发明提供的对单元载荷数据在结构节点上进行分载计算的流程图；
图4为本发明提供的物理坐标系下的四边形单元结构示意图；
图5为图4等参坐标系下的四边形单元结构示意图；
图6为本发明提供的对节点力载荷进行修正的方法流程图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步地详细描述，显然，所描述的实施例仅仅是本发明一部份实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
下面结合说明书附图对本发明实施例作进一步详细描述。
如图1所示，本发明实施例公开了一种气动载荷转换为结构载荷的方法，该方法包括：
气动面载荷预处理，形成气动面Cp数据；结构数据预处理，形成结构节点空间位置数据；
根据形成的气动面Cp数据和结构节点空间位置数据，插值计算得到结构节点Cp数据；
获得结构节点Cp数据后，在单元面上积分形成单元载荷数据，对单元载荷数据在结构节点上进行分载计算，得到结构节点力载荷。
优选地，本发明实施例结构数据的预处理可通过节点文件和单元文件读入结构数据。
本发明采用Python语言进行开发，提供方便的文本输入和文本内容显示窗口，通过读取气动载荷数据，继而形成压力场，通过读取结构有限元数据，并在形成的压力场内进行插值，形成单元面压力载荷，并最终转化为节点力载荷。本发明能够快速有效地将气动载荷转化为结构载荷，并保证转换前后的总载荷大小和压心相同。
优选地，本发明实施例气动面类型包括扁平面和曲面。本发明实施中，气动部门会根据气动面的几何特点给出气动载荷，若气动面接近平面，即气动面为扁平面时，会给出气动面轮廓数据和分布载荷Cp，通过对气动数据的处理，会得到加载面上的Cp值；若气动面是大曲率结构，即气动面为曲面时，会给出详细的气动面网格，可以直接读取预存的数据。
具体地，所述气动面为扁平面时，气动面载荷预处理流程包括：读取输入的气动面轮廓数据，形成分载网格面；读取气动分布载荷Cp，根据双线性插值方法，形成分载网格面Cp数据。所述气动面为曲面时，气动面载荷预处理流程包括：读取预存的气动面节点Cp数据。
本发明分载网格面Cp数据是一种特殊的气动面Cp数据，其数据量相对气动面Cp数据少，一般的气动面Cp数据由CFD计算得出，数据点较多，扁平面Cp数据一般由工程方法得出，数据点较少。
优选地，本发明所述根据双线性插值方法，形成分载网格面Cp数据的流程如下：
假设分载网格面的四个角点的坐标和Cp数据分别为：(x1，y1，cp1)、(x2，y2，cp2)、(x3，y3，cp3)和(x4，y4，cp4)；
由四个角点插值生成双线性曲面，则该面内任一点的参数坐标和Cp数据为：
x(u，w)＝x1+u(x2-x1)+[x3+u(x4-x3)-x1-u(x2-x1)]w
y(u，w)＝y1+u(y2-y1)+[y3+u(y4-y3)-y1-u(y2-y1)]w
cp(u，w)＝cp1(1-u)(1-w)+cp3(1-u)w+cp2u(1-w)+cp4uw
式中：u∈[0，1]；w∈[0，1]，u、w均为参数坐标。
优选地，如图2所示，本发明实施例在预处理完成后，就可以开始对结构节点的Cp数据进行插值计算，具体地，本发明所述插值计算得到结构节点Cp数据的具体流程包括：
4.1)根据结构节点的空间位置数据，确定在设定范围内与所有结构节点关联的气动Cp节点；
4.2)根据步骤4.1)确定的气动Cp节点的空间位置数据和Cp数据，通过样条基函数法或径向基函数法，计算得到局部气动Cp数据分布场的插值系数；
4.3)根据结构节点的空间位置数据和步骤4.2)中得到的局部气动Cp数据分布场的插值系数，计算得到结构节点Cp数据。
本发明需要根据结构节点的坐标位置，找到在一定范围内与之关联的气动Cp节点。具体实施中，重复步骤4.1)直到所有结构节点均找到与之关联的气动Cp节点。
本发明实施例通过径向奇函数法计算得到局部气动Cp数据分布场的插值系数的方式如下：
在定义域x∈Rd上，所有形如φ(x-c)＝φ(||x-c||)及其线性组合而成的函数空间称为由函数φ导出的径向基空间。当xj两两不同时，φ(x-xj)线性无关，故可选其作为某函数子空间的一组基，当xj几乎充满R时，由φ(x-xj)生成的函数空间几乎可以逼近所有函数。常用的基函数有：Gauss函数
Figure BDA0002818723570000081
Markoff函数
Figure BDA0002818723570000082
Multi-Quadric函数φ(r)＝(c2+r2)-β等等。
本发明实施例采用Multi-Quadric函数插值，即为对于给定数据点集{(xj,fj)}n，寻找函数
Figure BDA0002818723570000083
使其满足插值条件
Figure BDA0002818723570000084
k＝1,L,n。针对Multi-Quadric函数插值的唯一性定理：对于MQ基函数
Figure BDA0002818723570000085
插值问题解存在且唯一，同时增加一个多项式基函数，因此取插值函数为：
Figure BDA0002818723570000086
考虑到求解的复杂性，取Pj(x)为常数项，同时已知参数c取0，将线性方程组简化，因此插值函数为：
Figure BDA0002818723570000087
根据无穷远边界条件进行约束：
Figure BDA0002818723570000088
以及已知条件：s(xj)＝fj，即可确定权重aj和b。
本发明实施例通过样条基函数法计算得到局部气动Cp数据分布场的插值系数的方式如下：
在N维欧式空间内域D上有n个向量(每个向量代表一个节点坐标)及其对应的函数值，即已知Xi＝[X1,X2,L,XN]T,(i＝1,2,L,n)。设向量X的函数W＝W(X)(称为真实函数)是单值连续的，根据已知条件构成一个逼近函数，公式如下：
Figure BDA0002818723570000089
式中，
Figure BDA00028187235700000810
c1,c2,L,cN+1+n为待定系数，N是X的维度；ε为给定常量，对于一般平坦曲面ε＝10-2:1，对于有奇性的曲面可取ε＝10-5:10-6；Xpi为第i个已知节点的第p维坐标；Xpk为待求节点Xk的第p维坐标。
根据无穷远处边界条件可得：
Figure BDA0002818723570000091
根据n个已知节点函数值可确定待定系数。
优选地，如图3所示，本发明实施例在获得结构节点的Cp数据后，就可以开始对结构节点进行分载计算，具体地，本发明所述对单元载荷数据在结构节点上进行分载计算的具体流程包括：
5.1)获取结构单元的节点编号及其空间坐标；
5.2)根据结构单元的类型，计算结构单元法向量；
5.3)设置结构单元积分点数目；
5.4)根据积分点参数坐标和结构单元节点的物理坐标，计算得到结构单元高斯点的物理坐标；
5.5)根据高斯点物理坐标和插值系数计算得到高斯点的Cp数据；
5.6)根据结构单元法向量计算任一空间平面旋转至XY平面的旋转矩阵，并将空间单元旋转到XY平面；
5.7)通过高斯积分计算得到空间单元旋转后的压心和压力；
5.8)在旋转后的单元上进行分载；
5.9)将分载得到的基于XY平面单元的节点力旋转回原空间；
5.10)对属于不同结构单元的同一节点的节点载荷进行叠加，得到结构节点力载荷；
5.11)判断所有结构单元是否均完成处理，如没完成处理，则继续重复步骤5.1)-5.10)。
本发明所述结构单元的类型包括三角形和四边形，所述结构单元的类型为三角形时，通过平面单元两个相连的单元边向量叉乘获得垂直于平面的结构单元法向量；所述结构单元的类型为四边形时，通过对角线向量叉乘获得结构单元法向量。
如图4所示为本发明实施例物理坐标系(x,y)下的四边形单元，图5为转换到等参坐标系(ξ,η)下的四边形单元。可见，在等参坐标下，单元形状非常规则，有利于数值积分运算，因此在有限元计算中，常常采用等参元的高斯积分方法对单元进行积分计算。即将一个函数的定积分转化为若干个点的函数值的加权组合。
高斯点根据分布函数多项式的次数选择高斯点的个数，本发明实施例在程序中内置了最多3X 3个高斯点的坐标和权系数。
假定在等参单元上具有分布函数f(ξ,η),对单元进行积分计算即需要计算如下所示定积分：
Figure BDA0002818723570000101
通过高斯积分，可以将上式写为：
Figure BDA0002818723570000102
其中，Hi,Hj是对应的高斯点权系数，可以通过查表获得，f(ξi,ηi)为积分点(ξi,ηi)的函数值，(ξi,ηi)为高斯点的参数坐标，也可以通过查表获得。
优选地，本发明所述在旋转后的单元上进行分载计算可以选择使用最小优化法，或非线性规划法，或多点排直接分配法等三种方法进行计算。
本发明实施例采用多点排直接分配法进行计算的方式如下：假设有限元节点与气动点之间有一个无形的梁，它是以气动点一端为固支的悬臂梁，其自由端上的有限元节点分配到载荷Pj时产生的变形能为：
Figure BDA0002818723570000103
于是全系统的变形能为：
Figure BDA0002818723570000111
载荷的转换要对系统影响最小，即在满足静力等效的条件下，系统的变形能最小，即：
Figure BDA0002818723570000112
其中，n是单元中的节点数。
采用拉格朗日(Lagrange)乘子法建立拉格朗日函数：
Figure BDA0002818723570000113
式中，
Figure BDA0002818723570000114
λ、λx、λy为拉格朗日乘子。为使得F(λλxλy)最小，有：
Figure BDA0002818723570000115
因此有：
Figure BDA0002818723570000116
将上式代入问题方程组，可以得到如下线性代数方程组：
Figure BDA0002818723570000117
对上式求解，即可得到λ、λx、λy，代入
Figure BDA0002818723570000118
即可求得分载后的载荷Pj。
本发明实施例采用该方法满足了静力等效原则和传力路线不变，静力等效原则保证了总载荷和总压心不变，而传力路线不变主要体现在保证载荷真实传递。
本发明实施例采用最小优化法和非线性规划法进行计算的方式如下：
需要求解的问题可以描述为：对于一个处于平面上的单元，找到一组作用在节点上的力载荷
Figure BDA0002818723570000121
满足与作用在压心上的合力
Figure BDA0002818723570000122
和合力矩
Figure BDA0002818723570000123
等效，并使得
Figure BDA0002818723570000124
和或者
Figure BDA0002818723570000125
最小。
Figure BDA0002818723570000126
具体实施中，对于上述优化问题，可以直接使用最小优化方法和非线性规划两种方法进行求解。
实施中，本发明通过高斯积分计算得到空间单元旋转后的压心和压力的公式如下：
Figure BDA0002818723570000127
其中，F为计算得到的压力，xF，yF，zF是压力作用点，即压心的坐标。
优选地，本发明气动载荷转换为结构载荷的方法中，还包括在得到结构节点力载荷后，判断载荷转换前后总压和压心是否一致，如一致，则输出节点载荷文件，如不一致，则对节点力载荷进行修正。
本发明实施中，在获得结构节点的节点力载荷后，需要依据总压和压心与未转换前的气动载荷进行比较，如果两者一致，则不用修正，如果差距较大则需要对节点力载荷进行修正。
如图1、6所示，本发明对节点力载荷进行修正的具体流程包括：
10.1)读取气动面提供的总压和压心数据；
10.2)基于步骤10.1)的总压和压心数据，计算各个方向上的合力矩信息；
10.3)针对单个方向建立线性方程组，并求解系数；
10.4)根据步骤10.3)获得的系数，结合结构节点坐标计算得到修正后的单方向分量载荷；
10.5)重复步骤10.1)-10.4)，直到X、Y、Z三个方向分量载荷均处理完毕；
10.6)输出修正后的节点力载荷。
本发明实施中在求解系数时，假设节点载荷在该方向的分量的修正量与节点其它两个方向的坐标满足线性关系，从而建立线性方程组，求解得到系数。
具体地，本发明实施例根据求解系数，结合结构节点坐标计算得到修正后的分量载荷的修正方式如下：对于节点载荷的修正，假设每个节点载荷修正量与节点坐标为线性关系，通过待定系数法确定。
假设已知总载为(Fx,Fy,Fz)，压心为(X,Y,Z)，可以在x、y、z三个方向对节点载荷的三个分量进行修正。
对于X方向，合力为Fx，合力矩为Mxz＝FxgZ、Mxy＝FxgY，因此处理后的节点载荷在X方向的分量Pix(i＝1,2,L,n)在X方向合力及合力矩需与其总合力及力矩平衡。即：
Figure BDA0002818723570000141
其中，yi、zi为节点yz坐标，假设节点载荷在X方向的分量的修正量与节点坐标满足线性关系，即：ΔPix＝a0+a1yi+a2zi，则可以通过解线性方程组求出系数a0、a1、a2。
同理，可按照上述方法求出节点载荷在Y方向、Z方向上分量的修正量。
综上，本发明提供的气动载荷转换为结构载荷的方法，提高了气动载荷处理的效率和准确性。
尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例做出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (9)

1.一种气动载荷转换为结构载荷的方法，该方法包括：
气动面载荷预处理，形成气动面Cp数据；结构数据预处理，形成结构节点空间位置数据；
根据形成的气动面Cp数据和结构节点空间位置数据，插值计算得到结构节点Cp数据；其特征在于，
还包括获得结构节点Cp数据后，在单元面上积分形成单元载荷数据，对单元载荷数据在结构节点上进行分载计算，得到结构节点力载荷，
所述对单元载荷数据在结构节点上进行分载计算的具体流程包括：
5.1)获取结构单元的节点编号及其空间坐标；
5.2)根据结构单元的类型，计算结构单元法向量；
5.3)设置结构单元积分点数目；
5.4)根据积分点参数坐标和结构单元节点的物理坐标，计算得到结构单元高斯点的物理坐标；
5.5)根据高斯点物理坐标和插值系数计算得到高斯点的Cp数据；
5.6)根据结构单元法向量计算任一空间平面旋转至XY平面的旋转矩阵，并将空间单元旋转到XY平面；
5.7)通过高斯积分计算得到空间单元旋转后的压心和压力；
5.8)在旋转后的单元上进行分载；
5.9)将分载得到的基于XY平面单元的节点力旋转回原空间；
5.10)对属于不同结构单元的同一节点的节点载荷进行叠加，得到结构节点力载荷；
5.11)判断所有结构单元是否均完成处理，如没完成处理，则继续重复步骤5.1)-5.10)。
2.根据权利要求1所述的气动载荷转换为结构载荷的方法，其特征在于，所述气动面类型包括扁平面和曲面；
所述气动面为扁平面时，气动面载荷预处理流程包括：
读取输入的气动面轮廓数据，形成分载网格面；
读取气动分布载荷Cp，根据双线性插值方法，形成分载网格面Cp数据；
所述气动面为曲面时，气动面载荷预处理流程包括：
读取预存的气动面节点Cp数据。
3.根据权利要求2所述的气动载荷转换为结构载荷的方法，其特征在于，所述根据双线性插值方法，形成分载网格面Cp数据的流程如下：
假设分载网格面的四个角点的坐标和Cp数据分别为：(x1，y1，cp1)、(x2，y2，cp2)、(x3，y3，cp3)和(x4，y4，cp4)；
由四个角点插值生成双线性曲面，则该面内任一点的参数坐标和Cp数据为：
x(u，w)＝x1+u(x2-x1)+[x3+u(x4-x3)-x1-u(x2-x1)]w
y(u，w)＝y1+u(y2-y1)+[y3+u(y4-y3)-y1-u(y2-y1)]w
cp(u，w)＝cp1(1-u)(1-w)+cp3(1-u)w+cp2u(1-w)+cp4uw
式中：u∈[0，1]；w∈[0，1]，u、w均为参数坐标。
4.根据权利要求1或2所述的气动载荷转换为结构载荷的方法，其特征在于，所述插值计算得到结构节点Cp数据的具体流程包括：
4.1)根据结构节点的空间位置数据，确定在设定范围内与所有结构节点关联的气动Cp节点；
4.2)根据步骤4.1)确定的气动Cp节点的空间位置数据和Cp 数据，通过样条基函数法或径向基函数法，计算得到局部气动Cp数据分布场的插值系数；
4.3)根据结构节点的空间位置数据和步骤4.2)中得到的局部气动Cp数据分布场的插值系数，计算得到结构节点Cp数据。
5.根据权利要求4所述的气动载荷转换为结构载荷的方法，其特征在于，所述结构单元的类型包括三角形和四边形，所述结构单元的类型为三角形时，通过平面单元两个相连的单元边向量叉乘获得垂直于平面的结构单元法向量；所述结构单元的类型为四边形时，通过对角线向量叉乘获得结构单元法向量。
6.根据权利要求5所述的气动载荷转换为结构载荷的方法，其特征在于，所述在旋转后的单元上进行分载计算具体采用最小优化法，或非线性规划法，或多点排直接分配法进行计算。
7.根据权利要求6所述的气动载荷转换为结构载荷的方法，其特征在于，所述通过高斯积分计算得到空间单元旋转后的压心和压力的公式如下：
Figure FDA0004024119880000031
其中，F为计算得到的压力，xF，yF，zF是压力作用点，即压心的坐标。
8.根据权利要求7所述的气动载荷转换为结构载荷的方法，其特征在于，还包括在得到结构节点力载荷后，判断载荷转换前后总压和压心是否一致，如一致，则输出节点载荷文件，如不一致，则对节点力载荷进行修正。
9.根据权利要求8所述的气动载荷转换为结构载荷的方法，其特征在于，所述对节点力载荷进行修正的具体流程包括：
10.1)读取气动面提供的总压和压心数据；
10.2)基于步骤10.1)的总压和压心数据，计算各个方向上的合力矩信息；
10.3)针对单个方向建立线性方程组，并求解系数；
10.4)根据步骤10.3)获得的系数，结合结构节点坐标计算得到修正后的单方向分量载荷；
10.5)重复步骤10.1)-10.4)，直到X、Y、Z三个方向分量载荷均处理完毕；
10.6)输出修正后的节点力载荷。
CN202011416150.2A 2020-12-04 2020-12-04 一种气动载荷转换为结构载荷的方法 Active CN112434457B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011416150.2A CN112434457B (zh) 2020-12-04 2020-12-04 一种气动载荷转换为结构载荷的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011416150.2A CN112434457B (zh) 2020-12-04 2020-12-04 一种气动载荷转换为结构载荷的方法

Publications (2)

Publication Number Publication Date
CN112434457A CN112434457A (zh) 2021-03-02
CN112434457B true CN112434457B (zh) 2023-03-14

Family

ID=74691670

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011416150.2A Active CN112434457B (zh) 2020-12-04 2020-12-04 一种气动载荷转换为结构载荷的方法

Country Status (1)

Country Link
CN (1) CN112434457B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113553736B (zh) * 2021-06-25 2024-06-07 江苏锐天智能科技股份有限公司 一种结构有限元快速加载方法
CN114626235B (zh) * 2022-03-21 2024-06-07 北京航空航天大学 一种面向复杂多层级结构产品的载荷分解改进方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106156444A (zh) * 2016-07-27 2016-11-23 中国航空工业集团公司西安飞机设计研究所 一种飞机气动载荷处理方法及飞机强度计算方法
CN107038296A (zh) * 2017-04-06 2017-08-11 深圳数设科技有限公司 一种确定飞机翼面中有限元点载荷分布的方法与装置
CN107194079A (zh) * 2017-05-25 2017-09-22 中国商用飞机有限责任公司 像素细分载荷传递方法和系统
CN109726437A (zh) * 2018-12-04 2019-05-07 中国航空工业集团公司西安飞机设计研究所 一种舱门气动载荷等效节点力处理方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106156444A (zh) * 2016-07-27 2016-11-23 中国航空工业集团公司西安飞机设计研究所 一种飞机气动载荷处理方法及飞机强度计算方法
CN107038296A (zh) * 2017-04-06 2017-08-11 深圳数设科技有限公司 一种确定飞机翼面中有限元点载荷分布的方法与装置
CN107194079A (zh) * 2017-05-25 2017-09-22 中国商用飞机有限责任公司 像素细分载荷传递方法和系统
CN109726437A (zh) * 2018-12-04 2019-05-07 中国航空工业集团公司西安飞机设计研究所 一种舱门气动载荷等效节点力处理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于径向点插值方法的柔性螺旋桨气动弹性模拟;张宇等;《上海交通大学学报》;20200928(第09期);全文 *
飞机翼面气动载荷压力分布有限元节点插值方法研究;张建刚;《航空科学技术》;20171215(第12期);全文 *

Also Published As

Publication number Publication date
CN112434457A (zh) 2021-03-02

Similar Documents

Publication Publication Date Title
CN112434457B (zh) 一种气动载荷转换为结构载荷的方法
Coirier et al. An accuracy assessment of Cartesian-mesh approaches for the Euler equations
Hughes et al. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement
Mavriplis An assessment of linear versus nonlinear multigrid methods for unstructured mesh solvers
Long et al. Shear‐flexible subdivision shells
Gofuku et al. Point-tangent/point-normal B-spline curve interpolation by geometric algorithms
Rendall et al. Parallel efficient mesh motion using radial basis functions with application to multi‐bladed rotors
Roh et al. The application of geometrically exact shell elements to B-spline surfaces
CN114722686B (zh) 一种基于有限元分析的大型设备吊耳设计及优化方法
Mavriplis Turbulent flow calculations using unstructured and adaptive meshes
Nishida Fully simultaneous coupling of the full potential equation and the integral boundary layer equations in three dimensions
CN109459206B (zh) 地面试验非定常气动力加载方法
CN111027152B (zh) 一种围绕水翼空化流固耦合数值模拟方法
Mi et al. Isogeometric MITC shell
Hashemian et al. Nonlinear bifurcation analysis of statically loaded free-form curved beams using isogeometric framework and pseudo-arclength continuation
WO2019137236A1 (zh) 优化有限元软件计算精度的非协调插值函数构造方法、系统及存储介质
CN115862771A (zh) 一种基于体细分的超弹性材料模型等几何分析仿真方法
CN110633556A (zh) 一种陶瓷基复合材料流固耦合响应计算方法
Babuška et al. Validation of recipes for the recovery of stresses and derivatives by a computer-based approach
Canuto et al. Numerical solution of partial differential equations in random domains: An application to wind engineering
CN113591356B (zh) 一种保留尖锐特征的非均匀不规则样条基函数的构造方法
Schmitt et al. Coupled Actuator Line and Finite Element Analysis Tool
CN115995279B (zh) 一种材料力学特性评估方法、装置、设备及可读存储介质
CN112560363B (zh) 一种基于映射过程的cfd计算中的网格变形质量评价方法
CN112632825B (zh) 一种基于有限元超收敛性的静电场光滑有限元数值算法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant