Large strain field near a crack tip in a rubber sheet

The distribution of stress-strain near a crack tip in a rubber sheet is investigated by employing the constitutive relation given by Gao (1997). It is shown that the crack tip field is composed of two shrinking sectors and one expanding sector. The stress state near the crack tip is in uniaxial tension. The analytical solutions are obtained for both expanding and shrinking sectors.     

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.     

Non-linear elastic non-uniform torsion of bars of arbitrary cross-section by BEM

In this paper an elastic non-uniform torsion analysis of simply or multiply connected cylindrical bars of arbitrary cross-section taking into account the effect of geometric non-linearity is presented employing the boundary-element(BE) method. The torque-rotation relationship is computed based on the finite-displacement (finite-rotation) theory, that is the transverse displacement components are expressed so as to be valid for large rotations and the longitudinal normal strain includes the second-order geometric non-linear term often described as the “Wagner strain”. The proposed formulation does not stand on the assumption of a thin-walled structure and therefore the cross-section's torsional rigidity is evaluated exactly without using the so-called Saint-Venant's torsional constant. The torsional rigidity of the cross-section is evaluated directly employing the primary warping function of the cross-section depending on its shape. Three boundary-value problems with respect to the variable along the beam axis angle of twist, to the primary and to the secondary warping functions are formulated. The first one, employing the Analog Equation Method (a BEM-based method), yields a system of non-linear equations from which the angle of twist is computed by an iterative process. The remaining two problems are solved employing a pure BE method. Numerical results are presented to illustrate the method and demonstrate its efficiency and accuracy. The developed procedure retains most of the advantages of a BEM solution over a pure domain discretization method, although it requires domain discretization.     

Analysis of a rubber cone tensioned by a concentrated force

Using the constitutive equation of a rubber-like materials given by Gao (1997), this paper investigates the problem of a cone under tension of a concentrated force at its apex. Under consideration is the axial-symmetry case and the large strain is taken into account. The stress strain fields near the apex are obtained by both asymptotic analysis and finite element calculation. The two results are consistent well. When the cone angle is 180°, the solution becomes that of non-linear Boussinesq's problem for tension case.     

基于对数应力率大应变本构关系的参数标定研究

在所有率型弹塑性本构模型中,只有对数应力率对应的本构模型能够满足自适应准则.基于对数应力率,采用实心圆轴扭转实验,对大应变弹塑性本构模型中的参数标定问题进行了讨论.推导出了考虑Swift效应时端部自由实心圆轴扭转变形的变形率、对数旋率、Kirchhoff应力及Kirchhoff应力的对数应力率.对于等向强化大应变弹塑性本构关系,给出了由实心圆轴扭转实验标定的、基于Kirchhhoff应力对数应力率的本构关系中塑性刚度函数的表达式.分析了扭转圆轴的Swift效应对塑性刚度函数的影响.结果表明,实心圆轴扭转的轴向伸长变形和径向变形对基于对数应力率大应变本构关系中的塑性刚度函数都有影响.当不考虑Swift效应时,所得塑性刚度函数表达式与不考虑Swift效应时基于Jaumann应力率的塑性刚度函数表达式相同.     

The application of graph theory to the limit analysis for torsion of thin-walled bars with complex cross-sections

Graph theory is employed in this paper as a means to establish the topological model of complex thin-walled cross-sections. On this basis, the upper and lower bound theorems of the plastic limit analysis are applied to the analysis of the plastic limit shear flows on the cross-section of thin-walled bars under St. Venant torsion. Corresponding mathematical programming problems are formulated and their duality is shown. After solving the linear programming problem corresponding to the lower bound theorem, the limit torsional moment of a thin-walled cross-section can be calculated according to the shear stress distribution in the limit state. The formula for calculating the limit torsional moment is given. Furthermore, the limit state of thin-walled cross-sections under St. Venant torsion is also discussed and the concept of the limit tree is introduced. A computer program has been developed by the author. Results calculated by the program for typical complex cross-sections are given.     

Finite viscoelasticity of filled rubber: experiments and numerical simulation

