GlassEpoxy复合材料叠层板面内剪切连续损伤模型

合集下载

复合材料层合板损伤模型的建构方法及其应用

复合材料层合板损伤模型的建构方法及其应用

复合材料层合板损伤模型的建构方法及其应用复合材料是一种具有良好性能和广泛应用的材料。

然而,在复合材料的使用过程中,由于不同载荷条件和作业环境的影响,往往会出现各种各样的损伤。

为了研究复合材料损伤的特征和演变规律,构建一个准确可靠的损伤模型显得尤为重要。

在复合材料的层合板结构中,由于材料层的叠加和粘结,导致了一个复杂的结构体系。

损伤的发生和演化是由于层与层之间的应力分布不均匀,导致局部应力过大而引起的。

因此,构建一个层合板的损伤模型,需要考虑到材料层的叠加效应和相互作用。

首先,建构层合板的损伤模型需要选取适当的力学模型。

一般来说,可以采用线性弹性模型或非线性模型进行分析。

线性弹性模型是一种较为简单的模型,但在材料的强度有限和变形较大时,其适用性较差。

非线性模型则可以更准确地描述材料在不同应变和应力水平下的力学性能,但计算复杂度较高。

根据具体研究目的和可行性,选择合适的模型进行建模。

其次,需要确定层合板内部的应力分布情况。

根据层合板的几何形状和受力情况,可以通过结构力学原理和有限元分析等方法计算得到不同位置的应力值。

考虑到复合材料的各向异性特性,需要对应力进行向量描述,同时考虑拉伸、剪切和压缩等各个方向的作用。

然后,在选定力学模型和确定应力分布后,需要考虑层合板中各层之间的力学联系。

层合板的损伤扩展是由于材料层间的界面错位、粘结破坏或层板剥落等原因引起的。

因此,需要在模型中加入合适的界面层,考虑不同层之间的相互作用和分离。

最后,将建构的层合板损伤模型应用于实际问题中。

通过对该模型的分析和计算,可以获得层合板在不同载荷下的应力分布、位移变化和损伤演化情况。

根据模型结果,可以预测层合板的强度、疲劳寿命和失效机制等,为实际应用提供参考依据。

综上所述,构建复合材料层合板的损伤模型是一项复杂而重要的工作。

通过选择合适的力学模型、确定应力分布、考虑层间界面和进行实际应用,可以使模型更准确,并为复合材料的设计和工程应用提供理论支持综合考虑应力水平下的力学性能、应力分布情况、层间力学联系以及实际应用,构建复合材料层合板的损伤模型是一项复杂而重要的工作。

复合材料破损模式及损伤评估方法研究

复合材料破损模式及损伤评估方法研究

复合材料破损模式及损伤评估方法研究复合材料作为目前发展最快的结构材料之一,其广泛应用于航空、航天、汽车、船舶、建筑和体育器材等领域。

其轻便、强度高、抗腐蚀等特性受到人们的青睐,但同时也存在着不易损伤检测和修复难的问题。

因此,破损模式和损伤评估方法的研究对复合材料的应用和发展有着重要的意义。

一、复合材料破损模式分析1. 纵向破裂纵向破裂指的是工程承受纵向载荷时,纤维沿纵向方向产生拉伸破裂的现象。

这种破损模式通常发生在直角或圆形截面的薄壁复合材料构件中。

2. 层间剥离层间剥离是指复合材料中不同层之间由于应力集中或者其它因素导致的剥离现象。

这种损伤可能在第一层或者在中间层发生,在非对称复合材料构件中表现得更为明显。

3. 面内剪切面内剪切是指在载荷作用下,双向纤维绕组或多向针织复合材料中发生的层间滑移现象。

这种损伤模式容易在双向和针织复合材料中出现。

4. 晶格破坏在一些高强度高温复合材料中,由于材料长时间的受热膨胀和收缩等因素,会产生晶格结构的临界结构转换和断裂,这种破损模式可能导致材料的劣化和性能降低。

二、复合材料损伤评估方法1. 红外成像红外成像技术是利用红外光谱的成像方式,来检测材料内部的热量分布、热辐射和热传输特性,以预测材料中的粘结、裂纹、孔洞和表面腐蚀等损伤程度。

2. 声发射检测声发射检测是利用材料的声学反射和回声特性,通过记录声波传播的时间、波形、频率和幅度来识别并评估材料损伤程度和剥离程度,以提供损伤评估和修复指导。

3. 超声波检测超声波检测是利用声波传播在材料内部的特性,通过检测声波的时间、波形、频率和幅度等参数,来避免材料表面的干扰,对材料内部的损伤程度进行评估和定位。

4. X射线检测X射线检测是利用电磁波的透射能力,对材料内部的非均匀结构进行检测和成像。

对材料中的异常区域、孔洞、裂缝、层间剥离、残余应力等缺陷进行识别和定位。

三、结论复合材料在应用过程中容易发生各种各样的损伤,这些损伤对于材料性能的损害和材料寿命的缩短都有着直接的影响。

一种复合材料层合板分层扩展的修正内聚力模型

一种复合材料层合板分层扩展的修正内聚力模型

