飞机结构损伤容限分析

1第3章飞机强度计算方法飞机结构损伤容限分析21基本概念损伤容限指在规定的未经维修的使用阶段内，结构抵抗由于存在瑕疵、裂纹或损伤导致破坏的能力。损伤容限承认结构在使用前就带有初始缺陷，在使用中不可避免受到外来物的损伤，但必须把这些缺陷和损伤在规定的未修使用期内的增长控制在一定的范围内，使得裂纹不发生不稳定（快速）扩展，并在此期间，结构应满足规定的剩余强度要求，以满足飞机结构的安全性和可靠性。3临界值可检测性增加裂纹扩展可检门槛值裂纹检测周期飞行次数裂纹长度设计极限剩余强度降低剩余强度破损安全飞行次数应力待修理结构1基本概念41基本概念剩余强度：把含裂纹结构的承载能力称为该结构的“剩余强度”。很明显，剩余强度随着裂纹尺寸的增加而降低。裂纹扩展寿命：在设计应力谱作用下由初始裂纹a0扩展到临界裂纹acr时所达到的寿命。初姑缺陷：假定装配后飞机结构预存的缺陷尺寸，它刚小于无损检测的最大不可检缺陷尺寸。它的大小和特征受无损检测能力制约，无损检测能力提高了，初始缺陷尺寸a0就可定得小些。检查间隔：确保安全并由可检查度类别确定的两次之间的期限。51基本概念应当强调指出：损伤容限设计概念的引入，并不意味着对工艺质量要求降低，对缺陷和裂纹等不加注意。损伤容限设计是要求采取措施以保证含有一定损伤的结构在使用中的安全性。这些措施是彻底地、定期地对裂纹检查。而为了保证裂纹不会扩展到临界尺寸，首先要保证初始缺陷不会超过某一最大尺寸，这种检验要在飞机投入使用之前进行。6飞机在使用期间或制造初期允许出现裂纹（损伤）,甚至允许主要受力构件发生裂纹（并无危及结构安全）。利用断裂力学理论与实验结果,设计使得结构裂纹在一定限度内（损伤容限设计）,保证结构有足够的剩余强度、刚度（能继续承载）,利用定期的检测、维修保证飞机结构使用的安全可靠,而不致发生灾难事故。把握好含裂结构的裂纹的基本特性－－断裂力学1基本概念72断裂力学基础断裂力学概念以含裂体的特征参数(几何、载荷)表征其内力、形变规律，研究含裂体剩余强度规律及破坏准则,包括交变载荷作用下的裂纹演变规律及寿命估算分析等。8按裂纹的几何特征分类裂纹9按裂纹的力学特征分类1型(张开型):承受与裂纹面垂直的正应力，裂纹面位移沿y方向，裂纹张开。2型(滑开型):承受xy平面内的剪应力，裂纹面位移沿x方向，裂纹面沿x方向滑开。3型(撕开型):承受是在yz平面内的剪应力，裂纹面位移沿z方向，裂纹沿z方向撕开。裂纹10几何特征：一个方向的尺寸比另两个方向的尺寸小得多，例如等厚度平面薄板。受力特征：外力和约束，仅平行于板面作用，沿z方向不变化。应力特征：平面应力问题只有三个应力分量，仅为x、y的函数，与z无关。xyzy(,)xxxy(,)yyxy(,)xyyxxyxy平面应力11zxy几何特征：一个方向的尺寸比另两个方向的尺寸大得多，且沿长度方向几何形状和尺寸不变化。受力特征：外力和约束平行于横截面作用，沿长度z方向不变化。应变特征：如图选取坐标系，以任一横截面为xy面，任一纵线为z轴，则任一横截面均可视为对称面。平面应变问题只有三个应变分量，仅为x、y的函数，与z无关。(,)xxxy(,)yyxy(,)xyyxxyxy平面应变12线弹性断裂力学基础无穷大板中心一椭圆(a,b)孔边的应力集中1tKrtbaK或，时，，tKba10应力集中增大ba，0变成裂纹tK1.裂纹尖端应力应变场和应力强度因子131.裂纹尖端应力应变场和应力强度因子3cos1sinsin22223cos1sinsin22223sincossin2222xyxyKrKrKr1957年(IRWIN)2ijijKfrKa注意这里的是远场应力。为应力强度因子。14122(1)321coscos4222IKrukE122(1)321sinsin4222IKrvkE()xywdzE0w平面应变平面应力3431k平面应变平面应力1.裂纹尖端应力应变场和应力强度因子15裂端正前方：裂纹表面：0,22,0xyIyIxrKrK，0,00,xyyx，()()2IIIijijKfr()()IIiIiruKg裂纹尖端附近的应力和应变场，可以简写成如下形式1.裂纹尖端应力应变场和应力强度因子16()()2IIIijijKfr应力场公式的特点：1）在裂纹尖端，即r＝0处，应力趋于无穷大，应力在裂纹尖端出现奇异性。2）应力强度因子KI在裂纹尖端是有限量。3）裂纹尖端附近区域的应力分布是r和θ的函数，与无穷远处的应力和裂纹长度无关。1.裂纹尖端应力应变场和应力强度因子可以看出，用应力作为参量来建立如传统的强度条件失去了意义。但应力强度因子是有限量，它不代表某一点的应力，而是表征裂端应力应变场强度的参量。所以KⅠ可作为参量建立破坏条件是恰当的。强调：因KⅠ由线弹性理论推出，所以一般只适用于线弹性材料的断裂。由此建立起来的理论称为线弹性断裂力学。17应力强度因子可定义如下000000lim2lim2lim2yrxyrxzrKrKrKrⅠⅡⅢ一般可通过解析法、数值法或叠加法，典型的应力强度因子可以从手册中查到。1.裂纹尖端应力应变场和应力强度因子181.裂纹尖端应力应变场和应力强度因子例如，对于无限大板的对于有限尺寸的平板，则应进行修正，即KaⅠ()aKaWⅠKaⅠβ称为修正系数，或形状因子（函数），其它文献用Y或F表示。231()(1.770.277.05102.7)0.7ab2bW192.断裂韧度和断裂准则由前文分析可知：不管载荷类型、大小、方向和分布情况如何，也不管裂纹构件的几何形状和裂纹位置、长短与形状如何，只要是同一类型的裂纹，则裂纹尖端区域的应力和位移场的结构就完全相同。载荷情况反映在应力强度因子的大小上，带裂纹构件的几何形状以及裂纹的位置、长度、形状等也反映在应力强度因子上。如果两个构件的裂纹均属于I型裂纹，而且应力强度因子的数值相等，则裂纹尖端区域的应力和位移场就相同了。综上分析，我们可以想象到裂纹是否扩展与裂纹尖端区域的应力场强弱程度有关的，也就是与裂纹尖端应力强度因子大小有关。试验表明，对于一定厚度的平板，不论外加载荷、板件的几何形状和尺寸以及裂纹的情况如何，只要应力强度因子达到某一数值时，裂纹就失稳扩展，从而引起构件断裂。202.断裂韧度和断裂准则使裂纹发生失稳扩展的应力强度因子值，称为材料的断裂韧性（或断裂韧度），用Kc表示。断裂韧度代表含裂纹构件抵抗断裂破坏的能力。BKIC平面应变状态KC实验表明，材料断裂韧度随试件厚度B的增加而下降。这时因为薄板的裂纹尖端处于平面应力状态，裂纹不易扩展，其断裂韧度较高；随着板厚度的增加，裂纹尖端处于平面应变状态的部分增加，裂纹较易扩展，其断裂韧度降低。当板厚度增加到一定值以后，断裂韧度降到最低值，称为平面应变断裂韧度，用KIC表示。综上，I型裂纹失稳扩展的判据为,ICIICKKKK212.断裂韧度和断裂准则强调：虽然断裂韧度KIC是应力强度因子KI的临界值，两者存在密切的关系，但物理意义完全不同。KI是裂纹尖端应力场强弱程度的度量，与载荷、构件几何及裂纹情况有关；而断裂韧度KIC却是材料阻止裂纹失稳扩展能力的度量，只和材料的成分、热处理状态及加工工艺有关，它是一个材料常数。我们知道，只有当构件足够厚，裂纹尖端区域处于平面应变状态，材料的断裂韧度才是一个与厚度无关的常数，即平面应变的断裂韧度KIC，为了使裂纹尖端的塑性区尺寸远小于构件和裂纹尺寸，构件的韧带宽度W-a（W是构件的宽度）和裂纹尺寸均应足够大。实验表明，构件厚度、裂纹尺寸和韧带宽度之值均应满足下式22.5ICysKBaWa、、222.断裂韧度和断裂准则以上讨论的均为I型裂纹，对于II裂纹的断裂判据是否也为，尚缺乏充分实验数据的证实，因而没有被人们普遍接受。对于III型裂纹，其断裂判据也为IIIICKKIIIIIICKK至于复合型裂纹断裂判据，可参见相关书籍，这里只给出工程断裂判据，如下I-II型I-III型I-II-III型IIIICKKK22112IIIIICICKKKK22112IIIIIIICICKKKKK23IICKaK2.断裂韧度和断裂准则ICCKa21ICCKa临界应力或剩余强度临界（容限）裂纹尺寸,ICIICKKKK适用于脆性断裂和小范围屈服情况（需要修正）！实验证实，该判据对相对脆的材料或截面面积相当大的构件是合适的，许多高强刚、硬铝和高强度铝合金制成的构件。243.裂纹尖端塑性区前面曾提到，根据弹性解，在裂纹尖端处应力趋于无穷大，而实际上这是不可能的。事实上，当应力超过屈服应力时，必然在裂纹尖端邻近区域产生塑性变形，从而使裂纹尖端区的应力松弛，不可能达到无限大。应该说，材料一旦屈服，就不遵从弹性规律，故线弹性断裂力学不适用于屈服区。但如果屈服区很小(高强度、低中韧性材料即如此)则其周围的广大区域仍是弹性区，经过必要的修正后，线弹性断裂力学分析仍然有效。那么如何就塑性区影响进行修正呢？253.裂纹尖端塑性区裂纹尖端邻近塑性区的存在，自然要引起周围弹性应力分布的改变。在工程应用中，常引进一个被视为圆形区域的塑性范围。欧文为了避免弹塑性分析，建议将包含塑性区范围在内的有效半裂纹长度a+Ry作为假想裂纹，仍采用线弹性分析，并认为这一塑性修正对应力场而言是等效的。这种修正说明，塑性区的存在相当于裂纹长度的增加。Rya+Ryax弹性解σyσys等效模型263.裂纹尖端塑性区Rya+Ryax弹性解σyσys屈服范围对于面应力状态和平面应变状态(由Mises屈服准则求得)是不同的，欧文建议近似取为下式2222121,642IysyIyssKRKK平面应力状态平面应变状态这是因为，在平面应变状态下，沿厚度方向约束所产生的是拉应力,在三向拉应力下材料不易屈服而变脆。312()273.裂纹尖端塑性区在求得Ry后，即可计算等效裂纹长度a＝a＋Ry，然后计算等效应力强度因子KI，用等效裂纹长度代替原裂纹，有()IyKaRRya+Ryax弹性解σyσys28Rya+Ryax弹性解σyσys()IyKaR3.裂纹尖端塑性区2222121,642IysyIyssKRKK2212IysaK22142IysaK平面应力状态平面应变状态Ry是KI的函数，而KI是Ry的函数，所以严格上要用迭代的方法求得。欧文的塑性区模型只适用于小范围屈服情况，对大范围屈服，即裂纹尺寸和塑性区尺寸是同一量级，应采用弹塑性断裂力学方法。29IICKaKICCKa21ICCKa临界应力或剩余强度临界（容限）裂纹尺寸算例22142IysaK考虑形状修正和塑性修正30算例1板宽W＝100mm，板厚度B＝10mm，板中心具有穿透裂纹，其长度2a＝20mm。该板在远离裂纹处受拉应力作用，σ＝696MPa，板的材料为30CrMnSiNi2A，其σs＝1471MPa，KIC＝83.4MPa.m0.5，问该板是否会发生断裂？231()(1.770.2770.5102.7)解：2min2.58.04I