CN103077286A - 一种飞机颤振模型的频率误差修正方法 - Google Patents
一种飞机颤振模型的频率误差修正方法 Download PDFInfo
- Publication number
- CN103077286A CN103077286A CN2013100193137A CN201310019313A CN103077286A CN 103077286 A CN103077286 A CN 103077286A CN 2013100193137 A CN2013100193137 A CN 2013100193137A CN 201310019313 A CN201310019313 A CN 201310019313A CN 103077286 A CN103077286 A CN 103077286A
- Authority
- CN
- China
- Prior art keywords
- model
- multiple material
- material model
- error correction
- composite material
- 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
Abstract
一种飞机颤振模型的频率误差修正方法,属于飞机风洞复合材料颤振试验模型的制造技术领域。其特征是该方法根据复合材料结构相似模型,以下称作复材模型,各主要固有频率与设计值的误差,通过有限元仿真分析与优化方法计算得到修正方案,并按照修正方案,采用二次成型工艺在复材模型蒙皮内侧的各个修正区域增加复合材料铺层,在改变复材模型刚度分布的同时,通过配重调整保持复材模型重量分布特征不变,从而实现对复材模型各主要频率误差的修正,使之处于设计允差范围内。本发明的效果和益处是该修正方法针对性强,工艺简便,可以提高复材模型制造的精确性,降低废品率,缩短制造周期。
Description
技术领域
本发明属于飞机风洞复合材料颤振试验模型的制造技术领域,具体涉及一种基于优化技术及复合材料二次成型技术的复合材料结构相似跨声速飞机颤振模型的频率误差修正方法。
背景技术
飞机研制过程中必须通过大量的风洞试验来得到飞机的气动弹性响应特征,研究和校核飞机结构的颤振性能,为结构设计提供参考依据。飞机颤振模型的精确制造是风洞试验能否成功的关键,所以必须保证颤振模型的各主要固有频率满足设计要求。
目前,飞机的风洞颤振试验最常用的是金属骨架模型,金属骨架模型是由维形泡沫、非承载蒙皮和金属骨架等结构构成,结构刚度由金属骨架提供。由于金属材料的材料特性稳定,因此这类模型仅需要控制金属骨架的几何加工精度,即可控制模型的频率误差。但金属模型与飞机的模态相似性低,不能很好满足跨声速风洞颤振试验需求。复合材料结构相似模型,以下称作复材模型,它是由变厚度蒙皮、变截面梁、框等典型构件组合,各构件排布以及承载特性与原机结构相似,是与原机结构具有部分相似特征的缩比模型,与飞机结构动力学特征的相似度高,是一种先进的模型形式。复材模型的各主要承力构件由树脂基纤维增强复合材料制成,其中增强纤维包括具有平面编织结构的玻璃纤维织物、碳纤维织物、芳纶纤维织物或上述三种纤维的混编织物或由玻璃纤维、碳纤维制成的短切纤维毡;基体树脂包括不饱和聚酯树脂、乙烯基树脂或环氧树脂等热固性树脂。但由于模型设计过程中工程假设带来的计算误差、复合材料的材料的分散误差以及模型制造中的几何误差与装配中的附加刚度误差等,造成制造出的模型模态频率容易超出模型设计的误差许可范围。
因而,当制造出的复材模型,其各主要固有频率低于设计值且超出误差许可范围时,需提供一种模态频率的制造误差修正方法,修正复材模型的模态频率使之满足设计要求。
发明内容
本发明的目的在于提供一种针对飞机复合材料结构相似跨声速颤振试验模型模态频率的制造误差修正方法,当制造出的复材模型其主要固有频率低于设计值且超出误差许可范围时,通过误差修正使制造出的复材模型的各主要固有频率精确地保证在设计值的误差范围内。
本发明的技术方案是:
一种飞机颤振模型的频率误差修正方法:通过测量复材模型几何尺寸、重量特征和各主要阶固有频率,计算得到实际材料性质,然后根据各主要固有频率与设计值的误差,通过有限元仿真分析与优化方法计算得到修正方案,并按照修正方案,采用二次成型工艺在复材模型蒙皮内侧的各个修正区域增加复合材料铺层,在改变复材模型刚度分布的同时,通过配重调整保持复材模型重量分布特征不变,最后按照设计装配制成合格的风洞试验模型,实现复材模型频率误差的修正。
具体内容是:
1.复材模型几何尺寸、重量特征和各主要固有频率的测量
使测量复材模型主要固有频率时,需要将复材模型蒙皮、梁和肋按照设计要求临时装配,临时装配时所用粘接剂为软化温度为50℃-80℃的热塑性树脂。测量完成后,将复材模型置于烘箱内加热至临时装配用热塑性树脂的软化温度,分离复材模型的蒙皮、梁和肋等构件,冷却后使用无水酒精和丙酮对粘接装配面进行清洗,去除热塑性树脂。
2.复材模型实际材料密度和弹性模量的计算
复材模型实际的材料性质是根据实际测量得到的复材模型几何参数、重量参数和固有频率,通过建立有限元分析模型计算得到,其中复材模型的实际材料性质包括蒙皮的材料密度和弹性模量梁的材料密度 和弹性模量肋的材料密度和弹性模量
3.根据复材模型各主要固有频率误差以及实际的几何、材料参数,优化得到误差修正方案
复材模型频率误差的修正方案包括需要增加复合材料铺层的修正区域以及各个修正区域上的复合材料铺层参数。修正区域在蒙皮内侧,并按照梁肋所处位置进行划分;复合材料铺层参数包括复合材料铺层的增强纤维类型、基体树脂类型、纤维树脂配比、纤维铺层厚度和纤维铺设角度,其中基体树脂类型与蒙皮的基体树脂类型相同,纤维类型和纤维树脂配比形成的复合材料的材料密度ρ'和弹性模量E'必须满足
通过得到的复材模型实际材料密度和弹性模量建立有限元分析模型,在此基础上展开优化,以最小化模型主要固有频率与设计频率的相对误差平方和W为优化目标,以模型蒙皮各个区域上待增加的修正铺层厚度和纤维铺设角度为优化设计变量,优化得到复材模型制造频率误差的修正方案。
优化数学模型如下:
Find Xi(di,θi) i=1,2,6,m
s.t.0≤di≤H i=1,2,6,m
0≤θi≤π i=1,2,6,m
其中i为修正区域编号,m为进行误差修正的铺层区域总数,j为进行误差修正的模型主要固有频率编号,N为进行误差修正的模型主要固有频率总数,di为第i修正区域上修正铺层厚度,θi为第i修正区域上修正铺层纤维铺设角度,H为允许最大修正铺层厚度。
4.根据误差修正方案,对复材模型实施频率误差修正
使用常温或中温固化的接触成形工艺,在蒙皮内侧按照修正方案增加用于误差修正的复合材料铺层。
5.根据修正增加的铺层分布调整复材模型重量分布
根据调整后的变厚度蒙皮和模型的梁肋等构件重量分布特征,在复材模型梁肋结合处增加配重,使模型整体的重心位置以及模型各个梁肋间区域的重心位置与模型设计要求相符。
6.按照设计装配制成合格的风洞试验模型
使用不饱和聚酯或环氧树脂等热固性树脂,按照复材模型设计要求,将调整过刚度分布、质量分布后的复材模型构件粘接装配,形成最终的复材模型成品。
综上所述,本发明根据制造出复材模型的频率误差,采用有限元分析与优化计算得到误差修正补偿方案,结合特定厚度与铺层方向的复合材料二次成型贴补工艺,实现对复合材料结构相似跨声速颤振模型频率的制造误差修正,使复材模型能够满足设计要求。
本发明的效果和益处是:该修正方法充分利用了复合材料的增材制造工艺,修正过程针对性强,工艺简便,可以提高模型制造的精确性,降低废品率,缩短制造周期。
附图说明
图1是某机翼整体结构示意图。
图中:1是机翼根部接头;2是机翼蒙皮。
图2是某机翼骨架结构示意图。
图中:3是翼肋;4是机翼梁。
图3是某机翼进行误差修正后蒙皮结构示意图。
图中:5复合材料铺层。
具体实施方式
以下结合技术方案和附图详细叙述本发明的具体实施方式。
图1是某机机翼复合材料结构相似跨声速颤振模型结构图,它由机翼根部接头1、上下机翼蒙皮2、多根竖向翼梁3和横向翼肋4组成。
图2是复材模型去除蒙皮后的骨架结构图,它由多根横向机翼肋板3和竖向机翼梁架4组成。
复材模型的机翼蒙皮2、翼梁3和翼肋4均由树脂基纤维增强复合材料制成。
使用游标卡尺和侧厚卡规测量模型各构件的长度和厚度,包括变厚度蒙皮2、变截面翼梁3和翼肋4,得到各构件的几何参数ai。
使用电子天平测量模型各构件的重量参数,包括变厚度蒙皮2、变截面翼梁3和翼肋4,得到各构件的质量mi。
采用软化温度为50℃-80℃的热塑性树脂将复材模型各构件按照设计要求临时粘接装配,包括变厚度蒙皮2、变截面翼梁3和翼肋4。按照设计方案将模型通过根部接头2夹持在测试平台上,利用锤击法对模型进行模态测试,获得复材模型的各主要固有频率ωi,其设计值是从而得到频率误差。
复材模型主要固有频率测量完成后,将复材模型置于烘箱内加热至临时装配用热塑性树脂的软化温度,分离复材模型的蒙皮2、翼梁3和翼肋4等构件,冷却后使用无水酒精和丙酮对粘接装配面进行清洗,去除热塑性树脂。
选择蒙皮内侧为修正区域,并按照梁肋所处位置进行划分,划分后区域如图3所示。
复合材料铺层参数包括复合材料铺层的增强纤维类型、基体树脂类型、纤维树脂配比、纤维铺层厚度和纤维铺设角度。
基体树脂为与蒙皮的基体树脂类型相同的热固性树脂;选纤维类型为碳纤维,使纤维类型和纤维树脂配比形成的复合材料的材料密度ρ′和弹性模量E′满足
通过得到的复材模型实际材料密度和弹性模量建立有限元分析模型,在此基础上展开优化,确定各修正区域待增加复合材料铺层的纤维铺设方向和铺层厚度。以最小化复材模型主要固有频率与设计频率的相对误差平方和W为优化目标,以复材模型蒙皮各区域上待增加的修正铺层厚度和纤维铺设角度为优化设计变量,优化得到复材模型制造频率误差的修正方案。
优化数学模型如下:
Find Xi(di,θi) i=1,2,6,m
s.t.0≤di≤H i=1,2,6,m
0≤θi≤π i=1,2,6,m
其中i为修正区域编号,m为进行误差修正的铺层区域总数,j为进行误差修正的模型主要固有频率编号,N为进行误差修正的模型主要固有频率总数,di为第i修正区域上修正铺层厚度,θi为第i修正区域上修正铺层纤维铺设角度,H为允许最大修正铺层厚度。
按照修正实施方案,在蒙皮2内表面进行贴补,对贴补区域,进行砂纸打磨和酒精清洗,使用常温或中温固化接触成形工艺,在蒙皮内侧的不同修正区域贴补相应铺层材料、相应铺层厚度和相应铺层角度的复合材料铺层5,修正完后的蒙皮如图3所示。
根据调整后的变厚度蒙皮2和模型的翼梁3和翼肋4等构件重量分布特征,在复材模型梁肋结合处增加配重,使模型整体的重心位置以及模型各个梁肋间区域的重心位置与模型设计要求相符。
使用不饱和聚酯或环氧树脂等热固性树脂,按照复材模型设计要求,将调整过刚度分布、质量分布后的复材模型构件粘接装配,形成最终的复材模型成品。
Claims (10)
1.一种飞机颤振模型的频率误差修正方法,其特征在于:复合材料结构相似跨声速颤振模型也称复材模型,通过测量复材模型几何尺寸、重量特征和各主要阶固有频率,计算得到实际材料性质,然后根据各主要固有频率与设计值的误差,通过有限元仿真分析与优化方法计算得到修正方案,并按照修正方案,采用二次成型工艺在复材模型蒙皮内侧的各个修正区域增加复合材料铺层,在改变复材模型刚度分布的同时,通过配重调整保持复材模型重量分布特征不变,最后按照设计装配制成合格的风洞试验模型,实现复材模型频率误差的修正。
2.根据权利要求1所述的一种飞机颤振模型的频率误差修正方法,其特征在于:测量复材模型主要固有频率时,需要将复材模型蒙皮、梁和肋按照设计要求临时装配,临时装配时所用粘接剂为软化温度为50℃-80℃的热塑性树脂。
3.根据权利要求1所述的一种飞机颤振模型的频率误差修正方法,其特征在于:复材模型主要固有频率测量完成后,需要将复材模型置于烘箱内加热至临时装配用热塑性树脂的软化温度,分离复材模型的蒙皮、梁和肋等构件,冷却后使用无水酒精和丙酮对粘接装配面进行清洗,去除热塑性树脂。
5.根据权利要求1所述的一种飞机颤振模型的频率误差修正方法,其特征在于:复材模型频率误差的修正方案包括需要增加复合材料铺层的修正区域以及各个修正区域上的复合材料铺层参数,其中修正区域为蒙皮内侧,并按照梁肋所处位置进行划分,铺层参数包括复合材料铺层的增强纤维类型、基体树脂类型、纤维树脂配比、纤维铺层厚度和纤维铺设角度。
7.根据权利要求1所述的一种飞机颤振模型的频率误差修正方法,其特征在于:修正方案的获取采用优化设计方法得到,通过得到的复材模型实际材料密度和弹性模量建立有限元分析模型,在此基础上展开优化,以最小化复材模型各主要固有频率相对误差平方和为优化目标,以各个修正区域中所选复合材料的铺层厚度与纤维铺设角度为优化设计变量。
8.根据权利要求1所述的一种飞机颤振模型的频率误差修正方法,其特征在于:使用常温或中温固化的接触成形工艺,在蒙皮内侧按照修正方案增加用于误差修正的复合材料铺层。
9.根据权利要求1所述的一种飞机颤振模型的频率误差修正方法,其特征在于:根据调整后的变厚度蒙皮和模型的梁肋等构件重量分布特征,在模型梁肋结合处增加配重,使模型整体的重心位置以及模型各个梁肋间区域的重心位置与模型设计要求相符。
10.根据权利要求1所述的一种飞机颤振模型的频率误差修正方法,其特征在于:使用不饱和聚酯或环氧树脂等热固性树脂,按照复材模型设计要求,将调整过刚度分布、质量分布后的复材模型构件粘接装配,形成最终的复材模型成品。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310019313.7A CN103077286B (zh) | 2013-01-18 | 2013-01-18 | 一种飞机颤振模型的频率误差修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310019313.7A CN103077286B (zh) | 2013-01-18 | 2013-01-18 | 一种飞机颤振模型的频率误差修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103077286A true CN103077286A (zh) | 2013-05-01 |
CN103077286B CN103077286B (zh) | 2016-01-13 |
Family
ID=48153815
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310019313.7A Active CN103077286B (zh) | 2013-01-18 | 2013-01-18 | 一种飞机颤振模型的频率误差修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103077286B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103342167A (zh) * | 2013-05-18 | 2013-10-09 | 大连理工大学 | 一种机翼复合材料缩比模型的制造方法 |
CN104143014A (zh) * | 2014-04-17 | 2014-11-12 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种无独立配重空间筒状物颤振模型设计方法 |
CN104182565A (zh) * | 2014-06-25 | 2014-12-03 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种二级轻气炮试验模型的设计方法 |
CN105630600A (zh) * | 2014-10-31 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种高精度的有限元模型质量特性自动调整方法 |
CN106096088A (zh) * | 2016-05-31 | 2016-11-09 | 中国航空工业集团公司西安飞机设计研究所 | 一种螺旋桨飞机螺旋颤振分析方法 |
CN107357992A (zh) * | 2017-07-13 | 2017-11-17 | 东南大学 | 基于聚类分析的复合材料结构有限元模型修正方法 |
CN111125829A (zh) * | 2019-12-04 | 2020-05-08 | 江西洪都航空工业集团有限责任公司 | 一种优化全动平尾静气动弹性和颤振的方法 |
CN111752454A (zh) * | 2020-06-23 | 2020-10-09 | 厦门美柚股份有限公司 | 应用程序界面的皮肤控制方法、装置、终端及介质 |
CN113639956A (zh) * | 2021-10-18 | 2021-11-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种模型倾角测量装置的校准装置及其校准方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101509825A (zh) * | 2009-03-13 | 2009-08-19 | 西安交通大学 | 飞机颤振风洞模型一体化设计及制造方法 |
CN202177502U (zh) * | 2011-08-15 | 2012-03-28 | 中国航空工业集团公司西安飞机设计研究所 | 一种大展弦比机翼跨声速颤振模型 |
-
2013
- 2013-01-18 CN CN201310019313.7A patent/CN103077286B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101509825A (zh) * | 2009-03-13 | 2009-08-19 | 西安交通大学 | 飞机颤振风洞模型一体化设计及制造方法 |
CN202177502U (zh) * | 2011-08-15 | 2012-03-28 | 中国航空工业集团公司西安飞机设计研究所 | 一种大展弦比机翼跨声速颤振模型 |
Non-Patent Citations (1)
Title |
---|
章怡宁,朴春雨: "满足固有频率和振型要求的符合材料翼面结构动力学设计", 《航空学报》, vol. 28, no. 2, 31 March 2007 (2007-03-31), pages 309 - 313 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103342167A (zh) * | 2013-05-18 | 2013-10-09 | 大连理工大学 | 一种机翼复合材料缩比模型的制造方法 |
CN103342167B (zh) * | 2013-05-18 | 2015-06-10 | 大连理工大学 | 一种机翼复合材料缩比模型的制造方法 |
CN104143014A (zh) * | 2014-04-17 | 2014-11-12 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种无独立配重空间筒状物颤振模型设计方法 |
CN104143014B (zh) * | 2014-04-17 | 2017-06-16 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种无独立配重空间筒状物颤振模型设计方法 |
CN104182565B (zh) * | 2014-06-25 | 2017-04-19 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种二级轻气炮试验模型的设计方法 |
CN104182565A (zh) * | 2014-06-25 | 2014-12-03 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种二级轻气炮试验模型的设计方法 |
CN105630600A (zh) * | 2014-10-31 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种高精度的有限元模型质量特性自动调整方法 |
CN106096088A (zh) * | 2016-05-31 | 2016-11-09 | 中国航空工业集团公司西安飞机设计研究所 | 一种螺旋桨飞机螺旋颤振分析方法 |
CN107357992A (zh) * | 2017-07-13 | 2017-11-17 | 东南大学 | 基于聚类分析的复合材料结构有限元模型修正方法 |
CN111125829A (zh) * | 2019-12-04 | 2020-05-08 | 江西洪都航空工业集团有限责任公司 | 一种优化全动平尾静气动弹性和颤振的方法 |
CN111125829B (zh) * | 2019-12-04 | 2022-05-06 | 江西洪都航空工业集团有限责任公司 | 一种优化全动平尾静气动弹性和颤振的方法 |
CN111752454A (zh) * | 2020-06-23 | 2020-10-09 | 厦门美柚股份有限公司 | 应用程序界面的皮肤控制方法、装置、终端及介质 |
CN113639956A (zh) * | 2021-10-18 | 2021-11-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种模型倾角测量装置的校准装置及其校准方法 |
CN113639956B (zh) * | 2021-10-18 | 2022-01-18 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种模型倾角测量装置的校准装置及其校准方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103077286B (zh) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103077286B (zh) | 一种飞机颤振模型的频率误差修正方法 | |
CN103342167B (zh) | 一种机翼复合材料缩比模型的制造方法 | |
CN101509825B (zh) | 飞机颤振风洞模型一体化设计及制造方法 | |
CN202177502U (zh) | 一种大展弦比机翼跨声速颤振模型 | |
US9011616B2 (en) | Optimizing the shape of a composite structure | |
CN102867097B (zh) | 计及静弹性变形影响的光固化快速成型风洞模型设计方法 | |
CN103646131A (zh) | 一种考虑气动弹性约束的复合材料机翼多目标优化设计方法 | |
CN107256320A (zh) | 一种基于复合材料铺层的小型手抛无人机的有限元分析方法 | |
CN105447269B (zh) | 一种飞机在定常过载状态下的非结构质量配平方法 | |
Brooks et al. | Undeflected common research model (uCRM): an aerostructural model for the study of high aspect ratio transport aircraft wings | |
CN112611538B (zh) | 一种针对跨音速颤振风洞模型加工的设计参数控制方法 | |
CN103890764A (zh) | 包括板层堆叠和树脂的复合部件的几何建模 | |
CN102567582A (zh) | 基于有限元分析的复合材料构件热压罐成型工装型面设计方法 | |
CN112611539B (zh) | 一种全机跨音速颤振风洞模型及其设计方法 | |
CN109543344A (zh) | 一种飞机门框区截面垂向弯曲刚度的计算方法 | |
CN112149322A (zh) | 一种复合材料加筋结构固化变形的有限元分析方法 | |
CN102354324B (zh) | 跨音速颤振模型复合材料单梁结构设计及其刚度计算方法 | |
CN106777693A (zh) | 一种基于catia数模的实时优化设计方法 | |
CN113111538B (zh) | 考虑湿热效应的复合材料加筋板临界剪切屈曲载荷的预测方法 | |
CN104598693A (zh) | 一种确定薄壁结构高刚度连接区载荷传递的方法 | |
CN109902323A (zh) | 一种带有薄膜蒙皮结构的机翼承载能力分析方法 | |
CN112528537A (zh) | 一种面向压缩稳定性的变刚度复合材料结构分析方法 | |
CN207292477U (zh) | 一种微桁架翼面颤振模型 | |
CN102722609B (zh) | 一种无铰式旋翼桨叶模型及其铺层设计方法 | |
Stanford | Aeroelastic analysis and optimization of membrane micro air vehicle wings |
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 |