Summary  Constitutive equations are derived for the viscoelastic behavior of particle-re-inforced elastomers at isothermal deformation with finite strain. A filled rubber is thought of as a composite medium where inclusions with high and low concentrations of junctions between chains are randomly distributed in the bulk material. The characteristic length of the inhomogeneities is assumed to be small compared to the size of the specimen and substantially exceed the radius of gyration for macromolecules. Inclusions with high concentration of junctions are associated with regions of suppressed mobility of chains that surround isolated clusters and/or the secondary network of filler. Regions with low concentration of junctions arise during the preparation process due to a heterogeneity in the spatial distribution of the cross-linker and the filler. With reference to the concept of transient networks, the time-dependent response of an elastomer is attribute d to thermally activated rearrangement of strands in the domains with low concentration of junctions. Stress–strain relations for particle-reinforced rubber are developed by using the laws of thermodynamics. Adjustable parameters in the constitutive equations are found by fitting experimental data in tensile relaxation tests for several grades of unfilled and carbon black-filled rubber. It is demonstrated that even at moderate finite deformations (with axial elongations up to 100%), the characteristic rate of relaxation is noticeably affected by strain. Unlike glassy polymers, where the rate of relaxation increases with longitudinal strain, the growth of the elongation ratio results in a decrease in the relaxation rate for natural rubber (unfilled or particle-reinforced). The latter may be explained by (partial) crystallization of chains in the regions with low concentration of junctions. Received 16 October 2001; accepted for publication 25 June 2002 Present address: A. D. Drozdov Department of Production, Aalborg University, Fibigerstraede 16, DK-9220 Aalborg, Denmark We would like to express our gratitude to Dr. K. Fuller (TARRC, UK) for providing us with rubber specimens and to Prof. P. Haupt and Dr. S. Hartmann (University of Kassel, Germany) for sending their experimental data. We are indebted to Mr. G. Seifritz for his assistance in performing mechanical tests. ADD acknowledges stimulating discussions with Prof. N. Aksel (University of Bayreuth, Germany).     

Stress-strain field near the notch tip of a rubber sheet

Analized in this paper is the elastostatic field near a notch tip in a rubber-like thin sheet. The asymptotic equations for the notch tip field are derived based on the constitutive relation given by Ref. [1]. Near field solutions are obtained in regions that decreases and increases in size as the notch tip is approached. Their singular character depends on the constitutive parameters as well as the angle of notch that is evaluated numerically. The project supported by the National Natural Science Foundation of China     

Torsional vibrational properties of rubber cylinders

This paper describes an experimental and theoretical study of the torsional vibrational properties of rubber cylinders. Experiments were carried out on long cylinders using a modified Weissenberg Rheogoniometer. A comparison is made with a theory based on a modified generalized Maxwell model.     

实心圆试件扭转试验确定大应变本构关系

提出并完成了通过实心圆轴扭转试验建立大应变本构关系的方法。它比单向拉伸试验所得到的本构关系更为精确。因拉伸实验变形较大时试件伸长和变细对测量结果有影响;尤其在“颈缩”后,很难对有关力学量作有效测量和分析。即扭转本构关系的描绘更为完整。以低碳钢为例,扭转本构关系所描述的有效区间比拉伸本构关系大十余倍。该方法将有利于探讨研究更大应变下的材料力学行为。     

Deflections of a rubber membrane

This paper deals with the problem of the transverse deflection of a natural rubber membrane that is fixed along a circular boundary. Uniaxial experiments were performed in order to characterize the constitutive behaviour of the rubber material in terms of several constitutive models available in the literature. These constitutive models were used to develop computational estimates for the quasi-static load-displacement response of a rigid spherical indentor that deflects the rubber membrane in a controlled fashion and to determine the deflected shape of the membrane at specified load levels. Both axisymmetric and asymmetric deflections of the rubber membrane were investigated. The paper provides a comparison of the experimental results for the membrane deflections with results derived from computational simulations.     

壳体结构弹塑性大应变动态响应分析

本文利用假设剪切应和薄膜应变场的新型退化壳单元，并考虑大应主影响对壳体结构弹塑性动态分析进行了研究，并编制了相应的程度，此程序采用修正的Ｌａｇｒａｎｇｉａｎ公式，它不仅考虑了大变形大应变的影响，还对单元的节点坐标和厚度在每个时间步后都给予修正，通过对一些算例计算表明，此程序数值精度很好，可用于壳体结构的工程分析，并得出结论，对于应变问题考虑和忽略大应变影响将有较大的不同。     

大直径SHPB装置的数值模拟及实验误差分析

介绍了用于混凝土冲击试验的Φ100SHPB装置。利用数值模拟和应变测试相结合的方法讨论了应力波在该装置中的传播特性。二维效应导致应力波在大直径SHPB装置传播过程的弥散现象。同时在实验的大部分时间内大尺寸试样的应力(应变)存在不均匀性。因此大直径SHPB装置进行混凝土冲击压缩应力-应变曲线的实验研究会产生误差。利用数值模拟分析了误差影响程度,同时验证了用混凝土试样应变和透射杆应力得到应力-应变曲线初始段的可行性。     

