THE INVARIANT REPRESENTATION OF SPINS WITH APPLICATIONS IN THE THEORY OF FINITE DEFORMATIONS

Based on the general solution given to a kind of linear tensor equations,the spin of a symmetric tensor is derived in an invariant form.The result is applied to find the spins of the left and the tight stretch tensors and the relation among different rotation rate tensors has been discussed.According to work conjugacy,the relations between Cauchy stress and the stresses conjugate to Hill's generalized strains are obtained.Particularly,the logarithmic strain,its time rate and the conjugate stress have been discussed in de- tail.These results are important in modeling the constitutive relations for finite deformations in continuum me- chanics.

Finite Deformation, Finite Strain Nonlinear Dynamics and Dynamic Bifurcation in TVE Solids with Rheology

This paper presents a mathematical model consisting of conservation and balance laws (CBL) of classical continuum mechanics (CCM) and ordered rate constitutive theories in Lagrangian description derived using entropy inequality and the representation theorem for thermoviscoelastic solids (TVES) with rheology. The CBL and the constitutive theories take into account finite deformation and finite strain deformation physics and are based on contravariant deviatoric second Piola-Kirchhoff stress tensor and its work conjugate covariant Green’s strain tensor and their material derivatives of up to order m and n respectively. All published works on nonlinear dynamics of TVES with rheology are mostly based on phenomenological mathematical models. In rare instances, some aspects of CBL are used but are incorrectly altered to obtain mass, stiffness and damping matrices using space-time decoupled approaches. In the work presented in this paper, we show that this is not possible using CBL of CCM for TVES with rheology. Thus, the mathematical models used currently in the published works are not the correct description of the physics of nonlinear dynamics of TVES with rheology. The mathematical model used in the present work is strictly based on the CBL of CCM and is thermodynamically and mathematically consistent and the space-time coupled finite element methodology used in this work is unconditionally stable and provides solutions with desired accuracy and is ideally suited for nonlinear dynamics of TVES with memory. The work in this paper is the first presentation of a mathematical model strictly based on CBL of CCM and the solution of the mathematical model is obtained using unconditionally stable space-time coupled computational methodology that provides control over the errors in the evolution. Both space-time coupled and space-time decoupled finite element formulations are considered for obtaining solutions of the IVPs described by the mathematical model and are presented in the paper. Factors or the physics influencing dynamic response and dynamic bifurcation for TVES with rheology are identified and are also demonstrated through model problem studies. A simple model problem consisting of a rod (1D) of TVES material with memory fixed at one end and subjected to harmonic excitation at the other end is considered to study nonlinear dynamics of TVES with rheology, frequency response as well as dynamic bifurcation phenomenon.

Finite Deformation, Finite Strain Nonlinear Dynamics and Dynamic Bifurcation in TVE Solids

This paper presents the mathematical model consisting of conservation and balance laws (CBL) of classical continuum mechanics (CCM) and the constitutive theories derived using entropy inequality and representation theorem for thermoviscoelastic solids (TVES) matter without memory. The CBL and the constitutive theories take into account finite deformation and finite strain deformation physics. This mathematical model is thermodynamically and mathematically consistent and is ideally suited to study nonlinear dynamics of TVES and dynamic bifurcation and is used in the work presented in this paper. The finite element formulations are constructed for obtaining the solution of the initial value problems (IVPs) described by the mathematical models. Both space-time coupled as well as space-time decoupled finite element methods are considered for obtaining solutions of the IVPs. Space-time coupled finite element formulations based on space-time residual functional (STRF) that yield space-time variationally consistent space-time integral forms are considered. This approach ensures unconditional stability of the computations during the entire evolution. In the space-time decoupled finite element method based on Galerkin method with weak form for spatial discretization, the solutions of nonlinear ODEs in time resulting from the decoupling of space and time are obtained using Newmark linear acceleration method. Newton’s linear method is used to obtain converged solution for the nonlinear system of algebraic equations at each time step in the Newmark method. The different aspects of the deformation physics leading to the factors that influence nonlinear dynamic response and dynamic bifurcation are established using the proposed mathematical model, the solution method and their validity is demonstrated through model problem studies presented in this paper. Energy methods and superposition techniques in any form including those used in obtaining solutions are neither advocated nor used in the present work as these are not supported by calculus of variations and mathematical classification of differential operators appearing in nonlinear dynamics. The primary focus of the paper is to address various aspects of the deformation physics in nonlinear dynamics and their influence on dynamic bifurcation phenomenon using mathematical models strictly based on CBL of CCM using reliable unconditionally stable space-time coupled solution methods, which ensure solution accuracy or errors in the calculated solution are always identified. Many model problem studies are presented to further substantiate the concepts presented and discussed in the paper. Investigations presented in this paper are also compared with published works when appropriate.

