摘要
本文简单总结了损伤力学产生的背景及其基本理论,介绍了几种混凝土的局部损伤模型和非局部损伤模型,并分析了这两类模型的不同点及应用前景。为了克服局部损伤模型在数值模拟中存在的网格敏感性与病态局部化问题,在推导非局部模型的基本公式和控制方程时考虑了材料的内部尺度值,并引入了高阶变形梯度项,推导出了隐式梯度模型的基本公式及其有限元离散的控制方程,并将其有效应用于有限元求解程序中。文章重点讨论了求解过程中插值形函数的选取原则,提出了含裂缝混凝土板裂尖应变场奇异性的处理方法,即用线性外插法求解裂尖应变场。编制了有限元程序,并用该程序对平面应力状态下含裂缝板件的损伤和受集中载荷简支梁的三维损伤进行了数值模拟。分析了上两类构件在三种网格下的载荷位移曲线、及损伤分布云图和材料的损伤演变曲线。结果表明,隐式梯度模型可以很好地模拟混凝土材料的三维损伤特性,并且其损伤分布和承载力峰值对于网格的细化不敏感,解决了网格依赖性和零能量损耗问题。
In this paper, the procreant background and basic theory of damage mechanics was summarized, several kinds of local damage model and nonlocal damage model of concrete were introduced, and their differential and applied foreground was analysed. The micro structure criterion was considered and the high-order distor- tion grads items was introduced when the basic formula and manipulative equation of nonlocal damage model was deducted so as to overcome mesh-dependence and sick localization in numerical simulation. The basic formula and manipulative equation of implicit gradient damage model were deducted, and were applied to finite element computation program. The defining principle of the interpolations functions was discussed emphatically in the computation process in this paper, the singularity for the strain field at the crack tip of concrete board with crack was studied,and the numerical extrapolated method was applied at the crack tip. The damage of concrete plate with crack of plane-stress and three-dimensional dam- age of three curved beems were simulated. The curve of load-displacement and nephogram of damage distributing with three mesh were gained, and the curve of damage evolvement was analysed. Numerical calculation show that the implicit gradient model can simulate commendably the speciality of damage for concrete,and its distributing of damage and carrying capacity is no sensitivity to mesh- diminished. The pathological mesh-dependence and zero energy consuming are overcome.
引文
[1].Kaplan F M.J.Amer.Concrete Inst. 1961(58):591-610.
[2].尹双增,断裂、损伤理论及应用[M].北京:清华大学出版社,1993.
[3].谢和平.岩石、混凝土损伤力学[M].徐州:中国矿业大学出版社,1990.
[4].Kachanov L.M.Time of the rapture process under creep conditions [J].TVZ Akad.Nank.S.S.R.Otd.Tech.Nuak.1958,8.
[5].Rabonov Y N.Creep rapture[C],Proc 12,lenter.Congress,AppI.Mech. Stanford,Springer berlin,1969.
[6].Leckie F A,Hayhurst D R.Creep rupture of structure[C]. Proc.R.Soc. A340,1974.
[7].Hult J.Damagein induced tensile instability[C].Trans.3rd AMIRT, London,1974.
[8].Lemaitre J.Application of damage concepts to predict creep-fatigue failure[C].J.Eng.Mat.Tech,ASME,1979,101(1):202-209.
[9].Kajcionovic D.Continuum damage theory of brittle materials[J]. J.App1.Mech.1981.(48):809-824.
[10].Sidoroff F.Discription of anistropic damage application to elasti- city[C].IUTAM Colloquium, physical Nonlinearities in structural analysis,1981:237-244.
[11].Lemaitre J.预估结构中的塑性破坏或蠕变疲劳破坏的损伤模型,余天庆译,固体力学学报,1981,No.4.
[12].宋玉普.多种混凝土材料的本构关系和破坏准则[M].中国水利水电出版社.2002.
[13].李灏.损伤力学基础[M].山东:山东科学技术出版社,1992.
[14].余天庆,钱济成.损伤理论及其应用[M].北京:国防工业出版社,1993.
[15].王军.损伤力学的理论与应用[M].北京:科学出版社,1997.
[16].Dougill J W.et al.J.Eng[J].Mech.Div,ASCE,1976(102):333.
[17].Mazars J. Application de la mecanigue de Pendommagement au comportement non lineaire et a la rupture du beton de structure. These de Doctorat d'Etat.Univ. paris VI,1984.
[18].Loland K E.Concrete damage model for load-response estimation of Concrete[J].Cement and Concrete Research,1980(10):392-492.
[19].余天庆.混凝土的分段线性损伤模型.岩石、混凝土断裂与强度,1985 (2): 14-16
[20].钱济成,周建方.混凝土两种损伤模型及其应用[J]河海大学学报.1989(3) :40-47.
[21].Manson S.S,et al.Re-examination of cumulative analysis-Anenginee- ring prospectinve[J].Engng.Frac.Mech,1986.25-26.
[22].Supartono F,Sidoroff F.Anisotropic damage modeling for brittle elastic materials[C].Symposium of Franc-poland,1984.
[23].何明,徐道远.混凝土的损伤模型.福州大学学报1994(4).
[24].李兆霞,Z.Mroz.应变率相关的粘塑性损伤本构模型.岩石,混凝土裂纹和强度,学术会议论文集,国防科技大学出版社,1993.
[25].李兆霞.脆性材料损伤应变率效应及其本构模拟[J].固体力学学报,1998,19(4).
[26].杨廷毅.混凝土损伤断裂过程研究[M],浙江大学出版社,1993.
[27].杨光松.损伤力学与复合材料损伤[M],国防工业出版社,1995.
[28].陈建兵,李杰.混凝土结构非线性随机损伤状态的演化分析[PI,东南大学学报,2002.32(5):756-759.
[29].任青文,徐道远等.复杂条件下的高拱坝(300米级)建设中的应用基础研究,国家自然科学基金委员会,1997:1-53.
[30].刘华,蔡正敏,杨菊生等.混凝土结构三维损伤开裂破坏全过程非线性有限元分析[J].工程力学,1999,4:45-51.
[31].Noh.Sam-Yong,Kratzig,wilfried.B,Meskourid.konsiantin.Numerical simulation of service ability,damage evolution and failure of reinforced concrete shells. Computer and structures,v81,n8-1l, p843-857 May,2003.
[32].张强勇,朱维申,金亚兵.弹塑性损伤模型在某地下厂房中的应用[J],岩石力学与工程学报,1999.18( 6):654-657.
[33].Carmeliet J and Hens H.Probabilistic nonlocal damage model for continua with random field properties.J Engineering Mechanics, 1994,120(10):2013~2027.
[34].赵吉东,周维垣,刘元高,杨若琼.岩石类材料应变梯度损伤模型研究及应用[J].水利学报,2002,(7):70—74.
[35].赵吉东,周维垣,刘元高.基于应变梯度的损伤局部化研究及应用[J].力学学报,2002,34(3):445—452.
[36].赵扬锋,潘一山.岩石单轴压缩应变梯度损伤模型局部化带宽[J].辽宁工程技术大学学报,2006,25(3):355—357.
[37].王学滨.考虑应变梯度及刚度劣化的剪切带局部变形分析[J].工程力学,2006,23(10):101—106.
[38].燕秀发,戴耀.功能梯度材料裂纹高阶渐近场研究[J].机械强度,2006,28(4):593—597.
[39].Aifantis EC. On the role of gradient in the localization of deformation and fracture.Int J Engng Sci,1992,30:1279—1299.
[40].Roy A.B.Engelen,Marc G.D.Geers,Frank P.T.Baaijens.International Journal of Plasticity 19(2003)403-433.
[41].沈新普,沈国晓,陈立新,杨璐.梯度增强的弹塑性损伤非局部本构模型研究[J].应用数学和力学,2005,26(2):201—214.
[42].杨璐,朱浮声,沈新普.弹塑性损伤模型在应变局部化分析中的应用[J].力学与实践,2006,28(2):9—15.
[43].杨璐,朱浮声,沈新普.梯度依赖的混凝土弹塑性非局部损伤的本构模型[J].材料研究学报,2007,21(5):477—481.
[44].杨璐,沈新普,孙光.混凝土弹塑性损伤本构理论的研究[J].沈阳工业大学学报,2005,27(3):321—324.
[45].邓守春,刘金兴,张晶,梁乃刚.两种弹性损伤模型的基本方程与色散关系讨论[J].固体力学学报,2008,29(1):29—34.
[46].Borst R de,Pamin Jerzy,Marc G,et al.On coupled gradient-dependent plasticity and damage theories with a view to localization analysis.Eur J Mech A/Solids,1999,18:939—962.
[47].de Borst,R.,Pamin,J.,1996.Some novel developments in finite ele- ment procedures for gradient-dependent plasticity. Int.J.Numer. Meth.Eng.39(14),2477-2505.
[48].李运军.基于非局部损伤模型的混凝土损伤的数值模拟[D].安徽:合肥工业大学,2006.
[49].张伟.基于隐式梯度模型混凝土损伤的数值模拟[D].安徽:合肥工业大学,2008.
[50].徐俊祥.混凝土结构的损伤力学分析[D].江苏:河海大学,2004.
[51].周苏波.混凝土损伤的定量分析[D].江苏:河海大学,1999.
[52].刘军.混凝土损伤分析及其工程应用[D].辽宁:大连理工大学,2004.
[53].王向东.混凝土损伤理论在水工结构仿真分析中的应用[D].江苏:河海大学,2004.
[54].邓爱民.混凝土损伤行为特性研究[D].江苏:河海大学,2001.
[55].Eringen E,Mcdowell DL,Talreja R.Gradient concepts for evolution of damage. Mech Mater,1987,31:831—860.
[56].Lena S,Ristinmaa M.FE-formulation of a nonlocal plasticity theory. Comp Meth Appl Mech Engng,1996,136:127—144.
[57].Aifantis E C. On the microstructural origin of certain inelastic models[J]. Journal of Engineering Material Technology, 1984, 106:326-330.
[58].Nedjar B.Elastoplastic-damage modeling including the gradient of damage: Formulation and computational aspects.International Journal of Solids and Structures,2001,38:5421—5451.
[59].Ramaswamy S,Aravas N.Finite element implementation of gradient plasticity models.Part I:Gradient-dependent yield functions.Comp- uter Methods in Applied Mechanics and Engineering,1998,163:11~32.
[60].李锡夔,S.Cescotto.梯度塑性的有限元分析及应变局部化模拟[J].力学学报,1996,28(5):575-584.
[61].潘一山、徐秉业、王明洋.岩石塑性应变梯度与Ⅱ类岩石变形行为的研究[J].岩土工程学报,1998,21(4):471-474.
[62].徐秉业,刘信声.应用弹塑性力学[M].北京:清华大学出版社,1995.
[63].杨伯源,汪忠明,巫绪涛,李和平.权函数在非局部损伤模型中的应用研究[J].合肥工业大学学报,2005,28(9):1178—1181.
[64].李和平.高精度求解应力强度因子的数值外插法[J].机械强度, 2001,23(2):209—212.
[65].李和平,巫绪涛,杨伯源.数值外插法求解K因子的插值基础及误差估计[J].合肥工业大学学报,2000,23(6):962—965.
[66].巫绪涛,杨伯源.数值外插法求解空间裂纹应力强度因子的研究[J].合肥工业大学学报,1999,22(4):26—31.
[67].巫绪涛,杨伯源,李和平.一种新的求解K因子的位移插值法[J].河海大学学报,2003,31(6):662—665.
[68].I.M.Smith,D.V.Griffiths.有限元方法编程[M].北京:电子工业出版社,2003.
[69].王焕定,吴德伦主编,林家骥主审.有限单元法及计算程序[M].北京:中国建筑工业出版社,2004.
[70].丁皓江,何福保,谢贻权,徐兴.弹性和塑性力学中的有限单元法[M].北京:机械工业出版社,1992.
[71].彭国伦.Fortran 95程序设计[M].北京:中国电力出版社,2004.
[72].虞吉林.材料强化、损伤、破坏中的尺度效应[J].安徽:中国科学技术大学.
[73].Bazant ZP,Cedolin L.Stability of Structures: Elastic, Inelastic, Fracture and Damage Theories.New York:Oxford University Press,1991.
[74].Bazant ZP,Jirasek Milan. Nonlocal integral formulations of plasticity and damage:survey of progress.Journal of Engineering Mechanics,ASCE,2002,128:1119—1149.
[75].Dugdale D S. Yielding of steel sheets containing slits[J]. Mech Phys Solids, 1960,8:100-108.
[76].Barenblatt G I. The mathematical theory of equilibrium cracks in brittle fracture[J]. Advances Appl Mech, 1962, 7:55-129.
[77].Hillerborg A, Modéer M and Petersson P–E. Analysis of crack formation and crack growth in concrete by means of fracture mechanics and finite elements. Cement Concr Res, 1976, 6:773-781.
[78].A C Eringen. Nonlocal polar elastic continua[J]. Int J Solids Struct. 1972, 10:233~248.
[79].Peerlings R H J.,Borst R D, Brekelmans W A M etal. Localization issues in local and nonlocal continuum approaches to fracture[J]. European Journal of Mechanics A/Solids, 2002,21:175-189.
[80].Pijaudier-Cabot G, Ba?ant Z P,Tabbara M.Comparison of various models for strain-softening[J].Engineering Computation,1988,5:141-150.
[81].Peerlings R H J, Geers M G D, Borst R D, et al. A critical comparison of nonlocal and gradient-enhanced softening continua[J].Internat- ional Journal of Solids and Structures, 2001,38:7723-7746.
[82].Pijaudier-Cabot, G and Ba?ant Z P. Nonlocal damage theory[J]. Eng Mech., 1987,113:1512~1533
[83].D A. Fulk and D. W. Quinn. An analysis of 1-D smoothed particle hydrodyna- mics kernels[J]. Comput. Phs. 1996,126:165~180
[84].Lasry D,Belytschko T. Localization limiters in transient problems[J]. International Journal Solids Structures,1988, 24: 581-597.
[85].Huerta A, Pijaudier-Cabot G. Discretization influence on Regulari- zation by two localization limiters[J]. ASCE Journal of Engineering Mechanics,1994, 120:1198-1218.
[86].虞吉林.考虑微结构的固体力学的进展和若干应用.力学进展, 1985,15(1):82~89.
[87].Bazant ZP,Pijaudier-Cabot.Nonlocal continuum damage,localization instability and convergence.ASME J Appl Mech,1988,55:287~293.