RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CONTINUUM THEORIES (Ⅹ) —MASTER BALANCE LAW

Through a comparison between the expressions of master balance laws and the conservation laws derived by Noether's theorem, a unified master balance law and six physically possible balance equations for micropolar continuum mechanics are naturally deduced. Among them, by extending the well-known conventional concept of energymomentum tensor, the rather general conservation laws and balance equations named after energy-momentum, energy-angular momentum and energy-energy are obtained. It is clear that the forms of the physical field quantities in the master balance law for the last three cases could not be assumed directly by perceiving through the intuition. Finally, some existing results are reduced immediately as special cases.     

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CONTINUUM THEORIES (Ⅹ)--MASTER BALANCE LAW

Through a comparison between the expressions of master balance laws and the conservation laws derived by Noether's theorem, a unified master balance law and six physically possible balance equations for micropolar continuum mechanics are naturally deduced. Among them, by extending the well-known conventional concept of energymomentum tensor, the rather general conservation laws and balance equations named after energy-momentum, energy-angular momentum and energy-energy are obtained. It is clear that the forms of the physical field quantities in the master balance law for the last three cases could not be assumed directly by perceiving through the intuition. Finally,some existing results are reduced immediately as special cases.     

Rheology of polymer blends with matrix-phase viscoelasticity and a narrow droplet size distribution

A Hamiltonian framework of non-equilibrium thermodynamics is adopted to construct a set of dynamical continuum equations for a polymer blend with matrix viscoelasticity and a narrow droplet size distribution that is assumed to obey a Weibull distribution function. The microstructure of the matrix is described in terms of a conformation tensor. The variable droplet distribution is described in terms of two thermodynamic variables: the droplet shape tensor and the number density of representative droplets. A Hamiltonian functional in terms of the thermodynamic variables is introduced and a set of time evolution equations for the system variables is derived. Sample calculations for homogenous flows and constant droplet distribution are compared with data of a PIB/PDMS blend and a HPC/PDMS blend with high viscoelastic contrast. For the PIB/PDMS blend, satisfactory predictions of the flow curves are obtained. Sample calculations for a blend with variable droplet distribution are performed and the effect of flow on the rheology, droplet morphology, and on the droplet distribution are discussed. It is found that deformation can increase or decrease the dispersity of the droplet morphology for the flows investigated herein.     

Second-harmonic resonance with Faraday excitation

Capillary-gravity waves in an inviscid liquid exhibit second- or sub-harmonic resonance at precise frequencies. When the container performs small periodic vertical vibrations, either wave may also experience Faraday (‘parametric’) excitation. Equations describing this situation are derived, incorporating slight detuning from two-wave and Faraday resonances. Similar equations arise in other physical contexts.With Faraday forcing of the wave with lower frequency, the evolution equations (without detuning) are transformable to the corresponding unforced equations, the general solution of which is known. With Faraday forcing of the wave with higher frequency no such simplification is possible. Here, various transformed equations are considered and numerical results elucidate their solutions. For some initial data, solutions remain bounded; but other initial values give unbounded solutions. We establish the form of the boundaries that separate these two classes.     

Constitutive equations of a tensorial model for ductile damage of metals

Based on a micromechanical concept of void growth and change in void shape, a dissipation potential and constitutive equations for ductile damage of metals are presented. Multiplicative decomposition of the metric transformation tensor and thermodynamic formulation of the constitutive equations lead to a symmetric second-order tensor of damage which is physically meaningful. Its first invariant defines the damage related to plastic dilatation of the material due to the void growth. The second invariant of the deviatoric tensor accounts for the damage associated with a change in the void shape. Two physically motivated normalized measures allow us to represent the kinetic process of strain-induced damage by using the equivalent parameter of damage including the limit conditions for the onset of void coalescence and ductile failure. An experimental analysis of the evolution of ductile damage is presented for the case of uniaxial tension of sheet steel specimens with artificial defects.     

RESTUDY OF THEORIES FOR ELASTIC SOLIDS WITH MICROSTRUCTURE

