CN114239152A - 一种飞机矩形气密舱变形计算方法 - Google Patents
一种飞机矩形气密舱变形计算方法 Download PDFInfo
- Publication number
- CN114239152A CN114239152A CN202111666757.0A CN202111666757A CN114239152A CN 114239152 A CN114239152 A CN 114239152A CN 202111666757 A CN202111666757 A CN 202111666757A CN 114239152 A CN114239152 A CN 114239152A
- Authority
- CN
- China
- Prior art keywords
- deformation
- rectangular
- point
- calculation method
- equation
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本申请属于航空结构强度分析领域,为一种飞机矩形气密舱变形计算方法,对于非常规布局矩形气密舱,先建立矩形气密舱力学模型,根据力学模型结构对称、载荷对称的特点再次简化力学模型为1/4结构模型,先计算1/4模型角点处的弯矩,再计算整个模型长边和短边的变形方程,能够得到力学模型在任意尺寸下的变形量,从而得到变形规律;通过本发明计算获得的变形方程能够得到准确且唯一的任意尺寸下的矩形气密舱变形情况,不需要采用有限元分析,计算简单,且能够进行强度分析自动化的集成,具有广泛的适用性。
Description
技术领域
本申请属于航空结构强度分析领域,特别涉及一种飞机矩形气密舱变形计算方法。
背景技术
随着航空技术的发展,现代飞行器可以在广泛的空间范围内飞行。飞行高度高,飞行速度大,给飞机设备和机上人员的工作效率带来不利的影响。为了保证飞机座舱内的合适环境需要对座舱进行空气调节。气密座舱是飞机座舱环境控制系统的调节对象,是对旅客和机组人员提供舒适而安全的生活和工作环境的基本技术保证,飞机的气密舱是飞机重要构件之一。
翼身融合布局飞机具有大升阻比、低阻力、低噪声等优点,是未来飞机最具潜力的发展方向之一。翼身融合式布局飞机结构与传统飞机结构最主要的区别就是压力舱由传统的圆筒形变成为扁平的矩形盒子形状,如图2所示,由此使得飞机结构的受力形式发生了巨大的变化。
变形计算是飞机强度设计中的一项重要内容,应用有限元软件可以得到较为精确的变形计算结果,但是对于网格密度的选取、边界条件和载荷的施加等均有一定要求,不同的设计人员可能计算得到差异较大的结果。同时,进行有限元分析比较繁琐,且当结构尺寸改变时,需要重新建立有限元模型进行计算,且不利于强度分析自动化的集成。因此,在飞机设计中,应用本发明进行强度设计是不可或缺的。而现阶段缺少非常规布局矩形气密舱变形计算的本发明。
发明内容
本申请的目的是提供了一种飞机矩形气密舱变形计算方法,以解决现有技术中采用有限元分析对气密舱进行变形计算时计算繁琐、不同设计人员差异较大的问题。
本申请的技术方案是:一种飞机矩形气密舱变形计算方法,包括,
建立矩形气密舱简化的力学模型;以左下角点为坐标原点o,建立坐标系,x轴向右为正,y轴向上为正,z轴符合右手坐标系;根据对称性,建立简化的1/4结构模型,计算1/4结构模型角点处的弯矩值;分析并计算力学模型中长边和短边的变形方程;根据变形方程,分别研究力学模型长边和短边在不同长宽比下的变形特点,得出变形规律。
优选地,所述1/4结构模型角点处的弯矩值Mo计算方法为在对称剖面位置K处,剪力Fq为0,设定弯矩以内表面拉伸为正、剪力以使微元顺时针旋转为正,以轴向拉伸为正;根据对称性,得到K点处的轴力FNK为FNK=qL/2,其中q为单位长度上气密舱的内压载荷,L为气密舱剖面的长边长度,H为气密舱剖面的短边长度,矩形气密舱剖面的弯矩刚度为EI;根据载荷平衡,得到角点处的弯矩值大小Mo为
根据结构对称性,力学模型内四个角点的弯矩值大小均相同。
优选地,所述长边变形方程的计算方法为
优选地,所述短边变形方程的计算方法为
优选地,所述长边变形规律的计算方法为
令BC边变形方程的一阶导数为0,得到
根据
优选地,设定矩形舱体剖面的长宽比为k,即L/H=k,令DC边变形方程的一阶导数为0得到变形规律,所述短边变形规律的计算方法为
u(y2)″=u(y3)″>0,y2、y3为两个极小值点;u(y1)″<0,y1为极大值点;
极大值记为um,且有
本申请的一种飞机矩形气密舱变形计算方法,对于非常规布局矩形气密舱,先建立矩形气密舱力学模型,根据力学模型结构对称、载荷对称的特点再次简化力学模型为1/4结构模型,先计算1/4模型角点处的弯矩,再计算整个模型长边和短边的变形方程,能够得到力学模型在任意尺寸下的变形量,从而得到变形规律;通过本发明计算获得的变形方程能够得到准确且唯一的任意尺寸下的矩形气密舱变形情况,不需要采用有限元分析,计算简单,且能够进行强度分析自动化的集成,具有广泛的适用性。
附图说明
为了更清楚地说明本申请提供的技术方案,下面将对附图作简单地介绍。显而易见地,下面描述的附图仅仅是本申请的一些实施例。
图1为本申请整体流程示意图;
图2为背景技术结构示意图;
图3为本申请力学模型结构示意图;
图4为本申请1/4结构模型示意图;
图5为本申请BC边受力示意图;
图6为本申请长边任意剖面处受力示意图;
图7为本申请DC边受力示意图;
图8为本申请短边任意剖面处受力示意图;
图9为本申请k=1时的舱体短边变形示意图;
图14为本申请代入具体数据的算例示意图;
图15为本申请代入具体数据的计算变形结果示意图;
图16为有限元长边和短边计算结果图;
图17为有限元短边计算结果放大图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。
一种飞机矩形气密舱变形计算方法,如图1所示,包括以下步骤:
步骤S100,建立矩形气密舱简化的力学模型;
其简化的力学模型图见图3。
图3中,q为单位长度上气密舱的内压载荷,L为气密舱剖面的长边长度,H为气密舱剖面的短边长度,矩形气密舱剖面的弯矩刚度为EI;
气密舱左下、左上、右上、右下四个点分别记为A、B、C、D;
以左下角点A为坐标原点o,建立坐标系,x轴向右为正,y轴向上为正,z轴符合右手坐标系。
如图4所示,步骤S200,根据对称性,建立简化的1/4结构模型,计算1/4结构模型角点处的弯矩值;
由于力学模型中结构对称,载荷对称,力学模型能够简化为1/4结构模型,1/4结构模型角点处的弯矩值计算方法为:
对称剖面位置K处,剪力Fq为0,在K点处有轴力FNK以及弯矩Mt两个未知量,设定弯矩定义内表面拉伸为正;剪力为使微元顺时针旋转为正;轴向拉伸为正;
根据对称性,得到K点处的轴力FNK为:
FNK=qL/2 (1)
根据内压载荷q、轴力FNK和弯矩Mt在K点处产生的转角为0的变形协调方程,得到:
根据载荷平衡,得到角点C处的弯矩值大小Mo为:
根据结构对称性,力学模型内四个角点的弯矩值大小均相同。
步骤S300,分析并计算力学模型中长边和短边的变形方程;
如图5、图6所示,长边变形方程的计算方法为:
设定FBy和FCy分别为B、C点的剪力,根据力系平衡,得到:
设置BC边的任意剖面处的剪力为FQy,BC边任意剖面的弯矩为M(x),根据力矩平衡,得到:
设定BC边沿x轴的变形为v(x),则:
两次积分后得到:
根据边界条件:
v(0)=0,v(L)=0 (8)
将式(8)带入式(7),求解方程组得到:
将式(3)、式(9)代入式(7),得到BC边的变形方程为:
如图7、图8示意图,短边变形方程的计算方法为:
设定FDx和FCx分别为D、C点的剪力,根据力系平衡,得到:
设定FQx为任意剖面y处的剪力;M(y)为DC边任意剖面的弯矩,根据力矩平衡,得到:
设定DC边沿y轴的变形为u(y),得到:
两次积分后得到:
根据边界条件:
u(0)=0,u(H)=0 (15)
将式(15)代入式(14)并求解方程组获得:
将式(3)、式(16)代入式(14)可得DC边的变形方程为:
步骤S400,根据变形方程,分别研究力学模型长边和短边在不同长宽比下的变形特点,得出变形规律。
变形方程的计算方法为:
令式(10)的一阶导数为0,即:
v(x)′=0 (18)
得到:
根据:
设定矩形舱体剖面的长宽比为k,即:
L/H=k (22)
令式(17)的一阶导数为0,即:
u(y)′=0 (23)。
如图9-13所示,得到以下变形规律:
u(y2)″=u(y3)″>0,y2、y3为两个极小值点;u(y1)″<0,y1为极大值点;
极大值记为um,且有:
设定极大值点为um≤0,得到:
不同长宽比的变形规律结论:
1)气密载荷作用下,飞机矩形气密舱长边的变形全部为正,即气密载荷下长边全部向外鼓;
2)当长宽比k=1时为正方形气密舱,正方形舱体的变形两边一致,为全部向外鼓;
6)随着长宽比值的增大,气密舱短边的变形变化趋势为:全部向外鼓——仅中间段鼓,两侧内凹——正中间点变形为0,两侧变形内凹——全部内凹。
对于非常规布局矩形气密舱,先建立矩形气密舱力学模型,根据力学模型结构对称、载荷对称的特点再次简化力学模型为1/4结构模型,先计算1/4模型角点处的弯矩,再计算整个模型长边和短边的变形方程,能够得到力学模型在任意尺寸下的变形量,从而得到变形规律;通过本发明计算获得的变形方程能够得到准确且唯一的任意尺寸下的矩形气密舱变形情况,不需要采用有限元分析,计算简单,且能够进行强度分析自动化的集成,具有广泛的适用性。
作为一种具体实施方式,还包括一种飞机矩形气密舱变形计算方法,以一个具体实例对该飞机矩形气密舱变形方法进行计算,具体为:
如图14所示,已知:某一矩形剖面舱体,承受气密载荷q=0.001MPa,设定剖面弯曲刚度为EI=1.29×109Nmm2,H=3000mm,L=3600mm。
计算剖面任意点的变形。
1)若计算1点处(与B点的距离为1000mm)的变形
将x=1000mm、以及q、H、L、EI的值代入下式:
可得:v(1)=0.38mm。
2)若计算2点处(BC边中点,与B点距离1800mm)的变形
将x=1800mm、以及q、H、L、EI的值代入下式:
可得:v(2)=0.54mm。
求BC边任意点的位移值,仅需代入对应点的坐标值x即可。
3)若计算3点处的变形(与D点距离为500mm)
将y=500mm、以及q、H、L、EI的值代入下式:
可得:u(3)=-0.0379mm。
求DC边任意点的位移值,仅需代入对应点的坐标值y即可。
AB、AD边的变形可由对称性求得。
如图15所示,可以看到力学模型长边和短边的变形情况,与图16、图17中的有限元计算结果中的长边和短边变形量对比可以明显看出,两种方法得到的变形量基本相同,同时本申请所用方法与有限元模型计算相比具有更强的适用性。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种飞机矩形气密舱变形计算方法,其特征在于:包括,
建立矩形气密舱简化的力学模型;以左下角点为坐标原点o,建立坐标系,x轴向右为正,y轴向上为正,z轴符合右手坐标系;
根据对称性,建立简化的1/4结构模型,计算1/4结构模型角点处的弯矩值;
分析并计算力学模型中长边和短边的变形方程;
根据变形方程,分别研究力学模型长边和短边在不同长宽比下的变形特点,得出变形规律。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111666757.0A CN114239152B (zh) | 2021-12-31 | 2021-12-31 | 一种飞机矩形气密舱变形计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111666757.0A CN114239152B (zh) | 2021-12-31 | 2021-12-31 | 一种飞机矩形气密舱变形计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114239152A true CN114239152A (zh) | 2022-03-25 |
CN114239152B CN114239152B (zh) | 2022-09-20 |
Family
ID=80745093
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111666757.0A Active CN114239152B (zh) | 2021-12-31 | 2021-12-31 | 一种飞机矩形气密舱变形计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114239152B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107016181A (zh) * | 2017-03-30 | 2017-08-04 | 中国航空工业集团公司西安飞机设计研究所 | 一种舱门结构静强度计算方法 |
CN107066703A (zh) * | 2017-03-23 | 2017-08-18 | 西安飞机工业(集团)有限责任公司 | 一种舱门结构静强度计算方法 |
-
2021
- 2021-12-31 CN CN202111666757.0A patent/CN114239152B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107066703A (zh) * | 2017-03-23 | 2017-08-18 | 西安飞机工业(集团)有限责任公司 | 一种舱门结构静强度计算方法 |
CN107016181A (zh) * | 2017-03-30 | 2017-08-04 | 中国航空工业集团公司西安飞机设计研究所 | 一种舱门结构静强度计算方法 |
Non-Patent Citations (2)
Title |
---|
陈庆强等: "22000m~3液化气船整船和舱段三维有限元强度分析", 《船舶工程》 * |
陈庆强等: "2万2千方液化气船整船和舱段三维有限元强度分析", 《船舶力学》 * |
Also Published As
Publication number | Publication date |
---|---|
CN114239152B (zh) | 2022-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109460596A (zh) | 一种飞翼无人机非线性载荷计算方法 | |
CN110210130B (zh) | 针对工字梁二维模型的形状优化方法 | |
CN108515518A (zh) | 柔性支撑工业机器人的一种工作空间求解方法 | |
CN112001034B (zh) | 一种多面锥构型飞行风洞标模布局及设计方法 | |
CN109858071B (zh) | 一种考虑剪力滞后作用的薄壁箱梁结构动力特性分析方法 | |
CN110398339A (zh) | 一种满足多级全自由飞风洞试验相似律的条件确定方法 | |
CN114239152B (zh) | 一种飞机矩形气密舱变形计算方法 | |
CN107016218B (zh) | 一种确定飞机翼尖小翼翼面中有限元点载荷分布的方法与装置 | |
CN104951595B (zh) | 一种蒙皮零件工艺耳片建模方法 | |
Ritter et al. | Comparison of nonlinear aeroelastic methods for maneuver simulation of very flexible aircraft | |
CN114154434A (zh) | 一种无尾飞翼布局的多约束精细化气动优化设计方法 | |
CN106773782B (zh) | 一种气动伺服弹性混合建模方法 | |
CN112580240B (zh) | 一种适用于复杂大柔性飞机建模的非线性子结构方法 | |
CN109657337A (zh) | 一种快速模拟单面搭接件螺栓连接的建模方法 | |
CN117494611A (zh) | 一种基于有粘当地流活塞理论的超声速颤振分析方法 | |
CN112229420A (zh) | 一种用于飞机蒙皮对缝测量的线激光标定方法 | |
CN106647270B (zh) | 针对空间密频结构的稳定自适应模糊主动振动控制方法 | |
CN114491796A (zh) | 一种异型剖面气密舱载荷确定方法 | |
Rubio et al. | Modal characterization and structural aerodynamic response of a crane fly forewing | |
CN113704881B (zh) | 一种结构有限元模型惯性载荷施加方法 | |
US20240104267A1 (en) | Method of designing aircraft part and method of producing aircraft part | |
CN117828768A (zh) | 一种异型剖面气密舱载荷确定方法 | |
CN107844643A (zh) | 一种导弹弹体弹性变形下导弹压心变化分析方法 | |
Simpson | Unsteady aerodynamics, reduced-order modelling, and predictive control in linear and nonlinear aeroelasticity with arbitrary kinematics | |
CN117313276A (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 |