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.

Numerical Analysis of Finite Deformation of Overbroken Rock Mass in Gob Area Based on Euler Model of Control Volume

The overbroken rock mass of gob areas is made up of broken and accumulated rock blocks compressed to some extent by the overlying strata. The beating pressure of the gob can directly affect the safety of mining fields, formarion of road retained along the next goaf and seepage of water and methane through the gob. In this paper, the software RFPA＇2000 is used to construct numerical models. Especially the Euler method of control volume is proposed to solve the simulation difficulty arising from plastically finite deformations. The results show that three characteristic regions occurred in the gob area： （1） a naturally accumulated region, 0-10 m away from unbroken surrounding rock walls, where the beating pressure is nearly zero; （2） an overcompacted region, 10-20 m away from unbroken walls, where the beating pressure results in the maximum value of the gob area; （3） a stable compaction region, more than 20 m away from unbroken walls and occupying absolutely most of the gob area, where the beating pressures show basically no differences. Such a characteristic can exolain the easy-seeoaged “O”-ring phenomena around mining fields very well.

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.

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.

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.

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.

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.

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.

FINITE DEFORMATION ELASTO-PLASTIC THEORY AND CONSISTENT ALGORITHM

By using the logarithmic strain, the finite deformation plastic theory, corresponding to the infinitesimal plastic theory, is established successively. The plastic consistent algorithm with first order accuracy for the finite element method (FEM) is developed. Numerical examples are presented to illustrate the validity of the theory and effectiveness of the algorithm.

EXACT SOLUTION FOR FINITE DEFORMATION PROBLEMS OF CANTILEVER BEAM WITH VARIABLE SECTION UNDER THE ACTION OF ARBITRARY TRANSVERSE LOADS

This paper deals with finite deformation problems of cantilever beam with variable sec- tion under the action of arbitrary transverse loads.By the use of a method of variable replacement, the nonlinear differential equation with varied coefficient for the problem can be transformed into an equation with variable separable.The exact solution can be obtained by the integration method. Some examples are given in the paper,and the results of these examples show that this exact solution includes the existing solutions in references as special cases.

THE RANDOM VARIATIONAL PRINCIPLE IN FINITE DEFORMATION OF ELASTICITY AND FINITE ELEMENT METHOD

In the present paper, we hare mtroduced the random materials. loads. geometricalshapes, force and displacement boundary condition directly. into the functionalvariational formula, by. use of a small parameter perturbation method, a unifiedrandom variational principle in finite defomation of elastieity and nonlinear randomfinite element method are esiablished, and used.for reliability, analysis of structures.Numerical examples showed that the methods have the advontages of simple andconvenient program implementation and are effective for the probabilistic problems inmechanics.

A COMPUTATIONAL METHOD ON PROCESS OPTIMAL CONTROL FOR FINITE DEFORMATION PROBLEMS

On the principle of non-incremental algorithm, some basic ideas of process optimal control iterative algorithm, based on the Optimal Control Variational Principle in Mechanics, is proposed in this paper. Then the essential governing equations are presented. This work provides a new method to achieve the numerical solutions of the mechanic of finite deformation.