ＩｎｔｒｏｄｕｃｔｉｏｎＵｐｔｏｎｏｗｔｈｅｒｅｈａｓｂｅｅｎｖｅｒｙｍｕｃｈｗｒｉｔｔｅｎｗｏｒｋｏｎｔｈｅｓｕｂｊｅｃｔｓｏｆｃｏｎｔｉｎｕｕｍｔｈｅｏｒｉｅｓｉｎｗｈｉｃｈｔｈｅｄｅｆｏｒｍａｔｉｏｎｉｓｄｅｓｃｒｉｂｅｄｎｏｔｏｎｌｙｂｙｔｈｅｕｓｕａｌｖｅｃｔｏｒｄｉｓｐｌａｃｅｍｅｎｔｆｉｅｌｄ ,ｂｕｔｂｙｏｔｈｅｒｖｅｃｔｏｒｏｒｔｅｎｓｏｒｆｉｅｌｄｓａｓｗｅｌｌ.Ｉｎａｆａｍｏｕｓｍｏｎｏｇｒａｐｈ ,Ｅ .ＣｏｓｓｅｒａｔａｎｄＦ .Ｃｏｓｓｅｒａｔ[1]ｇａｖｅａｓｙｓｔｅｍａｔｉｃ…     

On the mechanics of solids with a growing mass

A general constitutive theory of the stress-modulated growth of biomaterials is presented with a particular accent given to pseudo-elastic soft living tissues. The governing equations of the mechanics of solids with a growing mass are revisited within the framework of finite deformation continuum thermodynamics. The multiplicative decomposition of the deformation gradient into its elastic and growth parts is employed to study the growth of isotropic, transversely isotropic, and orthotropic biomaterials. An explicit representation of the growth part of the deformation gradient is given in each case, which leads to an effective incremental formulation in the analysis of the stress-modulated growth process. The rectangular components of the instantaneous elastic moduli tensor are derived corresponding to selected forms of the elastic strain energy function. Physically appealing structures of the stress-dependent evolution equations for the growth induced stretch ratios are proposed.     

Integrity bases for a symmetric tensor and a vector — The crystal classes

Integrity bases are derived for a symmetric second-order tensor and a vector for the transformation groups corresponding to each of the crystal classes. In each case the irreducibility of the integrity basis is proven.     

Two Theories of Complex Bodies, one Lagrangean, the other Kinetic, Have a Common Ground

I have found in previous works that most special models proposed to represent bodies with some type of microstructure can be classified easily under the general umbrella of a theory where each element of the continuum is thought of as a Lagrangian system. To study phenomena in ‘kinetic’ continua I proposed an apparently different approach; the outcome is again a set of evolution equations. They mimic equations familiar in continua with affine microstructure: a Cauchy’s equation and an equation of balance of tensor moment of momentum, with the addition, however, of an equation of balance for a ‘Reynolds’ tensor’, an equation which, in a sense, shifts the boundary between kinetic and thermal properties of matter. I will show that there is no contrast between the two approaches. The latter one is based on an adequate and appropriately justified expression of the kinetic energy of the continuum, comprising the trace of the quoted Reynolds’ tensor and thus importing into the mechanical energy a term usually accounted by additional heat.     

Feedback minimization of first-passage failure of quasi integrable Hamiltonian systems

A nonlinear stochastic optimal control strategy for minimizing the first-passage failure of quasi integrable Hamiltonian systems (multi-degree-of-freedom integrable Hamiltonian systems subject to light dampings and weakly random excitations) is proposed. The equations of motion for a controlled quasi integrable Hamiltonian system are reduced to a set of averaged Itô stochastic differential equations by using the stochastic averaging method. Then, the dynamical programming equations and their associated boundary and final time conditions for the control problems of maximization of reliability and mean first-passage time are formulated. The optimal control law is derived from the dynamical programming equations and the control constraints. The final dynamical programming equations for these control problems are determined and their relationships to the backward Kolmogorov equation governing the conditional reliability function and the Pontryagin equation governing the mean first-passage time are separately established. The conditional reliability function and the mean first-passage time of the controlled system are obtained by solving the final dynamical programming equations or their equivalent Kolmogorov and Pontryagin equations. An example is presented to illustrate the application and effectiveness of the proposed control strategy.     

General invariant representations of the constitutive equations for isotropic nonlinearly elastic materials

This paper develops general invariant representations of the constitutive equations for isotropic nonlinearly elastic materials. Different sets of mutually orthogonal unit tensor bases are constructed from the strain argument tensor by using the representation theorem and corresponding irreducible invariants are defined. Their relations and geometrical interpretations are established in three dimensional principal space. It is shown that the constitutive law linking the stress and strain tensors is revealed to be a simple relationship between two vectors in the principal space. Relative to two different sets of the basis tensors, the constitutive equations are transformed according to the transformation rule of vectors. When a potential function is assumed to exist, the vector associated with the stress tensor is expressed in terms of its gradient with respect to the vector associated with the strain tensor. The Hill’s stability condition is shown to be that the scalar product of the increment of those two vectors must be positive. When potential function exists, it becomes to be that the 3 × 3 constitutive matrix derived from its second order derivative with respect to the vector associated with the strain must be positive definite. By decomposing the second order symmetric tensor space into the direct sum of a coaxial tensor subspace and another one orthogonal to it, the closed form representations for the fourth order tangent operator and its inversion are derived in an extremely simple way.     