JournalofMechanicalStrength2021ꎬ43(5):1262 ̄1266DOI:10 16579/j.issn.1001 9669 2021 05 035∗20200713收到初稿ꎬ20200905收到修改稿ꎮ∗∗贾云龙ꎬ男ꎬ1992年5月生ꎬ河南周口人ꎬ汉族ꎬ航天行云科技有限公司主管设计师ꎬ硕士ꎬ主要从事振动与强度研究ꎮ一种复合材料层合板分层扩展的修正内聚力模型∗ACORRECTEDCOHESIVEZONEMODELINTHEINTERLAMINARFRACTUREOFCOMPOSITE贾云龙∗∗1㊀㊀徐㊀琳2㊀㊀项㊀斌1(1.航天行云科技有限公司ꎬ武汉430040)(2.武汉理工大学高性能船舶技术教育部重点实验室ꎬ武汉430063)JIAYunLong1㊀XULin2㊀XIANGBin1(1.LeobitTechnologyCo.ꎬLtd.ꎬWuhan430040ꎬChina)(2.KeyLaboratoryofHighPerformanceShipTechnologyofMinistryofEducationꎬWuhanUniversityofTechnologyꎬWuhan430063ꎬChina)摘要㊀针对复合材料的分层现象ꎬ考虑初始断裂韧性和裂尖材料损伤的影响ꎬ提出了一种修正内聚力模型ꎬ并给出了内聚力模型参数确定的方法ꎮ根据相关试验ꎬ采用该方法ꎬ获得了考虑初始断裂韧性和裂尖材料损伤的内聚力模型ꎬ基于该内聚力模型对试验进行了仿真分析ꎮ结果表明ꎬ与仅考虑初始断裂韧性获取的内聚力模型相比ꎬ考虑初始断裂韧性和裂尖材料损伤的内聚力模型ꎬ可以更加准确地模拟复合材料的分层现象ꎮ关键词㊀复合材料㊀内聚力模型㊀裂纹扩展㊀数值仿真中图分类号㊀TB332Abstract㊀Inordertoconsiderbothmaterialdamageofthecracktipandtheinitialfracturetoughnessꎬacorrectedcohesivezonemodel(CZM)thatconsiderstheinfluenceofmaterialdamageoftheinitialcrackpropagationhasbeenproposed.Thecorrectedcohesivezonemodeobtainedfromtheendnotchflexure(ENF)testhasbeenimplementedintheABAQUS/Standard.Theresultsshowthatitisnecessarytoconsiderbothmaterialdamageofthecracktipandtheinitialfracturetoughnessꎬandthecorrectedcohesivezonemodecanaccuratelysimulatethecrackpropagation.Keywords㊀CompositeꎻCohesivezonemodeꎻCrackpropagationꎻFEMCorrespondingauthor:JIAYunLongꎬE ̄mail:whutjiayunlong@foxmail.comꎬTel:+86 ̄27 ̄59393300Manuscriptreceived20200713ꎬinrevisedform20200905.㊀㊀引言复合材料结构由于具有高比模量㊁高比强度㊁优良的能量吸收性能㊁耐腐蚀性ꎬ越来越广泛地应用于船舶㊁航空航天等领域ꎮ层合复合材料结构㊁夹芯复合材料结构㊁胶接复合材料结构等受到了广泛关注[1 ̄8]ꎮ在复合材料结构的服役期间ꎬ材料界面的初始缺陷或应力奇异点ꎬ使复合材料结构的界面处极易产生分层破坏ꎬ导致复合材料结构的失效[9 ̄13]ꎮ目前ꎬ内聚力模型(CZM)广泛地用于研究复合材料分层的现象ꎮ内聚力模型使用界面应力与界面相对位移(Traction ̄Separation)描述复合材料的界面分层现象ꎮ通过选择合理的参数ꎬ内聚力模型可以很好地模拟复合材料的分层现象[14 ̄15]ꎮ对于复合材料II型界面分层现象ꎬ1982年RussellAJ和StreetKN利用端部缺口弯曲(ENF)试验进行了研究[16]ꎮ根据端部缺口弯曲试验ꎬ基于线弹性断裂力学的分析表明ꎬ随着预置的裂纹长度的变化ꎬII型裂纹的初始断裂韧性也随着变化[17]ꎮ基于内聚力模型ꎬ对端部缺口弯曲试验的分析表明ꎬII型裂纹的初始断裂韧性和预置裂纹的长度无关[18]ꎮ当裂纹扩展时ꎬ裂纹尖端附近材料产生损伤ꎬ这种损伤对于复合材料的分层现象不可忽视[19][20]2443 ̄2455ꎮ因此ꎬ没有考虑裂纹尖端附近材料损伤的线弹性断裂力学就不适用了ꎮ在研究复合材料II型分层现象的端部缺口弯曲试验中ꎬ断裂韧性的准确计算ꎬ须考虑裂尖处的材料损伤ꎮ内聚力模型是一种简单有效考虑裂尖材料损伤的方法[21]ꎮ基于端部缺口弯曲试验ꎬ文献[22]473 ̄487对II型裂纹进行了研究ꎮ文献[23]174 ̄181采用内聚力模型对该试验进行了有限元模拟ꎬ指出双线性内聚力模型㊀第43卷第5期贾云龙等:一种复合材料层合板分层扩展的修正内聚力模型1263㊀㊀不能模拟II型裂纹的扩展ꎮ考虑到裂尖材料损伤和初始断裂韧性ꎬ提出了一种确定内聚力模型参数的方法ꎮ基于相关试验结果ꎬ采用该方法得到了考虑初始断裂韧性和裂尖材料损伤的内聚力模型ꎮ基于该内聚力模型ꎬ采用商业有限元软件Abaqusꎬ对试验现象进行了仿真ꎮ结果表明ꎬ较其他方法得到的内聚力模型ꎬ该内聚力模型可以更准确地模拟复合材料结构的分层现象ꎮ1㊀理论基础1 1㊀内聚力模型在使用内聚力模型进行复合材料分层现象的研究中ꎬ双线性内聚力模型和梯形内聚力模型是比较常用的ꎬ如图1所示ꎮ图1㊀常用内聚力模型Fig.1㊀ThecommonformsofCZM在初始加载阶段ꎬ内聚力τ与界面相对位移δ∗成线性关系ꎻ当内聚力τ等于内聚力强度τC时ꎬ界面材料的刚度随着界面相对位移的增加开始衰减ꎻ当界面的相对位移等于界面相对位移临界值δf时ꎬ界面材料彻底失效ꎬ刚度变为零ꎮ应力位移曲线所围成的面积等于断裂韧性GIIcꎮ对于复合材料分层现象ꎬ断裂韧性通常随裂纹扩展而逐渐增加ꎬ可以用R曲线表示ꎬ如图2所示ꎮ结合内聚力模型ꎬ描述复合材料分层现象的R曲线可做如下[23]174 ̄181GIIδ∗()=GIIcꎬini+ʏδ∗0τδ()dδ(1)式中ꎬτ为裂尖材料损伤区的剪应力ꎻGIIcꎬini为界面材料的初始断裂韧性ꎻ裂尖位置的界面相对位移为δ∗ꎻδ是裂尖材料损伤区的界面相对位移ꎮ根据公式(1)即可得到内聚力模型[24]τδ∗()=∂GII∂δ∗(2)㊀㊀当试验获取GII ̄δ∗曲线后ꎬ根据公式(2)就可得到描述复合材料分层现象的内聚力模型ꎮ图2㊀R曲线示意图Fig.2㊀AschematicshapeofR ̄curvebehavior1 2㊀ENF试验端部缺口弯曲试件如图3所示ꎮ试件预置裂纹的长度为aꎬ宽度为bꎬ长度为2Lꎮ试件材料为等厚度单向层合板ꎬ纤维方向沿试件长度方向ꎮP为试验加载载荷ꎬΔ为加载点处的位移ꎮ图3㊀端部缺口弯曲试件Fig.3㊀Theendnotchedflexure(ENF)specimen为了考虑裂尖材料的损伤ꎬ文献[20]2443 ̄2455[25]877 ̄897提出了基于等效裂纹长度的修正梁理论(CBTE)ꎮ等效裂纹长度的公式如下[26]ae=38E1bh3CC3-2L33(3)CC=ΔP-3L10G13bh(4)式中ꎬae为等效裂纹长度ꎬ其考虑到了裂尖材料损伤以及剪切变形ꎮE1㊁G13是材料的弹性模量ꎮ基于等效裂纹长度ꎬ可得G=3P2ae264bD1(5)㊀㊀根据试验数据ꎬ采用上述公式(5)ꎬ可得界面材料的断裂韧性ꎮ同时ꎬ根据试验中记录的界面相对位移δ∗ꎬ即可得到GII ̄δ∗曲线ꎮ从而ꎬ根据公式(2)得到内聚力模型ꎮ但是ꎬ直接根据公式(2)得到的内聚力模型忽略了初始断裂韧性GIIcꎬiniꎬ导致内聚力模型τ ̄δ∗曲线所围区域的面积为GIIcꎬss-GIIcꎬini()ꎮ为了考虑初始断裂韧性GIIcꎬiniꎬ文献[23]174 ̄181使内聚力强度τC保持不变ꎬ通过增加δfꎬ使τ ̄δ∗曲线所围成的面积为GIIcꎬssꎬ以此对内聚力模型进行了修正ꎮ㊀1264㊀机㊀㊀械㊀㊀强㊀㊀度2021年㊀2㊀修正内聚力模型2 1㊀修正方法断裂韧性与界面相对位移的关系为FGIIcꎬδ∗()=0(6)式中ꎬGIIc为断裂韧性ꎻδ∗为界面相对位移ꎮ考虑到裂尖材料损伤λꎬ将式(6)改写如下FGIIcꎬδ∗ꎬλ()=0(7)㊀㊀获取内聚力模型的难点在于如何考虑裂尖材料损伤的影响ꎮ基于等效裂纹长度概念ꎬ文献[20]2443 ̄2455[25]877 ̄897提出了修正梁理论(CBTE)ꎬ以考虑裂尖材料损伤λ对断裂韧性GIIc的影响ꎮ借用上述思路ꎬ假设可用等效界面相对位移c考虑裂尖材料损伤λ对FGIIcꎬδ∗()=0的影响ꎬ如下FGIIcꎬδ∗-c()=0(8)㊀㊀根据公式(7)ꎬ将试验实测能量释放率与界面相对位移GII ̄δ∗曲线向右进行平移距离cꎬ即得到考虑裂尖材料损伤的GII ̄δ∗曲线ꎬ如图4所示ꎮ图4㊀GII ̄δ∗曲线的修正示意图Fig.4㊀TheschematicofthecorrectedGII ̄δ∗curve根据公式(2)对平移后的GII ̄δ∗曲线进行微分ꎬ就可以得到考虑裂尖材料损伤的内聚力模型τ ̄δ∗曲线ꎮ但是ꎬ根据公式(2)微分得到的内聚力模型τ ̄δ∗曲线含有未知量cꎮ考虑到τ ̄δ∗曲线所围成的面积须为GIIcꎬssꎬ依此得到c的大小ꎬ从而确定内聚力模型ꎮ2 2㊀修正结果对于端部缺口弯曲试验ꎬ文献[22]473 ̄487中获取的GII ̄δ∗试验拟合曲线如下GIIδ∗()=-56 2δ∗2+15 3δ∗+0 69(9)㊀㊀由式(8)可得到的未修正τ ̄δ∗曲线为τδ∗()=-112 4δ∗+15 3(10)㊀㊀根据本文的方法得到的τ ̄δ∗曲线为τδ∗()=-112 4δ∗+19 7(11)㊀㊀根据文献[23]174 ̄181的方法得到τ ̄δ∗曲线为τδ∗()=-67 7δ∗+15 3(12)㊀㊀上述τ ̄δ∗曲线的对比ꎬ如图5所示ꎮ根据图5ꎬ当直接采用这些内聚力模型进行有限元图5㊀修正前后τ ̄δ∗曲线对比Fig.5㊀Comparisonofτ ̄δ∗curves仿真时ꎬ在δ∗=0时ꎬ其应力位移关系的跳跃变化将导致数值求解的奇异性ꎮ为了克服在δ∗=0处的数值求解困难ꎬ引入一个很大的初始刚度KIIꎬ如图6所示ꎮ图6㊀初始刚度示意图Fig.6㊀SchematicoftheInitialstiffness3㊀数值分析3 1㊀网格和边界条件基于2 2节的内聚力模型ꎬ采用有限元软件Abaqusꎬ对端部缺口弯曲试验进行仿真模拟ꎮ试件的材料参数如表1所示ꎮ表1㊀试件材料参数Tab.1㊀Materialpropertiesofunidirectionalglass/epoxy参数Parameter值ValueE1/MPa18000E2/MPa2200E3/MPa2200G12/MPa1800G13/MPa1800G23/MPa1600ν120 29ν130 29ν230 38试件尺寸如图7所示ꎬ宽度为20mmꎬ试验支座和加载压头的半径均为2mmꎬ试验加载载荷为PꎬΔ为试验加载点处位移ꎮ仿真模型使用六面体单元(C3DR8)进行网格划分ꎻ胶层模型使用0厚度的内聚力单元(COH3D8)模拟ꎬ预置裂纹模型使用Seam裂纹模型模拟ꎬ两个试验㊀第43卷第5期贾云龙等:一种复合材料层合板分层扩展的修正内聚力模型1265㊀㊀图7㊀试件尺寸示意图Fig.7㊀TheENFspecimendimensions支座与一个加载压头使用解析刚体模拟ꎮ在裂纹尖端处㊁支座区域与加载点区域ꎬ取0 5mm网格尺寸ꎮ对于距离裂纹尖端㊁支座与加载点较远的区域ꎬ取1mm网格尺寸ꎮ沿试件的宽度方向ꎬ其网格的尺寸取1mmꎮ网格的划分情况如图8所示ꎮ仿真时ꎬ开启几何非线性ꎮ图8㊀网格划分示意图Fig.8㊀FiniteelementmodelofthetestedENFspecimens3 2㊀仿真结果及分析将仿真结果与试验结果进行对比ꎬ如图9所示ꎮ由图9可知ꎬ相较于文献[23]174 ̄181中的内聚力模型ꎬ本文方法得到的内聚力模型ꎬ其仿真结果与试验结果的吻合度更高ꎮ文献[23]174 ̄181的内聚力模型对初始断裂韧性GIIcꎬini进行了考虑ꎬ修正前后ꎬ内聚力强度τC=15 30MPaꎮ采用考虑裂纹尖端材料损伤和初始断裂韧性的方法获取的内聚力模型ꎬ其内聚力强度τC=19 72MPaꎮ从仿真结果P ̄Δ来看ꎬ考虑初始断裂韧性和裂尖材料损伤的内聚力模型可以更准确地对试验现象进行模拟ꎮ仿真GII ̄δ∗曲线与试验GII ̄δ∗曲线的对比ꎬ如图10所示ꎮ对于仿真所得断裂韧性GII都是从零开始ꎬ随着界面相对位移δ∗的增加逐渐而增加ꎬ并最终趋于稳定值ꎮ试验实测断裂韧性GII是从0 69N/mm开始逐渐增加ꎬ这是由于在裂纹开始扩展之前ꎬ裂尖材料产生损伤ꎬ耗散了能量ꎮ因此ꎬ在使用内聚力模型进行仿真时ꎬ需在内聚力模型中考虑裂尖材料损伤的影响ꎮ根据图10ꎬ采用本文方法得到GII ̄δ∗曲线与文献[23]174 ̄181的GII ̄δ∗曲线存在较大差异ꎮ其原因是本文的内聚力模型考虑了初始断裂韧性和裂尖材料损图9㊀有限元计算值与试验值的比较Fig.9㊀ComparisonbetweentheexperimentalandnumericalP ̄Δ图10㊀GII ̄δ∗曲线的比较Fig.10㊀ComparisonofGII ̄δ∗curves伤ꎬ文献[23]174 ̄181的内聚力模型仅考虑了初始断裂韧性ꎮ综上ꎬ较于仅考虑初始断裂韧性的内聚力模型ꎬ考虑初始断裂韧性和裂尖材料损伤影响的内聚力模型ꎬ可以更准确地模拟复合材料的分层现象ꎮ4㊀结论针对复合材料的分层现象ꎬ考虑初始断裂韧性和裂尖处材料的损伤ꎬ提出了一种修正内聚力模型ꎬ并对文献[22]473 ̄487中的试验进行了仿真ꎮ主要结论如下:1)内聚力模型的能否准确模拟复合材料裂纹扩展现象取决于内聚力模型的获取方法ꎮ2)考虑初始断裂韧性和裂尖材料损伤ꎬ给出了一种内聚力模型的获取方法ꎮ3)仿真结果表明ꎬ计及初始断裂韧性和裂尖材料损伤的内聚力模型可以更准确地模拟复合材料分层现象ꎮ参考文献(References)[1]㊀蔡忠云ꎬ唐文勇ꎬ陈念众ꎬ等.复合材料船体层合板的极限强度分析[J].船舶力学ꎬ2009(1):72 ̄81.CAIZhongYunꎬTANGWenYongꎬCHENNianZhongꎬetal.Ultimatestrengthanalysisofcompositelaminatedshippanels[J].JournalofShipMechanicsꎬ2009(1):72 ̄81(InChinese).[2]㊀孙枭雄ꎬ任慧龙ꎬ唐㊀宇.轻质夹芯复合材料结构强度评估方法㊀1266㊀机㊀㊀械㊀㊀强㊀㊀度2021年㊀研究[C].纪念船舶力学创刊二十周年学术会议ꎬ2017:435 ̄450.SUNXiaoXiongꎬRENHuiLongꎬTANGYu.Investigationonstrengthassessmentmethodsforlightensandwichcompositestructure[C].AcademicConferencetoCommemoratethe20thAnniversaryofthePublicationofShipMechanicsꎬ2017:435 ̄450(InChinese). [3]㊀曾海艳ꎬ严仁军ꎬ徐㊀琳ꎬ等.弯曲载荷下复合材料夹芯 L 型接头强度和疲劳试验研究[J].船舶力学ꎬ2017ꎬ21(12):1540 ̄1550.ZENGHaiYanꎬYANRenJunꎬXULinꎬetal.Experimentalstudyonstrengthandfatigueofsandwichcompositel ̄jointunderbending[J].JournalofShipMechanicsꎬ2017ꎬ21(12):1540 ̄1550(InChinese).[4]㊀王纬波ꎬ李永胜.复合材料板与钢板胶接㊁螺接与混合对接接头的力学特性研究[J].船舶力学ꎬ2011ꎬ15(9):1052 ̄1064.WANGWeiBoꎬLIYongSheng.Mechanicalbehaviorsofadhesively ̄bondedꎬboltedandhybridcomposite ̄to ̄steeljoints[J].JournalofShipMechanicsꎬ2011ꎬ15(9):1052 ̄1064.[5]㊀胡明勇ꎬ王安稳ꎬ姜㊀伟ꎬ等.复合材料层合板的动力响应和横向应力分析[J].船舶力学ꎬ2008(5):778 ̄784.HUMingYongꎬWANGAnWenꎬJIANGWeiꎬetal.Dynamicresponseandtransversestressesanalysisofcompositelaminatedplates[J].JournalofShipMechanicsꎬ2008(5):778 ̄784(InChinese).[6]㊀陈㊀悦ꎬ朱㊀锡ꎬ李华东ꎬ等.含分层缺陷复合材料夹芯梁力学特性及失效模式的试验研究[J].海军工程大学学报ꎬ2016(6):65 ̄70.CHENYueꎬZHUXiꎬLIHuaDongꎬetal.Mechanicalbehaviorandfailuremodesofcompositesandwichcolumnswithface/coredebond[J].JournalofNavalUniversityofEngineeringꎬ2016(6):65 ̄70(InChinese).[7]㊀赵㊀洁ꎬ陆㊀华ꎬ王明春.筋条分层损伤复材加筋壁板的稳定性分析及修理[J].机械强度ꎬ2016ꎬ38(4):871 ̄874.ZHAOJieꎬLUHuaꎬWANGMingChun.Stabilityanalysisandrepairofstiffenedcompositewallslabwithadelaminatongatthestringerinterface[J].JournalofMechanicalStrengthꎬ2016ꎬ38(4):871 ̄874(InChinese).[8]㊀刘㊀佶ꎬ许希武.Z ̄pin增强对自动铺丝复合材料T型接头拉脱㊁剪切性能影响的试验研究及数值模拟[J].机械强度ꎬ2017ꎬ39(3):527 ̄533.LIUJiꎬXUXiWu.Experimentandsimulationstudyontheeffectofz ̄pinreinforcedonthepull ̄offandshearcarryingcapacityoffiberplacementcompositeT ̄joints[J].JournalofMechanicalStrengthꎬ2017ꎬ39(3):527 ̄533(InChinese).[9]㊀马存旺ꎬ李志国ꎬ鲁国富ꎬ等.整体复合材料结构分层特性研究进展(一)[J].飞机设计ꎬ2014(5):31 ̄39.MACunWangꎬLIZhiGuoꎬLUGuoFuꎬetal.Developmentondelamiantioncharacteristicsofintegratedcompositestructures:PartI[J].AircraftDesignꎬ2014(5):31 ̄39(InChinese).[10]㊀OdessaIꎬFrostigYꎬRabinovitchO.Modelingofinterfacialdebondingpropagationinsandwichpanels[J].InternationalJournalofSolidsandStructuresꎬ2017(148/149):67 ̄78.[11]㊀YazdaniSꎬRustWJHꎬWriggersP.Delaminationonsetandgrowthincompositeshells[J].Computers&Structuresꎬ2018(195)(SupplementC):1 ̄15.[12]㊀TsaiMYꎬMortonJ.Anexperimentalinvestigationofnonlineardeformationsinsingle ̄lapjoints[J].MechanicsofMaterialsꎬ1995ꎬ20(3):183 ̄194.[13]㊀曹㊀勇ꎬ冯蕴雯ꎬ薛小锋ꎬ等.夹持力对填充孔复合材料层合板层间应力的影响分析[J].机械强度ꎬ2015ꎬ37(5):904 ̄909.CAOYongꎬFENGYunWenꎬXUEXiaoFengꎬetal.Impactanalysisofbolt ̄clampingforcetofilled ̄holelaminatesinterlaminarstress[J].JournalofMechanicalStrengthꎬ2015ꎬ37(5):904 ̄909(InChinese).[14]㊀CavalliMNꎬThoulessMD.Theeffectofdamagenucleationonthetoughnessofanadhesivejoint[J].JournalofAdhesionꎬ2001ꎬ16(1):75 ̄92.[15]㊀胡波涛ꎬ柴亚南ꎬ陈向明ꎬ等.后屈曲复合材料加筋板筋条 ̄蒙皮界面失效表征[J].机械强度ꎬ2019ꎬ41(6):1473 ̄1479.HUBoTaoꎬCHAIYaNanꎬCHENXiangMingꎬetal.Stiffener ̄skininterfacefailurecharacterizationofthestiffenedcompositepanelatpost ̄bucklingstage[J].JournalofMechanicalStrengthꎬ2019ꎬ41(6):1473 ̄1479(InChinese).[16]㊀RussellAJꎬStreetKN.Factorsaffectingtheinterlaminarfractureenergyofgraphite/epoxylaminates[J].ProgressinScienceandEngineeringofComposite(ICCM ̄IV)ꎬ1982:279 ̄286. [17]㊀BachrachWEꎬHicksTRꎬHabasZSꎬetal.Mixedexperimental ̄structuralmodelforinterlaminarshearfracturetoughness[J].JournalofAerospaceEngineeringꎬ1991ꎬ4(1):108 ̄125. [18]㊀OuyangZꎬLiG.Nonlinearinterfaceshearfractureofendnotchedflexurespecimens[J].InternationalJournalofSolidsandStructuresꎬ2009ꎬ46(13):2659 ̄2668.[19]㊀MouraMꎬMoraisA.EquivalentcrackbasedanalysesofENFandELStests[J].EngineeringFractureMechanicsꎬ2008ꎬ75(9):2584 ̄2596.[20]㊀BlackmanBRKꎬBrunnerAJꎬWilliamsJG.ModeIIfracturetestingofcomposites:Anewlookatanoldproblem[J].EngineeringFractureMechanicsꎬ2006ꎬ73(16):2443 ̄2455.[21]㊀MoSNꎬBelytschkoT.Extendedfiniteelementmethodforcohesivecrackgrowth[J].EngineeringFractureMechanicsꎬ2002ꎬ69(7):813 ̄833.[22]㊀AnyfantisKNꎬTsouvalisNG.ExperimentalandnumericalinvestigationofModeIIfractureinfibrousreinforcedcomposites[J].JournalofReinforcedPlastics&Compositesꎬ2011ꎬ30(6):473 ̄487.[23]㊀Heidari ̄RaraniMꎬGhasemiAR.AppropriateshapeofcohesivezonemodelfordelaminationpropagationinENFspecimenswithR ̄curveeffects[J].TheoreticalandAppliedFractureMechanicsꎬ2017(90):174 ̄181.[24]㊀SuoZꎬBaoGꎬFanB.DelaminationR ̄curvephenomenaduetodamage[J].JournaloftheMechanicsandPhysicsofSolidsꎬ1992ꎬ40(1):1 ̄16.[25]㊀BlackmanBRKꎬKinlochAJꎬParaschiM.ThedeterminationofthemodeIIadhesivefractureresistanceꎬGIICofstructuraladhesivejoints:Aneffectivecracklengthapproach[J].EngineeringFractureMechanicsꎬ2005ꎬ72(6):877 ̄897.[26]㊀DeMoraisAB.NovelcohesivebeammodelfortheEnd ̄NotchedFlexure(ENF)specimen[J].EngineeringFractureMechanicsꎬ2011ꎬ78(17):3017 ̄3029.。

