基于界面层模型的三维界面裂纹分析

合集下载

基于xfem的垂直于双材料界面的裂纹扩展问题

基于xfem的垂直于双材料界面的裂纹扩展问题

文章主题:基于xfem的垂直于双材料界面的裂纹扩展问题在材料科学和工程领域,裂纹扩展问题一直是一个备受关注的研究课题。

特别是在双材料界面上的裂纹扩展问题,由于双材料特性的不均匀性和复杂性,增加了研究和分析的难度。

在本文中,我们将从基于xfem(扩展有限元法)的角度出发,探讨垂直于双材料界面的裂纹扩展问题,以期为这一领域的研究和实际应用提供新的思路和方法。

一、概述垂直于双材料界面的裂纹扩展问题是指在两种材料的交界面上,裂纹在垂直方向上的扩展行为。

这种情况下,裂纹扩展的受力和受约束条件都受到了双材料特性的影响,需要深入分析和研究。

传统的有限元法在模拟和分析这种问题时存在一定的局限性,而xfem则能够有效地刻画裂纹的扩展路径和受力情况,因此成为了研究这一问题的有力工具。

二、裂纹模型的建立在进行垂直于双材料界面的裂纹扩展问题建模时,需要考虑双材料界面的影响,分析裂纹在材料间传播的受力情况和速度。

利用xfem,可以方便地将裂纹扩展路径和扩展速度等参数纳入模型中,通过数值计算得到裂纹扩展的演化规律和裂纹尖端的受力情况。

这有助于更准确地理解和分析垂直双材料界面上的裂纹扩展问题。

三、影响因素分析垂直于双材料界面的裂纹扩展受到诸多因素的影响,其中包括材料性质、裂纹尺寸、应力状态等。

通过xfem的数值模拟,可以分析不同因素对裂纹扩展行为的影响程度,揭示裂纹扩展过程中的关键因素和规律。

这有助于为材料设计和工程应用提供更可靠的参考依据。

四、工程应用与展望垂直于双材料界面的裂纹扩展问题在工程应用中具有重要意义,例如在复合材料结构的设计和评估中。

通过对裂纹扩展行为的深入研究和分析,可以为工程实践提供更有效的裂纹控制和结构安全性评估方法。

未来,基于xfem的裂纹扩展问题研究还可以结合人工智能算法和大数据分析等方法,进一步提高模型精度和计算效率,拓展应用领域和深化理论研究。

个人观点与总结垂直于双材料界面的裂纹扩展问题是一个复杂而又具有挑战性的研究领域,需要运用先进的数值模拟方法和理论分析手段来探讨和解决。

用界面元分析层状陶瓷的三点弯曲断裂性能

用界面元分析层状陶瓷的三点弯曲断裂性能

