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

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

The Payne effect in finite viscoelasticity: constitutive modelling based on fractional derivatives and intrinsic time scales

As we know from experimental testing, the stiffness behaviour of carbon black-filled elastomers under dynamic deformations is weakly dependent on the frequency of deformation but strongly dependent on the amplitude. Increasing strain amplitudes lead to a decrease in the dynamic stiffness, which is known as the Payne effect. In this essay, we develop a constitutive approach of finite viscoelasticity to represent the Payne effect in the context of continuum mechanics. The starting point for the constitutive model resulting from this development is the theory of finite linear viscoelasticity for incompressible materials, where the free energy is assumed to be a linear functional of the relative Piola strain tensor. Motivated by the weak frequency dependence of the dynamic stiffness of reinforced rubber, the memory kernel of the free energy functional is of the Mittag Leffler type. We demonstrate that the model is compatible with the Second Law of Thermodynamics and equal to a fractional differential equation between the overstress of the Second Piola Kirchhoff type and the Piola strain tensor. In order to represent the dependence of the dynamic stiffness on the amplitude of strain, we replace the physical time by an intrinsic time variable. The temporal evolution of the intrinsic time is driven by an internal variable, which is a measure for the current state of the material's microstructure. The material constants of the model are estimated using a stochastic identification algorithm of the Monte Carlo type. We demonstrate that the constitutive approach pursued here represents the combined frequency and amplitude dependence of filler-reinforced rubber. In comparison with the micromechanical Kraus model developed for sinusoidal strains, the theory set out in this essay allows the representation of the stress response under arbitrary loading histories.     

短纤维增强三元乙丙橡胶横观各向同性黏-超弹性本构模型

短纤维增强三元乙丙橡胶包覆薄膜,是一种应用于固体火箭发动机缠绕包覆装药的新型复合材料.为了描述其在工作过程中受振动、冲击等载荷作用时的力学行为,基于黏弹性理论和纤维增强连续介质力学理论,提出了一种考虑应变率强化效应的横观各向同性黏-超弹本构模型.模型中应变能函数被分解为超弹性应变能和黏性应变能,其中超弹性应变能包括表征各向同性的橡胶基体应变能和表征各向异性的纤维拉伸应变能,黏性应变能采用表征橡胶和纤维黏性响应的宏观唯象模型.选取表征各应变能的函数形式,经过数学变换、替代、叠加,求解确定最终的应力应变形式,明确模型参数获取的具体步骤,将预测结果与实验结果对比分析,准确性较高.研究表明:该模型能有效预测材料在低应变率下纤维方向为0?～45?的非线性率相关力学特性;模型形式易于实现有限元开发,对固体火箭发动机装药结构完整性分析具有参考价值.     

平面应变不可压缩橡胶圆柱的大变形

针对一类新的橡胶材料应变能函数，推导了受内压作用下不可压缩橡胶圆柱的大变形公式，给出了位移、应力的解析表达式（用积分形式表示）。建立了适用于分析非线性不可压缩橡胶材料的罚有限元列式，算例表明：位移与应力能很好地与理论解吻合，并提出了控制计算稳定的方法，特别是详细地讨论罚因子的选取及对计算结果的影响。     

Large deformation of incompressible rubber cylinder under plane strain

The large deformation of incompressible rubber cylinder under inner pressure is analyzed by a kind of new rubber materials strain energy function. The theory formulation for the displacement and stress is presented. The penalty finite element formulation is established in order to analyze nonlinear rubber materials, and the results of finite element method agree well with theoretical ones. A new method for controlling the calculating stability and convergence rates is put forward. The selection of the appropriate penalty factor and its influence on calculated results are discussed.     

A constitutive model for the Mullins effect with changes in material symmetry