NONLINEAR WAVES AND PERIODIC SOLUTION IN FINITE DEFORMATION ELASTIC ROD

A nonlinear wave equation of elastic rod taking account of finite deformation, transverse inertia and shearing strain is derived by means of the Hamilton principle in this paper. Nonlinear wave equation and truncated nonlinear wave equation are solved by the Jacobi elliptic sine function expansion and the third kind of Jacobi elliptic function expansion method. The exact periodic solutions of these nonlinear equations are obtained, including the shock wave solution and the solitary wave solution. The necessary condition of exact periodic solutions, shock solution and solitary solution existence is discussed.

Three kinds of nonlinear dispersive waves in elastic rods with finite deformation

On the basis of classical linear theory on longitudinal, torsional and flexural waves in thin elastic rods, and taking finite deformation and dispersive effects into consideration, three kinds of nonlinear evolution equations are derived. Qualitative analysis of three kinds of nonlinear equations are presented. It is shown that these equations have homoclinic or heteroclinic orbits on the phase plane, corresponding to solitary wave or shock wave solutions, respectively. Based on the principle of homogeneous balance, these equations are solved with the Jacobi elliptic function expansion method. Results show that existence of solitary wave solution and shock wave solution is possible under certain conditions. These conclusions are consistent with qualitative analysis.

FINITE ELEMENT ANALYSIS OF THERMOELASTIC BEHAVIOR OF PIEZOELECTRIC STRUCTURES UNDER FINITE DEFORMATION

It is noted that the behavior of most piezoelectric materials is temperaturedependent and such piezo-thermo-elastic coupling phenomenon has become even more pronounced in thecase of finite deformation. On the other hand, for the purpose of precise shape and vibrationcontrol of piezoelectric smart structures, their deformation under external excitation must beideally modeled. This demands a thorough study of the coupled piezo-thermo-elastic response underfinite deformation. In this study, the governing equations of piezoelectric structures areformulated through the theory of virtual displacement principle and a finite element method isdeveloped. It should be emphasized that in the finite element method the fully coupledpiezo-thermo-elastic behavior and the geometric non-linearity are considered. The method developedis then applied to simulate the dynamic and steady response of a clamped plate to heat flux actingon one side of the plate to mimic the behavior of a battery plate of satellite irradiated under thesun. The results obtained are compared against classical solutions, whereby the thermal conductivityis assumed to be independent of deformation. It is found that the full-coupled theory predicts lesstransient response of the temperature compared to the classic analysis. In the steady state limit,the predicted temperature distribution within the plate for small heat flux is almost the same forboth analyses. However, it is noted that increasing the heat flux will increase the deviationbetween the predictions of the temperature distribution by the full coupled theory and by theclassic analysis. It is concluded from the present study that, in order to precisely predict thedeformation of smart structures, the piezo-thermo-elastic coupling, geometric non-linearity and thedeformation dependent thermal conductivity should be taken into account.

GEOMETRICAL NONLINEAR WAVES IN FINITEDEFORMATION ELASTIC RODS

By using Hamilton-type variation principle in non-conservation system, the nonlinear equation of wave motion of a elastic thin rod was derived according to Lagrange description of finite deformation theory. The dissipation caused due to viscous effect and the dispersion introduced by transverse inertia were taken into consideration so that steady traveling wave solution can be obtained. Using multi-scale method the nonlinear equation is reduced to a KdV-Burgers equation which corresponds with saddle-spiral heteroclinic orbit on phase plane. Its solution is called the oscillating-solitary wave or saddle-spiral shock wave. If viscous effect or transverse inertia is neglected, the equation is degraded to classical KdV or Burgers equation. The former implies a propagating solitary wave with homoclinic on phase plane, the latter means shock wave and heteroclinic orbit.