文章编号:100520930(2000)20420398207 中图分类号:TG 174.453 文献标识码:A ①用界面元分析层状陶瓷的三点弯曲断裂性能曾 东, 方岱宁(清华大学力学系,北京100084)摘要 采用无厚度界面单元对层状陶瓷的断裂性能进行了数值研究.在三点弯曲断裂试验中,可以观察到层状陶瓷裂纹扩展同时具有贯穿裂纹和层间裂纹两种形式.在本文的有限元计算中,对两种不同形式的裂纹扩展采用了不同的计算方法.对于纵向扩展的贯穿裂纹,采用最大拉应力破坏准则,而对于横向的层间裂纹扩展,则用无厚度的界面单元的破坏进行模拟.对层状陶瓷试件的层数与断裂性能关系的计算可以得到以下结论:随着层状陶瓷的层数的增加,其断裂性能逐步提高,但是增幅逐渐减少,到20层左右时,再增加层状陶瓷的层数,对提高断裂性能并无太大的帮助.关键词:层状陶瓷;三点弯曲;内聚力模型;界面单元;裂纹扩展众所周知,陶瓷是一种具有高刚度、高强度和高耐磨性的优质材料,但是其韧性较差,尤其是其脆性断裂,具有在裂纹扩展瞬间即导致材料整体的灾难性破坏这一致命弱点,从而大大限制了其应用的广泛性与可靠性.在1986年,美国材料研究室、空军W eigh t 研究室等多家单位的科学家由H .M .Burte 牵头,在全美复合材料会议上发文提倡将生物技术应用于材料的研究中,以推动新材料的研制.这是因为,在自然界中,存在着许多生物复合材料,如骨骼和木材,在许多方面,它们具有传统人工复合材料至今很难达到的程度.这其中,贝壳珍珠层的强韧化机理给予材料和力学科学工作者以很大的启示.珍珠层是贝壳类软体动物壳体最内部的一层,它具有保护贝壳不受侵害的作用.它的强度和韧度都非常好,而且是一种轻质材料.尤其重要的是,它的特殊结构对于裂纹扩展具有较大的抑制作用.珍珠层基本组成为钙盐、有机物和及极少量的水,其中钙盐以薄片状霰石的异构体形式存在,含量高达95%以上,中间夹杂有机物.在对贝壳珍珠层的断裂研究发现,珍珠层的断裂功比纯霰石大3000倍,可见有机物的含量虽然极小,但对增韧作用极大.Carry ,Jack s on [1]等人对贝壳珍珠层的细观结构和力学行为进行了开创性的研究,发现了珍珠层复杂的破坏机理.通过对贝壳珍珠层的分析,我们可以得到以下重要启示,由于将弱介质引入强介质中,使断裂过程发生质的变化,裂纹不仅在强介质中扩展,而且也在弱介质(水和有机物)内扩展,从而大大提高了断裂功.第8卷4期2000年12月 应用基础与工程科学学报JOU RNAL O F BA S I C SC IEN CE AND EN G I N EER I N G V ol .8,N o .4D ece m ber 2000①收稿日期:2000206227;修订日期:2000210223基金项目:国家自然科学基金重大资助项目(19891180)作者简介:曾东(1972—),男,博士研究生图1 层状陶瓷的裂纹扩展路径F ig .1 T he fracture surface of a lam inateceram ic tested in a th ree 2po in t test 基于以上的分析,有人提出了层状复合陶瓷的设计思路.1990年,英国剑桥大学的C legg [2]等人在N ature 杂志上发表文章说,他们首先通过设计,制造出了碳化硅 石墨叠层陶瓷.C legg 使用了一些厚度为0.2mm 的Si C 基片,将它们压迭起来,并采用石墨作为这些薄片的分割材料,从而构成一种具有简单的层状结构的复合材料.但是它所带来的优越的抗断裂、抗冲击性能却是以前的陶瓷基复合材料所无法比拟的.而且由于弱介质石墨的引入的比例较少,材料的整体强度降低不多.C legg 对层状陶瓷的断裂行为进行了三点弯曲试验观察,发现了如图1所示的断裂破坏形式和裂纹扩展路径,其压头载荷与位移的曲线如图2所示.图2 载荷与位移曲线F ig .2 T he l oad deflecti on behavi o r of a lam inated ceram ic tested in a th ree 2po in t bend test在图2中,我们可以比较清晰地看到,断裂裂纹几乎每扩展到一个石墨层时就会发生一次偏折,结果是增加了裂纹扩展的面积,并且消除了一次脆断,使得材料还可以承受载荷.而在传统的脆性陶瓷材料中,一次宏观的裂纹扩展就意味着材料的整体破坏.而且,其载荷2位移曲线也与普通的脆性陶瓷有很大差别.C legg 还估算了他所制备的Si C C 层状复合陶瓷的断裂韧性与断裂功.若按传统的计算方法,其断裂韧性K I C 可以达到15M Pa m 1 2远远高于传统的Si C 陶瓷块体材料(约3~4.5M Pa m 1 2).但由于层状复合陶瓷是混合断裂问题,而且是多裂纹扩展,因而用传统的计算方法所得的K I C 只是一个名义断裂韧性.但是材料韧性的增加是无疑的,而断裂功更可高达Si C 陶瓷料的十倍乃至几十倍以上,高达4625J m 2,这是其他的增韧方式所无法比拟的.由于研究开展的时间不长,因而对层状复合陶瓷材料特殊断裂性能的研究目前还只处于摸索阶段.开发材料性能,确定材料组成、结构、性能之间的关系,而这些材料学的问题又为力学研究提供了很多的研究方向.目前,在力学研究上,建立层状复合陶瓷断裂模型,建立界面力学模型研究界面应力对材料断裂行为及对材料整体性能的影响,对层状复993N o .4 曾东等: 用界面元分析层状陶瓷的三点弯曲断裂性能 合陶瓷的增韧机制作理论上的探讨与研究,为材料设计提出预测与指导,将是以后的研究重点,本文就是利用界面单元对层状陶瓷的断裂性能进行数值研究的初步尝试.1 界面单元本文计算中采用的界面单元是基于内聚力模型[3,4].内聚力实质上是物质原子或分子之间的相互作用力,复合材料的界面是由基体和增强相通过化学反应生成的一层结构复杂的物质.所以内聚力模型只是一种唯象简化,但通过适当地选取参数,可以反映出界面层物质的模量、强度、韧度等力学性质.界面层粘结力T 随其上下表面的相对位移∃而变化,其分量之间的关系,可以通过下式给定:T n =274Ρn (1-Κ)2∃n ∆n , T t =274Ρt (1-Κ)2∃t ∆t(1)Κ=∃m ax n ∆n 2+∃m ax t ∆t 21 2(2)下标n 、t 分别表示界面的法向和切向.其中∆n 为界面材料常数,表示可维系界面连接作用的最大法向分离距离.当∃n Ε∆n 时,T n =0,此时认为界面已被完全拉开,失去传递荷载的能力.Ρn 为另一界面材料常数,表示界面法向强度.相应的还有切向强度Ρt 及最大分离距离∆t .Κ是一个不可逆的量,取值范围从0至1,Κ=1时,对应界面完全脱粘,Κ可视为界面的损伤变量.界面法向负的相对位移不引起界面损伤.图3 归一化的法向内聚力模型F ig .3 Cohesive fo rce curve duringno r m al separati on 图3给出了法向的归一化内聚力曲线.随界面层被拉开的相对距离的增大,界面之间的粘结力经历了一个上升、到达最高点之后下降、最终为零的过程.曲线的这种软化段,可以模拟界面层的损伤过程,而且法向损伤和切向损伤是相互耦合的.只要给定法向和切向的界面强度Ρm ax 和最大张开位移∆,则上述的内聚力模型就完全确定了.图3曲线下的面积代表界面在正应力下破坏时所消耗的能量.利用以上的内聚力模型,可以在有限元中建立无厚度的界面元模拟界面的脱粘.图4 四节点界面单元F ig .4 In terface ele m en t 界面单元中以位移间断替代应变建立与应力的关系,可以得到无厚度的单元,是解决很薄夹层位移非连续问题的有效方法.本文以二维问题为例,给出了界面单元的推导.这里所用的界面单元,沿界面法向是初始无厚度的.图4所示为四节点界面单元,采用线性的位移插值函数.设单元节点位移向量为:a e =(u 1t ,u 1n ,…,u 4t ,u 4n )T (3)004 应用基础与工程科学学报 V o l .8单元上下表面的连续位移场为:u =(u -t ,u -n ,u +t ,u +n )T (4)式中上标“+”、“-”分别表示界面单元的上、下表面.线性位移插值函数为:Α1=12(1+Ν), Α2=12(1-Ν)(5)Ν为沿单元切向的局部坐标,界面相对位移场为:∃=(∃t ,∃n )T =(u +t -u -t ,u +n -u-n )T =B a e (6)其中:B =Α10Α20-Α20-Α100Α10Α20-Α20-Α1(7)则可以得到单元刚度矩阵为:K e =∫1-1B T D B d Ν(8)其中,这里的界面单元的切线刚度阵为:D =5T t 5∃t5T n 5∃n 5T n 5∃t 5T n 5∃n (9)当界面单元受挤压时,要人为加大界面单元的受压刚度,从而使得界面单元的嵌入变形不会过于显著.2 计算模型以下以Si 3N 4 BN 体系层状陶瓷为例进行三点弯曲断裂曲线的计算.为了节省计算量,取三点弯曲梁的一半作为计算模型,梁的半长是15mm ,高为3mm ,考虑梁处于平面应变状态.其中,硬层是各向同性的Si 3N 4陶瓷基片,软层为BN 夹层.在计算中,忽略了软层的厚度,布以界面元模拟裂纹在软层中的扩展,用到的Si 3N 4陶瓷基片的材料参数如表1所列.表1 Si 3N 4基体陶瓷的材料参数T able 1 Standard m aterial data used in the si m ulati on sE (GPa )ΜΡb (M Pa )3000.24700对于界面元所采用的内聚力模型的四个待定参数,通过实验测定是非常困难的.由于界面单元破坏实质是采用能量准则,所以在计算中,在保证界面分离能不变的前提下,一般选取法向和切向的模型参数一致,则待定参数减少至两个.T vergaard 和N eedle m an [3,4]在用界面元在计算纤维增强复合材料时,取界面的最大张开位移为纤维直径的1%,类似,在本文中可以取∆n 为硬层厚度的1%.从Ph ili p s [5,6]的论文可知,由BN 组成的界面层的能量释放率在10J m 2左右,所以,按照T vergaard 和N eedle m an[3,4]所提出的内聚力模104N o .4 曾东等: 用界面元分析层状陶瓷的三点弯曲断裂性能 型,有以下关系式:G I C =916Ρn ∆n ,这样就可以定下另一个界面参数Ρn 了.在层状陶瓷的三点弯曲断裂中,纵向的贯穿裂纹和水平的层间裂纹交替扩展,在计算中,沿高度方向的贯穿裂纹扩展采用最大拉应力破坏准则,而层间裂纹扩展则采用界面单元的破坏来模拟.对于各向同性、单一裂纹扩展的材料,一般采用断裂韧性K I C 作为断裂性能评估参数,而面对各向异性层状陶瓷的复杂断裂形态,采用断裂韧性作为材料断裂性能的评估参数是不太适宜的,这时通常采用载荷2位移曲线下的面积,即断裂功作为材料的断裂性能的评估参数.如果知道了断裂功W ,通过以下关系式,也可以确定材料的名义断裂韧性K C [2]:K C =EW1-Μ2(10)这里称其为名义断裂韧性,这是因为借用了 型裂纹扩展的结果,将实验或者计算得到的断裂功折算成为具有单一裂纹扩展性能的块体陶瓷所表现的断裂韧性,从而可以与基体陶瓷的断裂性能进行对比,比较直观地判定层状陶瓷断裂性能的改进.3 计算结果在计算中,分别对开缺口与不开缺口的层状陶瓷梁进行的三点弯曲断裂计算,若梁开口则缺口深度为梁高度的一半.当层状陶瓷的层数为1层时,块状基体陶瓷计算主要考虑层数对层状陶瓷试件断裂韧性的影响.图5 计算与试验的载荷位移曲线对比F ig .5 Computati onal and test l oad dis p lace m en t curves of the lam inated ceram icin th ree 2po in t bending test图5中实线表示计算得到的无缺口梁三点弯曲断裂载荷位移曲线,层状陶瓷的层数为15,虚线是选自Ph ili p s [6]论文中的实验曲线.可以看出,由于界面层的引入,层状陶瓷的断裂试验的载荷位移曲线呈现出阶梯状,从而提高了断裂功.对比两条曲线可以看出,由于界面层采用的是无厚度界面单元,没有将BN 软层对整体模量和强度的弱化考虑进来,所以计算曲线的斜率和峰值均大于实验曲线.由于层状陶瓷的制备具有许多不确定因素,无法做到象计算模型考虑的那样各层性质完全一致,因而其载荷位移曲线比之计算曲线缺乏规律性,不过在总体形貌上,二者还是具有相似性.曲线下的所围面积就是断裂功.从方程(10)可以看出,计算出断裂功,就可以得到层状陶瓷的名义断裂韧性,为了更方便的进行对比说明,以下的图均采用无量纲的数值,即W W 0和K C K C0,其中W 和K C 表示计算出的层状陶瓷的断裂功和名义断裂韧性,W 0和K C0是基体块状陶瓷的断裂功和断1裂韧性.以下是层数对断裂性能影响的计算结果.204 应用基础与工程科学学报 V o l .8从图6可以看出,层状陶瓷的断裂性能随着层数的增加是逐步提高的,但是当层数达到20层左右时,增幅并不显著.上图计算的是完整梁,即块体陶瓷试件和层状陶瓷实试件均不开口,断裂功增加倍数最大可达2倍左右.图6 完整梁的计算结果F ig .6 T he curve of the no r m alized toughen ingvs.the num ber of layers fo r the unnotched s peci m en 图7 开缺口梁的计算结果F ig .7 T he curve of the no r m alized toughen ing vs .the num ber of layers fo r the no tched s peci m en 图7显示的是开缺口梁的计算结果,基本的结论与图6一致,但断裂性能的增加倍数显著提高,这是由于对于块体陶瓷,开缺口与不开缺口,所计算出的断裂功相差甚大,而对于层状陶瓷,由于横向裂纹的存在,开缺口影响并不是很大.4 结论从本文的有限元计算中可以得到以下结论:1)层状陶瓷三点弯曲的载荷位移曲线具有不同于普通块体陶瓷的阶梯状,其中每一个台阶代表一层硬层的破坏和其后发生的裂纹在层间的扩展,正是由于层状陶瓷具有如此特殊的断裂形态,因而大大增加了断裂功.2)层状陶瓷的断裂性能相比于其基体块体陶瓷有较为显著的提高,考虑层数的影响,可以看出,随着层数的增加,层状陶瓷的断裂性能逐步提高,但是试件到了20层左右后,再增加层数,其断裂性能的改善并不显著.3)在对层状陶瓷断裂性能进行评估时,试件开口与不开口得到的评估效果相差很大.其主要的原因在于,开口与否对于基体块状陶瓷的断裂功计算影响很大,而对于层状陶瓷,其层间裂纹扩展会很快钝化开口裂纹,因而影响不大.4)另外,界面层的引进还可以增加陶瓷材料的柔度,计算中可以看出,对块状陶瓷进行三点弯曲,压头的最大位移略大于0.1mm ,而对于层状陶瓷,压头位移可以增大4~5倍.参考文献[1] Jacks on A P ,V incent J F V ,Turner R M .TheM echanical design of nacre [J ].P roc R Soc Lond ,1988,B 234:415~440[2] C legg W J ,Kendall K ,A lford N M .A si m p le w ay to m ake tough cera m ics [J ].N ature ,1990,347(10):455~457304N o .4 曾东等: 用界面元分析层状陶瓷的三点弯曲断裂性能 404 应用基础与工程科学学报 V o l.8[3] N eedle m an A.A continuum model for void nucleati on by inclusi on debonding[J].Journal of A pp lied M echanics,1987,54(9):525~531[4] Tvergaard V.Effect of fiber debonding in a w hisker2reinforced m etal[J].M aterials Science and Enginering A,1990,125:202~213[5] Philli p s A J,C legg W J,C lyne T W.F racture behavi or of cera m ic la m inates in bending( ).modeling of crackp ropagati on[J].A cta M etalM ater,1993,41(3):805~818[6] Philli p s A J,C legg W J,C lyne T W.F racture behavi or of cera m ic la m inates in bending( ).comparisi on ofmodelling w ith experi m ental data[J].A cta M etalM ater,1993,41(3):819~827Nu mer i cal Analysis of Fracture Behavi or ofLam i n ated Ceram i c i n Three-poi n t Bendi n gZEN G Dong, FAN G D ain ing(D epartm ent of Engineering M echanics,T singhua U niversity,Beijing100084)Abstract:This paper concerns the fracture behavi or of s peci m ens m ade up of cera m ic sheets.The num ericalanalysis of fracture behavi or of la m inated cera m ic in three2point bending has perfor m ed,based on boththrough2thickness crack p ropagati on w hen a critical stress is reached and interfacial crack advancing dictatedby the fracture of interface ele m ents.The essential conclusi ons related to fracture m echanis m can be dra w nfrom this num erical work.L ittle effect on the num ber of layer is p redicted,as l ong as the la m inated s peci m enhas more than about20layers.Keywords:la m inated cera m ic,fracture,interface ele m ent model,crack p ropagati on,three2point bending, finite ele m ent m ethod(FE M)。

