颗粒增强橡胶细观力学性能二维数值模拟

在细观层次上建立了具有随机分布形态的代表性体积单元,通过细观力学有限元方法对炭黑颗粒填充橡胶复合材料的宏观力学性能进行了研究。采用二维平面应变模型进行单轴压缩模拟仿真,通过施加周期边界条件保证了代表性体积单元变形场的协调性。重点研究讨论了颗粒的随机分布形态和粒径大小、刚度、体积分数对复合材料宏观应力-应变关系曲线和有效弹性模量的影响。结果表明:炭黑颗粒的填充显著提升了橡胶材料的刚度,在炭黑含量Vf=0.2513时,复合材料的有效弹性模量值高于橡胶初始模量值的2倍;复合材料的有效弹性模量随颗粒所占体积分数的增加而增大。     

炭黑颗粒增强橡胶界面脱粘单向拉伸力学行为研究

利用岛津万能试验机对天然橡胶与丁苯橡胶共混的胶料(NSBR)进行常温下多步松弛循环加卸载试验,拟合材料参数。通过电镜实验观察复合材料截面,发现炭黑体积分数较大时,炭黑颗粒存在聚集现象。建立了圆形炭黑颗粒增强橡胶复合材料二维多颗粒随机分布RVE模型,并运用ABAQUS软件对炭黑颗粒增强橡胶界面脱粘和颗粒聚集行为进行仿真。对比不同体积分数和不同粒径分布形态等情况下单轴拉伸的宏观力学响应,结果表明:炭黑体积分数越大,炭黑增强橡胶界面开始发生脱粘的应变越小,越容易发生脱粘。针对炭黑的不同聚集程度,并考虑颗粒的界面脱粘,发现考虑颗粒聚集的模型应力集中程度更高,更容易发生脱粘。因此,炭黑体积分数越高,越不能忽略颗粒与基体的界面脱粘以及炭黑颗粒团聚。     

炭黑填充橡胶材料的修正三链模型和Mullins模型

本文探讨了炭黑颗粒填充橡胶材料的本构模型。考虑到橡胶单个分子链与周围分子网络的约束作用和炭黑颗粒对橡胶的补强作用,提出了一种修正三链模型,用Edwards管模型描述分子链之间的相互作用和约束,采用应变放大因子来考虑炭黑含量的影响。并在修正三链模型的基础上,利用橡胶分子网络重构理论,提出了一种适合表征橡胶Mullins现象的本构模型。通过与实验数据比较分析,修正三链模型可较准确地表征未填充橡胶材料不同变形模式的力学性能和炭黑颗粒填充橡胶材料的单向拉伸力学行为,Mullins模型也可较好地描述橡胶材料的Mullins现象。     

基于扩展有限单元法的短纤维/橡胶复合材料裂纹扩展分析

为了研究短纤维/橡胶复合材料裂纹在拉伸载荷作用下的扩展演化规律,应用扩展有限元法对预制裂纹的短纤维/橡胶复合材料的裂纹扩展进行了数值模拟.采用随机顺序吸附算法在ABAQUS软件中生成短纤维分布模型,分析了失效准则参数(最大许用主应力和裂纹表面能)、短纤维细观参数(体积含量、长度和取向角)对裂纹扩展行为的影响规律,探究了...     

基于有限元压缩方法的复合材料RVE创建

文章提出一种有限元压缩方法,可以简单、高效地创建具有较高增强体体积分数的复合材料代表性体胞单元(RVE),其具体步骤如下:(1)基于随机顺序吸附(RSA)法生成具有较低增强体体积分数的复合材料周期性RVE;(2)在周期性边界条件约束下,采用有限元方法压缩前述创建的低增强体体积分数复合材料周期性RVE,得到有限元网格格式、具有较高增强体体积分数的复合材料周期性RVE;(3)通过后处理提取得到的高增强体体积分数复合材料周期性RVE中所有增强体的位置(和取向),进而创建CAD格式的复合材料周期性RVE.采用提出的有限元压缩方法,成功创建了体积分数达50.0%的球形增强体复合材料的周期性RVE.采用最近邻距离的概率分布函数、最近邻取向角的累积概率分布函数、Ripleys-K函数和对关联函数对创建的复合材料周期性RVE中球形增强体分布规律进行了统计分析,发现创建的复合材料周期性RVE中球形增强体空间随机分布.基于创建的复合材料周期性RVE和有限元均质法预测了不同类型的球形增强体复合材料的弹性性能,并与实验测试和双夹杂模型预测的结果进行了对比,验证了创建的复合材料周期性RVE及提出的有限元压缩方法...     

