The Plane Stress Crack-Tip Field for an Incompressible Rubber Material

The crack-tip field under plane stress condition for an incompressible rubbermaterial￣[1] is investigated by. the use of the fully nonlinear equilibrium theory. It isfound thai the crack-tip field is composed of two shrink sectors and one expansion se-ctor. At the crack-tip, stress and strain possess the singularity of R￣(-1) and R￣(-1n), respec-tively, (R is the distance to the crack-tip before deformation, n is the material const-ant). When the crack-tip is approached, the thickness of the sheet shrinks to zerowith the order of R￣(1.4n). The results obtained in this paper are consistent with that ob-tained in [8] when s→∞ .

CAVITY FORMATION AND ITS VIBRATION FOR A CLASS OF GENERALIZED INCOMPRESSIBLE HYPER-ELASTIC MATERIALS

The problem of radial symmetric motion for a solid sphere composed of a class of generalized incompressible neo-Hookean materials, subjected to a suddenly applied surface tensile dead load, is examined.The analytic solutions for this problem and the motion equation of cavity that describes cavity formation and growth with time are obtained. The e?ect of radial perturbation of the materials on cavity formation and its motion is discussed. The plane of the perturbation parameters of the materials is divided into four regions. The existential conditions and qualitative properties of solutions of the motion equation of the cavity are studied in di?erent parameters’ regions in detail. It is proved that the cavity motion with time is a nonlinear periodic vibration. The vibration center is then determined.

QUALITATIVE ANALYSIS OF SPHERICAL CAVITY NUCLEATION AND GROWTH FOR INCOMPRESSIBLE GENERALIZED VALANIS-LANDEL HYPERELASTIC MATERIALS

A cavitated bifurcation problem is examined for a sphere composed of a class of generalized Valanis-Landel materials subjected to a uniform radial tensile dead-load. A cavitated bifurcation equation is obtained. An explicit formula for the critical value associated with the vari- ation of the imperfection parameters is presented. The distinguishing between the left-bifurcation and right-bifurcation of the nontrivial solution of the cavitated bifurcation equation at the critical point is made. It is proved that there exists a secondary turning bifurcation point on the nontrivial solution branch, which bifurcates locally to the left. It is shown that the dimensionless cavitated bifurcation equation is equivalent to normal forms with single-sided constraint conditions at the critical point by using the singularity theory. The stability and catastrophe of the solutions of the cavitated bifurcation equation are discussed.

Cavitation for Incompressible Anisotropic Hyper-Elastic Materials

The cavitation problem in a solid sphere composed of an incompressible anisotropic hyper elastic material under a uniform radial tensile dead load was examined. A new analytical solution was obtained. The stress contributions were given and the jumping and concentration of stresses were discussed. The stability of solutions and the effect of the degree of anisotropy of the material were analyzed.

CAVITATED BIFURCATION FOR INCOMPRESSIBLE HYPERELASTIC MATERIAL

The spherical cavitated bifurcation for a hyperelastic solid sphere made of the incompressible Valanis-Landel material under boundary dead-loading is examined. The analytic solution for the bifurcation problem is obtained. The catastrophe and concentration of stresses are discussed. The stability of solutions is discussed through the energy comparison. And the growth of a pre-existing micro-void is also observed.

An improved dynamic model for a silicone material beam with large deformation

