Mechanical response of fiber-reinforced incompressible non-linearly elastic solids

The mechanical response of some fiber-reinforced incompressible non-linearly elastic solids is examined under homogeneous deformation. In particular, the materials under consideration are neo-Hookean models augmented with a function that accounts for the existence of a unidirectional reinforcement. This function endows the material with its anisotropic character and is referred to as a reinforcing model. The nature of the anisotropy considered has a particular influence on the shear response of the material, in contrast to previous analyses in which the reinforcing model was taken to depend only on the stretch in the fiber direction.     

Critical stretch for formation of a cylindrical void in a compressible hyperelastic material

Formation of a cylindrical void in an infinitely long compressible hyperelastic cylinder under axial and radial stretches is examined. The cavitation phenomenon is viewed here as a bifurcation of a solution with a cavity from the homogeneously deformed configuration, taking place when the applied radial stretch reaches a certain critical value. This amounts to modeling the underlying phenomenon as a kind of elastic instability. A formulation of shooting-method type is presented to derive an equation which gives the critical radial stretch for a prescribed axial stretch. For a special class of hyperelastic solids, called modified Blatz-Ko material, the obtained equation leads to an explicit expression for the critical stretches and stresses. Some analytical as well as numerical calculations are carried out to explicitly obtain the critical values for cavitation. The results are summarized in the form of cavitation curves in the two-dimensional space of axial and radial stretches or stresses. Influence of the axial stretch on the critical radial stretch is discussed. Throughout the paper, the corresponding results for a spherically symmetric void formation are referred to when appropriate and compared with the cylindrical case of the present interest. It is then indicated that in the state of equitriaxial stretch, cavitation into a cylindrical shape is likely to occur at lower stretch than into a spherical one.     

Mechanical response of neo-Hookean fiber reinforced incompressible nonlinearly elastic solids

The mechanical behavior of an incompressible neo-Hookean material, directionally reinforced by neo-Hookean fibers, is examined under homogeneous deformations. A composite model for this transversely isotropic material is developed based on a multiplicative decomposition of the deformation gradient which considers interaction between the fiber and the matrix. The so-called standard reinforcing model exhibits non-monotonic behavior in compression. The present composites-based approach leads to a modification of the standard reinforcing model in which monotonic behavior in compression is observed. This stems from the micromechanical basis of the model in which the fiber is treated as a neo-Hookean material. The conditions for loss of monotonicity and positivity in the stress-shear behavior in off-axis simple 2D shear are also obtained.     

Kink Surfaces in a Directionally Reinforced neo-Hookean Material under Plane Deformation: I. Mechanical Equilibrium

Stationary kinks (elastostatic shocks) are examined in the context of a base neo-Hookean response augmented with unidirectional reinforcing that is characterized by a single additional constitutive parameter for the additional fiber reinforcing stiffness. Previous work has shown that such a transversely isotropic material can lose ellipticity in plane deformation if the reinforcing is sufficiently large and the fiber direction is sufficiently compressed. Here we show that the same reinforcing levels can give rise to piecewise smooth plane deformations separated by a plane stationary kink. Attention is restricted to deformations in which, on one side of the kink, the load axis is aligned with the fiber axis. Then the fiber stretch on this side of the kink is a natural load parameter. It is found that such a deformation can support a planar kink for a certain range of this load parameter. This range is dependent on the reinforcing parameter, and can even involve fiber extension if the reinforcing is sufficiently large. The set of all deformation states on the other side of the kink is precisely characterized in terms of a one-parameter family of (kink orientation, kink strength)-pairs. The results are interpreted in terms of the associated fiber alignment discontinuity and fiber stretch discontinuity.     

Instabilities and loss of ellipticity in fiber-reinforced compressible non-linearly elastic solids under plane deformation

In a recent paper we examined the loss of ellipticity and its interpretation in terms of fiber kinking and other instability phenomena in respect of a fiber-reinforced incompressible elastic material. Here we provide a corresponding analysis for fiber-reinforced compressible elastic materials. The analysis concerns a material model which consists of an isotropic base material augmented by a reinforcement dependent on the fiber direction. The assessment of loss of ellipticity can be cast in terms of the eigenvalues of the acoustic tensors associated with the isotropic and anisotropic parts of the strain-energy function. For the anisotropic part, two different reinforcing models are examined and it is shown that, depending on the choice of model and whether the fiber is under compression or extension, loss of ellipticity may be associated with, in particular, a weak surface of discontinuity normal to or parallel to the deformed fiber direction or at an intermediate angle. Under compression the associated failure interpretations include fiber kinking and fiber splitting, while under extension fiber de-bonding and matrix failure are included.     

A three-dimensional nonlinear model for dissipative response of soft tissue