DIFFERENTIAL QUADRATURE METHOD FOR BENDING OF ORTHOTROPIC PLATES WITH FINITE DEFORMATION AND TRANSVERSE SHEAR EFFECTS

Based on the Reddy's theory of plates with the effect of higher-order shear deformations, the governing equations for bending of orthotropic plates with finite deformations were established. The differential quadrature (DQ) method of nonlinear analysis to the problem was presented. New DQ approach, presented by Wang and Bert (DQWB), is extended to handle the multiple boundary conditions of plates. The techniques were also further extended to simplify nonlinear computations. The numerical convergence and comparison of solutions were studied. The results show that the DQ method presented is very reliable and valid. Moreover, the influences of geometric and material parameters as well as the transverse shear deformations on nonlinear bending were investigated. Numerical results show the influence of the shear deformation on the static bending of orthotropic moderately thick plate is significant.

VARIATIONAL PRINCIPLES OF ASYMMETRIC ELASTICITY THEORY OF FINITE DEFORMATION

In this paper, based on the finite deformation S-R decomposition theorem, the definition of the body moment is renewed as the stem of its internal and external. The expression of the increment rate of the deformation energy is derived and the physical meaning is clarified. The power variational principle and the complementary power variational principle for finite deformation mechanics are supplemented and perfected.

New finite strain elastoplastic equations for accurately and explicitly simulating pseudoelastic-to-plastic transition effects of shape memory alloys

A new finite strain elatoplastic J2-flow model with coupling effects of both isotropic and anisotropic hardening is proposed with the co-rotational logarithmic rate.In terms of certain single-variable shape functions representing uniaxial loading and unloading curves,explicit multi-axial expressions for the three hardening quantities incorporated in the new model proposed are derived in unified forms for the purpose of automatically and accurately simulating complex pseudoelastic-to-plastic transition effects of shape memory alloys(SMAs)under multiple loading-unloading cycles.Numerical examples show that with only a single parameter of direct physical meaning for each cycle,accurate and explicit simulations may be achieved for extensive data from multiple cycle tests.

An efficient Galerkin meshfree formulation for shear deformable beam under finite deformation

This paper proposes a geometrically nonlinear total Lagrangian Galerkin meshfree formulation based on the stabilized conforming nodal integration for efficient analysis of shear deformable beam.The present nonlinear analysis encompasses the fully geometric nonlinearities due to large deflection,large deformation as well as finite rotation.The incremental equilibrium equation is obtained by the consistent linearization of the nonlinear variational equation.The Lagrangian meshfree shape function is utilized to discretize the variational equation.Subsequently to resolve the shear and membrane locking issues and accelerate the computation,the method of stabilized conforming nodal integration is systematically implemented through the Lagrangian gradient smoothing operation.Numerical results reveal that the present formulation is very effective.

Thermo-hydro-mechanical-air coupling finite element method and its application to multi-phase problems

In this paper, a finite element method （FEM）-based multi-phase problem based on a newly proposed thermal elastoplastic constitutive model for saturated/unsaturated geomaterial is discussed. A program of FEM named as SOFT, adopting unified field equations for thermo-hydro-mechanical-air （THMA） behavior of geomaterial and using finite element-finite difference （FE-FD） scheme for so/l-water-air three-phase coupling problem, is used in the numerical simulation. As an application of the newly proposed numerical method, two engineering problems, one for slope failure in unsaturated model ground and another for in situ heating test related to deep geological repository of high-level radioactive waste （HLRW）, are simulated. The model tests on slope failure in unsaturated Shirasu ground, carried out by Kitamura et al. （2007）, is simulated in the framework of soil-water-air three-phase coupling under the condition of constant temperature. While the in situ heating test reported by Munoz （2006） is simulated in the same framework under the conditions of variable temperature hut constant air pressure.

STUDY ON DYNAMIC CONSTITUTIVE RELATIONS FOR CONCRETE WITH FINITE DEFORMATION

The differences between finite deformation and infinitesimal deformation are discussed. They are exercised on elasto-viscoplastic constitutive relations of concrete. Then, a rate-dependent mechanics model was presented on the basis of Ottosen's four-parameter yield criterion, where different loading surface transferring laws were taken into account, when material was in hardening stage or in softening stage, respectively. The model is well established, so that it can be applied to simulate the response of concrete subject to impact loading. Green-Naghdi stress rate was introduced as objective stress rate. Appropriate hypothesis was postulated in accordance with many experimental results, which could reflect the mechanical behaviour of concrete with large deformation. Available thoughts as well as effective methods are also provided for the research on related engineering problems.

