橡胶楔体与刚性缺口接触大变形分析

利用Ｋｕｏｗｌｅｓ与Ｓｔｅｒｎｂｅｒｇ提出的非线性弹性大变形应变能函数,对橡胶楔体与刚性缺口接触问题进行大变形渐近分析,推导了楔体尖端场的渐近方程,得到楔体尖端附近的应力应变场及应力的奇异性指数与橡胶楔体角度、刚性缺口角度及材料常数有关的表达式;楔尖附近同一半径上应力分量为常数,同时,利用非线性有限元理论编制了大变形有限元程序,考虑楔体尖端与缺口接触边界条件,计算得到了与分析解一致的结论,当缺口角     

Kinetic relations and the propagation of phase boundaries in solids

This paper treats the dynamics of phase transformations in elastic bars. The specific issue studied is the compatibility of the field equations and jump conditions of the one-dimensional theory of such bars with two additional constitutive requirements: a kinetic relation controlling the rate at which the phase transition takes place and a nucleation criterion for the initiation of the phase transition. A special elastic material with a piecewise-linear, non-monotonic stress-strain relation is considered, and the Riemann problem for this material is analyzed. For a large class of initial data, it is found that the kinetic relation and the nucleation criterion together single out a unique solution to this problem from among the infinitely many solutions that satisfy the entropy jump condition at all strain discontinuities.     

Radial deflections of thin precompressed cylindrical rubber bush mountings

For an isotropic incompressible hyperelastic Varga material the plane stress (membrane) theory of thin sheets is employed to formulate the load-deflection relation for small superimposed radial deflections of a cylindrical rubber bush which is precompressed by a large uniform radial inflation. The Varga material is a prototype for rubber over a limited range of deformation and the load-deflection relation obtained provides an extreme lower bound to the practical situation. Moreover this relation complements existing results for cylindrical rubber bushes so that now at least some assessment can be made of the effect of precompression on the radial mode of deflection for bushes of finite length. Typical numerical values are given and are contrasted with corresponding values obtained from existing plane strain radial load-deflection relations for long precompressed cylindrical rubber bush mountings.     

Phase mixtures in dynamics of pseudoelasticity

We give a numerical treatment of phase mixtures in pseudoelasticity from a purely mathematical point of view. It is based on a surprising result that the approximate solution may consist of persistent oscillations in strain which resemble the experimentally observed interface patterns. Such a solution is obtained from a sequence of solutions for a rate-type viscoelastic problem with a non-monotone equilibrium stress-strain relation, for which in the limit as the viscosity tends to infinity the viscoelastic problem reduces to the rate-independent elastic problem describing phase transitions. In this manner, it seems to give yet another perspective for the phase mixture from dynamic point of view as the evolution of an unstable state, in contrast to the traditional treatment from stability analysis for phase equilibrium.     

Stress response of a rubber block with frictional contact under axial loading

Rubber blocks (or springs) are structural components widely used in many applications. Design characteristics of a rubber block under axial loading are an apparent modulus (or compression modulus) and normal and shear stresses on contact surfaces. These are affected by a contact condition of the rubber block in contact with two rigid plates and the shape of the block. The problem of a rubber block bonded to two rigid plates has been solved using various approaches. In contrast, for a rubber block whose one surface is bonded to a rigid plate and the other surface in contact with a frictional surface, there is little work in spite of practical applications. For this contact condition, approximate solutions for rectangular blocks in plane strain and for axi-symmetric discs are derived under the assumption of Coulomb frictional contact. The problem is treated as an extension of the problem of an incompressible rubber block bonded to two rigid plates with one of the plates having a frictional interaction with the rubber block. In the linear range of deformation, finite element analysis and experimental results for rubber blocks with shape factors ranging from 1 to 6 are compared for the validation of analytic results. It is found that friction coefficients play important roles in the design characteristics of the rubber block.     

壁面定常波纹状吹吸槽道流中湍流特性的研究

在非平衡湍流中,如具有周期性边界条件的流动,由于雷诺应力与平均流速的变形率有着 不同的性质,当周期性边界条件发生变化时,雷诺应力和平均流速变形率的相位对边界条件 的响应也不同,但是二者的相位差在相当大的范围内是稳定的. 这一特性加深了对雷诺应力 的认识,并对非平衡湍流中的模式理论及大涡模拟中亚格子雷诺应力模式的建立提出了许多 需要注意的问题. 利用层流模型,把空间周期性边界条件作为某种扰动,研究了扰动 及其非线性项的分布以及相位间的关系,得到了一些有益的结果.     