A set of three-dimensional constitutive equations is proposed for modeling the nonlinear dissipative response of soft tissue. These constitutive equations are phenomenological in nature and they model a number of physical features that have been observed in soft tissue. The equations model the tissue as a composite of a purely elastic component and a dissipative component, both of which experience the same total dilatation and distortion. The stress response of the purely elastic component depends on dilatation, distortion and the stretch of material fibers, whereas the stress response of the dissipative component depends on distortional deformation only. The equations are hyperelastic in the sense that the stress is obtained by derivatives of a strain energy function, and they are properly invariant under superposed rigid body motions. In contrast with standard viscoelastic models of tissues, the proposed constitutive model includes the total deformation rate in evolution equations that can reproduce the observed physical feature that the hysteresis loops of most biological soft tissues are nearly independent of strain rate (Biomechanics, Mechanical Properties of Living Tissues, second ed. (1993)). Material constants are determined which produce good agreement with uniaxial stress experiments on superficial musculoaponeurotic system and facial skin.     

Finite elastic deformations of transversely isotropic circular cylindrical tubes

We consider the finite radially symmetric deformation of a circular cylindrical tube of a homogeneous transversely isotropic elastic material subject to axial stretch, radial deformation and torsion, supported by axial load, internal pressure and end moment. Two different directions of transverse isotropy are considered: the radial direction and an arbitrary direction in planes normal locally to the radial direction, the only directions for which the considered deformation is admissible in general. In the absence of body forces, formulas are obtained for the internal pressure, and the resultant axial load and torsional moment on the ends of the tube in respect of a general strain-energy function. For a specific material model of transversely isotropic elasticity, and material and geometrical parameters, numerical results are used to illustrate the dependence of the pressure, (reduced) axial load and moment on the radial stretch and a measure of the torsional deformation for a fixed value of the axial stretch.     

Remarks on the Behavior of Simple Directionally Reinforced Incompressible Nonlinearly Elastic Solids

The effect of directional reinforcing in generating qualitative changes in the mechanical response of a base neo-Hookean material is examined in the context of homogenous deformation. Single axis reinforcing giving transverse isotropy is the major focus, in which case a standard reinforcing model is characterized by a single constitutive reinforcing parameter. Various qualitative changes in the mechanical response ensue as the reinforcing parameter increases from the zero-value associated with neo-Hookean response. These include (in order): the existence of a limiting contractive stretch for transverse-axis tensile load; loss of monotonicity in off-axis simple shear; loss of monotonicity in on-axis compression; loss of positivity in the stress-shear product in off-axis simple shear; and loss of monotonicity for plane strain in on-axis compression. The qualitative changes in the simple shear response are associated with stretch relaxation in the reinforcing direction due to finite rotation. This revised version was published online in August 2006 with corrections to the Cover Date.     

Large deformation of anisotropic austenitic stainless steel sheets at room temperature: Multi-axial experiments and phenomenological modeling

A newly developed multi-axial testing technique for sheet materials is employed to investigate the inelastic response of a temper-rolled stainless steel 301LN under isothermal quasi-static loading conditions at room temperature. The experimental technique consists of a flat sheet specimen, which is subject to combinations of shear and normal loading using a custom-made dual-actuator system. The large deformation behavior under monotonic loading is determined along more than 20 distinct radial paths in the stress space. The experimental results indicate that Hill's quadratic yield function along with an associated flow rule provides a good approximation of the initial yield behavior of this anisotropic two-phase FCC/BCC sheet material. Based on the experimental data for radial monotonic loading, it is concluded that conventional isotropic-kinematic hardening models cannot successfully describe the strain hardening of this austenitic steel. Instead, a non-associated anisotropic hardening model is proposed that relates the increase in yield strength to an isothermal martensitic transformation kinetics law. The comparison of the model predictions with the experimental results shows very good agreement for all biaxial and uniaxial experiments.     

钢筋混凝土靶侵彻的可压缩弹-塑性动态空腔膨胀阻力模型

在Forrestal素混凝土靶侵彻的可压缩弹-塑性球形动态空腔膨胀理论模型基础上,考虑粉碎区以内钢筋对混凝土的环向约束作用,提出了一个适用于刚性弹侵彻钢筋混凝土靶的阻力模型。论文通过体积配筋率的引入,获得了钢筋混凝土靶空腔表面径向应力的理论解,并讨论了配筋率对空腔壁面径向应力及各分区大小的影响。结果表明:钢筋对混凝土的环向约束效应影响了空腔膨胀过程中混凝土各区域的大小分布,并提高了空腔表面的径向应力。     

Gas bubble oscillations in a fluid in the presence of 2:1 resonance between the radial and deformation oscillation frequencies