基于粘结界面模型的三维裂纹扩展研究

基于粘结界面模型的三维裂纹扩展研究

1 粘 结 单 元
1 1 几 何 模 型 .
材料 的断 开是通 过粘结 单元来 描述 的 , 而粘 结单元 的张 开量是通 过上 下界 面 厂 的位移 变化 l 表示 l l l u 的. 图 1 如 所示 ,I I: 一u , 中 “ I “l 一其 表示 变形 后上下 界 面 的位 移. X ( =1 2 3 表示 三 维空 间 的笛 以 i , ,)
下 材 料 界 面 断裂 问题 的三 维 数 值 模 拟 .用 所 编 写 的 程 序 ( V E C F M) 分 别 对 I, Ⅱ, Ⅲ 型裂 纹 扩 展 问 题 进 行 数 值
模 拟 ,并 且 与 A a u . b q s 7计算 结果 进 行 对 比. 6
关 键 词 : 面 断 裂 ;粘 结 单 元 ;显 式 动 力 学算 法 界
中图 分 类 号 : 3 6 1 0 4. 文献 标 识 码 : A
O 引 言
工 程 实 际 中 的许 多 结 构 往 往 存 在 材 料 界 面 , 不 同 程 度 地 影 响 结 构 的整 体 力 学 行 为 , 是 造 成 事 故 发 生 会 这
的重要原 因之一 . 现有 的用 于解决 断裂力 学 问题 的数值 计算 方 法主要 包 括扩 展有 限元法 、 界元 法 、 网格 边 无 法 、 值 流 行 法 、 波 数 值 法 等. 中扩 展 有 限 元 法 可 以避 免 重 新 划 分 网格 , 需 要 相 当复 杂 的 插 值 函 数 小 其 却 数 ; 边界 元法必 须求 问题 的基本解 , 对于 非线性 问题 , 求解 非 常困难 ; 网格 法 中的 E G法在研究 断裂 无 F ; 波 数值 法 目前 用 于 断裂 力 学 问 题研 究 还 处 于初 级 阶 小