复合材料层合板面内剪切实验方法的评价_王言磊

复合材料层合板面内剪切实验方法的评价_王言磊

① /②
93% 96% 92% 106% 96%
注 :表中数值为 5个试件的平均值 。
5种层合板在 2种面内剪切实验方法下典型的 剪应力 -应变曲线如图 2 ~ 6。从图中可以看到 , 大多 数层合板的剪应力 -应变曲线在剪应变为 0. 6%以内 吻合地非常好 , 在此之后 , ±45°拉伸实验获得的曲 线渐渐落在 V 形开口轨道剪切实验曲线下面 , 只有 BX600层合板的剪应力 -应变曲线整体吻合地较好 。 除剪切强度偏低外 , ±45°拉伸实验获得的极限剪应 变也偏低 。
图 2 典型的 L700层合板的剪应力 -应变 F ig. 2 Typica l shea r stress-strain curve s of L700 lam ina te
图 3 典型的 LT600层合板的剪应力 -应变 F ig. 3 T ypical shear stre ss-strain cu rves o f LT600 lam ina te
3 试件参数及实验
为使 ±45°拉伸实验和 V形开口轨道剪切实验的 对比更具有一般性 , 本文选用 5种 E 玻璃纤维织物 / 乙烯基酯树脂手糊层合板 , 详细参数如表 1所示 。这 5种层合板均为最常见的层合板类型 , 具有一定的代 表性 。乙烯基酯树脂采用上纬公司的 SWANCOR-901 产品 。 5种层合板的纤维体积含量均在 30%左右 。 按照 ASTM D3518和 ASTM D7078制作的标准试件如 图 1所示 , 每组试件各 5个 。
4 结果与分析
5种复合材料层 合板的 增强纤 维铺 层各不 相 同 , 因此它们在 2种面内剪切实验方法下的破坏特
征也具有较大 的差异 。 其面 内剪切性能参 数见表 2。从表 2中可以发现 , ±45°拉伸实验在测试面内 剪切模量时 具有相对较好的精度 。与 V 形开口轨 道剪切实验相比 , 误差均在 7%以内 。 然而在测试 面内剪切强度方 面 , ±45°拉伸实验所得值 往往偏 低 , 最大相差 39%, 最小相差 12%。 以往研究结果 也表明 ±45°拉伸实验在测试剪切强度方面最大误 差可达 30 ~ 40%[ 3] 。 ±45°拉伸实验为测定层合板 面内剪切模量提供了一种简单的方法 , 且测试精度 一般在可接受范围内 。但在使用剪切强度和剪应变 结果时须十分小心 , 这是因为 ±45°拉伸试件是处于 双轴应力状态而不是纯剪切状态 , 且试件自由边附

缝合复合材料层板面内局部纤维弯曲模型及剪切强度预报

缝合复合材料层板面内局部纤维弯曲模型及剪切强度预报

缝合复合材料层板面内局部纤维弯曲模型及剪切强度预报李晨;许希武
【期刊名称】《固体力学学报》
【年(卷),期】2007(28)3
【摘要】提出了面内局部纤维弯曲模型,基于有限元法和周期性边界条件建立了缝合层板面内剪切强度分析方法,采用桥联模型和最大应力判据分析损伤扩展并获得面内剪切强度,预报结果与试验吻合较好,探讨了缝合参数对层合板面内剪切强度的影响规律,结果表明缝合削弱了层合板的面内剪切强度,缝合针距和行距越大对面内剪切强度越有利,较细的缝合线对面内剪切强度有利.
【总页数】6页(P287-292)
【关键词】缝合复合材料层板;面内剪切强度;纤维弯曲;有限元;桥联模型
【作者】李晨;许希武
【作者单位】南京航空航天大学航空宇航学院
【正文语种】中文
【中图分类】TB33
【相关文献】
1.缝合复合材料层合板面内拉伸强度的分析 [J], 徐建新;王春水;李强;蔡宇
2.缝合复合材料层板三维纤维弯曲模型及压缩强度预报 [J], 李晨;许希武
3.T700碳纤维复合材料层合板面内剪切过程的声发射特性 [J], 李伟;王骁;蒋鹏;黄远航;孙贺
4.基于Ls-Dyna软件2种材料模型的碳纤维复合材料层合板面内剪切有限元仿真[J], 孟宪明;钟正;程从前;曹铁山;赵杰;黄亚烽;吴瑶
5.缝合复合材料层板面内局部纤维弯曲模型及刚度预报 [J], 李晨;许希武
因版权原因,仅展示原文概要,查看原文内容请购买。

复合材料层压板疲劳损伤演化模型的综述与评估

复合材料层压板疲劳损伤演化模型的综述与评估

复合材料层压板疲劳损伤演化模型的综述与评估
冯炜森;杨成鹏;贾斐
【期刊名称】《材料导报》
【年(卷),期】2024(38)9
【摘要】疲劳损伤演化模型是材料疲劳性能退化与寿命预测的重要基础。

本文对
复合材料层压板的疲劳损伤演化模型进行了全面综述,针对材料损伤演化三阶段(Ⅰ、Ⅱ、Ⅲ)特征将其磨损模型分为两类:Ⅰ-Ⅱ-Ⅲ三阶段模型和Ⅰ-Ⅱ两阶段模型。

提出了模型评估的五个基本标准,选取了八组不同加载条件下的疲劳损伤演化实验数据,
对表征能力较强的疲劳磨损模型进行了分析评估。

现有研究表明:绝大多数疲劳损
伤演化模型属于宏观唯象模型,模型中引入的疲劳影响因素少,对实验的依赖程度高;对于Ⅰ-Ⅱ两阶段损伤,两阶段磨损模型具有很高的拟合精度;相比而言,Ⅰ-Ⅱ-Ⅲ三
阶段磨损模型可准确模拟疲劳损伤演化全过程,比Ⅰ-Ⅱ两阶段磨损模型更具普适性,其中Mao、吴富强和Shiri的模型应用较为广泛。

为了克服唯象模型的不足,须考
虑不同影响因素下材料的微细观损伤形式和机理,进而发展基于物理机制的疲劳损
伤演化理论和模型。

【总页数】11页(P226-236)
【作者】冯炜森;杨成鹏;贾斐
【作者单位】西北工业大学力学与土木建筑学院;西安电子科技大学机电工程学院【正文语种】中文
【中图分类】TB332
【相关文献】
1.含分层损伤的复合材料层压板后屈曲及低周疲劳分层扩展有限元模拟研究
2.复合材料层压板损伤修复与评估
3.含冲击损伤的T300复合材料层压板的疲劳特性
4.低周疲劳过程损伤变量的复合分析法和三阶段损伤演化模型
5.复合材料层压板冲击损伤的无损检测与评估
因版权原因,仅展示原文概要,查看原文内容请购买。

GLARE层板低速冲击的实验与模拟

GLARE层板低速冲击的实验与模拟

GLARE层板低速冲击的实验与模拟崔旭;李斌;王朔;薛红前;张震宇【摘要】研究玻璃纤维增强铝合金层合板(glass fiber reinforced aluminum laminates,GLARE)在落锤低速冲击下的材料行为,建立ABAQUS有限元模型进行模拟并对其进行实验验证。

针对纤维金属基体材料的特点,采用连续损伤模型(continuous damage model,CDM)分别给予落锤6.22 J、12.38 J和14.46 J的冲击能量,在ABAQUS中对模型设置相应的边界条件和载荷,得出落锤下落方向的速率-时间曲线和能量损耗曲线图。