Small nonlinear oscillations of an ellipsoidal bubble in a fluid in the presence of 2:1 frequency resonance between the radial and ellipsoidal modes are considered. The equations of motion are reduced to Hamiltonian form. The quadratic and cubic terms are taken into account in the expansion of the Hamiltonian. The Hamilton function is transformed to the normal form using the invariant normalization method in the first approximation. This makes it possible to construct an analogy between the system considered and the well-known problem of a pendulous spring. The radial and ellipsoidal bubble oscillation modes correspond to the vertical and horizontal coordinates of a material point, respectively. In the absence of resonance the solution of the nonlinear equations differs from the solution of the linear equations by only a small (quadratic in the amplitude) change in the oscillation frequency. In the resonance case the radial and ellipsoidal oscillation modes periodically change places and the energy of one mode is converted into that of the other. The interest in the system in resonance is associated with precisely this fact. The question of the dissipation effect in real media is considered. The decay rate depends significantly on the physical properties of the material and, in certain special cases, can be small enough for the energy transfer effect to manifest itself.     

A composites-based hyperelastic constitutive model for soft tissue with application to the human annulus fibrosus

This paper presents a composites-based hyperelastic constitutive model for soft tissue. Well organized soft tissue is treated as a composite in which the matrix material is embedded with a single family of aligned fibers. The fiber is modeled as a generalized neo-Hookean material in which the stiffness depends on fiber stretch. The deformation gradient is decomposed multiplicatively into two parts: a uniaxial deformation along the fiber direction and a subsequent shear deformation. This permits the fiber-matrix interaction caused by inhomogeneous deformation to be estimated by using effective properties from conventional composites theory based on small strain linear elasticity and suitably generalized to the present large deformation case. A transversely isotropic hyperelastic model is proposed to describe the mechanical behavior of fiber-reinforced soft tissue. This model is then applied to the human annulus fibrosus. Because of the layered anatomical structure of the annulus fibrosus, an orthotropic hyperelastic model of the annulus fibrosus is developed. Simulations show that the model reproduces the stress-strain response of the human annulus fibrosus accurately. We also show that the expression for the fiber-matrix shear interaction energy used in a previous phenomenological model is compatible with that derived in the present paper.     

Numerical study of strength properties for a composite material with short reinforcing fibers

A numerical approach to the determination of strength properties for concrete with short reinforcing fibers on the basis of the finite element method is proposed. The mathematical model takes into account various processes of nonlinear deformation of the concrete matrix under compressive and tensile loads, the possibility of developing the inelastic strains in the concrete matrix and reinforcing fibers and the nonlinear interaction between them. The effect of fiber concentration, various loading surfaces for the material matrix, and the bonding type on the deformation of a composite material is analyzed. Numerical examples of strength analysis are given.     

Stress field of orthotropic cylinder subjected to axial compression

Based on the constancy hypothesis of material volume, the circumferential and radial stresses of a cylinder specimen are analyzed when the cylinder is subject to a loading along the axial direction. The circumferential and radial stress distribution is a power function of radius parameter when the constitutive relation of specimen material is orthotropic. The stress distribution is a quadratic function of radius parameter for transversely isotropic material. Along the cylinder axial line, the circumferential and radial stresses are maximum and equal to each other. In the circumference boundary surface, the radial stress is zero and the circumferential stress value is minimal. The failure theory of maximum tensile circumferential strain is applied to calculate the critical axial loading. The circumference-boundary-layer failure criterion of orthotropic cylinders is described with the Hill-Tsai strength theory. The obtained strength theory is related to axial stress and mechanical properties of specimen material and to the specimen axialdeformation strain rate and the change rate of strain rate.     

Cavity formation and singular periodic oscillations in isotropic incompressible hyperelastic materials

In this paper, a dynamical problem is considered for an incompressible hyperelastic solid sphere composed of the classical isotropic neo-Hookean material, where the sphere is subjected to a class of periodic step radial tensile loads on its surface. A second-order non-linear ordinary differential equation that describes cavity formation and motion is proposed. The qualitative properties of the solutions of the equation are examined. Correspondingly, under a prescribed constant dead-load, it is proved that a cavity forms in the sphere as the dead-load exceeds a certain critical value and the motion of the formed cavity presents a class of singular periodic oscillations. Under periodic step loads, the existence conditions for periodic oscillation of the formed cavity are determined by using the phase diagrams of the motion equation of cavity. In each section, numerical examples are also carried out.     

翻转后的不可压缩neo-Hookean圆柱管的有限变形

研究了由径向横观各向同性不可压缩的neo-Hookean材料组成的圆柱形管在翻转后的有限变形问题。利用材料的不可压缩条件和半逆解法对相应的数学模型进行求解,并根据边界条件得到了翻转后的圆柱形管的内半径以及轴向伸长率应满足的非线性方程组。通过数值算例讨论了材料参数和结构参数对翻转后圆柱形管的内半径以及轴向伸长率变化的影响。结果表明：初始厚度对翻转后圆柱形管的内半径与轴向伸长率没有本质上的影响;而径向各向异性参数却有本质上的影响,特别是在轴向伸长率方面。     