When an unfilled or particle reinforced rubber is subjected to cyclic loading–unloading with a fixed amplitude from its natural reference configuration, the stress required on reloading is less than on the initial loading for a deformation up to the maximum value of the stretches achieved. The stress differences in successive loading cycles are largest during the first and second cycles and become negligible after about 4–6 cycles. This phenomenon is known as the Mullins effect. In this paper new experimental data are reported showing the change in material symmetry for an initially undamaged and isotropic material subjected to uniaxial and biaxial extension tests. The effect of preconditioning in one direction on the mechanical response when loaded in a perpendicular direction is discussed. A simple phenomenological model is derived to account for stress softening and changes in material symmetry. The formulation is based on the theory of pseudo-elasticity, the basis of which is the inclusion of scalar variables in the energy function. When active, these variables modify the form of the energy function during the deformation process and therefore change the material response. The general formulation is specialized to pure homogeneous deformation in order to fit the new data. The numerical results are in very good agreement with the experimental data.     

The influence of dewetting in filled elastomers on the changes of their mechanical properties

Summary The changes of mechanical properties of filled polymers and their dependence on deformation history are the subject of this paper. For most of filled polymers in practical use, the theory of linear viscoelasticity cannot be applied, even at small deformations. In this work, samples of glass bead filled polybutadiene rubber with different filler levels were investigated at small strains (<10%). The evolution of the relaxation modulus and Poisson's ratio was observed in cyclic experiments, which were also applied in inducing a defined deformation history for the succeeding relaxation experiments. In these experiments, the relaxation modulus and Poisson's ratio were measured as functions of time, with strain, strain rate, filler level and the preceeding deformation history as parameters. The results indicate dewetting as the main reason for the changes of the mechanical properties of the filled materials.The financial support of this work by the Deutsche Forschungsgemeinschaft (DFG) is gratefully acknowledged.     

Multi-cycle deformation of silicone elastomer: observations and constitutive modeling with finite strains

Observations are reported on a medical grade of silicone elastomer in uniaxial tensile tests up to breakage of specimens, short-term relaxation tests, and cyclic tests with monotonically increasing maximum elongation ratios. Experimental data in cyclic tests demonstrate the fading memory phenomenon: stress–strain diagrams for two specimens with different deformation histories along the first n?1 cycles and coinciding loading programs for the other cycles become identical starting from the nth cycle. A constitutive model is developed in cyclic viscoplasticity of elastomers with finite strains, and its adjustable parameters are found by fitting the experimental data. Ability of the stress–strain relations to predict the mechanical response in cyclic tests with various deformation programs is confirmed by numerical simulation.     

生长对超弹性球壳变形和稳定性的影响

应用连续介质力学有限变形理论建立受内压作用不可压超弹性球壳大变形问题的力学模型, 且运用基于变形梯度张量极分解的弹性体积生长理论分析生长对不可压超弹性球壳变形和稳定性的影响. 通过对球壳变形与内压关系式的数值计算得到不同生长条件下球壳的变形曲线和应力分布曲线及由生长引起的残余变形和残余应力分布. 计算结果表明生长对球壳变形特性有明显的影响, 生长影响球壳可产生不稳定变形的临界壁厚和临界内压, 且在某些情况下可改变球壳的稳定性. 生长在球壳中可产生一定的残余变形和残余应力, 对球壳中的应力分布有一定的影响. 另外当生长的程度足够大时, 即便没有外力作用,球壳仅在生长引起的残余应力作用下也可产生不稳定变形.      

三轴试验土样剪切破坏过程的判别

本文在土工三轴试验中应用土样全表面变形数字图像测量系统,通过跟踪包裹土样的橡皮膜表面的标记点,得到加载过程中其全表面的变形和应变分布。根据土样表面的应变分布和实测的应力应变关系曲线,可以得到每一时刻土样表面的应力分布。于是,土样表面剪切破坏开始出现的时刻和剪切破坏带贯穿的时刻可分别依据应变和应力判定。依据应变,将橡皮膜上各点的轴向应变的拐点最晚出现的时刻作为剪切带完全贯穿的时刻;依据应力,将全部点中某一个或几个点最早达到峰值应力的时刻作为剪切带开始出现的时刻;土样各点相继达到应力峰值的过程就是剪切破坏的过程。以剪切破坏的起始点和剪切带贯穿的点为分界,土样剪切破坏过程可以划分成三个阶段:破坏前阶段、破坏阶段和破坏后阶段,在不同阶段土样具有不同的变形特点和变形机理。     

多孔介质的流变模型研究