炭黑颗粒对丁苯橡胶材料力学性能的影响

利用电子万能材料试验机,对不同体积分数的炭黑填充丁苯橡胶复合材料进行了单向拉伸和循环拉伸加卸载等准静态力学试验,研究了炭黑填充橡胶材料单向拉伸应力应变关系、Mullins效应、调制应变相关性和拉伸断裂力学行为等.试验结果表明:调制应变越大,橡胶材料的刚度越小,橡胶的非线性应力应变关系越明显;炭黑含量越高,橡胶材料的初始模量和刚度越大,应力反翘现象和Mullins效应越显著;同时,随着炭黑体积分数的增加,炭黑填充丁苯橡胶材料的拉伸强度和断裂伸长率将呈现先增大后减小的趋势.     

改进的二维高精度通用单胞模型

基于复合材料细观力学周期性假设,利用高阶理论的改进算法,对高精度通用单胞模型的计算方法进行了改进．模型中用界面的平均量代替细观位移函数中解系数,并利用细观单元力学方程的分析与求解,建立了细观平均量与复合材料宏观平均量间的联系．改进高精度通用单胞模型的求解方程数目减少了大约60％,求解时间大大缩短;并且消除了亚子胞的概念,同时解耦了横向与纵向的方程．该模型的计算结果与试验结果及理论计算结果具有较好的一致性．     

炭黑增强对橡胶复合材料温度相关力学行为的影响

利用含温度箱的自动网格法测试系统,对不同炭黑含量的橡胶复合材料开展不同温度下的超弹性性能试验,考察了炭黑增强对其温度相关力学行为的影响。结果表明,温升过程中炭黑颗粒对橡胶基体的增强作用在逐渐减弱,与之相关的内能弹性效应也逐渐减小,加上橡胶基体本身的熵弹性效应,这两者共同决定了炭黑增强胶料会随着温度的升高"先变软后变硬",而非如纯胶料般随着温度的升高直接变"硬"。换言之,温度变化对颗粒与基体间相互作用造成的影响,使得炭黑增强橡胶的温度相关力学行为比不含炭黑的纯橡胶更为复杂。另外,依据试验研究的结果,本文还推导了修正八链模型的温度相关显式表达式,它能较好地表征胶料超弹性力学行为与温度的相关性。     

基于多尺度方法的三维编织C/C复合材料弯曲失效特性分析

为了预测三维编织C/C复合材料的弯曲失效行为,基于多尺度渐进展开理论,结合细观渐进损伤模型,建立了三维编织C/C复合材料宏细观多尺度分析模型。通过商业有限元软件ABAQUS用户子程序UMAT的二次开发,在宏观结构有限元分析中实时调用细观单胞模型进行细观渐进损伤分析,实现了宏细观尺度之间交互式信息传递和多尺度损伤模拟。利用上述模型对三点弯曲载荷下三维编织C/C复合材料梁的渐进损伤和失效过程进行了模拟,预测了梁的载荷-挠度曲线和弯曲强度,并与实验结果进行了对比分析,验证了基于多尺度方法的三维编织C/C复合材料弯曲强度预测模型的有效性,为此类材料及结构失效分析提供了一种手段。     

混凝土材料宏观力学特性分析的细观单元等效化模型

提出了一种混凝土材料宏观力学特性分析的新方法—细观单元等效化模型。该方法从描述混凝土材料的细观尺度入手,采用Monte Carlo法生成由骨料颗粒及砂浆基质组成的混凝土试件的随机骨料模型;然后,依据混凝土材料特征单元尺度来剖分有限元网格并投影到建立的随机骨料模型上,各细观单元的有效力学特性则采用复合材料等效化方法来确定。本文方法体现了材料非线性宏观力学特性源于其内在的不均匀性这一认识,而对不均匀性的描述则是以网格剖分是否影响其宏观力学特性为准则。因此,本文方法较其他细观力学方法最大的优点在于极大地减小了体系自由度数目（特别是对于三维问题）,从而提高了计算效率。算例分析初步验证了本文方法的高效性。     