Fatigue characterization of tire rubber

An investigation of the fatigue behavior of tire rubber was undertaken. The tire rubber used in this work is a blend of natural rubber and polybutadiene. It is vulcanized and filled with carbon black. Fatigue tests were conducted on edge-notched prismatic specimens subjected to tension–tension loading in a displacement-controlled machine. The maximum strain was 8% and the minimum strain was 0.8%. The waveform of the fatigue cycle was sinusoidal at a frequency of 10 Hz. Initial cracks of various lengths were cut into the edge of the rubber specimens with a sharp razor blade. The tearing energy was calculated during the test. A power law relation between the crack growth rate and tearing energy was obtained. The fatigue life of tire rubber was estimated based on the integration of the power law relations and the tearing energy versus fatigue cycles curves in conjunction with the value of the critical tearing energy at crack growth initiation.     

On hysteretic response and stationary phase fronts in rubber

We consider the dynamic response of natural, latex and synthetic, nitrile rubbers under non-monotonic dynamic loading conditions; in particular, we recreate an experiment first considered by Kolsky (Nature 224:1301, 1969) in which two segments of a long rubber specimen are initially maintained at different strain levels by external force and then allowed to evolve dynamically towards equilibrium. We show that as a result of the hysteretic behavior, a phase boundary that is stationary with respect to the material points can be established in both these materials. We also show that this phase boundary persists indefinitely in strain-crystallizing natural, latex rubber, but disappears quickly in the non-crystallizing nitrile rubber.     

Integral-based constitutive equation for rubber at high strain rates

High-speed experiments were conducted to characterize the deformation and failure of Styrene Butadiene Rubber at impact rates. Dynamic tensile stress–strain curves of uniaxial strip specimens and force–extension curves of thin sheets were obtained from a Charpy tensile impact apparatus. Results from the uniaxial tension tests indicated that although the rubber became stiffer with increasing strain rates, the stress–strain curves remained virtually the same above 280 s⁻¹. Above this critical strain rate, strength, fracture strain and toughness decreased with increasing strain rates. When strain rates were below 180 s⁻¹, the initial modulus, tensile strength and breaking extension increased as the strain rate increased. Between strain rates of 180 and 280 s⁻¹, the initial modulus and tensile strength increased with increasing strain rates but the extension at break decreased with increasing strain rates. A hyper-viscoelastic constitutive relation of integral form was used to describe the rate-dependent material behavior of the rubber. Two characteristic relaxation times, 5 ms and 0.25 ms, were needed to fit the proposed constitutive equation to the data. The proposed constitutive equation was implemented in ABAQUS Explicit via a user-defined subroutine and used to predict the dynamic response of the rubber sheets in the experiments. Numerical predictions for the transient deformation and failure of the rubber sheet were within 10% of experimental results.     

Residual strain evaluation of curved surface by grating-transferring technique and GPA

This paper investigates an advanced grating-transferring technique combined with geometric phase analysis (GPA) for residual strain evaluation of curved surface. A standard holographic grating is first transferred to a pre-produced epoxy resin film and then consolidated to a test region of curved surface. With a rubber mold and silicone rubber the deformed grating is replicated to a sheet metal after hole-drilling for release of residual stress. After that the grating is transferred from the sheet metal to the glass plate, which would be served as an analyzer grating (specimen grating). By GPA the local strain distributions related to the phase difference between the reference grating and analyzer grating for the released stress can be evaluated. A validation test has been conducted on the weld joint of a stainless steel tube and the obtained results demonstrate the ability of the method in measuring the residual strain of curved surface.     

Cyclic stress-softening model for the Mullins effect in compression