SOLITARY WAVES IN FINITE DEFORMATION ELASTIC CIRCULAR ROD

A new nonlinear wave equation of a finite deformation elastic circular rod simultaneously introducing transverse inertia and shearing strain was derived by means of Hamilton principle. The nonlinear equation includes two nonlinear terms caused by finite deformation and double geometric dispersion effects caused by transverse inertia and transverse shearing strain. Nonlinear wave equation and corresponding truncated nonlinear wave equation were solved by the hyperbolic secant function finite expansion method. The solitary wave solutions of these nonlinear equations were obtained. The necessary condition of these solutions existence was given also.

EFFECTIVE STRESS AND STRAIN IN FINITE DEFORMATION

Whether the concept of effective stress and strain in elastic-plastic theory is still valid under the condition of finite deformation was mainly discussed. The uni-axial compression experiments in plane stress and plane strain states were chosen for study. In the two kinds of stress states, the stress-strain curve described by logarithm strain and rotated Kirchhoff stress matches the experiments data better than the curves defined by other stress-strain description.

A new integration algorithm for finite deformation of thermo-elasto-viscoplastic single crystals

An algorithm for integrating the constitutive equations in thermal framework is presented, in which the plastic deformation gradient is chosen as the integration variable. Compared with the classic algorithm, a key feature of this new approach is that it can describe the finite deformation of crystals under thermal conditions. The obtained plastic deformation gradient contains not only plastic defor- mation but also thermal effects. The governing equation for the plastic deformation gradient is obtained based on ther- mal multiplicative decomposition of the total deformation gradient. An implicit method is used to integrate this evo- lution equation to ensure stability. Single crystal 1 100 aluminum is investigated to demonstrate practical applications of the model. The effects of anisotropic properties, time step, strain rate and temperature are calculated using this integration model.

A FINITE DEFORMATION ANALYSIS OF POLYCRYSTAL BY CONSIDERING THE INFLUENCE OF GRAIN BOUNDARY

By combining grain boundary (GB) and its influence zone, a micromechanic model for polycrystal is established for considering the influence of GB. By using the crystal plasticity theory and the finite element method for finite deformation, numerical simulation is carried out by the model. Calculated results display the microscopic characteristic of deformation fields of grains and are in qualitative agreement with experimental results.

INVERSE ASYMPTOTIC SOLUTION METHOD FOR FINITEDEFORMATION ELASTO-PLASTICITY

The development of modern mechanics in recent years has made many importantprogresse in the concepts and methods for nonlinear large deforntation mechanics([1],[2],[3]etc.). The presenl paper is aimed to show how the natural co-moving systemmethod and Stokes-Chen,s decomposition theorem can be effectively appliedasymptotically to solving problems of finite defomation elasto-plasticity by inverseasymptotic method for enyineering design purpose. Rigid punch problem is exanipfifiedin the paper.

ANALYSIS OF SIMPLE SHEAR ENDOCHRONIC EQUATIONSFOR FINITE PLASTIC DEFORMATION

Jaumann rate, generalized Jaumann rate,Fu rate and Wu rate were incorporated into endochronic equations for finite plastic deformation to analyze simple shear finite deformation. The results show that an oscillatory shear stress and normal stress response to a monotonically increasing shear strain occurs when Jaumann rate objective model is adopted for hypoelastic or endochronic materials. The oscillatory response is dependent on objective rate adopted,independent on elastoplastic models. Normal stress is unequal to zero during simple shear finite deformation.

Dynamic Plastic Response of Circular Plates Resting on Fluid with Finite Deformation

This paper concerns the dynamic plastic response of a circular plate resting on fluid subjected to a uniformly distributed rectangular load pulse with finite deformation. It is assumed that the fluid is incompressible and inviscous, and the plate is made of rigid-plastic material and simply supported along its edge. By using the method of the Hankel integral transformation, the nonuniform fluid resistance is derived as the plate and the fluid is coupled. Finally, an analytic solution for a circular plate under a medium load is obtained according to the equations of motion of the plate with finite deformation.