基于界面层模型的三维界面裂纹分析
基于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(图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区,使推力的分布更加均匀,减少推力与焊接点之间的应力集中。
当然,也可以对衬层材料进行必要的改性,以增加材料的可塑性和强度,以减少脱粘裂纹的发生。
在使用三区域界面模型改善固体推进剂衬层界面脱粘裂纹问题时,还需要注意其他因素。
例如,如果衬层材料的温度太高,可能导致固体推进剂内部物理性质的变化,进而影响衬层材料的机械性能。
此外,其它影响脱粘裂纹的因素还包括衬层材料的成分及其在表面的厚度,以及衬层材料的受力方式等。
因此,在应用三区域界面模型时,还要根据实际情况考虑其他影响因素,识别各个影响因素对脱粘裂纹产生的影响,并采取相应的措施。
三相压电复合本构模型中的弧形界面裂纹
三相压电复合本构模型中的弧形界
面裂纹
三相压电复合本构模型中的弧形界面裂纹是一种特殊的材料缺陷,通常发生在金属/复合材料表面之间。
这种缺陷形式也被称为“外壳式”或“拉伸式”缺陷。
弧形界面裂纹具有半径R和深度d这两个参数,它们是三相复合材料结构中最重要的缺陷参数,因为它们决定了材料宏观性能的大小。
弧形界面裂纹的存在受到复合材料的材料结构、环境因素和加工工艺的影响。
它们的形成通常是由于材料的结构不均匀导致的腐蚀,例如氧化物的形成。
此外,由于金属材料和复合材料的热膨胀系数不一致,也会导致弧形界面裂纹的形成。
弧形界面裂纹的形成往往会降低复合材料的整体性能,使得材料的力学性能受到影响。
弧形界面裂纹可能会导致材料出现拉伸,压缩,剪切和剪切剪切等失效模式,并可能导致材料的破坏。
此外,弧形界面裂纹也会影响复合材料的热传导性能。
由于弧形界面裂纹的存在,热量的传递会受到阻碍,从而使得复合材料的散热性能受到影响。
为了减少弧形界面裂纹的影响,有必要采取一些有效的措施来改善复合材料的材料结构和加工工艺,以减少弧形界面裂纹的形成。
针对复合材料的材料结构,可以采用抗氧化剂、抗腐蚀剂等技术;针对复合材料的加工工艺,可以尝试采用低温加工、热处理等工艺来改善材料的性能,以减少弧形界面裂纹的形成。
在三相压电复合本构模型中,弧形界面裂纹是一种特殊的材料缺陷,它影响着材料的宏观性能,因此,在设计和制造过程中必须注意控制和减少弧形界面裂纹的形成,以保证复合材料的可靠性和安全性。
基于APDL和UIDL联合技术计算三维表面裂纹应力强度因子
Journal of Mechanical Strength2020, 42(5):1223-1229DOI :10. 16579/j.issn.l001. 9669. 2020. 05. 030基于APDL和UIDL联合技术计算三维表面裂纹应力强度因子< CALCULATION OF 3-D SURFACE CRACK STRESS INTENSITYFACTOR BASED ON APDL AND UIDL车福炎"耿黎明…(中国舰船研究设计中心,武汉430064)CHE FuYan GENG LiMing(China Ship Development and Design Center,Wuhan 430064, China)摘要工程上常将裂纹看作是半椭圆表面裂纹。
由于建立半椭圆表面裂纹的全裂纹有限元模型并非易事,在应用有限元法求解裂纹的应力强度因子时,常取裂纹模型的1/4进行分析。
采用模型分块处理方法,利用AnS yS软件,联合APD1.和UIDL技术进行二次开发,实现通过简单界面操作建立半椭圆表面裂纹的全裂纹模型,并可将全裂纹模型方便地嵌人到实体结构中。
使用该建模方法分别求解了中心半椭圆表面裂纹和斜置半椭圆表面裂纹的应力强度因子,与理论解和文献中的求解结果基本一致。
表明使用该二次开发建立的全裂纹模型求得的应力强度因子是准确可靠的。
该联合开发的建模功能可以用来求解复杂外载和复杂结构的应力强度因子,对工程实际应用有重要的应用价值。
关键词半椭圆表面裂纹应力强度因子Ansys二次开发全裂纹模型中图分类号0346. 1Abstract In engineering, cracks are usually regarded as semi-elliptical surface cracks. It is not easy to establish a complete finite element model of a semi-elliptical surface crack, thus a quarter of the crack model is generally used to calculate the stress intensity factor by finite element method. This paper adopts APDL and UIDL technologies for secondary development of ANSYS software, establishes a complete model of semi-elliptical surface crack by simple interface operation with the method of model block processing. Using this method, the stress intensity factors of the central semi-elliptical crack and the oblique semielliptical crack were solved respectively. The results are basically consistent with the theoretical solution and the results in the literature, which shows that it is accurate and reliable to use the complete crack model established by this method to calculate the stress intensity factor. The combined development modeling function can be used to solve the stress intensity factors in complex structures under complex loads, and has important value for practical application of engineering.Key words Semi-elliptical surface crack;Stress intensity factor;ANSYS secondary development;Complete crack modelCorresponding author:CHE FuYan, E-mail:u20\3\****************y Fax:+86-27-88071042Manuscript received 20190507,in revised form 20191022.引言准确评估疲劳裂纹扩展寿命的一个关键因素是精 确计算裂纹的应力强度因子虽然应力强度因子手册中收录了多种裂纹应力强度因子的计算公式,但 是对于实际结构,应力强度因子A' —般不存在闭合 解。
基于ANSYS的穿透裂纹三维板数值分析
( 3 )
…
板中心裂纹
式中: 厂—— 围绕裂纹尖 端 的一条任 意 反时 针 回路 , 起端始 于裂 卜 弹塑性条件下 , 回路 ,上任 意点 ( Y ,)的应变 能密
图 1 函数 t) 0 取值 图 表 1 中 心裂 纹 计 算 结 果
・
5 6・
第3 8卷 第 2 9期
2 2 月 1 年1 0 0
山 西 建 筑
裂纹 的力学特性 。板表面均匀压力荷载及约束如图 8 图 9所示 。 阶梯荷载分别计算 。模拟结果与解 析解对 比, , 得到三 种不 同位 置
3 计 算方 法及 理论依 据
3 1 J积 分 .
3 10 . 1 46 5 .6 62 0 .2 7 75 .7 9 30 .3
误 差/ %
18 .3 18 .1 18 .2 18 .2 18 .2
K AJ Cl C
319 .2 4 64 .9 628 .5 783 .2 938 . 8
误差/ %
图0三 不 位 裂 情 的()值 线 1 种 同置 纹 况 , 取 曲。 为 ÷
[ 6
篓 研 院 力 度 子 册M 北 :学 版 究 ・ 强 因 手 [ ・京科 出 应 ]
On n m e i n l ss o e e r t d c a k t r e d m e so l t a e n ANS u rc a a y i fp n t a e r c h e . i n i n p a e b s d o YS
, ,
1 ia .2r 2
爹哼. 卧 参考文献 : x . [ ] 康颖安. 1 断裂力学的发展 与研 究现 状 [ ] 湖 南工程 学 院学 J.
各向异性双材料III型界面裂纹应力分析
第 5期
谢秀峰 等: 各向异性双材料 I I I 型界面裂纹应力分析
7 3 5
坐 由弹性力学可知, 位移函数的控制方程组为: 标 表 ( +2 ( , 一 不
极
( s )
+ 2 ( s ) ,
:
+ ( )
) a a w j
- 0 ,( = , 2 ). )・
作者简介 :谢秀 峰( 1 9 8 0 . ) , 女, 讲 师, 硕士 , 从 事偏 微分方程理论 及其应 用研究. E ma i l : x x f 2 1 @l 2 6 . c o n. r 基 金项 目: 山西省 自然科 学基金 ( No . 2 0 1 1 0 1 1 0 2 1 . 3 , N o . 2 0 1 2 0 1 1 0 1 9 . 2 ) ;山西省高校科技 项 目( N o . 2 0 1 1 1 0 9 0 1 2 3 0 1 0 ) 的资助.
为各向异性复合材料I , 其坐标轴与弹性主轴夹角为 , 材料工程常数为剪切模量( G 2 ) , ( G 3 。 ) 。 , 而Y< 0 为各向异 性复合材料I I , 其坐标轴与弹性主轴夹角为 , 材料工程常数为剪切模量( G 2 ) : , ( ) .
收稿 1 5 t 期 :2 0 1 3 - 0 5 . 2 4
第3 9卷第 5 期
J
o u r n a l o f S o u t h w es tUn i v er si t yf o r Na t i l i t i e s Na 1 on a t u r a l S c i e n c e E d i t i o n
.
1
( 太原科 技 大学应 用科 学学院 太原 0 3 0 0 2 4 )
裂缝三维地质建模的难点与对策
裂缝三维地质建模的难点与对策地质建模是石油勘探开发过程中的重要环节,其目的是通过对地质模型的构建和分析,了解油气藏的分布、形态、储量等关键信息,为油气勘探开发提供科学依据。
而在地质建模的过程中,裂缝是一个非常重要的地质构造,其对储层的物性、渗透率等有着重要的影响。
因此,裂缝三维地质建模成为了当前石油勘探开发领域中的热点问题。
本文将就裂缝三维地质建模的难点及对策进行探讨。
一、裂缝三维地质建模的难点1. 裂缝的复杂性裂缝是一种非常复杂的地质构造,其形态、大小、分布等都具有高度的不确定性。
裂缝的形态有多种,包括水平、斜向、垂直等不同类型,且同一类型的裂缝在不同地质条件下也会有不同的形态。
此外,裂缝的大小也存在很大的差异,从毫米级的微裂缝到数十米的大裂缝都有可能存在,这就对建模的精度和复杂度提出了很高的要求。
2. 数据获取的困难性裂缝三维地质建模需要大量的数据支持,包括地震数据、井位数据、钻井数据、岩心数据等。
然而,这些数据的获取都存在一定的困难性。
例如,地震数据需要通过地震勘探技术获取,但这种技术成本高昂,且受到地质条件等因素的限制;井位数据和钻井数据需要通过实地勘探获取,但勘探难度大,成本高,且数据质量不稳定;岩心数据需要通过钻井获取,但钻井成本高,取样难度大,且数据量有限。
因此,数据获取的困难性也成为了裂缝三维地质建模的一个难点。
3. 建模精度的要求裂缝三维地质建模的精度要求非常高,因为裂缝对储层的影响非常显著。
建模精度不足会导致模型中存在大量的误差,从而影响到后续的油气勘探开发工作。
此外,裂缝的分布也具有很大的随机性,需要根据实际情况进行合理的概率分布建模,这也增加了建模的难度。
二、裂缝三维地质建模的对策1. 数据采集技术的改进为了解决数据获取的困难性,可以采用一些新的技术手段来获取数据。
例如,可以利用地震反演技术获取地震数据,这种技术可以大大降低地震勘探的成本,同时也可以获得更为精确的数据;可以利用遥感技术获取地表的地质信息,这种技术可以大大缩短勘探周期,同时也可以提高数据的精度和质量。
双材料界面断裂力学模型与实验方法
3、在FRP混凝土界面断裂研究中 的应用优势
双材料梁界面力学模型在FRP混凝土界面断裂研究中具有以下优势: (1) 能够考虑界面两侧材料的异质性和不均匀性; (2)能够预测界面的破坏模式和 承载能力; (3)能够分析界面断裂的机理和影响因素,为优化FRP混凝土结构 设计提供理论支持。
FRP混凝土界面断裂研究
3、通过有限元模拟,我们发现界面裂纹萌生的临界应力与实验结果相符, 但在裂纹扩展过程中存在局部化现象,且对界面的整体承载能力产生影响。
4、根据断裂力学理论分析,我们发现界面的断裂准则可归纳为应力强度因 子与临界应力之间的比较,这也与实验和数值模拟结果相吻合。
结论与展望
本研究通过实验、数值模拟和理论分析相结合的方法,深入探讨了基于断裂 力学的钢筋、FRP与混凝土界面力学特性。虽然取得了一定的研究成果,但仍存 在以下不足之处:
研究方法
本研究采用实验、数值模拟和理论分析相结合的方法,以探讨基于断裂力学 的钢筋、FRP与混凝土界面力学特性。具体研究方法如下:
1、实验研究:通过准静态拉伸试验和疲劳试验,测定界面的粘结性能和断 裂行为。
2、数值模拟:采用有限元方法对界面进行建模和分析,模拟界面裂纹的萌 生、扩展和断裂过程。
材料选择
双材料的选择应遵循以下原则:
1、强度:选用的两种材料应具有较高的强度,以确保在承受载荷时不会发 生过早的破坏。
2、韧性:材料应具有一定的韧性,以吸收冲击能量,提高材料的抗断裂能 力。
3、加工性能:材料的加工性能应良好,以便于制作成各种形状和尺寸的构 件。
4、相容性:两种材料的相容性应良好,以避免在界面处产生化学腐蚀或电 化学腐蚀。
3、理论分析:基于断裂力学理论,对界面裂纹的扩展规律进行系统分析, 推导相应的断裂准则。
三维介质界面裂纹的动力分析
三维介质界面裂纹的动力分析
陶伟明;曹志远
【期刊名称】《地震工程与工程振动》
【年(卷),期】1995(15)1
【摘要】本文通过区域的连续映射推导了一种映射半解析边界元法,将其应用于界面裂纹问题的动力分析,可大大减小计算工作量。
典型算例的结果表明,本文方法是非常有效的,并可推广应用于三维复杂形状裂纹的计算问题。
【总页数】8页(P111-118)
【关键词】边界元法;动力分析;半解析法;界面裂纹
【作者】陶伟明;曹志远
【作者单位】同济大学
【正文语种】中文
【中图分类】O346.1
【相关文献】
1.含垂直界面有限裂纹层状介质的断裂分析 [J], 彭达仁;张起森
2.基于界面层模型的三维界面裂纹分析 [J], 王硕;张铮;
3.三维压电介质界面裂纹的边界积分—微分方程 [J], 赵明皞;郭长江;方平治;沈亚鹏
4.不同功能梯度压电压磁层状介质中共线界面裂纹的动态性能分析 [J], 张培伟;周振功;王彪
5.三维横观各向同性介质界面裂纹的边界积分方程方法 [J], 赵明皞;李冬霞;沈亚鹏
因版权原因,仅展示原文概要,查看原文内容请购买。
三维断裂分析软件FRANC3D教程浅析
三维断裂分析软件FRANC3D教程浅析张涛;杨玉玲【摘要】介绍了三维断裂分析软件FRANC3D(FRacture Analysis Code in 3 Dimensions)的各个组成部分以及在学习过程中应该注意的问题和常见的错误,就其初级教程的内容及原理做简单说明,以减少在以后应用软件的过程中产生错误.【期刊名称】《山西建筑》【年(卷),期】2010(036)025【总页数】3页(P366-368)【关键词】断裂分析;FRANC3D;初级教程;简单说明【作者】张涛;杨玉玲【作者单位】贵州大学土木建筑工程学院,贵州贵阳,550003;贵州大学土木建筑工程学院,贵州贵阳,550003【正文语种】中文【中图分类】TP391.30 引言三维断裂分析软件FRANC3D是美国康奈尔大学开发的一款专门用于断裂分析的自由软件,它和OSM(Object Solid Modeler)立体对象建模器、BES(Boundary Element System)边界元系统共同构成一套完整的断裂分析系统。
OSM主要用于前处理的工作,比如坐标的输入、线面体的形成等,BES主要用于后处理的工作(有限元的分析)。
初级教程里介绍了一个铝质材料的简单裂纹分析及扩展模拟。
以下就此教程做简要的概述和分析。
1 OSM(立体对象建模器)1.1 对象界面及按钮介绍运行OSM后,最先看见的就是如图1所示的界面,它主要由主菜单、信息栏以及一些显示控制按钮组成。
对于先前建好的OSM模型,此时可以直接选择“Read OSM Modle File”读入模型文件,第一次建模型则需选择“Create New Modle”,按照教程所给坐标点,输入坐标。
连点成线、连线成面,其中模型中的四分之一圆弧是由两个弧线连接而成的(如图2所示)。
另外,需注意每个面的外法线的方向(教程上要求每个面法线的方向都必须向外),查看法线的方向可以点选显示菜单(SELECT,DISPLAY)中的NORMALS。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
王硕,张铮
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 型复合的状态。 至于一般三维界面裂纹,其应力状态可能会表现得更加复杂。因此,考虑到界面裂纹的多型耦合现象,