This paper models the cyclic stress softening of an elastomer in compression. After the initial compression the material is described as being transversely isotropic. We derive non-linear transversely isotropic constitutive equations for the elastic response, stress relaxation, residual strain, and creep of residual strain in order to model accurately the inelastic features associated with cyclic stress softening. These equations are combined with a transversely isotropic version of the Arruda–Boyce eight-chain model to develop a constitutive relation that is capable of accurately representing the Mullins effect during cyclic stress softening for a transversely isotropic, hyperelastic material, in particular a carbon-filled rubber vulcanizate. To establish the validity of the model we compare it with two test samples, one for filled vulcanized styrene–butadiene rubber and the other for filled vulcanized natural rubber. The model is found to fit this experimental data extremely well.     

多孔硅橡胶有限变形的弹性行为

针对孔隙度较大 (孔隙度大于 5 0 % )的硅橡胶材料在压缩情况下的大变形 ,提出了可描述此类可压橡胶材料力学行为的应变能密度函数 ,推导了硅橡胶材料的本构方程。利用硅橡胶材料的单轴压缩实验进行了材料参数拟合 ,讨论了多孔硅橡胶的孔隙度和体积变形对压缩性能的影响     

Fully implicit formulation of elastoplastic homogenization problem for two-scale analysis

In this study, using a virtual work equation, a micro-/macro-kinematic relation and a linearized constitutive relation, a boundary value problem is fully implicitly formulated to determine perturbed displacement increment fields in elastoplastic unit cells for two-scale analysis. It is shown that this implicit homogenization problem can be iteratively solved with quadratic convergences by successively updating strain increment fields in unit cells, and that the boundary value problem formulated provides a computational algorithm which is versatile for initial setting of strain increment fields. The computational algorithm developed is then examined by performing a two-scale analysis of a holed plate with an elastoplastic micro-structure, subjected to tensile loading. This demonstrates that the convergence in iteratively solving the implicit homogenization problem strongly depends on the initial setting of strain increment fields in unit cells.     

Constitutive modeling of isotropic hyperelastic materials in an exponential framework using a self-contained approach

In this paper, an exponential framework for strain energy density functions of elastomers and soft biological tissues is proposed. Based on this framework and using a self-contained approach that is different from a guesswork or combination viewpoint, a set strain energy density functions in terms of the first and second strain invariants is rebuilt. Among the constructed options for strain energy density, a new exponential and mathematically justified model is examined. This model benefits from the existence of second strain invariant, simplicity, stability of parameters, and the state of being accurate. This model can capture strain softening, strain hardening and is able to differentiate between various deformation-state dependent responses of elastomers and soft tissues undergoing finite deformation. The model has two material parameters and the mathematical formulation is simple to render the possibility of numerical implementations. In order to investigate the appropriateness of the proposed model in comparison to other hyperelastic models, several experimental data for incompressible isotropic materials (elastomers) such as VHB 4905 (polyacrylate rubber), two various silicone rubbers, synthetic rubber neoprene, two different natural rubbers, b186 rubber (a carbon black-filled rubber), Yeoh vulcanizate rubber, and finally porcine liver tissue (a very soft biological tissue) are examined. The results demonstrate that the proposed model provides an acceptable prediction of the behavior of elastomers and soft tissues under large deformation for different applied loading states.     

导弹适配器受压及弹射时摩擦力研究

针对导弹适配器的结构特征,分别建立了可压缩橡胶泡沫和不可压橡胶圆筒轴对称平面应变问题有限变形的平衡方程,基于Blatz-Ko应变能函数和三次缩减多项式应变能函数,得到了相应的位移和应力模式;在此基础上求解了适配器受压问题的非线性方程组和导弹发射时适配器所受到的摩擦力.算例分析与有限元数值模拟比较表明:解析解与数值解非常吻合,径向应力在发射筒内外表面误差最大为0.558%,周向应力在粘合面误差最大为0.246%,导弹发射时的最大量纲为一的摩擦力为1.0228.适配器径向应力在材料粘合交界面上最小,在适配器外表面最大,均为压应力;橡胶泡沫和不可压橡胶的周向应力均为压应力,橡胶泡沫的周向应力由内向外变大,不可压橡胶的周向应力由内向外变小.橡胶泡沫的径向受压大于周向受压,不可压橡胶的周向受压大于径向受压.研究不同过盈量对应力和摩擦力的影响表明:过盈量每增加0.0013,橡胶泡沫 层和不可压橡胶层的径向应力约增加0.13,不可压橡胶层的周向应力约增加2.14,而摩擦力约增加0.22.过盈量对不可压橡胶层的周向应力和导弹所受到的摩擦力影响非常大,对橡胶泡沫的径向应力有一定的影响,周向应力变化很小.     

