CN104933259B - 一种飞机高升力系统载荷计算方法 - Google Patents

一种飞机高升力系统载荷计算方法 Download PDF

Info

Publication number
CN104933259B
CN104933259B CN201510367111.0A CN201510367111A CN104933259B CN 104933259 B CN104933259 B CN 104933259B CN 201510367111 A CN201510367111 A CN 201510367111A CN 104933259 B CN104933259 B CN 104933259B
Authority
CN
China
Prior art keywords
msub
mrow
mfrac
suspension mechanism
centerdot
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
CN201510367111.0A
Other languages
English (en)
Other versions
CN104933259A (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.)
Xian Aircraft Design and Research Institute of AVIC
Original Assignee
Xian Aircraft Design and Research Institute of AVIC
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 Xian Aircraft Design and Research Institute of AVIC filed Critical Xian Aircraft Design and Research Institute of AVIC
Priority to CN201510367111.0A priority Critical patent/CN104933259B/zh
Publication of CN104933259A publication Critical patent/CN104933259A/zh
Application granted granted Critical
Publication of CN104933259B publication Critical patent/CN104933259B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明涉及航空系统设计领域,具体涉及一种飞机高升力系统载荷计算方法,以解决现有计算方法效率低,无法满足设计要求问题。高升力系统中的增升装置与悬挂机构多连接点连接,悬挂机构与高升力系统中的多个运动机构点接触,依据传力路线上各环节的受力与约束情况,可以把传力路线分解为有机联系的增升装置‑悬挂机构、悬挂机构‑运动机构子系统,再基于各子系统设计理念和力学特征,分步解决复杂系统的传力计算问题,计算速度快、计算结果准确可靠,可以完全避免有限元软件所需的复杂而又冗长的过程。

Description