可压超弹性-塑性材料中的空穴生成

本文研究了材料的弹塑性性质对球体中空穴生成问题的影响，材料的弹性用一种可压超弹性材料的本构关系来描述，材料的塑性用满足材料的不可压条件和Tresca屈服条件的理想塑性材料的本构关系来描述。这类超弹性．塑性材料中可以发生空穴的生成现象，得到了在表面拉伸作用下球体中空穴生成时空穴半径与临界拉伸之间的关系式和临界拉伸。球体的变形可分为弹-塑性变形阶段和完全塑性变形阶段，球体中心首先形成塑性变形区域，并有空穴的突然生成；塑性变形区域能够快速增长，并且使球体很快进入完全塑性变形阶段；空穴在弹-塑性变形阶段迅速增长，但进入完全塑性变形阶段后增长较慢。同时给出了不同变形阶段球体中的应力分布。数值结果表明材料的塑性性质对材料中的空穴生成有明显的影响。     

Homogenization of long fiber reinforced composites including fiber bending effects

This paper presents a homogenization method, which accounts for intrinsic size effects related to the fiber diameter in long fiber reinforced composite materials with two independent constitutive models for the matrix and fiber materials. A new choice of internal kinematic variables allows to maintain the kinematics of the two material phases independent from the assumed constitutive models, so that stress–deformation relationships, can be expressed in the framework of hyper-elasticity and hyper-elastoplasticity for the fiber and the matrix materials respectively. The bending stiffness of the reinforcing fibers is captured by higher order strain terms, resulting in an accurate representation of the micro-mechanical behavior of the composite. Numerical examples show that the accuracy of the proposed model is very close to a non-homogenized finite-element model with an explicit discretization of the matrix and the fibers.     

Hydraulic Fracture Propagation with 3-D Leak-off

Hydraulic fracture models typically couple a fracture elasticity model with a geological reservoir model to forecast the rate of fluid leak-off from the propagating fracture. The most commonly used leak-off model is that originally specified by Carter, which involves the assumption that the fracture is embedded within an infinite homogenous porous medium where flow only occurs perpendicular to the fracture plane. The objectives of this paper are: (1) to show that assuming one-dimensional leak-off can lead to erroneous conclusions, (2) to present a robust numerical methodology for simulating three-dimensional leak-off from propagating hydraulic fractures, and (3) to present and compare a new analytical method based on assuming three-dimensional flow of an incompressible fluid through an incompressible porous formation from a circular planar fracture. Provided the fluid and formation compressibility can be ignored within the reservoir flow model, the three-dimensional leak-off from a circular planar fracture can be written in closed-form as a function, which depends linearly on fracture pressure and radial extent. This simple expression for leak-off can be easily coupled to a range of circular fracture elasticity models. As a comparison example, the Carter model, our new function and a three-dimensional numerical model of the full problem are coupled to the PK-radial fracture model. Comparison with the numerical model shows that our new function overestimates fracture growth during intermediate times but accurately predicts both the early and late-time asymptotic behavior. In contrast, the Carter model fails to replicate both the early and late-time asymptotic behavior. Our new function additionally improves on the Carter model by not requiring the evaluation of convolution integrals and allowing easy evaluation of both the spatial leakage flux distribution across the fracture face and the three-dimensional pressure distribution within the porous formation.     

A direct approach to fiber and membrane reinforced bodies. Part I. Stress concentrated on curves for modelling fiber reinforced materials

An approach is outlined to the equilibrium in fiber-reinforced materials in which the fibers are modeled as curves or lines with concentrated material properties. The system of forces representing the interaction of the fibers with the bulk matter is analyzed, and equilibrium of forces is derived from global laws. The displacements of the bulk matter are assumed to have continuous extension to the fibers. This forces the set of admissible deformations superquadratically integrable. This in turn forces the energy of the bulk of superquadratic growth. The material of the bulk matrix therefore cannot be linearly elastic. The energy of fibers can have a slower growth and can be quadratic. A formal set of assumptions is given under which an equilibrium state of minimum energy exists in the given external conditions. A weak form of equilibrium equations is derived for this equilibrium state. An explicitly calculable axisymmetric example is presented with an isotropic and quadratic energy of the matrix (linear elasticity) and linearly stretchable fiber. Since the superquadratic growth assumption is not satisfied, some peculiar features of the solution arise, such as the infinite limit of the radial displacement near the fiber. Nevertheless, from the obtained solution, we can compute the normal force in the fiber and the shear stress at the interface.