Problems in determining the elastic strain energy function for rubber

Lightly crosslinked natural rubber can be stretched by 600% or more, and recovers almost completely. It is often regarded as a model highly elastic material and characterized by a strain energy function to describe its stress-strain behavior under various types of deformation. A number of such functions have been proposed; some of them appear in current finite element programs. They are usually validated by comparison with measured stress-strain relations by Treloar [7] [L.R.G. Treloar, Stress-strain data for vulcanized rubber under various types of deformation, Trans. Faraday Soc. 40 (1944) 59-70] and Jones and Treloar [15] [D.F. Jones, L.R.G. Treloar, The properties of rubber in pure homogeneous strain, J. Phys. D Appl. Phys. 8 (1975) 1285-1304]. But Treloar pointed out that the relations at high strains became markedly irreversible, and he did not assign a strain energy function for strains greater than about 300%. Rivlin's universal relation between torsional stiffness and tensile stress [14] [R.S. Rivlin, Large elastic deformations of isotropic materials. Part V1: further results in the theory of torsion, shear and flexure, Philos. Trans. R. Soc. A 243 (1949) 251-288] is applied here to show that a typical elastic solid cannot be described by any strain energy function at strains greater than about 300%. Elastic strain energy functions for higher strains, or for other rubbery materials, are thus of doubtful value unless evidence for reversibility of stress-strain relations is adduced or the applicability of a strain energy function is demonstrated.     

Semi-analytic inverse method for rubber testing at high strain rates

Because of their low mechanical wave speed, high strain rate testing of rubber is highly difficult. Indeed, stress and strain homogeneity is hard to achieve. In this paper, a semi-analytic inverse solution is proposed. This solution is based on a uni-axial stress state assumption in the specimen. Moreover, a new-Hookean law is assumed for rubber. The new method is successfully applied to a high strain rate test on a synthetic rubber.     

ANALYSIS OF AN INCOMPRESSIBLE RUBBER WEDGE CONTACTING WITH A RIGID NOTCH

I. INTRODUCTION It is well known that rubber materials are much more important in modern society. It is, therefore,not surprising that a sizable number of investigations have been concerned with the rubber materials.Considerable attention has been paid …     

On the response of rubbers at high strain rates—I. Simple waves

In this series of papers, we examine the propagation of waves of finite deformation in rubbers through experiments and analysis; in the present paper, Part I, attention is focused on the propagation of one-dimensional waves in strips of natural, latex and synthetic, nitrile rubber. Tensile wave propagation experiments were conducted at high strain rates by holding one end fixed and displacing the other end at a constant velocity. A high-speed video camera was used to monitor the motion and to determine the evolution of strain and particle velocity in rubber strips. Analysis of the response through the theory of finite waves indicated a need for an appropriate constitutive model for rubber; by quantitative matching between the experimental observations and analytical predictions, an appropriate instantaneous elastic response for the rubbers was obtained. This matching process suggested that a simple power-law constitutive model was capable of representing the high strain-rate response for both rubbers used.     

A constitutive equation for the Mullins effect in stress controlled uniaxial extension experiments

A one-dimensional constitutive equation for the Mullins effect in rubberlike materials, which is motivated by the two phase microstructural material model proposed by Mullins and Tobin [I], is developed in [2]. The constitutive equation is used in [2] to predict the effect of stress softening on the small transverse vibration of a rubber string loaded in uniaxial extension. The two functions which comprise the constitutive equation were assumed to be monotone, but no further analysis of the actual nature of these functions was necessary.In this paper, we examine more closely how the results of a stress controlled uniaxial extension experiment can be used to gain insight into the specific nature of the microstructural strain and the strain amplification functions which comprise the constitutive equation. We examine experimental representations of the two functions which are independent of any special microstructural interpretations. Stress controlled uniaxial extension experiments with buna-n, neoprene, and silicone rubber cords are examined. We demonstrate how the experimental data can be applied to yield representations of the functions of interest to within a multiplicative constant; but no attempt is made to find specific analytical representations of these functions. For buna-n and neoprene samples, we observe behavior consistent with our monotone assumptions, while anomalous behavior is observed with silicone rubber.