CN110162826A - 薄壁结构热气动弹性动响应分析方法 - Google Patents
薄壁结构热气动弹性动响应分析方法 Download PDFInfo
- Publication number
- CN110162826A CN110162826A CN201910212466.0A CN201910212466A CN110162826A CN 110162826 A CN110162826 A CN 110162826A CN 201910212466 A CN201910212466 A CN 201910212466A CN 110162826 A CN110162826 A CN 110162826A
- Authority
- CN
- China
- Prior art keywords
- thin
- wall construction
- matrix
- dynamic response
- finite element
- 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.)
- Granted
Links
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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提出一种薄壁结构热气动弹性动响应分析方法,根据薄壁结构和边界条件,将薄壁结构离散为七个自由度壳单元,将加肋离散为七个自由度梁单元,建立薄壁的有限元模型,得到结构有限元网格;将薄壁结构温度场施加到结构有限元网格上,对有限元模型进行线性颤振分析,得到每个减缩频率k对应的非定常气动力矩阵Qk,拟合出非定常气动力矩阵的时域表达式;对得到的非定常气动力矩阵的时域表达式进行迭代计算,进行薄壁的非线性瞬态响应分析,得到薄壁结构非线性热气动弹性动响应。本发明的方法可以解决复杂薄壁结构的非线性热气动弹性动响应的计算问题。
Description
技术领域
本发明涉及热气动弹性响应分析,属于飞行器气动弹性技术技术领域。
背景技术
在超声速、高超声速飞行器上广泛使用大尺寸薄壁结构,如热防护结构、发动机保护罩和大尺寸翼面蒙皮等,这类薄壁结构容易发生弹性力、惯性力、气动力和热应力四者之间相互耦合引发的热气动弹性问题,这种耦合往往导致飞行器的性能下降甚至结构破坏。另外,由于薄壁结构横向挠度在板中面引起面内应力,面内应力在振动中产生一种硬弹性的效应,这是一种典型的结构几何非线性问题。
对于壁板颤振的经典方法如Galerkin方法、Rayleigh-Ritz方法、谐波平衡方法等仅均限与形状规则的矩形薄板,且边界条件为简支边界或者固支边界条件,对于形状不规则、边界条件复杂的情况,经典方法不适用。
采用有限元方法不仅适用于形状不规则、边界条件复杂的薄壁结构,而且能够考虑气动加热、面内内力及不同来流方向等复杂问题。航空航天领域常用的有限元分析软件MSC.Nastran提供了结构线性气动弹性响应的计算功能,但无法进行结构非线性的气动弹性动响应分析。
发明内容
本发明的目的在于克服现有技术不足,提供一种薄壁结构非线性气动弹性动响应分析方法,可以获得薄壁结构的响应和非定常气动力随时间变化历程。
本发明的技术解决方案:
一种薄壁结构热气动弹性动响应分析方法,包括以下步骤:
步骤一,根据薄壁结构和边界条件,将薄壁结构离散为空间壳单元,将加肋离散为梁单元,建立薄壁的有限元模型,所述的壳单元每个节点有七个自由度ΔC=[u v w w,x w,yw,xy θz]T,其中,u、v表示壳单元两个方向的面内位移,w、w,x、w,y、w,xy分别表示壳单元四个弯曲位移,θz表示壳单元扭转角;
所述的梁单元每个节点有七个自由度,ΔL=[u′ v′ w′ θx θy θ′z θx,x]
其中u′、v′表示梁单元两个方向的面内位移,w'表示梁单元弯曲位移,θx,θy,θ′z分别代表三个方向的角位移,θx,x表示梁单元的扭转角;
步骤二,根据薄壁结构建立气动面网格,并建立气动面网格和结构有限元模型之间的插值关系,得到结构有限元网格;
步骤三,将薄壁结构温度场施加到步骤二中得到的结构有限元网格上,对有限元模型进行线性颤振分析,得到每个减缩频率k对应的非定常气动力矩阵Qk并拟合出Qk的时域表达式;
步骤四,通过对步骤三得到的非定常气动力矩阵的时域表达式进行迭代计算,进行薄壁的非线性瞬态响应分析,得到薄壁结构非线性热气动弹性动响应;
所述的步骤三中非定常气动力矩阵的时域表达式为:其中Q0-Q2+j是实待定矩阵,s是拉氏复变量,rj是由经验给定的正实数,r1=0.1~0.25,r2=0.3~0.45,r3=0.5~0.65,r2=0.7~0.9,b为薄壁结构沿来流方向的长度,V为飞行器的速度。
所述的薄壁结构非线性热气动弹性动响应包括薄壁结构非线性热气动弹性气动力和薄壁结构的位移。
所述的薄壁结构的位移求解表达式为:
w=HwTbwb,u=HuTmwm,v=HvTmwm
其中
其中,wm={u,v}T,wb={w,w,x,w,y,w,xy,θz}T,为薄壁结构沿来流方向的宽度,为薄壁结构沿来流方向的长度;
Hw=[1 x y x2 xy y2 x3 x2y xy2 y3 x3y x2y2 xy3 x3y2 x2y3 x3y3],Hu=[1 x yxy 0 0 0 0],Hv=[0 0 0 0 1 x y xy],x,y,z分别表示三个坐标方向。
所述的薄壁结构非线性热气动弹性气动力表达式为:
式中,qa为气流动压,M∞为飞行器马赫数,γ为比热比,w,t为弯曲运动速度,Hw,x为矩阵Hw对x的导数。
本发明与现有技术相比的有益效果:
(1)本发明通过构建加肋壁板结构的壳单元和空间梁单元,并在壳单元和空间梁单元上增加了表示扭转自由度的扭转角,可以更真实模拟非线性结构的运动效应,通过扭转自由度的加入,将扭转角引入到了单位刚度矩阵中,应用在有限元模型的建立中,由此方法建立的有限元模型对复杂结构形式和混合边界条件的薄壁具有很好适应性,适于对飞行器薄壁结构的气动弹性动响应分析;
(2)本发明提出了非线性壁板气动弹性动响应瞬态过程的非定常气动力求解方法,可直接用于评估壁板颤振对局部气动特性的影响;
(3)本发明将时域气动力矩阵代入到非线性瞬态响应分析程序中,兼顾了非线性瞬态响应,可以实现薄壁结构几何非线性特性的模拟、同时引入了非定常气动力,实现了薄壁结构的非线性热气动弹性动响应分析。
附图说明
图1为本发明非线性热气动弹性动响应分析方法流程图;
图2为本发明薄壁结构示意图;
图3为本发明薄壁结构壳单元;
图4为本发明薄壁结构空间梁单元;
图5为本发明薄壁结构随时间变化的位移响应图;
图6为本发明薄壁结构气动弹性动响应相位图。
具体实施方式
下面结合具体实例及附图对本发明进行详细说明。
根据图1所示,一种薄壁结构热气动弹性动响应分析方法,包括以下步骤:
步骤一,根据薄壁结构和边界条件建立薄壁的有限元模型,薄壁结构如图2所示,薄壁结构离散为空间壳单元,如图3所示,对于加筋薄壁结构,还需要将加筋离散为空间梁单元,如图4所示。
对于薄壁结构,每个节点有七个自由度ΔC=[u v w w,x w,y w,xy θz]T,包括两个面内位移wm={u,v}T和四个弯曲位移wb={w,w,x,w,y,w,xy,θz}T,单元如图5所示;空间梁单元每个节点有七个自由度ΔL=[u′ v′ w′ θx θy θ′z θx,x],如图4所示其中θz和θx,x均代表扭转角。
对于壳单元,位移函数可以表示为:
w=α1+α2x+α3y+α4x2+α5xy+α6y2+α7x3+α8x2y+α9xy2+
α10y3+α11x3y+α12x2y2+α13xy3+α14x3y2+α15x2y3+α16x3y3
=Hwα
u=ζ1+ζ2x+ζ3y+ζ4xy=Huζ
v=ζ5+ζ6x+ζ7y+ζ8xy=Hvζ
其中Hw=[1 x y x2 xy y2 x3 x2y xy2 y3 x3y x2y2 xy3 x3y2 x2y3 x3y3],α=[α1α2 … α16]T,Hu=[1 x y xy 0 0 0 0],Hv=[0 0 0 0 1 x y xy],ζ=[ζ1 ζ2 … ζ8]T。
为了确定α和ζ,可将图3中结点1,2,3,4的坐标代入下面的表达式
其中,wi、wi,x、wi,y、wi,xy表示结点i的弯曲位移,Hwi、Hwi,x、Hwi,y、Hwi,xy表示结点i的弯曲位移的计算系数,ui、vi表示结点i的面内位移,Hui、Hvi表示结点i的面内位移的计算系数。
对于空间梁单元,采用线性插值计算,公式如下:
u'=u1(1-ξ)+u2ξ
v'=v1F1+θz1lF2+v2F3+θz2lF4
w'=w1F1+θy1lF2+w2F3+θy2lF4
θx=θx1F1+(θx,x)1lF2+θx2F3+(θx,x)2lF4
其中l为梁单元长度,ξ=x/l,F1=1-3ξ2+2ξ3,F2=ξ-2ξ2+ξ3,F3=3ξ2-2ξ3,F4=ξ3-ξ2,u1、u2、v1、v2为x和y方向的位移、F1、F2、F3、F4为系数、θz′1、θz′2、θy1、θy2、θx1、θx2为两个节点在三方向的角位移、(θx,x)1、(θx,x)2为两个节点的扭转角,下标1和2分别代表空间梁单元的两个节点。
设F=(F1 F2 F3 F4),则梁单元的弯曲刚度矩阵为:
其中,EI为弯曲刚度。
则梁单元的扭转刚度矩阵:
其中,GJ为扭转刚度。
因为拉压、扭转和双向弯曲互不耦联,从而可由拉压、扭转和两个方向的弯曲单元根据位移编号顺序通过组合来形成单元刚度矩阵,单元刚度矩阵为对称矩阵,表示如下:
其中aL=EA/l,bL=12EIz/l3,cL=6EIz/l2,dL=12EIy/l3,eL=6EIy/l2,fL=6GJ/5l,gL=GJ/10,hL=2EIy/l,iL=2EIz/l,jL=2GJl/15,kL=-GJl/30,EA为单元轴向拉压刚度,EIz为xoy平面内的弯曲刚度,EIy为xoz平面内的弯曲刚度,Iy是筋条对面板中面的惯性矩,设A1为筋条在yoz平面的截面面积,GJ为扭转刚度,J=Iy+Iz,通过单元刚度矩阵可以得到结构的整体刚度矩阵Ks,在有限元模型的建立时使用;
步骤二,根据薄壁结构建立气动面网格,并建立气动面网格和结构有限元模型之间的插值关系,插值矩阵为Gka,在MSC.Nastran的前后处理软件MSC.Patran将壁板结构温度场施加到结构有限元网格上,采用MSC.Nastran求解序列SOL103对有限元模型进行线性颤振分析,得到不同减缩频率k下的非定常气动力矩阵Qk,利用Roger有理函数拟合的方法拟合出非定常气动力的时域表达式,非定常气动力的有理函数表达式为:
其中Q0-Q2+j是实待定矩阵,s是拉氏复变量,rj是由经验给定的正实数,r1=0.2,r2=0.4,r3=0.6,r4=0.8;
在一个实施例中,减缩频率k取为颤振减缩频率附近的若干k值,每个k值对应一个非定常气动力矩阵,矩阵包含实部和虚部,计算得到的一组k下的非定常气动力矩阵Qk,将简谐条件带入式,并分解为实部、虚部,其中是拉氏复变量、eiωk是简谐条件的表达方式。在一个实施例中,取7个为一组,则Qk中每个元素q的待定系数为7个,其中实部为虚部为
采用最小二乘法可求得系数q0-q6,由此完成了气动力从频域到时域的转化,得到时域非定常气动力。
壁板气动弹性动响应系统质量矩阵为M=Ms-Ma,阻尼矩阵为C=Cs(t)-Ca,刚度矩阵为K=Ks(t)-Ka,其中Ms、Cs(t)、Ks(t)为结构刚度矩阵、阻尼矩阵、刚度矩阵,为气动力产生的附加质量矩阵,为气动阻尼矩阵,为气动刚度矩阵;
步骤三,在MSC.Nastran的前后处理软件MSC.Patran将壁板结构温度场施加到结构有限元网格上,并进行非线性瞬态响应分析的设置,生成求解序列为SOL129的计算文件.bdf;利用MSC.Nastran的二次开发语言DMAP将壁板系统的质量矩阵、阻尼矩阵和刚度矩阵写入非线性瞬态响应分析的计算文件.bdf中,替换有限元模型自身的质量矩阵、阻尼矩阵和刚度矩阵,作为薄壁气动弹性动响应分析的初始刚度矩阵,其中阻尼矩阵和刚度矩阵均为响应的函数,因此在每一个非线性迭代步进行更新;
步骤四,采用MSC.Nastran求解序列SOL129进行壁板非线性瞬态响应分析,可完成薄壁结构非线性热气动弹性动响应分析。
通过求解可以得到结构任意点的位移:
w=HwTbwb,u=HuTmwm,v=HvTmwm
其中
非定常气动力表示为:
结合结构响应的表达式,可得到耦合非线性结构响应的非定常气动力求解方法为:
式中,qa=ρaV2/2为气流动压,M∞为飞行器马赫数,γ为比热比,w,t为弯曲运动速度,Hw,x为矩阵Hw对x的导数。
本发明未详细说明部分为本领域技术人员公知技术。
Claims (4)
1.一种薄壁结构热气动弹性动响应分析方法,其特征在于:包括以下步骤,
步骤一,根据薄壁结构和边界条件,将薄壁结构离散为空间壳单元,将加肋也离散为壳单元,建立薄壁的有限元模型,所述的壳单元每个节点有七个自由度ΔC=[u v w w,x w,yw,xy θz]T,其中,u、v表示壳单元两个方向的面内位移,w、w,x、w,y、w,xy分别表示壳单元四个弯曲位移,θz表示壳单元扭转角;
所述的梁单元每个节点有七个自由度,ΔL=[u′ v′ w′ θx θy θz′ θx,x]
其中u′、v′表示梁单元两个方向的面内位移,w'表示梁单元弯曲位移,θx,θy,θz′分别代表在x,y,z方向的角位移,θx,x表示梁单元的扭转角;
步骤二,根据薄壁结构建立气动面网格,建立气动面网格和结构有限元模型之间的插值关系,得到结构有限元网格;
步骤三,将薄壁结构温度场施加到步骤二中得到的结构有限元网格上,对有限元模型进行线性颤振分析,得到每个减缩频率k对应的非定常气动力矩阵Qk,拟合出非定常气动力矩阵的时域表达式;
步骤四,通过对步骤三得到的非定常气动力矩阵的时域表达式进行迭代计算,得到薄壁结构非线性热气动弹性动响应。
2.根据权利要求1所述的薄壁结构热气动弹性动响应分析方法,其特征在于:非定常气动力矩阵的时域表达式为:
其中Q0-Q2+j是实待定矩阵,s是拉氏复变量,rj是由经验给定的正实数,r1=0.1~0.25,r2=0.3~0.45,r3=0.5~0.65,r2=0.7~0.9,b为薄壁结构沿来流方向的长度,V为飞行器的速度。
3.根据权利要求1所述的薄壁结构热气动弹性动响应分析方法,其特征在于:所述的薄壁结构的位移求解表达式为:
w=HwTbwb,u=HuTmwm,v=HvTmwm
其中
其中,wm={u,v}T,wb={w,w,x,w,y,w,xy,θz}T,为薄壁结构沿来流方向的宽度,为薄壁结构沿来流方向的长度,
Hw=[1 x y x2 xy y2 x3 x2y xy2 y3 x3y x2y2 xy3 x3y2 x2y3 x3y3],Hu=[1 x y xy 00 0 0],Hv=[0 0 0 0 1 x y xy],x,y,z分别表示三个坐标。
4.根据权利要求1所述的薄壁结构热气动弹性动响应分析方法,其特征在于:所述的薄壁结构非线性热气动弹性气动力表达式为:
式中,qa为气流动压,M∞为马赫数,γ为比热比,w,t为弯曲运动速度,Hw,x为矩阵Hw对x的导数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910212466.0A CN110162826B (zh) | 2019-03-20 | 2019-03-20 | 薄壁结构热气动弹性动响应分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910212466.0A CN110162826B (zh) | 2019-03-20 | 2019-03-20 | 薄壁结构热气动弹性动响应分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110162826A true CN110162826A (zh) | 2019-08-23 |
CN110162826B CN110162826B (zh) | 2021-05-11 |
Family
ID=67639012
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910212466.0A Active CN110162826B (zh) | 2019-03-20 | 2019-03-20 | 薄壁结构热气动弹性动响应分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110162826B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112364544A (zh) * | 2020-11-19 | 2021-02-12 | 中国空气动力研究与发展中心超高速空气动力研究所 | 再入气动环境致结构热力响应有限元求解优化方法 |
CN112487577A (zh) * | 2020-11-27 | 2021-03-12 | 重庆大学 | 一种结构稳态非线性动力学响应快速计算方法及系统和存储介质 |
CN113361145A (zh) * | 2021-07-19 | 2021-09-07 | 重庆大学 | 一种拉压不同模量矩形薄板热应力的确定方法 |
CN113935116A (zh) * | 2021-12-16 | 2022-01-14 | 清华大学 | 航空发动机波纹管低速冲击响应分析方法 |
CN114970033A (zh) * | 2022-06-07 | 2022-08-30 | 中国核工业二三建设有限公司 | 一种大型薄壁设备吊装过程快速有限元求解方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102866637A (zh) * | 2012-10-07 | 2013-01-09 | 西北工业大学 | 一种基于二次降阶的带操纵面机翼非定常气动力模拟方法 |
CN104133933A (zh) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | 一种高超声速飞行器热环境下气动弹性力学特性分析方法 |
CN105046021A (zh) * | 2015-08-25 | 2015-11-11 | 西北工业大学 | 非定常气动力最小状态有理近似的非线性优化算法 |
CN106096088A (zh) * | 2016-05-31 | 2016-11-09 | 中国航空工业集团公司西安飞机设计研究所 | 一种螺旋桨飞机螺旋颤振分析方法 |
CN107103103A (zh) * | 2016-02-22 | 2017-08-29 | 上海机电工程研究所 | 二维间隙结构非线性气动弹性模型建模方法 |
CN108052772A (zh) * | 2017-12-30 | 2018-05-18 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性分析方法 |
CN108363843A (zh) * | 2018-01-25 | 2018-08-03 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性全机配平方法 |
CN109359418A (zh) * | 2018-11-09 | 2019-02-19 | 深圳大学 | 一种悬索桥非线性风致响应的预测方法、系统及存储介质 |
CN109492234A (zh) * | 2017-09-11 | 2019-03-19 | 北京机电工程研究所 | 一种改进的流固耦合插值方法 |
-
2019
- 2019-03-20 CN CN201910212466.0A patent/CN110162826B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102866637A (zh) * | 2012-10-07 | 2013-01-09 | 西北工业大学 | 一种基于二次降阶的带操纵面机翼非定常气动力模拟方法 |
CN104133933A (zh) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | 一种高超声速飞行器热环境下气动弹性力学特性分析方法 |
CN105046021A (zh) * | 2015-08-25 | 2015-11-11 | 西北工业大学 | 非定常气动力最小状态有理近似的非线性优化算法 |
CN107103103A (zh) * | 2016-02-22 | 2017-08-29 | 上海机电工程研究所 | 二维间隙结构非线性气动弹性模型建模方法 |
CN106096088A (zh) * | 2016-05-31 | 2016-11-09 | 中国航空工业集团公司西安飞机设计研究所 | 一种螺旋桨飞机螺旋颤振分析方法 |
CN109492234A (zh) * | 2017-09-11 | 2019-03-19 | 北京机电工程研究所 | 一种改进的流固耦合插值方法 |
CN108052772A (zh) * | 2017-12-30 | 2018-05-18 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性分析方法 |
CN108363843A (zh) * | 2018-01-25 | 2018-08-03 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性全机配平方法 |
CN109359418A (zh) * | 2018-11-09 | 2019-02-19 | 深圳大学 | 一种悬索桥非线性风致响应的预测方法、系统及存储介质 |
Non-Patent Citations (4)
Title |
---|
夏巍等: "超声速飞行器壁板非线性颤振响应分析的时域法与频域法对比研究", 《固体力学学报》 * |
彭家忠: "某大长径比弹箭飞行器的气弹特性研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
聂国隽等: "考虑约束扭转的薄壁梁单元刚度矩阵", 《计算力学学报》 * |
高怡宁等: "基于最小状态法的气动力有理函数拟合精度影响因素研究", 《航空科学技术》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112364544A (zh) * | 2020-11-19 | 2021-02-12 | 中国空气动力研究与发展中心超高速空气动力研究所 | 再入气动环境致结构热力响应有限元求解优化方法 |
CN112364544B (zh) * | 2020-11-19 | 2022-04-12 | 中国空气动力研究与发展中心超高速空气动力研究所 | 再入气动环境致结构热力响应有限元求解优化方法 |
CN112487577A (zh) * | 2020-11-27 | 2021-03-12 | 重庆大学 | 一种结构稳态非线性动力学响应快速计算方法及系统和存储介质 |
CN113361145A (zh) * | 2021-07-19 | 2021-09-07 | 重庆大学 | 一种拉压不同模量矩形薄板热应力的确定方法 |
CN113361145B (zh) * | 2021-07-19 | 2022-12-23 | 重庆大学 | 一种拉压不同模量矩形薄板热应力的确定方法 |
CN113935116A (zh) * | 2021-12-16 | 2022-01-14 | 清华大学 | 航空发动机波纹管低速冲击响应分析方法 |
CN114970033A (zh) * | 2022-06-07 | 2022-08-30 | 中国核工业二三建设有限公司 | 一种大型薄壁设备吊装过程快速有限元求解方法及系统 |
CN114970033B (zh) * | 2022-06-07 | 2023-09-19 | 中国核工业二三建设有限公司 | 一种大型薄壁设备吊装过程快速有限元求解方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110162826B (zh) | 2021-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110162826A (zh) | 薄壁结构热气动弹性动响应分析方法 | |
Palacios et al. | Structural and aerodynamic models in nonlinear flight dynamics of very flexible aircraft | |
Li et al. | Aeroelastic flutter analysis for 2D Kirchhoff and Mindlin panels with different boundary conditions in supersonic airflow | |
Smith et al. | CFD-based analysis of nonlinear aeroelastic behavior of high-aspect ratio wings | |
Zhao et al. | Modal analysis of a rotating thin plate via absolute nodal coordinate formulation | |
Song et al. | Investigations on the flutter properties of supersonic panels with different boundary conditions | |
CN109933876A (zh) | 一种基于广义气动力的非定常气动力降阶方法 | |
Kantor et al. | Nonlinear structural, nonlinear aerodynamic model for static aeroelastic problems | |
Attar et al. | Modeling delta wing limit-cycle oscillations using a high-fidelity structural model | |
Zhang et al. | Linearized Euler solver for rapid frequency-domain aeroelastic analysis | |
Tian et al. | Approximate global mode method for flutter analysis of folding wings | |
CN108182330A (zh) | 一种基于b样条计算柔性矩形薄板刚柔耦合动力学响应的方法 | |
Natarajan et al. | Supersonic flutter analysis of functionally graded material plates with cracks | |
CN107766686B (zh) | 基于matlab计算fgm薄板刚柔耦合动力学响应的仿真方法 | |
Abbas et al. | Panel flutter analysis of plate element based on the absolute nodal coordinate formulation | |
Shin et al. | Nonlinear aeroelastic analysis for a control fin with an actuator | |
Malik et al. | Aeroelastic analysis of a slender wing | |
Lamberson et al. | Aeroelastic simulations with modal and finite-element structural solvers using CREATE-AV/Kestrel v5 | |
Palacios et al. | Structural models for flight dynamic analysis of very flexible aircraft | |
Kandil et al. | Buffet responses of a vertical tail in vortex breakdown flows | |
CN109726454B (zh) | 管路系统的流固耦合建模方法及装置 | |
De Marqui Junior et al. | Identification of flutter parameters for a wing model | |
Levin et al. | Limit cycle oscillations of plate-type fins using increased-order models | |
Huang et al. | Correlation studies of geometrically nonlinear aeroelastic formulations with beam and shell elements | |
Gang et al. | Active control law design for Flutter/LCO suppression based on reduced order model method |
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 |