一种飞机高升力系统载荷计算方法
技术领域
本发明涉及航空系统设计领域,具体涉及一种飞机高升力系统载荷计算方法。
背景技术
现有技术中,采用商用有限元软件NASTRAN或MARC进行飞机高升力系统载荷计算,但是,对于设计载荷情况多而复杂的高升力系统,以往的方法计算效率低,甚至不能完成高升力系统的载荷计算工作,制约着我国飞机系统设计。
发明内容
本发明的目的是提供一种飞机高升力系统载荷计算方法,以解决现有计算方法效率低,无法满足设计要求问题。
本发明的技术方案是:
一种飞机高升力系统载荷计算方法,所述高升力系统中的增升装置与悬挂机构多连接点连接,所述悬挂机构与所述高升力系统中的多个运动机构点接触,包括如下步骤:
步骤一、将所述高升力系统的传力路线分解为增升装置-悬挂机构、悬挂机构-运动机构;
步骤二、为所述高升力系统施加预定载荷,使所述增升装置受力,并将该力通过所述多连接点传递至所述悬挂机构;
步骤三、获取所述多连接点的刚度,根据所述预定载荷和所述多连接点的刚度得到各连接点的载荷,从而得到所述悬挂机构的总载荷;
步骤四、根据所述悬挂机构与所述多个运动机构的力学平衡方程,得到所述悬挂机构与每一个所述运动机构的接触点载荷。
可选地,在所述步骤三中:
首先根据最少势能原理公式(1):
其中,a1...an为增升装置悬挂点待定位移;
得到公式(2):
其中,a为悬挂点待定位移向量,P为所述预定载荷,K为各所述连接点的刚度矩阵;
再取位移泛函(3)为:
其中,N为形函数,u为位移泛函,为各悬挂点的近似位移场函数,i为连接点个数,ai为连接点数为i时的增升装置悬挂点待定位移,a为悬挂点待定位移向量;
再根据公式(4)得到各连接点的载荷∏:
其中,u为位移泛函,uT为u的转置矩阵,K为各所述连接点的刚度矩阵。
可选地,在所述步骤四中,所述悬挂机构与所述运动机构的接触点为三个,对应的载荷分别为滑轨上前滚轮载荷Fq、为滑轨上后前滚轮载荷Fh、啮合齿廓之间的正压力Ft,所述悬挂机构与所述运动机构的力学平衡方程为:
其中,FY、FZ、MX为所述悬挂机构的总载荷在X、Y、Z三个方向上的分载荷,α是Fq的方向角,β为Fh的方向角,θ为Ft的方向角,Rhg为滑轨外径,L为Ft到滑轨轴线距离,L1为齿轮与齿条间摩擦力到滑轨轴线距离,μ为齿轮摩擦系统。
本发明的有益效果:
本发明的高升力系统载荷计算方法,依据传力路线上各环节的受力与约束情况,可以把传力路线分解为有机联系的增升装置-悬挂机构、悬挂机构-运动机构子系统,再基于各子系统设计理念和力学特征,分步解决复杂系统的传力计算问题,计算速度快、计算结果准确可靠,可以完全避免有限元软件所需的复杂而又冗长的过程。
附图说明
图1为本发明所述实施例的流程图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。
本发明提供了一种飞机高升力系统载荷计算方法;其中,高升力系统中的增升装置与悬挂机构通过多连接点连接,悬挂机构与高升力系统中的多个运动机构点接触。
本发明的方法包括如下步骤:
步骤S101、将所述高升力系统的传力路线分解为增升装置-悬挂机构、悬挂机构-运动机构。
增升装置-高升力系统是复杂流场下的一个复杂传力系统,通常是气动外载-翼面(缝翼或前中后襟翼)-悬挂机构-非线性环节(运动机构)和驱动装置构成了一条可分级研究的传力路线。依据传力路线上各环节的受力与约束情况,可以把传力路线分解为有机联系的至少两个子系统,即增升装置-悬挂机构、悬挂机构-运动机构子系统。
后续再基于各子系统设计理念和力学特征,应用力学变分原理或二次规划进行力学分析,从而分步解决复杂系统的传力计算问题,计算速度快、计算结果准确可靠,可以完全避免有限元软件所需的复杂而又冗长的过程。。
步骤S102、为所述高升力系统施加预定载荷,使所述增升装置受力,并将该力通过所述多连接点传递至所述悬挂机构。
步骤S103、获取所述多连接点的刚度,根据所述预定载荷和所述多连接点的刚度得到各连接点的载荷,从而得到所述悬挂机构的总载荷。
由热力学第一定律,如令物体的动能为Ek,应变能为U,则在微小的δt时间间隔内,物体从一种状态过渡到另一种状态时,总能量的变化为:
δEk+δU=δW+δQ,
其中,δW为体力Fbi与面力pi所完成的功,δQ为物体由其周围介质所吸收(或向外散发)的热量。
假定弹性变形过程是绝热的,则对于静力平衡问题有δQ=0,δEK=0,δU=δW。
由虚位移原理:在外力作用下处于平衡状态的变形体,当给予物体微小虚位移时,外力的总虚功等于物体的总虚应变能,即δU=δW。
由变分原理,微分方程的边值问题可化为一泛函的极值问题:
由里兹法,取位移泛函为:
其中,N为形函数,u为位移泛函,为各悬挂点的近似位移场函数,i为连接点个数,ai为连接点数为i时的增升装置悬挂点待定位移,a为悬挂点待定位移向量。
泛函的变分为零相当于将泛函对所包含的待定参数进行全微分,并令所得的方程等于零,即:
由于δa1,δa2,…是任意的,满足上式时必然有…都等于零。因此可以得到一组方程为:
这是与待定参数a的个数相等的方程组,用以求解a;对于二次泛函,上式退化为一组线性方程为:
其中,a1...an为增升装置悬挂点待定位移,a为悬挂点待定位移向量,P为所述预定载荷,K为各所述连接点的刚度矩阵。
再根据如下公式得到各连接点的载荷∏:
其中,u为位移泛函,uT为u的转置矩阵,K为各所述连接点的刚度矩阵。
对于基于二次规划的传力计算方法,有以下原理:
分析表明,悬挂机构的变形能为正定二次型,由等式约束二次规划的一般形式:
s.t.Ax=b,
H∈Rn×m为对称阵。特别当H正定时,目标函数为凸函数,线性约束下可行域又是凸集,二次规划为凸二次规划,此时,二次规划具有如下良好的性质:
K-T条件不仅是最优解的必要条件而且是充分条件;局部最优解就是全局最优解。
基于以上原理及分析(对于程序接口、方程组和二次规划的解法和策略等不再赘述),针对该传力计算问题,对于增生装置-悬挂机构子系统传力计算,基于变分法与二次规划原理实现复杂翼面上气动分布载荷到挂点集中载荷的计算。通过该方法,可以计算得到增升装置各悬挂机构(如:主襟翼各挂点、缝翼各滑轨挂点)载荷,每个运动机构各挂点(连接点)载荷合成即得到其总载荷。
步骤S104、根据所述悬挂机构与所述多个运动机构的力学平衡方程,得到所述悬挂机构与每一个所述运动机构的接触点载荷。
悬挂机构-非线性环节(运动机构传力问题)子系统传力计算,在已知总载情况下,通过简化分析、隔离体分析,建立滑轨或滑轮架平衡方程即可完成非线性环节的传力计算。构建可考虑变摩擦情况的高升力系统运动机构全参数方程组:
其中,FY、FZ、MX为所述悬挂机构的总载荷在X、Y、Z三个方向上的分载荷,α是Fq的方向角,β为Fh的方向角,θ为Ft的方向角,Rhg为滑轨外径,L为Ft到滑轨轴线距离,L1为齿轮与齿条间摩擦力到滑轨轴线距离,μ为齿轮摩擦系统。通过以上计算即可完成:悬挂机构传力、运动机构各接触体载荷、高升力系统载荷计算。
实施实例
某飞机高升力系统设计载荷计算。
在采用商用有限元软件MSC/nastran进行大运飞机高升力系统设计载荷计算时,出现了计算速度慢、结果奇异等问题,不能满足设计要求。这些问题是由于nastran软件的非线性模块的缺陷造成的,其他的非线性商用有限元软件,无不需要通过有限元建模、载荷分配、非线性分析等过程来完成。对于飞机设计多载荷情况、复杂系统构型、多缝道、多舵面结构等来说,这些过程太过复杂冗长。
经过深入研究,本专利给出了基于变分法与二次规划的高升力系统载荷计算方法,计算速度快,可以完全避免有限元软件所需的复杂而又冗长的过程,完成了大运飞机的高升力系统设计载荷计算。
飞行试验与地面验证试验表明:本专利给出的计算方法准确可靠。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (3)

