本科毕业论文复合材料结构的概率损伤容限设计方法研究.docx
- 文档编号:2146796
- 上传时间:2023-05-02
- 格式:DOCX
- 页数:52
- 大小:513.54KB
本科毕业论文复合材料结构的概率损伤容限设计方法研究.docx
《本科毕业论文复合材料结构的概率损伤容限设计方法研究.docx》由会员分享,可在线阅读,更多相关《本科毕业论文复合材料结构的概率损伤容限设计方法研究.docx(52页珍藏版)》请在冰点文库上搜索。
本科毕业论文复合材料结构的概率损伤容限设计方法研究
编号
南京航空航天大学
毕业设计
题目复合材料结构的概率损伤容限设计方法研究
学生姓名
学号
学院
专业
班级
指导老师
二OO九年六月
南京航空航天大学
本科毕业设计(论文)诚信承诺书
本人郑重声明:
所呈交的毕业设计(论文)(题目:
复合材料结构的概率损伤容限设计方法研究)是本人在导师的指导下独立进行研究所取得的成果。
尽本人所知,除了毕业设计(论文)中特别加以标注引用的内容外,本毕业设计(论文)不包含任何其他个人或集体已经发表或撰写的成果作品。
作者签名:
时间:
2009年6月4日
复合材料结构的概率损伤容限设计方法研究
摘要
随着先进复合材料在飞机上的广泛应用,由于变量的随机性和安全系数的保守性,确定性方法已经无法满足设计要求,从而我们采用概率设计的方法,将结构安全程度定量化。
本文是建立在国内外的参考文献的数据基础之上,对四大概率设计方法之一的概率损伤容限设计方法(TsAGI方法)进行研究。
使用MATLAB语言编程模拟概率损伤容限设计,通过蒙特卡罗方法计算复合材料飞机结构的失效概率。
以Mig-29垂直尾翼为例进行详细的分析和编程计算。
接着在通过使用商业有限元软件ANSYS,对一个带孔的复合材料层合板,使用参数化设计语言(APDL)和概率设计模块(PDS),在随机输入变量的情况下,计算此板的失效概率。
其目的是为了说明用ANSYS有限元软件分析结构的可靠性是可行的。
最后还阐述了此方法的优缺点和以后的研究工作内容和方向。
关键词:
概率损伤容限设计,PDS,TsAGI方法,蒙特卡罗
Abstract
Asadvancedcompositematerialsarewidelyusedintheaerospacestructures,becauseoftherandomofvariablesandtheconservativeofsafefactor,thedeterministicdesignanalysiscannotmeetthedesignrequirements.Weadoptthemethodofprobabilitydesign,whichquantifiesstructuralsafetydegree.
Thisarticleisonthebaseofthedomesticandforeignrelateddata,oneoffourprobabilitymethodsisProbabilisticDesignofDamageTolerantmethod,thatisthemethodofTsAGI.Firstly,usingMATLABsimulatesprobabilisticdesignofdamagetolerant,andthefailureprobabilityofthecompositesofaircraftstructureiscalculatedwiththeMonteCarlomethod,forexampleMig-29fin.Secondly,throughtheuseofcommercialfiniteelementsoftwareANSYS,incompositeskinwithhole,theuseofparametricdesignlanguage(APDL)andprobabilisticdesignmodule(PDS)withrandominputvariables,thefailureprobabilityiscalculated.ItspurposeisthatitisfeasibleusingANSYSfiniteelementanalysissoftwaretoanalyzereliabilityofstructure.Finally,theadvantagesanddisadvantagesofthemethodareexpoundedandthecontentandthedirectionofresearchinfuturearegiven.
Keywords:
ProbabilisticDesignofDamageTolerant,PDS,themethodofTsAGI,MENTECARLO
第一章绪论
1.1研究背景
现代飞机由于其结构的复杂性以及飞行环境和载荷的随机性,出于安全和经济性的考虑,飞机结构需要具有耐久性、高可靠性、重量轻和低成本。
复合材料就符合了这些要求,从而被广泛地应用于现代飞机结构设计中。
但是,复合材料在制造过程中的特点引起复合材料结构的性能分布在一定范围内,从而导致了飞机结构设计内部变量的不确定性。
内部设计变量就包括了在制造过程中纤维和基体材料的性能、纤维的含量、铺层方向以及铺层的厚度。
设计变量还包括了结构尺寸、飞行的载荷情况以及飞行的环境。
如果仅由传统的安全系数来考虑结构设计是很危险的,安全系数可能太大,或某些情况下却太小,无法确定结构的可靠性,最严重极端的情况就是整个方法可能导致混合及低效的设计。
这样就需要概率设计方法来分析复合材料结构的可靠性,使安全程度定量化。
1.2国内外研究现状
由于复合材料在飞机结构上的使用和它的特殊性能,概率分析的方法在飞机结构的认证和设计中的优点很明显。
随着长期的发展,概率设计的方法已经逐渐成熟,并且开始运用于工程实践。
最近的几十年里,国外很多研究中心都致力于对概率分析方法的研究,还有一些公司将研究成果应用到了实际的飞机设计分析中,使得概率分析的方法在实践中有进一步的发展。
在国内从事概率设计的研究还是很少的,特别是有关复合材料的概率设计。
国外有关复合材料结构主要有四种概率设计方法,这四种方法分别是NASA刘易斯中心的IPACS方法,安全水平的方法,TsAGI的方法和NGCAD的方法。
各种方法有各自的优缺点。
我们要合理将它们应用在工程上。
本文主要研究的TsAGI方法,主要关于复合材料概率损伤容限方法。
这种方法是由俄罗斯航空联邦局下的航空流体力学研究中心,简称TsAGI,提出。
这是一种计算复合材料飞机结构可靠性的方法,复合材料飞机结构的损伤容限概率设计方法。
为设计者、工程师和分析人员提供了一个自动进行复合材料飞机结构的损伤容限可靠性分析的方法。
由于复合材料具有高的强度和刚度,对于在服役期间的循环载荷下裂纹的产生和扩展具有高的阻止能力。
经验表明在复合材料飞机结构产生灾难性破坏的最大的原因是在制造过程和服役期间的力学冲击损伤中产生的缺陷,这些缺陷目视不易发现。
因此就需要应用使用中的损伤容限准则预测复合材料结构在损伤发生后的情况。
1.3本文的研究工作
本文是建立在国内外的参考文献的数据基础之上,对四大概率设计方法的之一的概率损伤容限设计方法(TsAGI方法)进行研究。
使用MATLAB语言编程模拟概率损伤容限设计,通过蒙特卡罗方法计算复合材料飞机结构的失效概率。
接着在通过使用商业软件ANSYS,使用参数化设计语言(APDL)和概率设计模块(PDS),对一个带孔的复合材料层合板,考虑各个随机变量情况下,计算此板的失效概率。
最后还阐述了以后的研究工作内容和方向。
本文的主要内容如下:
第一章主要介绍了本文的研究背景,国内外关于复合材料概率设计的方法的研究现状,以及本文的研究工作内容。
第二章主要介绍了本文研究的概率方法TsAGI方法考虑主要随机变量,计算方法,还有损伤容限设计的整个模拟过程。
第三章主要对Mig-29飞机的主要随机变量的数据统计和处理。
得到各个变量的概率分布图和其他变量的样条曲线图。
第四章主要用MATLAB语言对整个损伤过程进行蒙特卡洛模拟。
绘制流程图并编写程序计算分析可靠度。
第五章主要介绍用有限元软件ANSYS中参数化设计语言APDL和概率设计模块PDS对带孔复合材料层合板进行可靠性分析。
第六章主要总结本文方法的优缺点和阐述了以后的研究方向。
第二章复合材料结构概率损伤容限设计方法(TsAGI方法)
2.1引言
俄罗斯航空联邦局下的航空流体力学研究中性,简称TsAGI,提出了一种计算复合材料飞机结构的可靠性的方法,复合材料飞机结构的损伤容限的概率设计方法。
为设计者、工程师和分析人员提供了一个自动进行复合材料飞机结构的损伤容限可靠性分析的方法。
由于复合材料具有高的强度和刚度,对于在服役期间的循环载荷下裂纹的产生和扩展具有高的阻止能力。
经验表明在复合材料飞机结构产生灾难性破坏的最大的原因是在制造过程和服役期间的力学冲击损伤中产生的缺陷,这些缺陷目视不易发现。
因此就需要应用使用中的损伤容限准则预测复合材料结构在损伤发生后的情况。
2.2主要的变量
这个问题的解决需要考虑到很多复合材料的损伤容限影响因素,包括一系列的工艺和服役因素、使用载荷、环境条件和材料的力学性质,这些都是随机变量。
这就需要在大量的数据下进行多参数的分析。
如图2-1所示,这些数据包括了:
(1)复合材料的力学性质;
(2)制造过程的类型和复合材料结构的设计;
(3)制造、测试和服役过程中非破坏性的检测结果;
(4)样本、组分和整体结构的经验数据。
要使损伤容限的设计有效,需要在整个过程中建立分析模型,分析模型中包括统计模型的结果和下面的参数:
(1)制造缺陷的类型和尺寸的统计分布;
(2)服役期间损伤的类型和尺寸的统计分布;
(3)典型的力学冲击下的损伤特点;
(4)估计制造缺陷和损伤对剩余强度和耐久性的影响;
(5)参数对剩余强度和耐久性的统计分布函数;
(6)设计条件、使用载荷、环境因子和它们的统计特性;
(7)检测和修复时间;
(8)估计失效概率的方法。
这些是结构安全性和可靠性预测的基础。
图2-1复合材料飞机结构概率损伤容限设计的主要输入变量
2.3计算方法
复合材料飞机结构概率损伤容限设计分析的方法利用了蒙特卡洛模拟,对应力和材料强度的分布考虑了材料在使用寿命内的性质、制造特性、操作结构损伤(包括修复引起的损伤)和操作环境(温度、吸湿、紫外线暴露)。
失效概率通过确定应力大于强度的概率计算出来。
在每一个时间间隔,将应力和强度相比较,这种模拟依赖于初始强度状态和操作过程中的随机变量。
结构载荷应用概率的方法模拟,确定在操作期间内的最大的载荷值。
应力和强度的比较确定局部结构的失效。
如果在N次载荷过程中,剩余强度的模拟发生了M次失效,失效概率就定义为:
(2-1)
根据现在的主要结构强度的实际情况,需要在以一系列的外在环境下进行强度的分析,就称为设计载荷情况。
有很多设计载荷情况,但是对于特定的结构就只考虑一些关键性的载荷情况。
则在时间t内结构的失效概率即为:
(2-2)
其中,
表示单块结构在
载荷情况下的失效概率。
选择的时间一般是服役期间,因此给出的是在飞机寿命期内的元件的失效概率。
可以从中得到一次飞行的失效概率。
损伤和未损伤结构的失效概率的计算方法不同。
因为,如果在所模拟的飞机寿命期内没有损伤产生,元件的强度没有退化,失效概率的公式即为:
(2-3)
是在结构位置在
载荷情况下,载荷承受能力的累积分布函数。
表示在时间t的最大载荷情况下的概率密度函数。
如果
与时间t无关,则
和
相独立。
在服役期间的缺陷的数量用泊松分布来模拟。
如果损伤的数量超越了一个,就应用蒙特卡洛模拟来确定这个位置上的失效概率。
如果没有损伤,应用近似的复合材料方法来分析。
2.4模拟过程
整个程序的模拟过程如图2-2所示:
图2-2模拟过程
a)对于每种设计载荷情况,每次迭代都开始于初始的强度值(
)。
结构的强度由很多设计变量确定,最后的强度值是一个概率统计的变量。
这些载荷情况一般是高斯分布或Weibull分布;
b)在相同的时间
下,产生了各种类型的缺陷,这些缺陷的数量和尺寸不同。
所有的设计载荷情况下的剩余强度
可以计算出来。
制造缺陷看成是在
时产生的损伤,每次的迭代开始时要随机选取一个初始的结构强度值;
c)产生的各种类型的操作损伤的数量;
d)如果在寿命期内没有制造缺陷和操作损伤,就产生了各种载荷情况下的最大的载荷。
在最大载荷出现的时候,结构的随机温度也产生。
剩余强度值依赖于温度和载荷。
如果载荷大于强度,结构失效就记录为M=M+1;
e)操作损伤的分布在结构寿命期内是相同的;
f)分析损伤尺寸的值所对应的合适的超越数曲线;
g)由于损伤的产生,使得结构的强度减弱,要计算强度减弱后的失效情况就需要计算在时间
的剩余强度。
每次的损伤,检测和修复的时间都是随机产生的,这里假设检测和修复的时间是相同的。
事实上,不同类型的损伤的检测间隔都是根据结构的构型,查找统一的检测间隔时间表得到的。
结构的检测时间是不同检测类型的检测时间中的最短的时间。
利用这一点可以确定时间常数
;
h)估算在时间间隔,
,所有载荷情况下的最大载荷。
在最大载荷产生的时候,结构的随机温度产生。
剩余强度值
依赖于温度和载荷。
如果载荷超过强度,结构失效就记录为M=M+1;
i)在时间
,由于结构的检测和修复,结构的强度值恢复到一定的水平,强度值的计算方法和步骤7相同。
重复上面的过程就可以计算出一次模拟过程中,结构在寿命期内的失效概率;
j)重复上面的迭代过程,就可以计算出多次模拟的结构的失效概率。
对于每个时间段内的失效概率的计算,都是通过计算每次飞行中的失效概率来确定的。
第三章数据处理
3.1引言
复合材料结构损伤损伤容限设计的变量有很多,这些变量主要是随机变量。
这些变量之间存在着一定的关系。
可以讲变量主要分为两大类:
一是和强度相关的;二是和应力相关的。
这些数据主要用于模拟失效概率的干涉模型。
对这些数据的处理有两种方法:
一是概率分布拟合;二是样条曲线拟合。
本章主要是用MATLAB语言对Mig-29飞机垂直尾翼的数据进行分析和拟合,为下面取随机数和建立数学模型做准备,最终为了计算其失效概率。
数据来自参考文献[5]。
3.2概率分布拟合
3.2.1载荷应力超越数以及分析拟合
在亚音速设计载荷发生的情况下,我们判断载荷大小的发生概率。
载荷大小从0-37的载荷应力是不发生的,因为载荷应力太小了。
根据超越数的定义,我们分别算出载荷数据的累积概率。
(3-1)
其中,P是某个载荷应力的累积概率,E是某个应力下的超越数,
是最大超越数。
亚音速设计载荷下的应力载荷超越数的数据来自参考文献[5],如表3-1所示。
处理后的累积概率如表3-2所示,对这些样本数据进行拟合分布曲线,通过描点法,从曲线形状可以看出偏移后的指数分布,MATLAB语言拟合后的结果如图3-1所示。
在亚音速设计载荷情况下,我们可以把在每个飞行小时最大超越数看作超出期望值,那样的话,我们就可以通过泊松分布产生每个飞行小时的载荷应力的发生次数。
表3-1亚音速设计载荷下应力载荷超越数
载荷应力MPa
每个飞行小时的超越数
每个寿命的超越数
0
54
153900
37
54
153900
55
12.4
35340
70
5.5
15675
85
1.7
4845
102
0.42
1197
120
0.057
162
137
0.02
57
154
0.004
11.4
表3-2亚音速设计载荷下应力载荷的累积概率
应力载荷MPa
累积概率P
0
0
37
0
55
0.7704
70
0.8981
85
0.9685
102
0.9922
120
0.9989
137
0.9996
154
0.9999
图3-1载荷应力累积概率分布
3.2.2温度超越数以及分析拟合
亚音速载荷情况下的温度超越数来自参考文献如表3-3所示。
按前面的处理方法,结果如表3-4所示。
对这些样本数据进行拟合分布曲线,通过描点法,从曲线形状可以看出是正态分布,MATLAB语言拟合后的结果如图3-2所示。
表3-3亚音速载荷情况下的温度超越数
温度
超越数
-73
1
-53
0.9998
-33
0.9974
-13
0.9772
7
0.8849
27
0.6554
47
0.3446
67
0.1151
87
0.02275
107
0.002555
119
5.191E-04
127
1.591E-04
147
5.417E-06
167
9.983E-08
187
9.901E-10
207
5.262E-12
表3-4亚音速载荷情况下的温度累积概率
温度
累积概率
-73
0
-53
0.0002
-33
0.0026
-13
0.0228
7
0.1151
27
0.3446
47
0.6554
67
0.8849
87
0.9772
107
0.9974
119
0.9995
127
0.9998
147
1.0000
167
1.0000
187
1.0000
207
1.0000
图3-2温度的累积概率分布
3.2.3损伤超越数以及分析拟合
(1)服役损伤
根据服役中垂直尾翼的检测数据,从而获得两类损伤超越数:
一是有关分层或刻痕的;二是有关孔+裂纹的。
刻痕的超越数的参考数据来自参考文献[5],如表3-5所示。
按前面的处理方法,刻痕累积概率如表3-6所示。
对这些样本数据进行拟合分布曲线,通过描点法,从曲线形状可以看出是偏移后的指数分布,MATLAB语言拟合后的刻痕累积概率分布如图3-3所示。
孔+裂纹的超越数如表3-7所示,处理后的孔+裂纹的累积概率如表3-8所示。
对这些样本数据进行拟合分布曲线,通过点描述法,从曲线形状可以看出是偏移后的指数分布,MATLAB语言拟合后的结果如图3-4所示。
表3-5服役损伤刻痕的超越数
尺寸
(mm)
每平方每小时
超越数10^-5
每0.31平方米每个寿命
超越数
0
10
0.093
25
10
0.093
30
7.8
0.07254
35
6
0.0558
40
4.6
0.04278
60
2.6
0.02418
100
0.65
0.00605
表3-6服役损伤刻痕的累积概率
刻痕尺寸
(mm)
每个寿命每0.31平方米
累积概率
0
0
25
0
30
0.2200
35
0.4000
40
0.5400
60
0.7400
100
0.9349
图3-3刻痕损伤的累积概率分布
表3-7服役损伤孔+裂纹的超越数
尺寸
(mm)
每平方米每小时
超越数10^-5
每0.31平方米每寿命
超越数
0
5
0.0465
10
5
0.0465
25
2
0.0186
30
1.3
0.01209
35
1
0.0093
40
0.8
0.00744
60
0.4
0.00372
100
0.1
0.00093
表3-8服役损伤孔+裂纹的累积概率
孔+裂纹
mm
每个寿命每0.31平方米
累积概率
0
0
10
0
25
0.6000
30
0.7400
35
0.8000
40
0.8400
60
0.9200
100
0.9800
图3-4服役损伤孔+裂纹的累积概率分布
整个过程中,损伤尺寸发生次数,我们可以把在每个寿命每0.31平方米中最大超越数看作超出期望值,那样的话,我们就可以通过泊松分布产生每个寿命每0.31平方米的发生次数。
(2)制造缺陷
根据制造中垂直尾翼的检测数据,从而获得两类损伤超越数:
一是有关分层或刻痕的;二是有关孔+裂纹的。
刻痕的超越数的参考数据来自参考文献[5],如表3-9所示。
按前面的处理方法,刻痕累积概率如表3-10所示。
对这些样本数据进行拟合分布曲线,通过描点法,从曲线形状可以看出是偏移后的指数分布,MATLAB语言拟合后的刻痕累积概率分布如图3-5所示。
孔+裂纹的超越数如表3-11所示,处理后的孔+裂纹的累积概率如表3-12所示。
对这些样本数据进行拟合分布曲线,通过点描述法,从曲线形状可以看出是偏移后的指数分布,MATLAB语言拟合后的结果如图3-6所示。
表3-9制造缺陷刻痕的超越数
尺寸
(mm)
每平方米每1000小时
超越数
每0.31平方米每寿命
超越数
0
0.027527
0.0256
37.5
0.014946
0.0139
75
0.005667
0.00527
表3-10制造缺陷刻痕的累积概率
刻痕尺寸
(mm)
每个寿命每0.31平方米
累积概率
0
0
37.5
0.457
75
0.7941
图3-5制造缺陷刻痕的累积概率分布
表3-11制造缺陷孔+裂纹的超越数
尺寸
(mm)
每平方米每1000小时
超越数
每0.31平方米每个寿命超越数
0
0.027527
0.0256
37.5
0.014946
0.0139
75
0.005667
0.00527
表3-12制造缺陷孔+裂纹的累积概率
损伤尺寸
mm
每个寿命每0.31平方米
累积概率
0
0
37.5
0.4554
75
0.7939
图3-6制造缺陷孔+裂纹的累积概率分布
初始损伤尺寸发生次数,我们可以把在每个寿命每0.31平方米中最大超越数看作超出期望值,那样的话,我们就可以通过泊松分布产生每个寿命每0.31平方米的发生次数。
3.3样条曲线拟合
3.3.1损伤表面的强度退化函数
考虑到冲击损伤的两种类型:
一是刻痕损伤类型;二是孔+裂纹的损伤类型。
在不同损伤类型的情况下,强度退化的函数是不一样的。
刻痕损伤尺寸对应的剩余强度数据来自参考文献[5],如表3-13所示。
经过样条曲线拟合后,刻痕损伤尺寸对应的剩余强度曲线如图3-7所示。
孔+裂纹损伤尺寸对应的剩余强度数据如表3-14所示。
经过样条曲线拟合后,孔+裂纹损伤尺寸对应的剩余强度曲线如图3-8所示。
表3-13刻痕损伤尺寸对应的剩余强度
刻痕尺寸(mm)
拉伸表面(MPa)
压缩表面(MPa)
0
470
-664
15
371
-311
20
359
-282
30
338
-244
40
319
-218
50
301
-199
60
285
-183
100
227
-144
125
193
-125
150
159
-108
图3-7刻痕损伤影响下的剩余强度
表3-14孔+裂纹损伤尺寸对应的剩余强度
孔+裂纹尺寸(mm)
拉伸载荷表面(MPa)
压缩载荷表面(MPa)
0
470
-664
5
282
-355
10
228
-287
15
196
-246
20
175
-220
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 本科毕业 论文 复合材料 结构 概率 损伤 容限 设计 方法 研究