多孔介质在应力作用下具有弹性变形和黏性变形两种完全不同的变形机制，多孔介质的弹性变形是由介质的本体有效力所致，而黏性变表则是由介质的结构有效应力所致。多孔介质的总变形为弹性变形和黏性变形的叠加，计算多孔介质总应变量的流变模型必须同时采用本体有效应力和结构有效应力（双重有效应力），而传统的流变模型仅采用Terzaghi有效应力是不妥当的，它无法正确描述多孔的应变行为，采用了双重有效应力之后的流变模型，通过调节介质特性参数，可以拟合介质的实际应变行为，并且把多孔介质与普通固体联系了起来。     

Torsional Swelling of a Hyperelastic Tube with Helically Wound Reinforcement

This paper examines the combination of radial deformation with torsion for a circular cylindrical tube composed of a transversely isotropic hyperelastic material subject to finite deformation swelling. The stored energy function involves separate matrix and fiber contributions such that the fiber contribution is minimized when the fiber direction is at a natural length. This natural length is not affected by the swelling. Hence swelling preferentially expands directions that are orthogonal to the fiber. The swelling itself is described via a swelling field that prescribes the local free volume at each location in the body. Such a treatment is a relatively simple generalization of the conventional incompressible theory. The direction of transverse isotropy associated with the fiber reinforcement is described by a helical winding about the tube axis. The swelling induced preferential expansion orthogonal to this direction induces the torsional aspect of the deformation. For a specific class of strain energy functions we find that the twist increases with swelling and approaches a limiting asymptotic value as the swelling becomes large. The fibers reorient such that fibers at the inner portion of the tube assume a more circumferential orientation whereas, at least for small and moderate swelling, the fibers in the outer portion of the tube assume a more axial orientation. For large swelling the fibers in the outer portion of the tube reorient beyond the axial orientation, and so are described by helices with orientation in the opposite sense to that in the reference configuration.      

Onset of failure in finitely strained layered composites subjected to combined normal and shear loading

A limiting factor in the design of fiber-reinforced composites is their failure under axial compression along the fiber direction. These critical axial stresses are significantly reduced in the presence of shear stresses. This investigation is motivated by the desire to study the onset of failure in fiber-reinforced composites under arbitrary multi-axial loading and in the absence of the experimentally inevitable imperfections and finite boundaries.By using a finite strain continuum mechanics formulation for the bifurcation (buckling) problem of a rate-independent, perfectly periodic (layered) solid of infinite extent, we are able to study the influence of load orientation, material properties and fiber volume fraction on the onset of instability in fiber-reinforced composites. Two applications of the general theory are presented in detail, one for a finitely strained elastic rubber composite and another for a graphite-epoxy composite, whose constitutive properties have been determined experimentally. For the latter case, extensive comparisons are made between the predictions of our general theory and the available experimental results as well as to the existing approximate structural theories. It is found that the load orientation, material properties and fiber volume fraction have substantial effects on the onset of failure stresses as well as on the type of the corresponding mode (local or global).     

Analytic Solution of the Problem of Additive Formation of an Inhomogeneous Elastic Spherical Body in an Arbitrary Nonstationary Central Force Field

We study the processes of additive formation of spherically shaped rigid bodies due to the uniform accretion of additional matter to their surface in an arbitrary centrally symmetric force field. A special case of such a field can be the gravitational or electrostatic force field. We consider the elastic deformation of the formed body. The body is assumed to be isotropic with elasticmoduli arbitrarily varying along the radial coordinate.We assume that arbitrary initial circular stresses can arise in the additional material added to the body in the process of its formation. In the framework of linear mechanics of growing bodies, the mathematical model of the processes under study is constructed in the quasistatic approximation. The boundary value problems describing the development of stress–strain state of the object under study before the beginning of the process and during the entire process of its formation are posed. The closed analytic solutions of the posed problems are constructed by quadratures for some general types of material inhomogeneity. Important typical characteristics of the mechanical behavior of spherical bodies additively formed in the central force field are revealed. These characteristics substantially distinguish such bodies from the already completely composed bodies similar in dimensions and properties which are placed in the force field and are described by problems of mechanics of deformable solids in the classical statement disregarding the mechanical aspects of additive processes.     

Mixed-mode fracture analysis of orthotropic functionally graded materials under mechanical and thermal loads

