CN102968526B - 一种大展弦比机翼颤振模型配重确定方法 - Google Patents
一种大展弦比机翼颤振模型配重确定方法 Download PDFInfo
- Publication number
- CN102968526B CN102968526B CN201210452011.4A CN201210452011A CN102968526B CN 102968526 B CN102968526 B CN 102968526B CN 201210452011 A CN201210452011 A CN 201210452011A CN 102968526 B CN102968526 B CN 102968526B
- Authority
- CN
- China
- Prior art keywords
- axis
- counterweight
- eta
- mass
- point
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000003014 reinforcing effect Effects 0.000 claims description 12
- 239000002023 wood Substances 0.000 claims description 10
- 230000007704 transition Effects 0.000 claims description 7
- 238000013461 design Methods 0.000 abstract description 16
- 238000011161 development Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 abstract description 2
- 238000007796 conventional method Methods 0.000 abstract 1
- 230000007423 decrease Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 8
- 238000005457 optimization Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于航空结构力学领域,涉及一种大展现比机翼颤振模型配重确定方法。主要如下:根据配重计算的输入数据,计算配重目标点特性;进一步结合配重目标点质量大小mP、绕刚轴的惯量大小IPx,以及坐标yP计算yA和yB;更进一步根据特征点G、H、I和J和质心目标的平面坐标,计算参数n;最后确定配重点C、D、E和F坐标和配重点质量大小。本发明思路简洁明晰,减少了模型设计的不确定性,缩短了确定配重的时间,提高了颤振模型的设计效率。相比于现有一般方法,本发明的计算效率可大幅提高,以往需要数小时甚至数天完成的工作,目前只需要数秒钟即可,大大缩短了颤振模型的设计周期,并为颤振模型试验工作和飞机研制带来便利。
Description
技术领域
本发明属于航空结构力学领域,特别是涉及到一种大展弦比机翼颤振模型配重确定方法。
背景技术
颤振模型可以用来获取飞机及其部件的颤振特性,对于颤振模型的设计,通过调整配重大小和位置,保证模型质量大小、质心位置及绕轴惯量和目标值相等,配重对于颤振模型的特性模拟具有举足轻重的作用;具体针对大展弦比机翼颤振模型,需要保证机翼分框段质量、质心及绕刚轴惯量和目标值相等;
现有方法一般根据框段内机翼模型的总质量、质心及绕刚轴惯量和框段内梁、木框、蒙皮、加强肋等部件的质量、质心及绕刚轴惯量获得配重目标的质量、质心及绕刚轴惯量;并进一步通过人工试凑与迭代优化的方法获得具体的配重位置和大小。
现有的方法没有引入模型的几何边界特性,存在配重位置超出模型边界的可能;而且其难做到逆向设计,试凑与迭代优化的方法给模型设计带来了很多的不确定性,且效率低下。
发明内容
本发明的目的是:提供一种大展弦比机翼颤振模型配重确定方法,能够获得大展弦比机翼颤振模型设计的模型配重大小及位置确定方法,考虑模型几何边界特性,提高设计效率。
本发明的技术方案是:
步骤一、建立坐标系:以机翼根部为原点o,以机翼的刚心轴为x轴,在机翼平面内建立y轴;
步骤二、确定配重目标点P质量的大小mP、惯量的大小IPx及坐标(xP,yP):
其中,
mP=mq-ml-mkd.................................[1]
IPx=Iqx-Ilx-Ikdx.................................[2]
xP=(mqxq-mlxl-mkdxkd)/mP.................................[3]
yP=(mqyq-mlyl-mkdykd)/mP.................................[4]
其中,
mq是框段内机翼模型的总质量,
ml是框段内梁的质量,
mkd是框段内木框、蒙皮、加强肋等部件的总质量,
Iqx是框段内机翼模型绕刚心轴x轴的总惯量,
Ilx是框段内梁绕刚心轴x轴的惯量,
Ikdx是框段内木框、蒙皮、加强肋等部件绕刚心轴x轴的总惯量,
(xq,yq)是框段内机翼模型质心,
(xl,yl)是框段内梁质心,
(xkd,ykd)是框段内木框、蒙皮、加强肋等部件的总质心,
步骤三、以配重目标点P为原点,以ξ轴为横坐标,以η为纵坐标,且ξ轴平行于x轴,η轴平行于y轴,建立局部坐标系;
步骤四、确定局部坐标下的特征点的位置
其中i分别表示四个特征点G、H、I、J;
步骤五、确定过渡点A,B的坐标
式[6]要求
其中,
步骤六、根据杠杆原理确定过渡点A和B质量大小
步骤七、确定配重点C、D、E和F坐标;
由于点C、D位于线段GK、KH上,点E、F位于线段IL、LJ上,可进一步求得xC、xD、xE、xF;直线CE和DF平行于ox轴;
步骤八、根据杠杆原理确定配重点C、D、E和F质量大小,完成单框配重;
步骤九、重复上述步骤直至完成其他框的配重。
本发明的优点是:本发明方法引入了模型几何边界特性,即在配重设计初始就引入了颤振模型框段的边界位置信息,因此在配重设计时减少了不确定性,在一定程度上避免了配重位置可能超出模型前后缘的情况,减小了重复设计的可能性;由于本发明使用了逆向设计方法,避免了试凑及人工优化可能带来的不确定性和效率低下问题;另外这种逆向设计方法原理简单,形式明晰流畅,非常适合于各种通用计算程序的编制运行,很容易为工程技术人员所掌握;相比于现有一般方法,本发明的计算效率可大幅提高,以往需要数小时甚至数天完成的工作,目前只需要数秒钟即可,大大缩短了颤振模型的设计周期,并为颤振模型试验工作和飞机研制带来便利。
附图说明
图1为本发明原理图;图中坐标系xoy原点o位于机翼根部,x轴为机翼的刚心轴,y轴在机翼平面内垂直于x轴;局部坐标系ζPη原点位于P点,ξ轴平行于x轴,η轴平行于y轴;线段GK、KH和线段IL、LJ为模型加强肋的位置,配重点C、D、E和F位于其上;点A、B为过渡点。
具体实施方式
下面结合附图对本发明做进一步详细描述,请参阅图1。
如图1所示,一种大展弦比机翼颤振模型配重确定方法,包括以下步骤:
步骤一、建立坐标系:以机翼根部为原点o,以机翼的刚心轴为x轴,机翼的刚心轴是拟合而成的一条直线。在机翼平面内建立y轴,并且y轴垂直于x轴向上;
步骤二、确定配重目标点P质量的大小mP、绕刚轴惯量的大小IPx以及坐标位置(xP,yP):
其中,
mP=mq-ml-mkd.................................[1]
IPx=Iqx-Ilx-Ikdx.................................[2]
xP=(mqxq-mlxl-mkdxkd)/mP.................................[3]
yP=(mqyq-mlyl-mkdykd)/mP.................................[4]
其中,
mq是框段内机翼模型的总质量,
ml是框段内梁的质量,
mkd是框段内木框、蒙皮、加强肋等部件的总质量,
Iqx是框段内机翼模型绕刚心轴x轴的总惯量,
Ilx是框段内梁绕刚心轴x轴的惯量,
Ikdx是框段内木框、蒙皮、加强肋等部件绕刚心轴x轴的总惯量,
(xq,yq)是框段内机翼模型质心,
(xl,yl)是框段内梁质心,
(xkd,ykd)是框段内木框、蒙皮、加强肋等部件的总质心,
步骤三、以配重目标点P为原点,以ξ轴为横坐标,以η为纵坐标,且ξ轴平行于x轴,η轴平行于y轴,建立局部坐标系;
步骤四、确定局部坐标系下的特征点的位置
其中i分别表示四个特征点G、H、I、J,而这四个特征点是框段的四个边界顶点;
步骤五、确定过渡点A,B的纵坐标
式[6]要求
其中,
步骤六、根据杠杆原理确定过渡点A和B质量大小,即将配重目标点P的质量mP,分为mA和mB,
步骤七、确定配重点C、D、E和F坐标;
由于点C、D位于线段GK、KH上,点E、F位于线段IL、LJ上,可进一步求得xC、xD、xE、xF;直线CE和DF平行于ox轴;
步骤八、根据杠杆原理确定配重点C、D、E和F质量大小,完成单框配重;
步骤九、重复上述步骤直至完成其他框的配重。
实施例
以某飞机大展弦比机翼颤振模型其中一个肋段的配重设计为例,进行计算验证。
(1)根据配重计算输入数据,计算配重目标点特性;
(2)进一步结合配重目标点质量大小mP、绕刚轴的惯量大小IPx,以及坐标yP计算yA和yB;
(3)进一步根据特征点G、H、I和J和质心目标的平面坐标,计算n;
(4)确定配重点C、D、E和F坐标和配重点质量大小。
对本发明所述方法,进行计算验证,根据表1给出的配重计算输入数据计算得到配重目标点特性,如表2所示;结合表3给出的特征点二维坐标,进行配重坐标和大小设计,设计结果如表4所示。由于本发明方法使用了逆向设计的方法,计算效率也大大提高,以往需要数小时甚至数天完成的工作,目前只需要数秒钟即可。
表1配重目标计算输入数据
表2配重目标点特性
表3特征点二维坐标
表4配重计算结果
Claims (1)
1.一种大展弦比机翼颤振模型配重确定方法,其特征在于,包括以下步骤:
步骤一、建立坐标系:以机翼根部为原点o,以机翼的刚心轴为x轴,在机翼平面内建立y轴;
步骤二、确定配重目标点P质量的大小mP、惯量的大小IPx及坐标(xP,yP):
其中,
mP=mq-ml-mkd……………………………[1]
IPx=Iqx-Ilx-Ikdx……………………………[2]
xP=(mqxq-mlxl-mkdxkd)/mP……………………………[3]
yP=(mqyq-mlyl-mkdykd)/mP……………………………[4]
其中,
mq是框段内机翼模型的总质量,
ml是框段内梁的质量,
mkd是框段内木框、蒙皮、加强肋部件的总质量,
Iqx是框段内机翼模型绕刚心轴x轴的总惯量,
Ilx是框段内梁绕刚心轴x轴的惯量,
Ikdx是框段内木框、蒙皮、加强肋部件绕刚心轴x轴的总惯量,
(xq,yq)是框段内机翼模型质心,
(xl,yl)是框段内梁质心,
(xkd,ykd)是框段内木框、蒙皮、加强肋部件的总质心,
步骤三、以配重目标点P为原点,以ξ轴为横坐标,以η为纵坐标,且ξ轴平行于x轴,η轴平行于y轴,建立局部坐标系;
步骤四、确定局部坐标系下的特征点的位置
其中i分别表示四个特征点G、H、I、J;
步骤五、确定过渡点A,B的坐标
式[6]要求
其中,
步骤六、根据杠杆原理确定过渡点A和B质量大小
步骤七、确定配重点C、D、E和F坐标;
由于点C、D位于线段GK、KH上,点E、F位于线段IL、LJ上,可进一步求得xC、xD、xE、xF;直线CE和DF平行于ox轴;
步骤八、根据杠杆原理确定配重点C、D、E和F质量大小,完成单框配重;
步骤九、重复上述步骤直至完成剩余框的配重。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210452011.4A CN102968526B (zh) | 2012-11-12 | 2012-11-12 | 一种大展弦比机翼颤振模型配重确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210452011.4A CN102968526B (zh) | 2012-11-12 | 2012-11-12 | 一种大展弦比机翼颤振模型配重确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102968526A CN102968526A (zh) | 2013-03-13 |
CN102968526B true CN102968526B (zh) | 2016-06-22 |
Family
ID=47798663
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210452011.4A Active CN102968526B (zh) | 2012-11-12 | 2012-11-12 | 一种大展弦比机翼颤振模型配重确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102968526B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104143014B (zh) * | 2014-04-17 | 2017-06-16 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种无独立配重空间筒状物颤振模型设计方法 |
CN105954000A (zh) * | 2016-05-31 | 2016-09-21 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞行颤振模型 |
CN106326664B (zh) * | 2016-08-30 | 2018-07-31 | 湖北航天技术研究院总体设计所 | 一种导弹分布机载载荷确定方法及系统 |
CN107665284B (zh) * | 2017-10-27 | 2021-03-23 | 北京航空航天大学 | 一种基于最小二乘法的飞机装配型架骨架优化布局方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN202177502U (zh) * | 2011-08-15 | 2012-03-28 | 中国航空工业集团公司西安飞机设计研究所 | 一种大展弦比机翼跨声速颤振模型 |
CN102661730A (zh) * | 2012-05-11 | 2012-09-12 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机高速颤振主梁模型空心梁截面尺寸的确定方法 |
-
2012
- 2012-11-12 CN CN201210452011.4A patent/CN102968526B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN202177502U (zh) * | 2011-08-15 | 2012-03-28 | 中国航空工业集团公司西安飞机设计研究所 | 一种大展弦比机翼跨声速颤振模型 |
CN102661730A (zh) * | 2012-05-11 | 2012-09-12 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机高速颤振主梁模型空心梁截面尺寸的确定方法 |
Non-Patent Citations (2)
Title |
---|
段斌修.《大型振动筛偏心配重质量优化设计探讨》.《武钢技术》.2011,第49卷(第6期),第53-56页. * |
罗务揆,霍应元.《弹性防护绳对颤振模型特性影响研究》.《第十一届全国空气弹性学术交流会》.2009,第107-110页. * |
Also Published As
Publication number | Publication date |
---|---|
CN102968526A (zh) | 2013-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107491616B (zh) | 一种适用于格栅构型舵面的结构有限元参数化建模方法 | |
CN104697462B (zh) | 一种基于中轴线的航空叶片型面特征参数提取方法 | |
CN102968526B (zh) | 一种大展弦比机翼颤振模型配重确定方法 | |
CN108052772A (zh) | 一种基于结构降阶模型的几何非线性静气动弹性分析方法 | |
CN109145471B (zh) | 基于cad与实测数据共融模型的虚拟装配系统及方法 | |
CN107391891A (zh) | 一种基于模型融合方法的大展弦比机翼优化设计方法 | |
CN106504329B (zh) | 一种基于牙体长轴的质点弹簧模型的牙龈变形仿真方法 | |
CN102867097B (zh) | 计及静弹性变形影响的光固化快速成型风洞模型设计方法 | |
CN105550388A (zh) | 钣金工艺的工艺模型生成方法及装置 | |
CN107972034A (zh) | 一种基于ros平台的复杂工件轨迹规划仿真系统 | |
CN110458174A (zh) | 一种无序点云关键特征点精确提取方法 | |
CN109726437B (zh) | 一种舱门气动载荷等效节点力处理方法 | |
CN106897501A (zh) | 面向自适应加工中基于叶片类零件变形的定位优化方法 | |
CN105109705A (zh) | 一种飞机翼面结构刚度计算方法 | |
CN112560365A (zh) | 一种基于全局映射变换的表面结构网格自动生成方法 | |
CN103760827A (zh) | 跃度约束的数控加工进给率的离线规划方法 | |
CN111159872A (zh) | 基于人机工程仿真分析的三维装配工艺示教方法及系统 | |
CN107145643B (zh) | 一种钣金雕塑型面钢架的智能设计系统及设计方法 | |
Lane et al. | A surface parameterization method for airfoil optimization and high lift 2D geometries utilizing the CST methodology | |
CN101276483B (zh) | 一种实现平移敏感的Laplacian网格编辑方法 | |
CN103942366B (zh) | 基于四段有理Bézier曲线表示的曲率连续的翼型及其生成方法 | |
CN108536932A (zh) | 基于互扭约束条件下的航空叶片积叠轴垂直度计算方法 | |
CN101763065A (zh) | 飞机复杂构件数控加工工装定位器自动设计方法 | |
CN106447781A (zh) | 一种基于Minkowski和面向自动装配的碰撞检测方法 | |
CN103043224A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |