有限元计算细观力学对复合材料力学行为的数值分析.pdf

1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/第28卷第2期1998年5月25日力学进展ADVANCES I N M ECHAN ICSVol128No12M ay 25,19983国家自然科学基金资助项目收稿日期:1996207213,修回日期:1997211221有限元计算细观力学对复合材料力学行为的数值分析3方岱宁周储伟清华大学工程力学系北京100084摘要有限元计算细观力学的发展是近十年来细观计算力学发展的主要特征和推动力.本文综述了有限元计算细观力学近十年来应用于复合材料力学行为分析研究方面的进展.介绍了基本的数值模型和计算方法,重点评述了强度和损伤等协同效应问题上的最新研究成果.最后对有限元计算细观力学应用于材料设计的前景做了展望.关键词细观力学,有限元,复合材料,刚度,强度,损伤,任意分布,材料设计1引言复合材料的就位特性、各向异性和呈层性所产生的各种复杂的力学现象,使得有限元计算技术对于求解复合材料及其结构的力学问题得到了相当广泛的应用1,2.在这一领域可分为两个分支:一是有限元法应用于复合材料结构(如板、壳等)力学问题1;二是有限元技术应用于复合材料细观力学行为的模拟分析3,4.前者追求真实工程环境下的工程结构问题的解决,后者侧重于材料细观结构与力学性能的关系分析.有限元法与细观力学和材料科学相结合产生了有限元计算细观力学3,4.作为细观计算力学5的最主要的组成部分,有限元计算细观力学的发展一直是近十年来细观计算力学发展的主要特征和推动力.它主要研究组分材料间力的相互作用和定量描述细观结构与性能间的关系.由于复合材料综合了不同单相材料的长处,对其材料力学行为的有意义的研究必须借助于细观力学进行.界面行为,损伤和动态行为对复合材料尤为重要.因此,有限元计算细观力学在求解复合材料细观力学问题中的应用正是在70年代随着细观力学的起飞而发展起来的57.但是,该领域却是在80年代末随着计算材料科学(Computationalmaterials Science)8,9或称计算机辅导材料设计(Computer-A idedM aterialsDesign)兴起而真正得到迅猛发展.这主要由于下述因素促成的:(1)细观力学理论解析的方法,至今还主要限于解决复合材料有效刚度混合效应的问题,尚不能解决与复杂损伤强度相关的协同效应、非比例加载响应和其有尖棱角(非旋转体)增强相的细观结构等问题3,10;(2)复合材料在力学加载下的细观结构信息不可能在实验中以系统的方法获得3,10;(3)超级计算机的发展和有限元计算371 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/软件的商业化,基本克服了有限元细观计算力学的最大缺点-输入数据工作量大和花费比较长的计算机时;(4)最重要的还是在于有限元细观计算力学方法能够描述复合材料的细观结构对宏观响应的影响的关系,使得特别设计的细观结构对载荷是如何响应和如何失效的问题可以进行数值模拟3.有限元细观计算力学的最大优点在于它能够获得纤维(或颗粒)直径尺度下的完整的应力2应变场来反映复合材料宏观应力-应变响应特征.这样,它能够分析宏观有效性能对细观结构的依赖关系.例如:能定量描述诸如纤维(或颗粒)的形状、尺寸、分布和体积含量等这些细观结构参量对宏观力学性能的影响.而这些优点正是计算材料科学在材料细观结构设计时所必需的11.在复合材料结构设计中,可以控制界面条件,纤维-基体的排列方式,颗粒(纤维)的形状和尺寸,这样就可能修改其强度和其他有关的力学性质,满足指定的功能要求.这种在计算机指导下设计具有特殊性能的复合材料细观结构的要求给有限元计算细观力学发展提供了机会和挑战.本文首先介绍了主要的材料模型和基本的计算方法,重点评述了有限元计算细观力学对复合材料力学行为中较为困难的塑性流动、强度及损伤协同效应问题和增强相任意分布问题数值模拟求解方面的进展,讨论了复合材料细观结构对宏观性能的影响,最后对有限元计算细观力学应用于材料设计做了展望.2材料模型和计算方法211材料模型有限元细观计算力学应用于复合材料力学行为数值模拟的本质,是将有限元计算技术与细观力学和材料学相结合,根据复合材料具体细观结构,建立代表性细观计算体元、界面条件和边界条件,求解受载下体元中具有夹杂的边值问题.从而建立起细观局部场量与宏观平均场量间的关系,最终获得复合材料的宏观力学响应.虽然现代计算技术能容易地获得全场解,但在有限元计算细观力学中,大多数数值解是应用了较理想的增强相周期分布的材料模型.代表性体元的材料模型应满足:(1)相对于细观分析的合适尺度,即基体中的增强相尺寸和增强相间的平均间距要大于细观结构的特征尺寸(如,晶界尺寸和位错运动距离);(2)反映细观结构的几何形状、分布和界面条件.根据不同复合材料的细观结构,设计出的有代表性的计算体元.其材料模型包括如下细观结构参量:(1)增强相2基体的分布:(i)纤维或颗粒的周期排列分布,分为立方、六角、长方等分布;(ii)纤维或颗粒的任意分布,又分为具有平移对称性的任意分布和一般任意分布.(2)增强相的形状(旋转体和非旋转体):方形、圆形、椭圆形或任意多边形.(3)增强相的不同取向:颗粒的取向,为0,?6,?4,?3,?2.(4)界面条件:理想粘合,或局部脱粘.(5)组分材料的特性:增强相和基体的刚性,弹性,粘性,塑性等.(6)增强相的尺寸:大小混杂分布.(7)增强相的含量(对于体积百分比,颗粒 50,纤维 80):增强相分布的不同可反映在体积百分比的计算上.例如,半径为r长度为2r的短圆柱纤维,其体积百分比为12,13f=l(r?R)3(1)这里2R是体元的长度.对立方体元,l=p?4;六方体元,l=23p?9.8)为分析上述增强相几何因素对材料性能的影响,对三维问题分别引入了颗粒长径比ap和体元长宽比acap=l?dac=H?w(2)471 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/这里l是颗粒在加载方向的长度,d是颗粒在垂直于加载方向上的直径或宽度,H是体元的长度,w是体元的宽度.对二维平面问题,引入了纤维横观长宽比 f和体元横观长宽比 c14,15f=l?wc=L?W(3)这里w是纤维在横观平面上垂直于加载方向上的直径或宽度,l是纤维在横观平面上沿加载方向上的长度,W,L分别为体元在横观平面上垂直于加载方向上的直径或宽度和沿加载方向上的长度.212计算方法这里仅简单地介绍根据细观力学建立的细观局部场量和宏观平均场量的关系,基本的数学公式和数值实现方法.这包括周期性条件,位移加载条件,本构关系,以及各种体单元的边界条件的实施.在有限元计算细观力学中,假设连续性本构关系能够合适地描述所有相的机械性能,在初值和边界条件下通过对控制方程和能量平衡的求解来确定宏观变形史.并认为加载速率相对较低,这样材料惯性和热梯度可忽略.对周期分布的增强相,其弹性和柔度张量C和S,是位置的周期函数,周期解可表达为(x)=(x+d),(x)=(x+d)(4)d=3i=12miaiei(5)mi(i=1,2,3)为任意整数,a=(aiei)是体元边矢.若平均应变、应力表达为 和,则在给定边界位移条件u09V=0 x下有1215:=1VVV(x)dV=1V9V12(nu0+u0n)ds=mnm nijdVm n(6)=1VV(x)dV=1VVC(x):(x)dV=mnm nijdVm n(7)n是对高斯点数求和,m对单元数求和.有效弹性模量可由各向异性弹性本构关系确定=C3:,=S3:(8)这里C3、S3分别是有效弹性和柔度张量.为计算方便,它们可由工程常数表示为66矩阵,如:S3ij=1E311-v321E322-v331E3330001E322-v332E3330001E333000symm1G323001G31301G312(9)571 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/E3为有效杨氏模量,v3为有效泊松比,G3为有效剪切模量.若基体材料为率无关的延性材料,如金属基复合材料,计算中多采用J2流动理论.为应力偏量,e=(3:?2)1?2为有效应力,应力增量和应变增量的关系为=E1+vI+1-2-32ag:(10)g=1?1+23(1+)EEt-1-12e(11)为Kronecker delta,I四阶等同张量,Et是当前应力状态下的割线模量,是基体的泊松比.在单轴加载条件下,应力2应变可用Ramberg and O sgood关系式=E+0E0n(12)0是单轴屈服应力,n为应变硬化指数,可取为3?7.图1两种不同分布的三维体元示意图在给定位移边界条件u09V=0 x下,周期单胞体元的边界条件必须满足连续性条件和周期条件.对二维平面问题和轴对称问题,N eedleman等16,17介绍了完整的边界条件,本文作者1215和其他人18,19给出了三维单胞体元的边界条件.所有的边界条件的数值实施都采用了671 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/限制单元技术.例如,对立方或长方体元而言,体元的外表面应满足对称性条件和周期条件,六方体元的边界条件较复杂,图1显示的是典型的四分之一的六方体元,该体元平面变形限制为变形后仍保持为平面,而斜面上的位移u1,u2被限制于反对称斜面的中线,这样斜面变形后仍保持为斜面1215.3刚度问题在70年代,有限元计算细观力学随着细观力学的起飞而得到了实质性的发展,主要应用于求解含玻璃增强相的树脂基复合材料的有效模量2022.当时的计算模型还需要应用叠加法来求解弹性边值问题21,22.随后十来年,由于细观力学的近似分析方法(如自洽方法、广义自洽法、积分法、M ori2Tanaka平均场法等)的发展和完善,使得有限元计算细观力学在求解复合材料刚度问题方面的进展有限.直到80年代末,由于先进的金属基和陶瓷基复合材料的广泛应用和计算机技术的发展,人们可利用有限元计算细观力学来获得复合材料的细观结构与宏观有效模量之间的定量关系12,23,例如,Brown等人23分析了纤维截面形状对金属基复合材料有效模量的影响,方岱宁等14计算了正交各向异性复合材料的九个柔度常数14,分析了长纤维形状、分布对有效剪切模量和杨氏模量的影响.方岱宁等12还数值模拟了颗粒形状、取向、分布对有效性能的影响.所获得的主要结论是:(1)连续纤维的分布和形状对纤维轴向刚度没有影响,而对横截面上杨氏模量有影响,所有的有效剪切模量和泊松比都依赖于纤维的形状和分布形式;(2)颗粒的取向、形状和分布对所有的有效模量都有影响;(3)增强相的取向、形状和分布对有效模图2不同模型计算的有效轴向弹性模量结果与所有结果的比较量的影响是随着增强相的增加而增大.图2是用轴对称单胞和三维体单胞模型对含不同分布的玻璃球颗粒增强的高分子基复合材料轴向刚度的计算结果和实验结果的对比12,结论是六方分布的计算结果和实验值吻合的最好.国内应用有限元计算细观力学求解刚度问题上有创意的工作还有杨庆生等人2426发展的自洽有限元法.对混合效应的刚度问题,求解的方法很多,也相对成熟,故在本文中不作为介绍的重点.需要指出是有限元计算细观力学能够获得材料细观结构与有效模量间的定量关系,这是细观力学理论分析方法所做不到的.771 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/4强度问题411塑性流动特征对纤维增强复合材料的屈服进行有限元计算细观模型分析也是起始于70年代2730.80年代N eedleman和Tvergaard等人将孔隙材料的轴对称细观模型运用到金属基复合材料3133.同时,三维单胞体元模型也开始应用于复合材料抗塑性流动的数值分析中3438.对塑性问题,考虑到复合材料韧性一般较低,认为其塑性流动仍在小变形范围.大多数工作都是基于连续性本构关系适用于各相组分材料这一假设,而应用J2流动理论.Devorak等人在纤维增强复合材料的塑性流动特征模拟方面做了许多工作3943.对所报道的复合材料抗塑性流动的有限元数值模拟工作,可总结为如下几个方面:(1)复合材料的单向拉伸变形行为和塑性流动特征模拟3138,获得宏观应力-应变曲线、应变硬化率以及流动强度.(2)复合材料在多轴载荷下屈服面的数值模拟2830,38,39,获得三维模型的屈服面变化形状,所得结论是:复合材料屈服对静水压力是敏感的29,30.(3)复合材料的循环塑性变形的数值模拟44,45,N eedleman等人44数值模拟了短纤维增强金属基复合材料在循环载荷下的Bauschinger效应,数值结果显示出:在正负加载过程中基体中受约束的塑性流动应力的水平不同是造成Bauschinger效应的主要原因.进一步地,N eedleman等人45还数值模拟了循环变形下的循环硬化和饱和过程.L evy等人36数值分析了金属基复合材料在热循环载荷下的塑性流动特征.图3 颗粒形状对金属基复合材料弹塑性变形的影响412细观结构的影响在复合材料的细观结构对宏观弹塑性变形影响的模拟和数值评价方面已做了许多的工作1215,3134,37,38,4648.分析的重点是增强相(颗粒或纤维)的几何参数(如形状、分布、取向、尺寸、聚集度等)对复合材料的宏观塑性流动的影响:(1)增强相的形状对复合材料塑性流动的影响的数值模拟13,31,38,4648.Bao等人38和方岱宁等人13的结果都表明颗粒或短纤维形状的影响可由其长径比来反映.图3显示了在给定增强871 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/相体积分率下含4种不同形状(球形、立方形、短柱、长方形、椭圆形)的陶瓷颗粒增强金属基复合材料的单向拉伸应力-应变曲线13.流动应力和硬化速率不仅和颗粒沿拉伸方向的长度(即方面比)有关,而且和颗粒形状有关,即具有尖棱角的颗粒有更高的基体流动约束而导致更高的流动应力和硬化率.对长纤维形状的计算结果38表明,长纤维形状不影响沿纤维轴向拉伸的塑性流动,而对横观截面上的拉伸强度有较大的影响.(2)颗粒或短纤维的取向对复合材料塑性流动影响的数值模拟必须用三维单胞体元模型13,47.Sorensen等人47研究了周期分布的短纤维取向(取向角:0?4)对单向拉伸应力2应变响应的影响,计算结果显示,当取向角较小时(即纤维偏离拉伸方向的角度较小时),拉伸流动应力水平要小于直线排列纤维的应力水平,而当取向角足够大时,拉伸应力2应变曲线不敏感于取向角.方岱宁13分析了非旋转颗粒(立方和长方形)的取向(取向角:0?2)对单向拉伸应力-应变响应的影响,注意到 的取值范围是Sorensen的两倍.图4为长方形陶瓷颗粒的取向角对金属基复合材料拉伸应力2应变响应的结果.当=90 时,沿拉伸方向的颗粒方面比最大,所以造成最大的流动应力水平和硬化率.图4长方形颗粒的取向对金属基复合材料弹塑性变形的影响(3)纤维和颗粒分布和排列方式对复合材料的横观截面的拉伸塑性流动的影响也是相当显著的14,4850.对不同分布的颗粒增强复合材料在横观截面上单向拉伸下的弹塑性变形数值模拟结果示于图5,包括了立方对角、长方分布、周期六方和立方分布.当沿颗粒轴向拉伸时,应力-应变曲线不敏感于颗粒的分布,而在横观截面单向拉伸时,图5显示不同分布之间的差别是明显的,并且这种差别是随着体积分率的增加而增加.(4)一些研究者3134对短纤维的聚集度和短纤维排列端部的影响进行了数值分析.通常是采用纤维周期对差排列的轴对称模型来模拟短纤维端部间距的影响31,32,也有用三维单胞体元模型来分析的34.获得的结论是:(i)纤维长径比越小,即纤维越粗,则纤维间距愈小,这样沿纤维方向的拉伸流动应力水平也愈小,而沿横向拉伸塑性流动则愈强;(ii)短纤维差排的端部影响是使得沿纤维方向的拉伸流动应力水平要比纤维对中直线排列的塑性流动应力低.971 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/图5不同分布对在横观截面上单向拉伸应力2应变曲线的影响413残余应力对塑性流动的影响在热机械过程中,延性基体的脆性增强相能够显著地影响到材料细观结构的演化.复合材料在制备过程中,从熔解温度冷却到室温,由于增强相和基体的热膨胀系数不同,而产生残余应力,造成复合材料基体局部塑性屈服,所产生的残余应力对材料受载的塑性流动将有较大的影响.有限元计算细观力学能够求解和模拟变化的温度场下的残余应力场35,37,49,5154.Davis37模拟了复合材料制备冷却过程中压缩残余应力的形成过程,并分析了该压缩残余应力所造成的塑性屈服现象.W eissenbek49分析了含不同分布的纤维周期复合材料的拉伸平均应变与温度之间的关系.Taggart51,Povirk52,L evy35,Tvergaard53和Davis37等分别数值分析了复合材料制备过程中所产生残余应力对复合材料抗塑性流动的影响.数值结果表明,压缩残余应力导致在拉伸和压缩过程中不同的载荷传递机理,压缩残余应力造成拉伸时的较高效率的载荷传递,这样拉伸流动应力大于压缩流动应力.另一方面,随着残余应力水平的增加,拉伸和压缩下的硬化速率都随着增加.Povirk52计算了不同温度下的拉伸应力2应变曲线,计算结果显示出,温度愈高,流动应力水平愈低.他还指出:纤维间距对残余应力分布和界面与纤维端部的应力影响较大.L evy35用三维模型分析了由于热处理而引起的残余应力对金属基复合材料的塑性流动的影响.从模型中所获得的结论是,残余应力是非均匀分布的,小的拉伸残余应力会提高压缩应力2应变曲线.5损伤问题511细观损伤模拟复合材料在制造、加工过程中,不可避免地会出现一些细观缺陷,例如纤维、颗粒等增强相的断裂、它们与基体的脱粘以及基体中产生微空洞微裂纹等.在工作过程中,伴随着新的缺陷不断地萌生,原有的缺陷进一步发展,使材料的性能劣化.这类细观损伤大体可分为三类:颗粒、纤维等增强相的断裂;增强相与基体材料脱粘甚至完全剥离或拔出;基体中空洞形核,聚合081 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/而产生韧性破坏.前两种损伤是空洞形核的主要诱因.由于上述复合材料细观损伤的产生及发展强烈地依赖于夹杂界面、空洞附近及裂纹尖端的局部应力应变场,而这些局部场的变化一般很复杂,必须应用有限元计算细观力学来对细观损伤进行数值分析,才可以克服求解析解时作出的过多的简化,使结果更加接近实际.纤维、颗粒等增强相的断裂,必然引起材料力学性质的改变.文献55引入一种三相损伤单元来研究颗粒断裂对材料整体力学性质的影响.这种单元计及了完好的与断裂的两种颗粒在复合材料中所起的作用,数值计算结果表明,断裂颗粒数较少时,模型是合理的,但断裂颗粒数较多或颗粒上有多重断裂时,怎样确定它们的空间分布有困难.而且这种模型没有涉及同其他损伤的相互作用,也没有涉及损伤的演化过程.颗粒、纤维等与基体在界面上的脱粘也是一种典型的损伤.有的工作将界面层假设为一层弹性介质56,或假设为沿界面法向和切向分布的两组弹簧57,58,这类模型数值计算较易实现,但需给出界面层的模量和厚度等参数,也未能很好地反映损伤演化.有些复合材料的界面层很薄,且无明显的分层,模量也难以测定,上述模型的应用受到了限制.文献59在有限元数值计算中利用打开相邻的基体单元和纤维单元之间的结点来模拟脱粘,可以将脱粘同其他类型的损伤用统一的数值方法处理,但这种模型尚不能较精确地反映界面损伤.文献16,60,6164利用了基于内聚力的唯象模型,此模型是无厚的,而且法向和切向的损伤是耦合的,能较好地反映一类复合材料的界面损伤.文献16利用此模型在统一的计算模式下数值模拟了脱粘萌生、发展直至完全剥离并引起空洞形核的全过程,并揭示出界面层的相对厚度与脱粘的韧性和脆性的关系.但此界面层模型是建立在连续介质基础上的,难于反映界面脱开后再接触时的一些接触性质.脱粘界面的裂纹尖端的高应力梯度也难以反映出来.界面脱粘实际上是一种断裂扩展,文献65建议了一种从断裂角度处理脱粘的界面离散单元,将起裂准则及裂纹扩展准则引入界面层,可以处理I型,II型及混合型断裂(界面脱粘不一定沿主应力方向,可能属于混合型断裂),并能解决上述脱开界面再接触的问题.纤维从基体中拔出也是一种界面损伤,一种基于率相关的内变量摩擦理论并同有限元结点相联系的模型66,67可以模拟实验观察到的结果68.实际纤维的拔出,应同脱粘及界面附近的基体材料延性破坏有关,这种模型没有考虑这些因素.颗粒、纤维的断裂与界面的脱粘这两种损伤是相互影响的.每一种损伤还同颗粒等夹杂的形状,体积比,远场应力状态等参数有关16,55,61.文献61讨论了上述各种因素对损伤演化的影响,并揭示空洞萌生、发展与上述损伤的联系.但文献61限于讨论颗粒界面在初始状态为完好的情况,若有初始缺陷存在,则演化过程将更加复杂.金属基复合材料的基体在高温下表现出很强的粘弹塑性质.文献62研究了热残余应力对起始脱粘的作用以及基体粘性对材料韧性的改善.空洞在复合材料延性的基体材料中萌生,扩展和聚合而导致破损是复合材料一种主要的延性损伤过程.一种基于GURSON理论的含有微空洞的损伤本构单元可以描述这种损伤过程6972.这种模型可以形象地描述延性损伤积累直至断裂的全过程.由于在有限元计算细观力学分析中引入了微空洞平均间距及微空洞体积比两个参数,能够反映出损伤的细观结构.值得指出的是,Xia和Shih6972发展了孔洞单胞破坏元来数值模拟裂纹前端的孔洞损伤演化.每一破坏单胞元含一个孔洞,随着孔洞损伤的演化,破坏单胞内的孔洞将逐步连通,并和主裂纹连通而造成最终的破断.国内在对复合材料损伤数值模拟方面有代表性的工作应为陈陆平、潘敬哲和钱令希等人的创造性工作57,58,73.他们利用参数拟规划化和非连续线弹性本构模型,构造了复合材料纤维?基体界面失效问题的细观力学模型并进行了有限元细观计算分析,得到了与实验相吻合的结果.细观损伤力学,正是由于其能形象、具体地反映材料细观损伤的物理涵义及过程而备受关注.有限181 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/元数值方法与细观损伤理论相结合,能够模拟更复杂的破坏过程,得到更精确的结果.512蠕变问题对复合材料损伤的数值研究的一个很重要的方面是蠕变问题.在航天航空领域,对材料的要求十分苛刻,除保证构件力学性能外,还要求重量轻,耐高温和抗冲击性好.加上构件的尺度要求,材料的结构效率问题显得更加突出.树脂基复合材料的性能已不能满足这方面的要求(其使用温度低于350).近年来国际正在迅速研究开发适用于350 1200使用的各种高性能增强体增强的各种耐高温复合材料.在高温工作环境下,复合材料损伤破坏的主要形式是蠕变破坏.进入90年代以来,有限元计算细观力学才开始应用于复合材料蠕变损伤问题求解7482,特别是最近两年,蠕变损伤数值模拟的研究不仅仅针对金属基复合材料,也针对陶瓷基复合材料.由于增强相和基体的蠕变速率不同而造成两相间发生应力传递77,增强相断裂和基体中裂纹扩展伴随着两相界面的滑移和脱粘75,76.这些现象形成复合材料蠕变的三种基本的材料模型74:(1)基体蠕变,但增强相为弹性;(2)增强相蠕变而基体为弹性;(3)基体和增强相均发生蠕变.对陶瓷基复合材料而言,大多数蠕变损伤现象是属于材料模型(2)74,而金属基复合材料则为材料模型(1)或模型(3)77,78.在复合材料蠕变问题的有限元计算细观力学分析中,大多数模型为周期分布下的一维蠕变本构和稳态蠕变77,78.对应材料模型(2)的陶瓷基复合材料,Holmes等人75,76的数值结果显示,最初蠕变使得应力由纤维传递给基体,造成基体中裂纹形成,而基体中这种裂纹扩展损伤又反过来使得应力集中回到纤维,最终导致进入瞬态蠕变阶段和断裂.Begley74的模型不仅考虑了纤维蠕变,而且用端部体元分析了两相界面的摩擦滑移和脱粘.他的结果表明,基体裂纹的非桥联区存在,使得在非桥联区和桥联区之间的过渡区中的纤维形成应变集中,而导致纤维蠕变断裂和基体裂纹的加速扩展.Sofronis81分析了颗粒增强金属基复合材料的蠕变抗力,指出只有当复合材料的约束保持住,高温下的复合材料的蠕变强度才能存在.他提出的模型是分析球颗粒界面扩散和滑移对蠕变强度的影响,结论是两相界面间的质量传递扩散和界面滑移能够减低复合材料的约束,并松弛基体应力,从而导致材料的蠕变抗力降低,数值结果显示,严重时可使复合材料的蠕变抗力甚至低于基体材料的蠕变强度.Kim等人80应用与Sofronis相同的模型研究了短纤维增强金属基复合材料的蠕变抗力,所示的结果是相同的.需要指出的是,在数值模拟中,Kim和Sofronis所用的都是基体服从一维指数蠕变率,而假设增强相是刚性的.而A ravas等人的工作82发展了三维蠕变本构关系,并应用于数值计算中.他们的数值模拟模型对应的是材料模型(3),即纤维和基体都发生蠕变,计算结果表明数值模拟结果和他们提出的三维理论模型结果吻合的较好.在复合材料蠕变分析中的另一值得注意的问题是延性基体中孔洞损伤演化8388.在这种情况下,考虑到孔洞损伤演化的本构模型应引入蠕变数值分析中85.6任意分布问题在有限元计算细观力学中,大多数的计算是基于理想化了的周期分布单胞体元模型.复合材料显微观察结果显示,增强相实际上是任意分布的48.显然,任意分布的增强相给细观数值求解边值问题带来了很大的困难.事实上,第二相粒子的空间分布问题长期以来一直是材料科学感兴趣的问题89.因为二相夹杂的分布是材料细观结构的一个重要特征.分布的定量描述对理解细观结构与宏观性能之间的关系是十分重要的.虽然许多研究者从理论和计算上已做了很大的努力,但在这一领域的研究进展是有限的.这里,对有限的计算技术成果进行简单的介绍.在对二相材料的任意分布进行有限元计算细观力学模拟中,一般可分为两种情况,即准任意分281 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/布和完全任意分布.(1)准任意分布又称为周期任意分布.图6是典型的周期任意分布例子,该分布应用了晶体学中的平移对称性原理.含任意个夹杂的代表性体元内的每一个夹杂的位置是可以任意指定的,但每一夹杂具有周期重复的中心,即代表性体元是周期重复的.利用这种准任意分布的模型,N akamura48,方岱宁15,N eedleman90等分别模拟了纤维和孔洞的不同的任意分布对复合材料有效模量和抗塑性流动能力的影响.数值分析的结果表明,对柱形纤维而言,准任意分布的复合材料流动应力水平比立方周期排列的分别地要低,而比六方周期排列的要高.另一方面,孔隙材料的应力2应变曲线依赖于孔洞的指定的任意分布91.图6典型的夹杂准任意分布示意图93图7典型的夹杂完全任意分布的网格离散化示意图93(2)对于增强相完全任意分布的问题,数值方法中必须引入统计技术来描述89,典型的方法是称为格形化技术(tessellation technique)91,这是一种具有统计意义的离散技术.W andow ski等92用这种方法与多边形技术相结合来模拟分析二维任意分布的问题和颗粒聚集度的影响.Cruz等人93也利用格形化技术发展了一种称为平行M onte2Carlo的有限元方法.图7是典型的格形化后的网格示意图.这里应指出两点:(i)为保证一定的计算精度,任意分布的计算量是较大的;(ii)每一多边形必须由格形化技术来保证满足边界连续性条件.Zhang等人94发展了一种杂交有限元方法来分析任意分布问题.类似地,Ghosh等95也提出了一种n边多边形网格有限元方法来分析多相非均匀材料.他们的共同点都是:每一多边形只含有一个夹杂,并看作为一个单元.可想而知,其计算精度是十分有限的.通过对上述报道的任意分布研究工作进行评述,我们可以获得初步的结论是:(1)所有的方法仍然只能解决二维任意分布问题.虽然从理论上讲格形化技术可以解决三381 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/维任意分布问题,但至今尚未见到有关对空间任意分布的数值模拟的成果报道.(2)为保证一定的计算精度,任意分布问题的计算量是比较可观的.(3)严格地说,复合材料的增强相任意分布数值模拟问题仍是有待解决和完善的一个重要问题,目前所报道的进展仍然是十分有限的.7应用于复合材料设计前景展望今天已进入材料设计时代,计算材料科学或称计算机辅助材料设计是近年来一个迅猛发展起来的多学科交叉的新兴研究领域8.它吸收和整合不同的学科知识和方法,对材料进行计算机辅助设计8,9.美国国家科学基金(N SF)近年来对计算材料科学研究给予了强力度的资助9.我国目前也重视到这一领域的发展96.复合材料的设计涉及:增强相和基材的组分,界面粘结,细观结构,成型工艺,性能测定等,与化学、物理、力学、材料学等多种学科有着广泛的内在联系,并互相渗透和互相推动,形成材料科学和工程.然而,上面评述的大多数的工作多是关于用有限元计算细观力学对复合材料的单类结构性能问题进行数值模拟分析,仅有的优化设计报道的是通过最小残余应力设计来得到夹杂的形状优化和界面性能优化56,97.由于高性能材料和技术属于高技术领域,各国对先进复合材料设计和制备技术进行严格保密,因此对结构复合材料在计算机辅助下进行结构优化设计制备的报道和信息却鲜见.在把物理、力学和材料科学的宏观定量化方法与复合材料细微观结构设计相结合的基础上,材料学工作者和力学工作者可不断修改复合材料设计以达到在复杂使用环境下最优的材料表现行为.这一进展将使复合材料细微观结构设计从定性走入定量,并得以从多层次角度研究损伤和断裂机理.结构复合材料设计可有几个层次:(1)传统的配方型经验设计;(2)强韧化原理指导下的定性设计;(3)应用专家系统和计算机模拟进行半定量设计;(4)材料、力学和计算机相结合的定量设计;(5)从颗粒甚至原子层次计算材料力学性质的定量设计.在这几种材料设计方法中,目前最优先和最重视发展的还是计算机辅助下材料定量化设计,原因在于计算机辅助材料设计的细观结构计算设计方法能够描述复合材料的细观结构对宏观响应的影响之间的定量关系,这使得对特别设计的细观结构对载荷是如何响应和如何失效问题可以进行数值模拟.在复合材料结构设计中,可以控制界面条件,纤维2基体的排列方式,颗粒(纤维)的形状和尺寸.这样就可能修改其强度和其他有关的力学性质,满足指定的功能要求.这种计算机指导设计具有特殊性能的复合材料细观结构的要求给材料定量设计发展提供了机会和挑战.在计算机辅助复合材料设计中,必须应用有限元计算细观力学来对复合材料的细观结构进行计算机模拟设计,即计算机辅助材料设计主要强调复合材料细观结构与宏观性能的定量关系.将材料细观结构参数(如组分材料和界面的刚度、强度、韧度等物理参数和增强相的分布、形状、尺寸、取向、体积分率、长径比等几何参数)作为设计变量,而将复合材料宏观性能作为设计目标.这样,建立细观结构与性能间关系;建立细观局部场和宏观平均场间关系;研究组分材料在复合状态下的结构力学行为以及它们对复合材料宏观性能的影响.设计原理是:用材料物理学和细观力学建立材料计算模型、边界条件和界面条件.用计算机进行二维和三维位移场和温度场的计算,再用复合材料细观力学进行后处理分析,获得所设计的目标和最优化有效机械性能.因此,有限元计算细观力学作为计算机辅助材料设计的基本方法是具有很大的应用和发展前景的.参考文献1Griffin O H.Evaluation of finite element software packages for stress analysis of lam inated composites.Compos T ech R ev,1982,4:1361412Noor A K,Burton W S.A ssessment of computationalmodels for multilayered composite shells.App l M ech R ev,1990,481 1994-2010 China Academic Journal Electronic Publishing House.All rights reserved.http:/43:67973Daw son P R,Needleman A.Issues in the finite elementmodeling of polyphases plasticity.M ater S ci Eng,A,1994,175:43484BaskesM I,Hoagland R G,Needleman A,M ater.Summary report:computational issues in the mechanical behavior ofmetals and intermetallics.S ci and Eng,A,1992,159:1345杨卫.细观力学和细观损伤力学.力学进展,1992,22:196Ko W L.Finite element m icroscopic stress analysis of cracked composite system s.J ComposM ater,1978,12:971157Haritons G K,Hager J W,Amos A K,Salkind M J.M esomechanic:the m icrostructure2mechanics connection.Int JS olids S truc,1988,24:108110968Landman U.On some issues in computationalmaterials science.ComputationalM ater S ci,1994,2:2092119Taggart G B.Computationalmaterials research at national science foundation.ComputationalM ater S ci,1994,2:14314810Brockenbrough J R,Suresh S,W ienecke,H A.Deformation ofmetal2matrix compositesw ith continuous fibers:geometri2cal effects of fiber distribution and shape.A cta M etallM ater,1991,39:73575211杨庆生.清华大学博士后出站