Dynamic modeling for incompressible hyperelastic materials with large deformation is an important issue in biomimetic applications. The previously proposed lower-order fully parameterized absolute nodal coordinate formulation（ANCF） beam element employs cubic interpolation in the longitudinal direction and linear interpolation in the transverse direction, whereas it cannot accurately describe the large bending deformation. On this account, a novel modeling method for studying the dynamic behavior of nonlinear materials is proposed in this paper. In this formulation, a higher-order beam element characterized by quadratic interpolation in the transverse directions is used in this investigation. Based on the Yeoh model and volumetric energy penalty function, the nonlinear elastic force matrices are derived within the ANCF framework. The feasibility and availability of the Yeoh model are verified through static experiment of nonlinear incompressible materials. Furthermore,dynamic simulation of a silicone cantilever beam under the gravity force is implemented to validate the superiority of the higher-order beam element. The simulation results obtained based on the Yeoh model by employing three different ANCF beam elements are compared with the result achieved from a commercial finite element package as the reference result. It is found that the results acquired utilizing a higher-order beam element are in good agreement with the reference results,while the results obtained using a lower-order beam element are different from the reference results. In addition, the stiffening problem caused by volumetric locking can be resolved effectively by applying a higher-order beam element. It is concluded that the proposed higher-order beam element formulation has satisfying accuracy in simulating dynamic motion process of the silicone beam.

QUALITATIVE ANALYSIS OF DYNAMICAL BEHAVIOR FOR AN IMPERFECT INCOMPRESSIBLE NEO-HOOKEAN SPHERICAL SHELL

The radial symmetric motion problem was examined for a spherical shell composed of a class of imperfect incompressible hyper-elastic materials, in which the materials may be viewed as the homogeneous incompressible isotropic neo-Hookean material with radial perturbations. A second-order nonlinear ordinary differential equation that describes the radial motion of the inner surface of the shell was obtained. And the first integral of the equation was then carded out. Via analyzing the dynamical properties of the solution of the differential equation, the effects of the prescribed imperfection parameter of the material and the ratio of the inner and the outer radii of the underformed shell on the motion of the inner surface of the shell were discussed, and the corresponding numerical examples were carded out simultaneously. In particular, for some given parameters, it was proved that, there exists a positive critical value, and the motion of the inner surface with respect to time will present a nonlinear periodic oscillation as the difference between the inner and the outer presses does not exceed the critical value. However, as the difference exceeds the critical value, the motion of the inner surface with respect to time will increase infinitely. That is to say, the shell will be destroyed ultimately.

Finite elastic torsional instability of incompressible thermo-hyperelastic cylinder

Torsional instability of an incompressible thermo-hyperelastic cylindrical rod, subjected to axial stretching and large torsions, is examined within the framework of finite elasticity. When the cylinder is stretched and twisted by a sufficiently large degree, a knot may form suddenly at one point. This inherent elastic instability is analyzed with the minimum potential energy principle and the critical values of torsion are obtained. The distribution of stresses as well as the tensile force and the torque are studied. Effect of tem- perature change is specifically discussed.

QUALITATIVE STUDY OF CAVITATED BIFURCATION FOR A CLASS OF INCOMPRESSIBLE GENERALIZED NEO-HOOKEAN SPHERES

The problem of spherical cavitated bifurcation was examined for a class of incompressible generalized neo-Hookean materials, in which the materials may be viewed as the homogeneous incompressible isotropic neo-Hookean material with radial perturbations. The condition of void nucleation for this problem was obtained. In contrast to the situation for a homogeneous isotropic neo-Hookean sphere, it is shown that not only there exists a secondary turning bifurcation point on the cavitated bifurcation solution which bifurcates locally to the left from trivial solution, and also the critical load is smaller than that for the material with no perturbations, as the parameters belong to some regions. It is proved that the cavitated bifurcation equation is equivalent to a class of normal forms with single-sided constraints near the critical point by using singularity theory. The stability of solutions and the actual stable equilibrium state were discussed respectively by using the minimal potential energy principle.

Highly accurate symplectic element based on two variational principles