考虑金属层与复材层间黏结层的作用,采用凝聚层(cohesive)将金属层和复合材料层粘接。

在仿真中观察层间的纤维和基体拉伸和压缩损伤状态及破坏情况,并与实验得出结果进行对比。

结果显示:有限元仿真可以准确模拟落锤冲击之后GLARE板背面的裂纹和鼓包的实效情况以及基体和纤维的损伤情况,很好地预测复合材料内部的损伤情况。

%The material behavior of glass fiber reinforced aluminum laminates (GLARE) under low-speed impact of a falling hammer was studied, and a finite element model of ABAQUS for simulation and experimental verification was built. According to the characteristics of fiber metal matrix materials, the continuous damage mechanics (CDM) model was adopted, and the impact energy values of the falling hammer were set as 6.22 J, 12.38 J and 14.46 J respectively. The corresponding boundary conditions and loads of the model in ABAQUS were confirmed to obtain the velocity-time curve and energy loss curve of the falling direction of the falling hammer. The cohesive layer between metal layer and composite layer was taken into consideration, and the ‘Cohesive Layer’ in ABAQUS was adopted tobond metal layer and composite material layer. The tensile and compressive damage states of fibers and matrix in simulation were observed and compared with the experimental results. The results show that the finite element simulation can accurately simulate the crack and bulge effect on the back of GLARE laminates after the impact, as well asthe damage situation of the matrix and fiber, so as to well predict the internal damage of the composite material.【期刊名称】《航空材料学报》【年(卷),期】2019(039)002【总页数】7页(P68-74)【关键词】纤维金属层板; GLARE层板; 低速冲击; 有限元模拟; 连续损伤模型;【作者】崔旭;李斌;王朔;薛红前;张震宇【作者单位】沈阳航空航天大学,沈阳 110000;沈阳航空航天大学,沈阳 110000;沈阳航空航天大学,沈阳 110000;西北工业大学,西安 710000;沈阳航空航天大学,沈阳 110000【正文语种】中文【中图分类】TB333纤维金属层板(fiber metal laminates,FML)是一种新型的金属基体复合材料,即以金属面板和纤维/树脂复合材料通过胶层黏合在一起[1]。

复合材料层合板拉压和面内剪切性能的分散性实验研究

复合材料层合板拉压和面内剪切性能的分散性实验研究