1.一种飞机高升力系统载荷计算方法,所述高升力系统中的增升装置与悬挂机构多连接点连接,所述悬挂机构与所述高升力系统中的多个运动机构点接触,其特征在于,包括如下步骤:
步骤一、将所述高升力系统的传力路线分解为增升装置-悬挂机构、悬挂机构-运动机构;
步骤二、为所述高升力系统施加预定载荷,使所述增升装置受力,并将该力通过所述多连接点传递至所述悬挂机构;
步骤三、获取所述多连接点的刚度,根据所述预定载荷和所述多连接点的刚度得到各连接点的载荷,从而得到所述悬挂机构的总载荷;
步骤四、根据所述悬挂机构与所述多个运动机构的力学平衡方程,得到所述悬挂机构与每一个所述运动机构的接触点载荷。
2.根据权利要求1所述的飞机高升力系统载荷计算方法,其特征在于,在所述步骤三中:
首先根据最少势能原理公式(1):
<mrow> <mi>&amp;delta;</mi> <mo>&amp;Pi;</mo> <mo>=</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mo>&amp;Pi;</mo> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>a</mi> <mn>1</mn> </msub> </mrow> </mfrac> <msub> <mi>&amp;delta;a</mi> <mn>1</mn> </msub> <mo>+</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mo>&amp;Pi;</mo> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>a</mi> <mn>2</mn> </msub> </mrow> </mfrac> <msub> <mi>&amp;delta;a</mi> <mn>2</mn> </msub> <mo>+</mo> <mo>...</mo> <mo>+</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mo>&amp;Pi;</mo> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>a</mi> <mi>n</mi> </msub> </mrow> </mfrac> <msub> <mi>&amp;delta;a</mi> <mi>n</mi> </msub> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中,a1...an为增升装置悬挂点待定位移;
得到公式(2):
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <mo>&amp;Pi;</mo> </mrow> <mrow> <mo>&amp;part;</mo> <mi>a</mi> </mrow> </mfrac> <mo>&amp;equiv;</mo> <mi>K</mi> <mi>a</mi> <mo>-</mo> <mi>P</mi> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
其中,a为悬挂点待定位移向量,P为所述预定载荷,K为各所述连接点的刚度矩阵;
再取位移泛函(3)为:
<mrow> <mi>u</mi> <mo>&amp;ap;</mo> <mover> <mi>u</mi> <mo>~</mo> </mover> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>N</mi> <mi>i</mi> </msub> <msub> <mi>a</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>N</mi> <mi>a</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
其中,N为形函数,u为位移泛函,为各悬挂点的近似位移场函数,i为连接点个数,ai为连接点数为i时的增升装置悬挂点待定位移,a为悬挂点待定位移向量;
再根据公式(4)得到各连接点的载荷∏:
<mrow> <mo>&amp;Pi;</mo> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>u</mi> <mi>T</mi> </msup> <mi>K</mi> <mi>u</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中,u为位移泛函,uT为u的转置矩阵,K为各所述连接点的刚度矩阵。
3.根据权利要求1或2所述的飞机高升力系统载荷计算方法,其特征在于,在所述步骤四中,所述悬挂机构与所述运动机构的接触点为三个,分别为滑轨上前滚轮载荷Fq、为滑轨上后前滚轮载荷Fh、啮合齿廓之间的正压力Ft,所述悬挂机构与所述运动机构的力学平衡方程为:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msub> <mi>F</mi> <mi>q</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <msub> <mi>&amp;mu;</mi> <mn>1</mn> </msub> <msub> <mi>F</mi> <mi>q</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>sin</mi> <mo>(</mo> <mfrac> <mi>&amp;pi;</mi> <mn>2</mn> </mfrac> <mo>+</mo> <mi>&amp;alpha;</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>sin</mi> <mi>&amp;beta;</mi> <mo>+</mo> <msub> <mi>&amp;mu;</mi> <mn>2</mn> </msub> <msub> <mi>F</mi> <mi>h</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>sin</mi> <mo>(</mo> <mfrac> <mi>&amp;pi;</mi> <mn>2</mn> </mfrac> <mo>+</mo> <mi>&amp;beta;</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>F</mi> <mi>t</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>sin</mi> <mi>&amp;theta;</mi> <mo>-</mo> <mi>&amp;mu;</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mo>&amp;CenterDot;</mo> <msub> <mi>F</mi> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>F</mi> <mi>Y</mi> </msub> <mo>=</mo> <mn>0.0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>F</mi> <mi>q</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <msub> <mi>&amp;mu;</mi> <mn>1</mn> </msub> <msub> <mi>F</mi> <mi>q</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>cos</mi> <mo>(</mo> <mfrac> <mi>&amp;pi;</mi> <mn>2</mn> </mfrac> <mo>+</mo> <mi>a</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;beta;</mi> <mo>+</mo> <msub> <mi>&amp;mu;</mi> <mn>2</mn> </msub> <msub> <mi>F</mi> <mi>h</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mo>(</mo> <mfrac> <mi>&amp;pi;</mi> <mn>2</mn> </mfrac> <mo>+</mo> <mi>&amp;beta;</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>F</mi> <mi>t</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mo>-</mo> <mi>&amp;mu;</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> <mo>&amp;CenterDot;</mo> <msub> <mi>F</mi> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>F</mi> <mi>Z</mi> </msub> <mo>=</mo> <mn>0.0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mrow> <mi>h</mi> <mi>g</mi> </mrow> </msub> <mo>&amp;times;</mo> <mo>(</mo> <mfrac> <mn>0.07</mn> <mn>2.55</mn> </mfrac> <msub> <mi>F</mi> <mi>q</mi> </msub> <mo>+</mo> <mfrac> <mn>0.07</mn> <mn>1.80</mn> </mfrac> <msub> <mi>F</mi> <mi>h</mi> </msub> <mo>)</mo> <mo>+</mo> <msub> <mi>F</mi> <mi>t</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>L</mi> <mo>-</mo> <mi>&amp;mu;</mi> <mo>&amp;CenterDot;</mo> <msub> <mi>F</mi> <mi>t</mi> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>L</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>M</mi> <mi>X</mi> </msub> <mo>=</mo> <mn>0.0</mn> </mtd> </mtr> </mtable> </mfenced>
其中,FY、FZ、MX为所述悬挂机构的总载荷在X、Y、Z三个方向上的分载荷,α是Fq的方向角,β为Fh的方向角,θ为Ft的方向角,Rhg为滑轨外径,L为Ft到滑轨轴线距离,L1为齿轮与齿条间摩擦力到滑轨轴线距离,μ为齿轮摩擦系统。
CN201510367111.0A 2015-06-29 2015-06-29 一种飞机高升力系统载荷计算方法 Active CN104933259B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510367111.0A CN104933259B (zh) 2015-06-29 2015-06-29 一种飞机高升力系统载荷计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510367111.0A CN104933259B (zh) 2015-06-29 2015-06-29 一种飞机高升力系统载荷计算方法

Publications (2)

Publication Number Publication Date
CN104933259A CN104933259A (zh) 2015-09-23
CN104933259B true CN104933259B (zh) 2018-04-13

Family

ID=54120425

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510367111.0A Active CN104933259B (zh) 2015-06-29 2015-06-29 一种飞机高升力系统载荷计算方法

Country Status (1)

Country Link
CN (1) CN104933259B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777743B (zh) * 2016-12-28 2017-10-24 石家庄铁道大学 液压悬挂结构的有限元建模方法
CN107515109B (zh) * 2017-08-04 2019-06-28 中国航空工业集团公司西安飞机设计研究所 一种试验件总载测量方法
CN108009334B (zh) * 2017-11-22 2021-03-30 中国航空工业集团公司西安飞机设计研究所 一种飞机货运系统吊轨接头载荷计算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1882474A (zh) * 2003-11-12 2006-12-20 德国空中客车公司 飞机高升力系统的驱动系统中的载荷限定方法
CN103502095A (zh) * 2011-04-28 2014-01-08 空中客车德国运营有限责任公司 用于飞行器的高升力系统和用于影响飞行器的高升力特性的方法
CN203756912U (zh) * 2013-12-04 2014-08-06 中国飞机强度研究所 一种载荷分解装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102004044961B4 (de) * 2004-09-16 2007-02-01 Liebherr-Aerospace Lindenberg Gmbh Vorrichtung zur Erfassung von Gleichlauffehlern von Hochauftriebsflächen an Flugzeugen

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1882474A (zh) * 2003-11-12 2006-12-20 德国空中客车公司 飞机高升力系统的驱动系统中的载荷限定方法
CN103502095A (zh) * 2011-04-28 2014-01-08 空中客车德国运营有限责任公司 用于飞行器的高升力系统和用于影响飞行器的高升力特性的方法
CN203756912U (zh) * 2013-12-04 2014-08-06 中国飞机强度研究所 一种载荷分解装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
民机高升力系统操作试验的符合性验证和设计方法研究;刘彦生 等;《航空工程进展》;20120229;第3卷(第1期);第125-130页 *

Also Published As

Publication number Publication date
CN104933259A (zh) 2015-09-23

Similar Documents

Publication Publication Date Title
CN104933259B (zh) 一种飞机高升力系统载荷计算方法
CN102804187B (zh) 物理量计算方法、数值解析方法、物理量计算装置及数值解析装置
CN104461677B (zh) 一种基于cfd和fem技术的虚拟热试验方法
CN103077259A (zh) 高超声速导弹多场耦合动力学一体化仿真分析方法
CN105205293B (zh) 用于获得飞机部件气动载荷的方法和系统
CN102508436B (zh) 机械手摩擦力动力学精确分析与控制应用方法
CN102865866B (zh) 基于双星敏感器的卫星姿态确定方法及定姿误差分析方法
CN106295146B (zh) 矿井瓦斯爆炸冲击波对人体损伤破坏的仿真评估方法
CN107391858A (zh) 一种获取风洞模型静气动弹性变形影响量的方法
CN105205267A (zh) 一种机翼整体油箱载荷计算方法
She et al. On the impact force of human-robot interaction: Joint compliance vs. link compliance
CN105468814A (zh) 一种考虑全机气动干扰的直升机飞行特性计算方法
CN105117524A (zh) 一种采用差动刹车控制飞机转弯过程的动态仿真方法
CN107451354A (zh) 一种鸭式布局飞行器气动参数的仿真方法及终端设备
Kiviaho et al. A robust and flexible coupling framework for aeroelastic analysis and optimization
Ghoreyshi et al. Computational aerodynamic modeling for flight dynamics simulation of ram-air parachutes
CN105260499A (zh) 一种三维柔性梁系统力学的运动仿真方法及系统
CN103123668B (zh) 一种基于混合单元法的空间绳系机器人系统的仿真方法
CN103678897B (zh) 一种基于凯恩方程的飞轮隔振平台专用动力学建模方法
CN110287507A (zh) 一种应用于恒压变量液压柱塞泵疲劳寿命分析方法
CN106773782A (zh) 一种气动伺服弹性混合建模方法
CN106339542A (zh) 一种轮胎的整体带束层安全倍数计算方法
CN111090908A (zh) 一种群伞系统设计参数选取方法
CN114218672A (zh) 跨超声速大迎角配平翼的抖振响应获取方法及相关装置
CN102673801B (zh) 井字型反扭矩舵结构布局方法

Legal Events

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