基于空心球代表性体积无构造孔隙介质非线性宏观本构势的合理性研讨.pdfVIP

甚于空心球代衰性体积元构造孔隙 介质非线性宏观本构势的合理性研究， 张 震 黄筑平 (北京大学力学与工程科学系，北京，100871) . 要 代表性体积元是细观力学中预侧材料宏观本构关系的基本分析单元.如果对该单元 施加不同的边界条件，则可能会得到不同的宏观本构关系.比如施加均匀应力边界条件，可 以得到宏观应变势的上界，而施加线性位移边界条件，则可得到宏观应变势的下界.但是， 在一定的徽结构下，一个合理的代表性体积元应该使在不同边界条件下得到的宏观本构关系 基本上相同.因此，研究代表性体积元的合理选取问鹰，对建立正确的宏观本构关系将是十 分IWK诊文以基体为非线性幕次材料的孔隙介质为研究对象，通过构造相应的宏观本构 势，攀点讨论了目前被国际上广泛采用的空心球壳代表性体积单元的合理q!(夙本·文通过 构造在线性位移边界条件下的空心球壳胞体内的细观位移许可场，得到了宏观应力势 W(E)，并由此求出了相应的宏观应变势TE(蜀.进而可计算文献[I〕中通过应力边界条 件得到的宏观应变势TI(x与本文得到的宏观应变势TE(助之间的差，并以此作为对代 表体积单元合理性的相应判据.计算结果表明:当荃体为线弹性材料时，以上两种边界条件 下的宏观应变势非常接近。其相对误差不超过15%;当基体的非线性不断增强时，其偏差也 随之不断增大卜但总的来看，二者还是比较接近的，由此可见，球壳模型具有较好的合理性， 特别是当基体材料的非线性不太强的时喊 关.一 代表性体积元，应力势，应变势，孔晾材料，非线性，幕次 0 引 言 采用代表体积单元 R〔VE,RepresentativeVolumeElement)研究复合材料的宏观性 能是细观力学常用的分析方法.作为理论分析模型，最为简单也最为常用的RVE是一个 球形胞体.较早运用球形RVE研究复合材料的工作可以追溯到Hashin(1962).他根据球 体问题的线性弹性精确解，得到了含球形夹杂复合材料的宏观线性本构关系的上下界 在有关多孔材料的研究中，空心球壳RVE模型一直被广泛地应用，如Gurson模型等. 如果对代表性体积元施加不同的边界条件，比如应力边界条件或位移边界条件，所得到 宏观本构关系一般来说是不同的.而一个合理的代表性体积元应该使不同边界条件下得 到的宏观本构关系基本上相同，所以代表性体积元选取的合理性问题值得讨论 *国家自然科学基金资助课题 256 对于孔隙材料，考虑球壳RVE胞体fZ(如图1)，其内半径为a，外半径为b，凡和 Vv分别是基体体积和孔洞体积，总体积为V二玖十V,.胞体上的平均应力和平均应变分 别为E和E，而基体中的微观应力场和微观应变场分别为U和￡ 田1珠充RVE 如果在RVE边界上给定位移边界条件 “=x.E (1) 式中:E是二阶对称张量，则材料的平均应变￡就是宏观应变E，并且可以通过构造 胞体兑 内的位移许可场，由最小势能原理获得相应的宏观应力势: OEE()--}-0E()=icyfglx(,e)dV· (2) 式中:47(x,E)为基体材料的徽观应力势:U是胞体内的位移许可场.当求出(2)式的 p4E(E)，我们便可利用如下的恒等式求得相应的宏观应变势: TE=va.E一40 (3) dE 如果在RVE的边界上给定面力边界条件t=n-E，则平均应力0就是E.通过构 造球壳胞体Q内的应力许可场，由最小余能原理可获得相应的宏观应变势: `FE(E)+y+(E)一infW(x,6(x)=inf典fV(x,v)dV. (4) 一 叮哟 --一 otfV 式中:w(x,Q)是基体材料的应变势:.S是胞体内的应力许可场. 一