三维断裂分析软件FRANC3D

三维断裂分析软件FRANC3D
万方数据
应变方程确定这些点的扩展方向,它们的扩展尺寸 则由各点相对应力强度因子的大小和用户确定的 最大扩展尺寸共同决定。这些扩展后的点连线后重 新组成新的裂纹面。用户可以重新分析扩展后裂纹 前缘的应力强度因子分布情况。FRANC3D还有一 个优点就是,在裂纹扩展后,只需对很小的一部分 区域重新划分网格,这大大减少了网格划分时间。
断裂分析器FRANC3D(图2)是整个系统的核 心部分。它具有以下功能: 2.2.1 几何模型完善
对于从OSM读入的几何模型,FRANC3D可 以对之进行修改,比如增加点、线、面等。同时,
第6期
贾学明,等:三维断裂分析软件FRANC3D
图1立体对象建模器0sM示意图 Fig.1 Sketch map of 0SM
蕊万K翥7一一丽丽1[lsSi1n旷z∥口十+I㈡了Jz cos2 0]J{ (1)
式中0为决定裂纹前缘点位置的方位角,E(Ⅳ)为 第二类椭圆积分
第6期
贾学明,等:三维断裂分析软件FRANC3D
E(Ⅳ)一I /1一Ⅳ2 sin2 0d0, 舻一L—F∑ (2)
J0
C。
对于给定的a和C值,E(托)一1.211096。
万方数据
在划分好网格以后,把结果写入.bes文件,就 可以进行求解了。求解使用边界元系统BES系统进 行。边界元计算是整个分析中最耗时的。随着单元数 的增加,求解时间也显著增加。用户可以根据需要选 择所需的求解器(QR分解,Gauss消元法、迭代 法),用户可以根据问题的规模选择使用内存(in core)或者外存(out of core)来计算;甚至可以使用 BES的并行计算程序在多台机器上同时计算。 2.2.3 应力强度因子计算和裂纹自动扩展
大的点裂纹的扩展速度相应的也要大。最终的扩展 趋势是裂纹各点的应力强度因子相等。对于无限大 体内埋裂纹,只有圆形裂纹的各点应力强度因子相 等。 3.2 CCNBD试样的应力强度因子计算

各向异性与正交异性双材料Ⅲ型界面裂纹分析

各向异性与正交异性双材料Ⅲ型界面裂纹分析

+2 ( Q 4 5 ) -
o y
+(
dy
_0
( 1 )
基 金项 目: 山西省 自然科学基金 ( 2 0 1 1 0 1 1 0 2 1— 3 ) ; 山西省高校科技项 目( 2 0 1 1 1 0 9 3 ) 作 者简 介 : 冀亚仙 ( 1 9 8 8一) 女, 硕士研究生 , 主要研 究方向为偏微分方程理论及应用 。
者的重视 。由于结合材料 的破坏多发 生在界面附
近, 因此 , 界 面 裂 纹 的 分 析 对 结 合 材 料 的 强度 和 可 靠性评 价 具有 重要 意义 。文 献 [ 1 ] 对 各 向同性 双 材 料界 面裂 纹进 行研 究 , 得 到 I型 和 Ⅱ型裂 纹 尖 端应 力 具有 振 荡奇 异 性 。文 献 [ 2] 采 用 函数 变 量 方 法 , 对 含有 边缘 界 面裂 纹 的 不 同正 交 各 向异 性 平 板 在

( 2
1. 2)
( 4)
其 中
= +s i Y =r ( c o s 0+s i s i n 0 ) ,
第3 4 卷
第1 期








V o 1 . 3 4 N o . 1
F e b . 2 0 1 3
2 0 1 3年 2月
J O U R N A L O F T A I Y U A N U N I V E R S I T Y O F S C I E N C E A N D T E C HN O L O G Y
随着复合材料 的工业应用范 围不 断扩大 , 由不
同材 料组 成 的界 面 力 学 行 为 越 来 越 受 到 国 内外 学

任意三维裂纹扩展分析

任意三维裂纹扩展分析

任意三维疲劳裂纹扩展分析1.前言在工程实际中,真实的构造总是存在众多缺陷或裂纹,对于一个含裂纹或缺陷的构件,多在其服役荷载远低于容许强度的情况下就发生了破坏。

实际工程构造在经受长时间多因素综合作用下,产生变形、裂纹等缺陷,从而导致整个构造的失效。

构造的失效主要由疲劳引起,其最终失效形式即为断裂,有大约80%以上的工程构造的断裂与疲劳有关,由疲劳引起的巨大经济损失及灾难性的后果不胜枚举。

我们通常不能仅仅因为某个构件出现了裂纹就简单的认为该构件不平安或不可靠,尤其是对于大型设备的重要构件,因为这将使企业消耗高昂的本钱。

对于出现的裂纹,以往多采用以下几种处理方法:一是对出现裂纹的构件进展更换,这对于含裂纹但仍能工作的构件是一个巨大的浪费。

二是强行停顿使用进展维修,这样会带来巨大的经济损失;三是冒险继续使用,但这样会带来巨大风险,甚至会造成人员伤亡。

所以,人们更想知道,出现的裂纹是否会在既定载荷〔包括疲劳载荷在内的任意载荷〕下扩展成不平安或失效的临界尺寸,因此,出现了疲劳裂纹扩展分析。

疲劳裂纹扩展分析是采用断裂力学的理论和方法对含裂纹等缺陷构件的失效过程进展分析,以评估产品的平安性和可靠性,可以进展损伤容限评估和剩余寿命预测等,已经在化工机械、飞行器、核工业等各个工程领域得到了广泛应用,并得到了世界各国政府及学术机构的重视。

2.疲劳裂纹扩展分析软件在工程实践中,疲劳裂纹扩展分析已成为评估产品性能、改进产品设计和提高服役寿命的一个重要工具。

目前,疲劳裂纹扩展分析主要有解析法和数值法这样两种方法,下面分别介绍这两种方法。

1〕解析法解析法主要依据相应的标准和经历公式,将复杂的三维问题简化为二维问题,并对复杂的裂纹形状和荷载状态进展简化,然后用经历的方法对裂纹平安性进展评估。

但对于大量构造复杂的工程实际问题却无能为力,况且其简化后的分析准确度及是否真实逼近服役情况也值得探讨。

目前,工程上有几款基于解析法而开发的裂纹扩展分析软件,它们主要应用于航空标准构造的裂纹扩展分析,包括DARWIN、NASGRO、AFGROW等。

叠层板基体开裂分析的三维半解析有限元法

叠层板基体开裂分析的三维半解析有限元法

( ) 忽略 开 裂 对 位 移 的 影 响 。这 样 ,层 板 3 内的位移场 为 :
M =M , ,) ( Yz 1

界 积分方 程—— 混合边 界元 法 ,该 方 法具有 很高 的 精度 ,但计算 过 程 中涉 及到边 界 积分方 程 的离散 和 超奇 异积分 的计算 ,所 以计算 工 作量较 大 ,过程 复
叠层板 由于具有 比强度 高,比刚度大等优点 , 因而在工程中获得 了广泛的应用。然而在制造和使
用过 程 中 ,此类结 构经 常会产 生种 种缺 陷 ,如基体 开裂 ,基 体 和 纤 维 界 面 的脱 离 ,纤 维 的断 裂 和 拔
杂。秦太验,汤任基 等 对无 限大或有 限体 三维 裂纹问题进行了数值分析 ,他们求解的问题限于内 埋裂纹问题 。Pn n1 a ̄ 【 引利用交替法研究 了三维裂
用位 移外 推 法 得 到 了计 算 三 维 问 题 SF的 计 算 公 I
适合于任意形状 的横 向基体 开裂问题 ,应用 范围 广。与 A s 软件数值结果的比较表 明,本文方法 ny s 具有较高的精度 。
l 空间半解析有 限元
如 图 1所示处 于单 向拉 伸状 态 ( 0 层 方 向 ) 沿 。 的正 交叠层 板 ,现假 设某 9 。 内出现 一段 平 行 于 0层
t :t J £ £ ) , , (
(, } ,) ,
() 1
收 稿 日期 :20 - —0 0 60 1 4
基金项 目:国家 自然科学基金资助项 目 ( 07 14 ;广东省 自然科学基金资助项 目 (3 5 2 16 29 ) 0 15 )
作者简介 :富明慧 ( 96年生 ) 16 ,男 ,教授 ,博士生导师 ;E m i tmh alss.d .a — a :s f @m i yu eu c l s .

固体推进剂-衬层界面脱粘裂纹的三区域界面层模型

固体推进剂-衬层界面脱粘裂纹的三区域界面层模型

固体推进剂-衬层界面脱粘裂纹的三区域界面层模型摘要:本文提出了一种三区域界面层模型,用于解释固体推进剂衬层界面脱粘裂纹的来源和形成机制。

关键词:固体推进剂,衬层界面,脱粘裂纹,三区域界面层模型。

正文:本文首先介绍了固体推进剂衬层界面脱粘裂纹的性质及其形成的机理,然后提出了一种新颖的三区域界面层模型,用以描述固体推进剂衬层界面脱粘裂纹的来源和形成机制。

此外,该模型还被应用到现实情况中,以说明脱粘裂纹的形成机理。

最后,该文提出了一些建议,以帮助解决该问题并改善固体推进剂性能。

三区域界面模型可以有效地解释固体推进剂衬层界面脱粘裂纹的形成机制。

首先,该模型提出了两个重点:首先,衬层界面以及衬层内部的内层-外层界面;其次,衬层的表面的三区域,即表面受到推力的区域(A区)、单元之间的焊接区域(B区)以及上游外侧和下游外侧的焊接区域(C区)。

这三个区域都存在不同的应力分布,由于衬层也可能存在材料弱点,使得这些弱点受到应力超过其承受能力,因此导致衬层材料发生脱粘裂纹。

此外,该模型可以用来提出针对脱粘裂纹的改善方案。

根据该模型,可以改变衬层材料的组成成分以改善耐受负载,并增加焊接点的长度以减小附近的应力集中,从而降低脱粘裂纹的数量和尺寸。

另外,还可以调整固体推进剂的结构,将推力的作用位置转移到B区或C区,使推力的分布更加均匀,减少推力与焊接点之间的应力集中。

当然,也可以对衬层材料进行必要的改性,以增加材料的可塑性和强度,以减少脱粘裂纹的发生。

在使用三区域界面模型改善固体推进剂衬层界面脱粘裂纹问题时,还需要注意其他因素。

例如,如果衬层材料的温度太高,可能导致固体推进剂内部物理性质的变化,进而影响衬层材料的机械性能。

此外,其它影响脱粘裂纹的因素还包括衬层材料的成分及其在表面的厚度,以及衬层材料的受力方式等。

因此,在应用三区域界面模型时,还要根据实际情况考虑其他影响因素,识别各个影响因素对脱粘裂纹产生的影响,并采取相应的措施。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
DOI: 10.12677/ijm.2018.72005 39 力学研究
王硕,张铮
a cos ϕ = + r cos θ cos ϕ x 1 − e 2 sin 2 ϕ a 2 − b2 ,其中 e 2 = 为椭圆离心率。 y = r sin θ a2 2 b sin ϕ z = + r cos θ sin ϕ a 1 − e 2 sin 2 ϕ
International Journal of Mechanics Research 力学研究, 2018, 7(2), 36-43 Published Online June 2018 in Hans. /journal/ijm https:///10.12677/ijm.2018.72005
2. 界面层理论
Williams [2]利用均质材料裂纹的理论解, 对界面裂纹尖端场进行研究, 由于界面两侧材料性质不同, 发现直接应用断裂力学的理论模型存在振荡奇异性等不合理的情况,即在理论上裂纹尖端上下表面互相 嵌入,如图 1 所示,这显然不符合物理常识。 据此, Erdogan 认为界面应当是具有细观结构的界面层, 无厚度的理想界面不能够体现出界面的性质。 经过实验测量与分析,他提出界面层是厚度一致的非均匀各向同性材料,其弹性模量沿厚度连续变化, 而泊松比则可以视为常数,模型如图 2 所示,结合材料中①区域和②区域的弹性模量为常数,而在两区 域之间材料属性随着纵坐标的变化而变化,这里,E(y)表示界面层的弹性模量。这样的模型与均质材料 裂纹模型相似,有效地消除了应力振荡的现象。
DOI: 10.12677/ijm.2018.72005 38 力学研究
王硕,张铮
材料1
z O h x
y
材料2
Figure 3. Diagram of three-dimensional interface layer model 图 3. 三维界面层模型图
z 裂纹前缘
φ O P0 b
a
x P0 n P t 切线
Open Access
1. 引言
随着各种新型功能性材料的日益发展,工程领域中对于轻质、强度、刚度的要求也逐渐提高。在这 些新型材料中,我们把两种不同的材料利用粘接、焊接等技术连接在一起所使用的结构称为结合材料。 界面问题在结合材料中广泛存在,因此对于结合材料界面的力学性能研究也需要更加深入[1]。目前关于 均质材料裂纹的研究已经趋于成熟[2], 但是如果将经典的裂纹模型引入界面问题中会出现两种特殊现象: 界面的应力奇异性和界面裂纹尖端的应力振荡奇异性, 因此许多学者从不同的角度出发提出了改进模型, 顺利地解决经典模型存在的问题, 比较经典的有 Comninou 的接触区理论、 Erdoga 的界面层理论以及 Rice 的小范围屈服理论[3],已经被广泛接受。 界面层理论的最大优势就在于,符合物理常识并且能够消除界面裂纹的应力振荡奇异性,该理论认 为,结合材料的界面是实际存在的非均匀各向同性材料,其厚度均匀并且很薄,弹性模量沿厚度连续变 化,与功能梯度材料(FGM)的性质契合[4] [5]。虽然到目前为止关于界面的研究很多,但主要还是以平面 应力和平面应变问题为研究对象[6] [7] [8],对于一般性的三维问题,相关的研究还远远不足。本文将界 面层理论扩展至三维问题,建立基于功能梯度材料的三维模型并建立相应的理论分析体系。
(4)
法平面
P0
θ
n
r y
P
Figure 5. Local coordinate system in the front of crack 图 5. 裂纹前缘局部坐标系 DOI: 10.12677/ijm.2018.72005 40 力学研究
x2 z 2 + ≤ 1, y = 0 a 2 b2
椭圆边界上的点可以用参数方程表示为(如图 4 所示)
(1)
= = x a cos φ , z b sin φ
(2)
为得到较简单的平衡方程并对裂纹及其邻域进行描述,过边界上一点 P0 (参看图 4 左图)建立局部坐 标系 P0tny (参看图 4 右图),P0 点为局部坐标原点,P0n 为 P0 点在裂纹扩展平面内的法线,与 Ox 轴(即椭 圆长轴)夹角为 φ,P0t 为 P0 点在裂纹扩展平面内的切线;在 P0 点的法平面 nP0y (即垂直于切线 P0t 的平


在工程领域中,界面问题广泛存在,界面的特殊性质导致结合材料的破坏往往始于该处,所以为了保证
文章引用: 王硕, 张铮. 基于界面层模型的三维界面裂纹分析[J]. 力学研究, 2018, 7(2): 36-43. DOI: 10.12677/ijm.2018.72005
王硕,张铮
整个材料的使用安全必须对界面裂纹进行研究。与均质材料裂纹不同,界面裂纹表现出了一些特殊的性 质,如裂纹尖端场应力的振荡奇异性、应力强度因子的量纲等问题,不少学者也提出了改进原有不足的 新理论。考虑到界面裂纹具有的某些特殊性质,本文依据界面层理论,将结合材料的界面视为具有厚度 的功能梯度材料(FGM)层,并建立针对界面层内椭圆裂纹的断裂力学理论体系进行分析,对三维界面裂 纹破坏机理的研究提供了理论方法。
4.1. 三维界面裂纹的弹性理论方程
如上所述,在曲线裂纹前缘的每一点,本文建立相应的局部柱坐标系,用于描述裂纹前缘的力学分 析体系,进而建立各种力学基本方程。 首先,结合柱坐标系下的平衡方程,可得到该局部坐标系下的平衡方程(无体力):
∂σ r ∂σ θ r σ r − σ θ 1 ∂σ rϕ + + + 0 + (σ r − σ ϕ ) cos θ − σ θ r sin θ = r ∂θ r Φ ∂ϕ ∂r ∂σ rθ ∂σ θ 2σ θ r 1 ∂σ θϕ + + + 0 + σ rθ cos θ − (σ ϕ − σ θ ) sin θ = r ∂θ r Φ ∂ϕ ∂r ∂σ σ ∂σ ∂σ ϕ rϕ + θϕ + rϕ + 1 0 + 2 (σ rϕ cos θ − σ θϕ sin θ ) = r ∂θ r Φ ∂ϕ ∂r
θ
r t y
n 法线
Figure 4. Coordinate system in three-dimensional model 图 4. 三维模型的坐标系
面)内引入极坐标系(r, θ),即 rP0θ 为过 P0 点的法平面上的极坐标系(这里,y 轴垂直于 P0 点的切线 P0t 和 法线 P0n)。 于是,可以利用上述平面坐标系与参数 φ (即坐标系(φ, r, θ))来描述空间内任一点,其中,整体坐标 系与局部坐标系的变换关系为
Figure 1. Oscillation singularity of interfacial crack stress field 图 1. 界面裂纹应力场的振荡奇异性
Figure 2. Diagram of interface layer model [4] 图 2. 界面层模型示意图[4]
依据界面层理论的相关内容,本文将现有的界面层理论推广到三维界面层裂纹(如图 3 所示)。 实际中常见的三维裂纹可大致认为是椭圆形。设在无限大材料中有一椭圆形深埋片状裂纹(位于界面 层内),基于该椭圆裂纹的中心建立坐标系(参看图 3),其方程为
th th th
Abstract
In the field of engineering, interface problems exist extensively, and the special nature of the interface leads to the destruction of bonded materials that often begins at this location. Therefore, in order to ensure the safety of the entire material, the interface crack must be studied. Different from the homogeneous material cracks, the interface cracks show some special properties, such as the oscillation singularity of the field stress at the crack tip, and the dimension of the stress intensity factor. Numerous scholars also proposed a new theory to improve the existing deficiencies. Considering the special properties of interfacial cracks, according to the interface layer theory, the interface of bonded materials is regarded as a functionally graded material (FGM) layer with thickness, and a theoretical system of fracture mechanics for the elliptic cracks in the interfacial layer is established for analysis, which provides a theoretical method for the study of three-dimensional interface crack failure mechanism.
3. 三维界面裂纹模型
与均质材料裂纹不同的是, 即使远场仅受单向应力作用, 界面裂纹仍然处于 I、 II、 III 型复合的状态。 至于一般三维界面裂纹,其应力状态可能会表现得更加复杂。因此,考虑到界面裂纹的多型耦合现象,
相关文档
最新文档