Mixed-mode fracture problems of orthotropic functionally graded materials (FGMs) are examined under mechanical and thermal loading conditions. In the case of mechanical loading, an embedded crack in an orthotropic FGM layer is considered. The crack is assumed to be loaded by arbitrary normal and shear tractions that are applied to its surfaces. An analytical solution based on the singular integral equations and a numerical approach based on the enriched finite elements are developed to evaluate the mixed-mode stress intensity factors and the energy release rate under the given mechanical loading conditions. The use of this dual approach methodology allowed the verifications of both methods leading to a highly accurate numerical predictive capability to assess the effects of material orthotropy and nonhomogeneity constants on the crack tip parameters. In the case of thermal loading, the response of periodic cracks in an orthotropic FGM layer subjected to transient thermal stresses is examined by means of the developed enriched finite element method. The results presented for the thermally loaded layer illustrate the influences of the material property gradation profiles and crack periodicity on the transient fracture mechanics parameters.     

Time-dependent thermoelastic creep analysis of rotating disk made of Al–SiC composite

Time-dependent creep stress redistribution analysis of rotating disk made of Al–SiC composite is investigated using Mendelson’s method of successive elastic solution. All mechanical and thermal properties except Poisson’s ratio are radial dependent based on volume fraction percent of SiC reinforcement. The material creep behavior is described by Sherby’s constitutive model using Pandey’s experimental results on Al–SiC composite. Loading is an inertia body force due to rotation and a distributed temperature field due to steady-state heat conduction from inner to outer surface of the disk. Using equations of equilibrium, stress strain, and strain displacement, a differential equation, containing creep strains, for displacement is obtained. History of stresses and deformations are calculated using method of successive elastic solution. It is concluded that the uniform distribution of SiC reinforcement does not considerably influence on stresses. However, the minimum and most uniform distribution of circumferential and effective thermoelastic stresses belongs to composite disk of aluminum with 0% SiC at inner surface and 40% SiC at outer surface. It has also been found that the stresses, displacement, and creep strains are changing with time at a decreasing rate so that after almost 50 years the solution approaches the steady-state condition.     

Transversely isotropic material: nonlinear Cosserat versus classical approach

We consider a specific case of unidirectional reinforced material under applied tensile load. The reinforcement of the material is inclined with 45° to the direction of the tensile resultant. Different approaches are discussed: one experiment and three computational models. Two models use the classical Cauchy continuum theory whereas the third computational model is based on a Cosserat continuum. It is well known that test specimen with inclination between unidirectional reinforcement and tensile direction show, besides Poissons effect, additional deformation perpendicular to the load direction. The classical transversely isotropic continuum theory predicts this deformation as typical S-shape. In the Cosserat continuum the orientation of the inner structure is incorporated. Thus, structural parameters influence the deformation. With the proposed geometrically non-linear Cosserat model classical and non-classical behaviour can be modelled. In the non-classical case, the transverse deformation is not described by one S-shape but by multiple S-shaped modes. The additional rotational parameters in the Cosserat continuum are responsible for the non-classical behaviour which is due to non-symmetric strain.     

短纤维增强EPDM包覆薄膜超弹性本构模型

短纤维增强三元乙丙橡胶(EPDM)包覆薄膜用于一种新型缠绕包覆工艺,主要解决复杂构型自由装填药柱外表面可靠性包覆问题.为了描述其在固体火箭发动机工作过程中产生的大变形、非线性和各向异性等力学行为,根据纤维增强复合材料连续介质力学理论,提出了各向异性超弹性本构模型.该模型中单位体积的应变能函数被解耦成两部分:表征各向同性的橡胶基体应变能和表征各向异性的纤维拉伸应变能,通过引入纤维方向对纤维应变能进行修正,给出了通过单轴拉伸、偏轴拉伸实验数据获取模型参数的具体方法.研究结果表明,该模型能够很好地预测材料在纤维方向0°~45°时的各向异性力学特性,并将预测结果与实验数据对比,误差在5%以下.所建立的各向异性超弹性本构模型准确性高、易于实现数值开发,在一定程度上能够为固体火箭发动机的装药结构完整性分析提供理论依据.