空心颗粒填充复合材料的建模与弹塑性仿真计算

空心玻璃微珠（hollow glass microballon,简称HGM）填充复合材料称为复合泡沫材料,近来已被广泛应用在工业中。利用RSA(Random Sequential Adsorption)方法生成了含不同体积比的HGMs填充代表体元模型,然后用有限元方法计算得出了材料的应力-应变关系,将材料属性简化为双线性随动强化模型,分析了HGM填充比、壁厚对材料的有效弹性常数、屈服极限的影响,并分析了材料内部细观应力场和塑性应变分布情况。结果发现,HGM壁厚比对材料有效弹性模量和屈服极限的增减起着决定性作用,而对于任意填充模式来说,其比模量和比强度总是大于纯树脂,这一点体现了该材料轻质的优良特性。材料基体的应力集中部位分布以及塑性应变区域的分布也取决于HGM的壁厚比。     

A numerical approximation to the elastic properties of sphere-reinforced composites

Three-dimensional cubic unit cells containing 30 non-overlapping identical spheres randomly distributed were generated using a new, modified random sequential adsortion algorithm suitable for particle volume fractions of up to 50%. The elastic constants of the ensemble of spheres embedded in a continuous and isotropic elastic matrix were computed through the finite element analysis of the three-dimensional periodic unit cells, whose size was chosen as a compromise between the minimum size required to obtain accurate results in the statistical sense and the maximum one imposed by the computational cost. Three types of materials were studied: rigid spheres and spherical voids in an elastic matrix and a typical composite made up of glass spheres in an epoxy resin. The moduli obtained for different unit cells showed very little scatter, and the average values obtained from the analysis of four unit cells could be considered very close to the “exact” solution to the problem, in agreement with the results of Drugan and Willis (J. Mech. Phys. Solids 44 (1996) 497) referring to the size of the representative volume element for elastic composites. They were used to assess the accuracy of three classical analytical models: the Mori-Tanaka mean-field analysis, the generalized self-consistent method, and Torquato's third-order approximation.     

Dependence of nonlinear elasticity on filler size in composite polymer systems

Nanosized filler particles enhance the mechanical properties of polymer composites in a size-dependent fashion. This is puzzling, because classical elasticity is inherently scale-free, and models for the elasticity of composite systems never predict a filler-size dependence. Here, we study the industrially important system of silica-filled rubbers, together with a well-characterized model-filled crosslinked gel and show that at high filler content both the linear and nonlinear elastic properties of these systems exhibit a unique scaling proportional to the cube of the volume fraction divided by the particle size. This remarkable behavior makes it possible to predict the full mechanical response of particle-filled rubbers for small but finite deformations based solely on the rheology of the matrix and the size and modulus of the filler particles.     

N330炭黑增强天然橡胶材料力学性能的实验研究

为了研究炭黑对橡胶材料力学性能的影响,对9种不同体积含量的炭黑填充橡胶材料进行了准静态力学实验研究。利用循环拉伸加卸载实验分析了炭黑对橡胶Mullins效应及能量损耗的影响,通过单轴拉伸实验研究了炭黑对橡胶材料刚度和起始模量的影响,采用多步松弛拉伸加卸载实验研究了炭黑对橡胶材料应力行为的应变历史相关性的影响。实验结果表明:炭黑填充量越高,橡胶材料的刚度越大,初始模量越大,Mullins效应也越明显;随着炭黑填充量的增加,橡胶在加卸载循环中所产生的迟滞损耗、Mullins效应相对能量损耗以及残余应变都呈现出非线性增长趋势;随着炭黑填充量的升高,橡胶在加卸载过程中的应力松弛现象越明显,其平衡态迟滞损耗以及与时间相关部分的迟滞损耗也越大。     