On the final boundary value problems in linear thermoelasticity

In the present study we derive some uniqueness criteria for solutions of the Cauchy problem for the standard equations of dynamical linear thermoelasticity backward in time. We use Lagrange-Brun identities combined with some differential inequalities in order to show that the final boundary value problem associated with the linear thermoelasticity backward in time has at most one solution in appropriate classes of displacement-temperature fields. The uniqueness results are obtained under the assumptions that the density mass and the specific heat are strictly positive and the conductivity tensor is positive definite.     

Optimal bounded control for minimizing the response of quasi-integrable Hamiltonian systems

A procedure for designing optimal bounded control to minimize the response of quasi-integrable Hamiltonian systems is proposed based on the stochastic averaging method for quasi-integrable Hamiltonian systems and the stochastic dynamical programming principle. The equations of motion of a controlled quasi-integrable Hamiltonian system are first reduced to a set of partially completed averaged Itô stochastic differential equations by using the stochastic averaging method for quasi-integrable Hamiltonian systems. Then, the dynamical programming equation for the control problems of minimizing the response of the averaged system is formulated based on the dynamical programming principle. The optimal control law is derived from the dynamical programming equation and control constraints without solving the dynamical programming equation. The response of optimally controlled systems is predicted through solving the Fokker-Planck-Kolmogrov equation associated with fully completed averaged Itô equations. Finally, two examples are worked out in detail to illustrate the application and effectiveness of the proposed control strategy.     

On kinematic,thermodynamic, and kinetic coupling of a damage theory for polycrystalline material

The paper proposes a new consistent formulation of polycrystalline finite-strain elasto-plasticity coupling kinematics and thermodynamics with damage using an extended multiplicative decomposition of the deformation gradient that accounts for temperature effects. The macroscopic deformation gradient comprises four terms: thermal deformation associated with the thermal expansion, the deviatoric plastic deformation attributed to the history of dislocation glide/movement, the volumetric deformation gradient associated with dissipative volume change of the material, and the elastic or recoverable deformation associated with the lattice rotation/stretch. Such a macroscopic decomposition of the deformation gradient is physically motivated by the mechanisms underlying lattice deformation, plastic flow, and evolution of damage in polycrystalline materials. It is shown that prescribing plasticity and damage evolution equations in their physical intermediate configurations leads to physically justified evolution equations in the current configuration. In the past, these equations have been modified in order to represent experimentally observed behavior with regard to damage evolution, whereas in this paper, these modifications appear naturally through mappings by the multiplicative decomposition of the deformation gradient. The prescribed kinematics captures precisely the damage deformation (of any rank) and does not require introducing a fictitious undamaged configuration or mechanically equivalent of the real damaged configuration as used in the past.     

Mathematical model of an incompressible viscoelastic Maxwell medium

Nonstationary motions of incompressible viscoelastic Maxwell continuum with a constant relaxation time are considered. Because in an incompressible continuous medium, pressure is not a thermodynamic variable but coincides with the stress-tensor trace to within a factor, it follows that, separating the spherical part from this tensor, one can assume that the remaining part of the stress tensor has zero trace. In the case of an incompressible medium, the equations for the velocity, pressure, and stress tensor form a closed system of first-order equations which has both real and complex characteristics, which complicates the formulation of the initial-boundary-value problem. Nevertheless, the resolvability of the Cauchy problem can be proved in the class of analytic functions. Unique resolvability of the linearized problem was established in the classes of functions of finite smoothness. The class of effectively one-dimensional motions for which the subsystem of three equations is a hyperbolic one was studied. The results of an asymptotic analysis of the latter imply the possible formation of discontinuities during the evolution of the solution. The general system of equations of motion admits an infinite-dimensional Lie pseudo-group which contains an extended Galilean group. The theorem of the invariance of the conditions on the a priori unknown free boundary was proved to obtain exact solutions of free-boundary problems. The problem of deformation of a viscoelastic strip subjected to tangential stresses applied to the free boundary is considered as an example of application of this theorem. In this problem, a scale effect of short-wave instability caused by the absence of diagonal dominance of the stress tensor deviator was found.     