s a i t a n l ss s o h t 3 p r m e e i u l d s rb to t t i la a y i h ws t a - a a t r we b l it i u i n, No ma it i u i n a d l g o ma s c r l d s rb to n o n r l d s rb to r o d f n to st e c i e t e d s rb to r p r iso h e h n c l h r c e s i it i u i n a e g o u c i n o d s rb h it i u i n p o e t f em c a ia a a t r .F — e t c
( p r m e fEng ne rn e h nis, o t De a t nto i e i g M c a c N r hwe t r se n
Po y e h ia n v r i Xi a 0 ,Ch n ) lt c n c l U i e st y, ’ n 7 7 1 0 2 ia
关键词 : 合材料 ; 学性能 ; 散性 ; 复 力 分 实验 研 究 中 图 分 类 号 : 3 TB 3 2 文献 标 识 码 :A 文 章 编 号 :1 O - 3 1 2 1 ) 60 2 — 6 O 14 8 ( 0 0 0 — 0 00
Ab ta t o i e en o c d T3 0 QY8 1 c m p st a ia e h e sl ,c mp e sa d i- ln sr c :F rf rr if r e 0 / b 9 o o i lm n t ,t e tn i l e e o rs n n pa e
变 曲 线 都 几 乎 为直 线 , 具有 脆 性 破 坏 的特 点 。并 且 利 用 统 计 学 原 理 对 实 验 结 果 进 行 统 计 分 析 , 究 各 力 学 性 能 参 数 的 分 研 布特性 , 果表明 : 参数威布尔分布 、 态分布 、 结 三 正 对数 正 态分 布 能 较 好 地 描 述 各 力 学 性 能 参 数 的分 布 规 律 , 后 给 出 复 最 合 材 料 层合 板 拉 、 、 切 性 能 的 分 散性 参 数 。 压 剪
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

收到初稿日期 : 2 0 1 2-1 2-1 8;收到修改稿日期 : 2 0 1 3-0 3-2 2;网络出版时间 : 2 0 1 3-0 5-1 4 1 2: 1 4: 5 4 / / / 网络出版地址 :www. c n k i . n e t k c m s d e t a i l 1 1. 1 8 0 1. T B. 2 0 1 3 0 5 1 4. 1 2 1 4. 0 0 9. h t m l ) 基金项目 :国家空间碎片专题研究基金项目 ( k s 0 6 2 1 1 j p : 通讯作者 :庞宝君 ,教授 ,博导 ,主要从事材料动力学行为 、空间环境效应及其对策研究 E-m a i l i t . e d u. c n a n b @h p g j
* , , , P ANG B a o u n CHE N Y o n Z HE NG W e i F E NG S h u o j g
( ,H ,H ) S c h o o l o f A s t r o n a u t i c s a r b i n I n s t i t u t e o f T e c h n o l o a r b i n 1 5 0 0 8 0, C h i n a g y : / , l a n e l a s s A b s t r a c t I n o r d e r t o i n v e s t i a t e t h e i n-p s h e a r b e h a v i o u r o f S 2-g e o x c o m o s i t e l a m i n a t e c o n t i n u m g p y p / d a m a e m e c h a n i c s( C DM)m o d e l i n v o l v i n e l a s t i c i t a n d f o r S 2 e o x l a m i n a t e u n d e r i n s h e a r l a s t i c i t l a s s l a n e -g -p g g y p y p y , w a s s t u d i e d i n t h i s i n t e n s i l e s h e a r t e s t s o f t h e c o m o s i t e w e r e c o n d u c t e d t w o t e s o f l o a d i n a e r . T i c a l l a n e -p p y p g p p y p C DM m o d e l s w e r e e s t a b l i s h e d b a s e d o n t h e s t r e s s a s s u m t i o n a n d t h e w e r e i d e n t i f i e d . T h e f o r m e r l a n e a r a m e t e r s p p p C DM m o d e l c o n s i d e r e d n o l a s t i c i t w h i l e t h e l a t e r o n e d i d . B o t h t w o C DM m o d e l s w e r e c o d e d t o VUMA T t h r o u h p y g / , A B AQU S E x l i c i t f i n i t e e l e m e n t a n a l s i s( F E A) f o r i n l a n e s h e a r b e h a v i o u r o f t h e c o m o s i t e w a s c o n d u c t e d b -p p y p y t h e m f o l l o w i n e x e r i m e n t a l t e s t s .T h e r e s u l t s o f F E A w e r e c o m a r e d w i t h e x e r i m e n t s a n d t h e i n f l u e n c e o f g p p p l a s t i c i t l a n e o n t h e i n s h e a r b e h a v i o u r o f t h e c o m o s i t e w a s d i s c u s s e d . E l e m e n t t e s o n t h e r e s u l t s o f F E A w e r e -p p y p y p l a s t i c i t r e d i c t a l s o a n a l z e d . T h e r e s u l t s s h o w t h a t t h e C DM m o d e l c o n s i d e r i n n o c a n e x a c t l t h e f a i l u r e s t r e n t h p y p y g y g / ,b l a s s l a n e r e d i c t o f S 2 e o x c o m o s i t e u n d e r i n-p s h e a r l o a d i n u t i t c a n n o t w e l l t h e n o n l i n e a r m e c h a n i c a l -g p p y p g ,w , l a s t i c i t e r f e c t l r e d i c t b e h a v i o u r .H o w e v e r h e n c o u l e d w i t h i s o t r o i c c r i t e r i a t h e C DM m o d e l c a n t h e p y p y p p p , n o n l i n e a r m e c h a n i c a l b e h a v i o r a n d t h e r e s u l t o f f a i l u r e s t r e n t h i s a l s o a c c e t a b l e c o m a r e d w i t h t h e e x e r i m e n t a l g p p p , r e s u l t s . T h e C DM m o d e l s c a n a l s o b e u s e d f o r F E A s i m u l a t i o n b s h e l l e l e m e n t s b u t t h e r e s u l t s o f l o a d-d e f l e c t i o n y c u r v e a r e a l i t t l e l o w e r t h a n t h a t u s i n l a n e s t r e s s e l e m e n t s b e c a u s e o f t h e t r a n s v e r s e s h e a r b e h a v i o u r o f s h e l l g p , e l e m e n t s .M o r e o v e r r e d u c e d i n t e r a t i o n m e t h o d i s f o u n d t o b e b e n e f i c i a l f o r F E A s i m u l a t i o n . g : / ; ; ; ;VUMA l a s s l a n e l a s t i c i t K e w o r d s e o x i n s h e a r h a r d e n i n c o n t i n u u m d a m a e T g -p p y p y g g y
/ C o n t i n u u m d a m a e m e c h a n i c s m o d e l o f l a s s e o x c o m o s i t e p y p g g
l a m i n a t e l a n e u n d e r i n s h e a r l o a d i n -p g
复 合 材 料 学 报 A c t a M a t e r i a e C o m o s i t a e S i n i c a p
( ) 文章编号 : 1 0 0 0-3 8 5 1 2 0 1 3 0 5-0 2 0 1-0 8
第3 0 卷 第 5 期 1 0月 2 0 1 3年 V o l . 3 0 N o . 5 O c t o b e r 0 1 3 2
/ G l a s s E o x p y 复合材料叠层板面内剪切连续损伤模型
庞宝君 * ,陈 勇 ,郑 伟 ,丰 硕
( ) 哈尔滨工业大学 航天学院 ,哈尔滨 1 5 0 0 8 0
/ 摘 要: 2 玻 璃 纤 维/环 氧 树 脂 ( S 2-G l a s s E o x 2- 为确定 S p y)叠 层 复 合 材 料 面 内 剪 切 应 力 - 应 变 关 系 ,对 S / G l a s s E o x p y 叠层复合材料面内剪切拉伸载荷下的弹 、塑性连续损伤本构模型及应用进 行 了 研 究 。 基 于 平 面 应 力 状态下的连续损伤力学模型 ,通过典型面内剪切拉伸实验 ,分别建立了忽略塑性应变和考虑塑性应 变 的 两 种 连 续 / 损伤力学 ( 模 型 ,并 确 定 相 关 参 数 。 通 过 A C DM) B AQU S E x l i c i t 用 户 子 程 序 VUMA T 接 口 ,分 别 采 用 两 种 p / C DM 模型对 S 2-G l a s s E o x p y 叠层复合材料 证 模 型 可靠性 ,并分析单元类型对有限元计算结果的影响 。 研究结果表 明 :忽 略 塑 性 应 变 的 C DM 模 型 可 以 很 好 地 预 测 复合材料面内剪切失效强度 ,但不能较好地预测其非线性力学响应 ;考虑塑性 应 变 ,将 塑 性 硬 化 与 损 伤 耦 合 后 的 C DM 模型则能较好的预测复合材料非线性力学响应 和 面 内 剪 切 失 效 强 度 ;该 平 面 应 力 状 态 下 建 立 的 C DM 模 型 可用于壳单元进行复合材料有限元数值计算 ,横向剪切作用导致传统壳单元数值计算的载荷位移曲线略低于平面 应力单元计算结果 ;减缩积分算法有利于提高有限元数值计算结果的准确性 。 / ;面内剪切 ;塑性硬化 ;连续损伤 ;用户材料子程序 关键词 : l a s s E o x G p y 中图分类号 : B 3 3 0. 1 文献标志码 : T A
相关文档
最新文档