含界面相Voronoi单元的等效弹性常数的计算

采用含界面相Voronoi单元有限元法,根据广义胡克定律,计算了在给定边界条件下,颗粒增强复合材料的等效弹性常数。建立了含多个随机分布的椭圆形夹杂及界面相的VCFEM模型,分析了夹杂体分比,界面相厚度和界面相弹性模量等因素对颗粒增强复合材料等效弹性常数的影响,并利用普通有限元方法对比验证。结果表明,当界面相弹性模量小于基体与夹杂时,材料的等效弹性模量会随着界面相厚度的增大而减小,随着夹杂体分比的增大而减小,并且界面过薄时,材料的等效弹性模量会随着夹杂体分比的增大而增大;当界面相弹性模量大于基体或夹杂时,材料的等效弹性模量会随着夹杂体分比和界面相厚度的增大而增大。而界面相的厚度和弹性模量对材料的等效泊松比的影响较小,材料的等效泊松比主要受夹杂体分比的影响,与其呈反比关系。     

Analysis of elastoplastic deformation in metal-matrix composites with particulate reinforcements

In this paper, elastoplastic stress-strain behavior during tensile deformation of an aluminum alloy matrix composite containing alumina circular and non-circular particles is analyzed. In terms of cell models in conjunction with continuum plasticity theory, various periodic arrays of particles are assumed in a three-dimensional finite element simulation. The geometrical effects of particle volume fraction, shape, aspect ratio, array and distribution, as well as non-circular particle orientation on the overall elastoplastic stress-strain behavior are examined in view to design optimum microstructures of the composites. The project supported by the National Natural Science Foundation of China and the State Education Commission of China     

A new formulation of the effective elastic–plastic response of two-phase particulate composite materials

In this paper the double-inclusion model, originally developed to determine effective linear elastic properties of composite materials, is reformulated and extended to predict the effective nonlinear elastic–plastic response of two-phase particulate composites reinforced with spherical particles. The resulting problem of elastic–plastic deformation of a double-inclusion embedded in an infinite reference medium subjected to an incrementally applied far-field strain is solved by the finite element method. The proposed double-inclusion model is evaluated by comparison of the model predictions to the available exact results obtained by the direct approach using representative volume elements containing many particles. It is found that the double-inclusion formulation is capable of providing accurate prediction of the effective elastic–plastic response of two-phase particulate composites at moderate particle volume fractions.     

Material instability-induced extreme damping in composites: A computational study

We investigate the effective viscoelastic performance of particle-reinforced composite materials whose particulate phase undergoes a material instability resulting in temporarily non-positive-definite elastic moduli. Recent experiments have shown that phase transitions in geometrically-constrained composite phases (such as in particles embedded in a stiff matrix) can lead to stable non-positive-definite elastic moduli, and they hinted at strong damping increases that can be achieved from such metastable composite phases. All previous theoretical efforts to explain such phenomena have used simplistic one-dimensional models or they were based on composite bounds and specific two-phase solids. Here, we study particle–matrix composites with periodic randomized particle dispersion. A finite element discretization is used in combination with a sophisticated nonlinear solver in order to perform the numerous calculations in a feasible amount of computing time. Our computational analysis shows that stable non-positive-definite inclusion moduli can indeed lead to extreme damping increases (i.e. greatly exceeding the intrinsic damping of each composite phase) and that such extreme damping arises from a shift in microstructural mechanisms.     

A thermodynamics based damage mechanics model for particulate composites

A micro-mechanical damage model is proposed to predict the overall viscoplastic behavior and damage evolution in a particle filled polymer matrix composite. Particulate composite consists of polymer matrix, particle fillers, and an interfacial transition interphase around the filler particles. Yet the composite is treated as a two distinct phase material, namely the matrix and the equivalent particle-interface assembly. The CTE mismatch between the matrix and the filler particles is introduced into the model. A damage evolution function based on irreversible thermodynamics is also introduced into the constitutive model to describe the degradation of the composite. The efficient general return-mapping algorithm is exploited to implement the proposed unified damage coupled viscoplastic model into finite element formulation. Furthermore, the model predictions for uniaxial loading conditions are compared with the experimental data.