For the stability requirement of numerical resultants, the mathematical theory of classical mixed methods are relatively complex. However, generalized mixed methods are automatically stable, and their building process is simple and straightforward. In this paper, based on the seminal idea of the generalized mixed methods, a simple, stable, and highly accurate 8-node noncompatible symplectic element(NCSE8) was developed by the combination of the modified Hellinger-Reissner mixed variational principle and the minimum energy principle. To ensure the accuracy of in-plane stress results, a simultaneous equation approach was also suggested. Numerical experimentation shows that the accuracy of stress results of NCSE8 are nearly the same as that of displacement methods, and they are in good agreement with the exact solutions when the mesh is relatively fine. NCSE8 has advantages of the clearing concept, easy calculation by a finite element computer program, higher accuracy and wide applicability for various linear elasticity compressible and nearly incompressible material problems. It is possible that NCSE8 becomes even more advantageous for the fracture problems due to its better accuracy of stresses.

DYNAMICAL FORMATION OF CAVITY IN A COMPOSED HYPER-ELASTIC SPHERE

The dynamical formation of cavity in a hyper_elastic sphere composed of two materials with the incompressible strain energy function, subjected to a suddenly applied uniform radial tensile boundary dead_load, was studied following the theory of finite deformation dynamics. Besides a trivial solution corresponding to the homogeneous static state, a cavity forms at the center of the sphere when the tensile load is larger than its critical value. An exact differential relation between the cavity radius and the tensile land was obtained. It is proved that the evolution of cavity radius with time displays nonlinear periodic oscillations. The phase diagram for oscillation, the maximum amplitude, the approximate period and the critical load were all discussed.

DYNAMIC CHARACTERISTICS IN INCOMPRESSIBLE HYPERELASTIC CYLINDRICAL MEMBRANES

In this paper, the dynamic characteristics are examined for a cylindrical membrane composed of a transversely isotropic incompressible hyperelastic material under an applied uniform radial constant pressure at its inner surface. A second-order nonlinear ordinary differential equation that approximately describes the radial oscillation of the inner surface of the membrane with respect to time is obtained. Some interesting conclusions are proposed for different materials, such as the neo-Hookean material, the Mooney-Rivlin material and the Rivlin-Saunders material. Firstly, the bifurcation conditions depending on the material parameters and the pressure loads are determined. Secondly, the conditions of periodic motion are presented in detail for membranes composed of different materials. Meanwhile, numerical simulations are also provided.

Comparison of two projection methods for modeling incompressible flows in MPM

Material point method（MPM）was originally introduced for large deformation problems in solid mechanics applications.Later,it has been successfully applied to solve a wide range of material behaviors.However,previous research has indicated that MPM exhibits numerical instabilities when resolving incompressible flow problems.We study Chorin＇s projection method in MPM algorithm to simulate material incompressibility.Two projection-type schemes,non-incremental projection and incremental projection,are investigated for their accuracy and stability within MPM.Numerical examples show that the non-incremental projection scheme provides stable results in single phase MPM framework.Further,it avoids artificial pressure oscillations and small time steps that are present in the explicit MPM approach.

Numerical Studies of Vanka-Type Smoothers in Computational Solid Mechanics

In this paper multigrid smoothers of Vanka-type are studied in the context of Computational Solid Mechanics(CSM).These smoothers were originally developed to solve saddle-point systems arising in the field of Computational Fluid Dynamics(CFD),particularly for incompressible flow problems.When treating(nearly)incompressible solids,similar equation systems arise so that it is reasonable to adopt the‘Vanka idea’for CSM.While there exist numerous studies about Vanka smoothers in the CFD literature,only few publications describe applications to solid mechanical problems.With this paper we want to contribute to close this gap.We depict and compare four different Vanka-like smoothers,two of them are oriented towards the stabilised equal-order Q_(1)/Q_(1)finite element pair.By means of different test configurations we assess how far the smoothers are able to handle the numerical difficulties that arise for nearly incompressible material and anisotropic meshes.On the one hand,we show that the efficiency of all Vanka-smoothers heavily depends on the proper parameter choice.On the other hand,we demonstrate that only some of them are able to robustly deal with more critical situations.Furthermore,we illustrate how the enclosure of the multigrid scheme by an outer Krylov space method influences the overall solver performance,and we extend all our examinations to the nonlinear finite deformation case.