Optimal control strategies for stochastically excited quasi partially integrable Hamiltonian systems

In this paper two different control strategies designed to alleviate the response of quasi partially integrable Hamiltonian systems subjected to stochastic excitation are proposed. First, by using the stochastic averaging method for quasi partially integrable Hamiltonian systems, an n-DOF controlled quasi partially integrable Hamiltonian system with stochastic excitation is converted into a set of partially averaged Itô stochastic differential equations. Then, the dynamical programming equation associated with the partially averaged Itô equations is formulated by applying the stochastic dynamical programming principle. In the first control strategy, the optimal control law is derived from the dynamical programming equation and the control constraints without solving the dynamical programming equation. In the second control strategy, the optimal control law is obtained by solving the dynamical programming equation. Finally, both the responses of controlled and uncontrolled systems are predicted through solving the Fokker-Plank-Kolmogorov equation associated with fully averaged Itô equations. An example is worked out to illustrate the application and effectiveness of the two proposed control strategies.     

A novel internal dissipation inequality by isotropy and its implication for inelastic constitutive characterization

In the present work a novel inelastic deformation caused internal dissipation inequality by isotropy is revealed. This inequality has the most concise form among a variety of internal dissipation inequalities, including the one widely used in constitutive characterization of isotropic finite strain elastoplasticity and viscoelasticiy. Further, the evolution term describing the difference between the rate of deformation tensor and the “principal rate” of the elastic logarithmic strain tensor is set, according to the standard practice by isotropy, to equal a rank-two isotropic tensor function of the corresponding branch stress, with the tensor function having an eigenspace identical to the eigenspace of the branch stress tensor. Through that a general form of evolution equation for the elastic logarithmic strain is formulated and some interesting and important results are derived. Namely, by isotropy the evolution of the elastic logarithmic strain tensor is embodied separately by the evolutions of its eigenvalues and eigenprojections, with the evolution of the eigenprojections driven by the rate of deformation tensor and the evolution of the eigenvalues connected to specific material behavior. It can be proved that by isotropy the evolution term in the present dissipation inequality stands for the essential form of the evolution term in the extensively applied dissipation inequality.     

Recursive-operator method in vibration problems for rod systems

Using linear differential equations with constant coefficients describing one-dimensional dynamical processes as an example, we show that the solutions of these equations and systems are related to the solution of the corresponding numerical recursion relations and one does not have to compute the roots of the corresponding characteristic equations. The arbitrary functions occurring in the general solution of the homogeneous equations are determined by the initial and boundary conditions or are chosen from various classes of analytic functions. The solutions of the inhomogeneous equations are constructed in the form of integro-differential series acting on the right-hand side of the equation, and the coefficients of the series are determined from the same recursion relations. The convergence of formal solutions as series of a more general recursive-operator construction was proved in [1]. In the special case where the solutions of the equation can be represented in separated variables, the power series can be effectively summed, i.e., expressed in terms of elementary functions, and coincide with the known solutions. In this case, to determine the natural vibration frequencies, one obtains algebraic rather than transcendental equations, which permits exactly determining the imaginary and complex roots of these equations without using the graphic method [2, pp. 448–449]. The correctness of the obtained formulas (differentiation formulas, explicit expressions for the series coefficients, etc.) can be verified directly by appropriate substitutions; therefore, we do not prove them here.     

Optimal Bounded Control of First-Passage Failure of Quasi-Integrable Hamiltonian Systems with Wide-Band Random Excitation

The optimal bounded control of quasi-integrable Hamiltonian systems with wide-band random excitation for minimizing their first-passage failure is investigated. First, a stochastic averaging method for multi-degrees-of-freedom (MDOF) strongly nonlinear quasi-integrable Hamiltonian systems with wide-band stationary random excitations using generalized harmonic functions is proposed. Then, the dynamical programming equations and their associated boundary and final time conditions for the control problems of maximizinig reliability and maximizing mean first-passage time are formulated based on the averaged Itô equations by applying the dynamical programming principle. The optimal control law is derived from the dynamical programming equations and control constraints. The relationship between the dynamical programming equations and the backward Kolmogorov equation for the conditional reliability function and the Pontryagin equation for the conditional mean first-passage time of optimally controlled system is discussed. Finally, the conditional reliability function, the conditional probability density and mean of first-passage time of an optimally controlled system are obtained by solving the backward Kolmogorov equation and Pontryagin equation. The application of the proposed procedure and effectiveness of control strategy are illustrated with an example.