Fatigue strength of metals and composites under repeated high-cycle torsion

The limiting stresses are determined and constant fatigue life diagrams for high-cycle torsion with repeated stress cycle are plotted using the limiting-state models obtained based on the hypothesis of a unified constant-life diagram, which is invariant to the number of cycles to failure. The unified constant-life diagram is given by a transcendental power function whose exponent is an additional parameter characterizing the sensitivity of the material to the asymmetry of the stress cycle. The calculated results and experimental data for carbon and alloyed steels and composite materials are in good agreement __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 2, pp. 19–28, February 2008.     

Establishing strain response model of vibrational structure by strain modal analysis approach

An analysis of the stress field under a dynamic load in strength design and fatigue life estimation is a topic of great importance. Based on the theory of vibrational displacement modal analysis, the principle and techniques of vibrational strain modal analysis are developed. The peculiarities of the strain transfer function matrix, the test approaches of it and the relation with the displacement transfer function are expounded. Resistance strain gauges are used to measure the strain transfer function and the modal parameters are identified on this basis. After that, the strain response expression and thereby the stress response expression are obtained.     

Nonuniform warping torsion of orthotropic composite beams

Summary  The proposed work deals with the behaviour of orthotropic composite beams subjected to torsion, taking into consideration the effects due to nonuniform warping. We start from a one-dimensional approach based on the integral parameter method and consider any complex cross-sectional shape. The basic relations obtained lead to numerical computation based on the finite element method. As a special case of the general results obtained, an example of a sandwich beam is worked out in detail. Received 11 May 1999; accepted for publication 30 November 1999     

Boundary element method to solve torsion problem of cracked cylinder

Separating the discontinuous solution by use of the single crack solution, together with the regular solution of harmonic function, the torsion problem of a cracked cylinder is reduced to solving a set of mixed-type integral equations and its numerical technique is then proposed by combining the numerical method of singular integral equation with the boundary element method. Several numerical examples are calculated which will be useful to engineering practice. The method proposed is characterized by its fine accuracy and convenience for using, which can be extended to the cases of multiple crack.The project supported by National Natural Science Foundation of China.     

HOMOGENIZATION—BASED TOPOLOGY DESIGN FOR PURE TORSION OF COMPOSITE SHAFTS

In conjunction with the homogenization theory and the finite element method, the mathematical models for designing the corss-section of composite shafts by maximizing the torsion rigidity are developed in this paper. To obtain the extremal torsion rigidity, both the cross-section of the macro scale shaft and the representative microstructure of the composite material are optimized using the new models. The micro scale computational model addresses the problem of finding the periodic microstructures with extreme shear moduli. The optimal microstructure obtained with the new model and the homogenization method can be used to improve and optimize natural or artificial materials. In order to be more practical for engineering applications, cellular materials rather than ranked materials are used in the optimal process in the existence of optimal bounds for the elastic properties. Moreover, the macro scale model is proposed to optimize the cross-section of the torsional shaft based on the tailared composites. The validating optimal results show that the models are very effective in obtaining composites with extreme elastic properties, and the cross-section of the composite shaft with the extremal torsion rigidity. The project supported by the National Natural Science Foundation of China (10172078 and 10102018)     

A Molecular-Statistical Basis for the Gent Constitutive Model of Rubber Elasticity

Molecular constitutive models for rubber based on non-Gaussian statistics generally involve the inverse Langevin function. Such models are widely used since they successfully capture the typical strain-hardening at large strains. Limiting chain extensibility constitutive models have also been developed on using phenomenological continuum mechanics approaches. One such model, the Gent model for incompressible isotropic hyperelastic materials, is particularly simple. The strain-energy density in the Gent model depends only on the first invariant I ₁ of the Cauchy–Green strain tensor, is a simple logarithmic function of I ₁ and involves just two material parameters, the shear modulus μ and a parameter J _m which measures a limiting value for I ₁−3 reflecting limiting chain extensibility. In this note, we show that the Gent phenomenological model is a very accurate approximation to a molecular based stretch averaged full-network model involving the inverse Langevin function. It is shown that the Gent model is closely related to that obtained by using a Padè approximant for this function. The constants μ and J _m in the Gent model are given in terms of microscopic properties. Since the Gent model is remarkably simple, and since analytic closed-form solutions to several benchmark boundary-value problems have been obtained recently on using this model, it is thus an attractive alternative to the comparatively complicated molecular models for incompressible rubber involving the inverse Langevin function. This revised version was published online in June 2006 with corrections to the Cover Date.