A bifurcation-based decohesion model for simulating the transition from localization to decohesion with the MPM

Based on recent experimental observations on the formation of localization before delamination, a bifurcation-based decohesion approach is proposed in this paper to simulate the transition from localization to decohesion involved in the delamination process of compressed films. The onset and orientation of discontinuous failure are identified from the discontinuous bifurcation analysis. A discrete constitutive model is then formulated based on the bifurcation analysis to predict the evolution of material failure as decohesion or separation of continuum. The Material Point Method, that does not employ fixed mesh-connectivity, is developed as a robust spatial discretization method to accommodate the multi-scale discontinuities involved in the film delamination. To demonstrate the potential of the proposed approach, a parametric study is conducted to explore the effects of aspect ratio and failure mode on the evolution of failure patterns under different boundary conditions, which provides a better understanding on the physics behind the film delamination process.Received: January 16, 2003; revised: August 11, 2003     

Closed-form solutions, extremality and nonsmoothness criteria in a large deformation elasticity problem

The pure azimuthal shear problem for a circular cylindrical tube of nonlinearly elastic material, both isotropic and anisotropic, is examined on the basis of a complementary energy principle. For particular choices of strain-energy function, one convex and one non-convex, closed-form solutions are obtained for this mixed boundary-value problem, for which the governing differential equation can be converted into an algebraic equation. The results for the non-convex strain energy function provide an illustration of a situation in which smooth analytic solutions of a nonlinear boundary-value problem are not global minimizers of the energy in the variational statement of the problem. Both the global minimizer and the local extrema are identified and the results are illustrated for particular values of the material parameters.      

Parabolic Type Decay Rates for 1-D-Thermoelastic Systems with Time-Dependent Coefficients

 The purpose of this paper is to derive L ^p − L ^q decay estimates for linear thermoelastic systems with time-dependent coefficient in one space variable. When all coefficients in the system have the same growth speed with small oscillations, we obtain a parabolic type decay estimate. For the system with time-dependent coefficients, we need to investigate the delicate asymptotic behaviour of characteristic roots and the remainder of diagonalization, which will be treated by dividing the phase space into three regions. Received September 15, 2001; in revised form April 20, 2002     

Weak solutions of the interior boundary value problems of plane Cosserat elasticity

In this paper we formulate interior Dirichlet and Neumann boundary value problems of plane Cosserat elasticity in Sobolev spaces, show that these problems are well-posed and find the corresponding weak solutions in the form of integral potentials. Received: April 7, 2005     

Thermomechanics of twinning from a parent to a twin with mirror image symmetry

In the paper we develop a modeling with multiple configurations and mirror image of parent crystal in the twinned structure, to describe the behavior of partially twinned structure. In the constitutive framework we take into account that: (1) the untwinned and twinned material have distinct natural configurations by virtue of their miscrostructure being different, (2) the material symmetry groups of the untwinned and twinned structures characterize the peculiar feature that the presence of the mirror image structure is related to the untwinned structure, but it can exist only as a counterpart of the previous one. The partially twinned structure is described by the evolution equations for the growth of twins, characterized by a pair of a deformation like tensorial variable and a scalar field with meaning of the volume fraction for the twins. The capability of the material to twin and untwin at a constant rate of strain in uniaxial compression has been analyzed and the oscillatory behavior predicted by the model reveals qualitative agreement with experimental evidences.      

A dual-parametric finite element method for cavitation in nonlinear elasticity

A dual-parametric finite element method is introduced in this paper for the computation of singular minimizers in the 2D cavitation problem in nonlinear elasticity. The method overcomes the difficulties, such as the mesh entanglement and material interpenetration, generally encountered in the finite element approximation of problems with extremely large expansionary deformation. Numerical experiments show that the method is highly efficient in the computation of cavitation problems. Numerical experiments are also conducted on discrete problems without the radial symmetry to show the validity of the method to more general settings and the potential of its application to the study of mechanism of cavity nucleation in nonlinear elastic materials.     

Survey of Modified Classical Bifurcations Due to Deformation-Sensitive Loading with Nonlinear,Nonsmooth and Nonconvex Effects

Some loading devices show certain deformation-sensitive but conservative characteristics leading to changes of the classical postbifurcation behavior of structures. To compare the effect of dead and configuration-dependent loading devices, classical bifurcation models are considered. Fully nonlinear behavior is considered, in sense of material, geometry and loading. Originally or approximately nonsmooth functions are included with nonconvexity in the material or loading behaviour. Global equilibrium paths are modified by deformation-sensitive loading devices, nonsmoothness and nonconvexity.     

Variable electrical and thermal conductivity in the theory of generalized thermoelastic diffusion

This paper deals with the problem of a thermoelastic half-space with a permeating substance in contact with the bounding plane in the context of the theory of generalized thermoelastic diffusion with one relaxation time and with variable electrical and thermal conductivity. The bounding surface of the half-space is taken to be traction free and subjected to a time dependent thermal shock. The solution is obtained in the Laplace transform domain by a direct approach. A numerical technique is employed to obtain the solution in the physical domain. It is found that there exist two coupled waves, one of which is elastic and the other is thermal, and a third wave affects diffusion mainly. A comparison is made with the results obtained in a thermoelastic medium with and without diffusion in the following cases : (a) the electrical and thermal conductivities have constant values, (b) the presence of magnetic field and (c) the generalized theory in thermoelasticity. Received: June 1, 2005     

Boundary conditions in intrinsic nonlinear elasticity

The intrinsic formulation of the displacement-traction problem of nonlinear elasticity is a system of partial differential equations and boundary conditions whose unknown is the Cauchy–Green strain tensor field instead of the deformation as is customary. We explicitly identify here the boundary conditions satisfied by the Cauchy–Green strain tensor field appearing in such intrinsic formulations.     

Concentrators of deformations and fracture of plastic bodies

An approach to the investigation of shape discontinuity regions as strain concentrators is proposed. The near-concentrator strain fields are determined on the basis of the theory of ideal rigid-plastic body; under the condition of plane deformation, their determination is reduced to integration of ordinary differential equations. The deformation as a function of the location of the plastic region and its shape evolution in the process of plastic flow is studied. The plastic flow is demonstrated to be not unique (within the framework of solution completeness). A deformation criterion for the choice of the preferred plastic flow is suggested. The fracture of a V-notched strip is considered. On the basis of the solutions obtained, an approach to the investigation of the fracture processes for more complicated models is formulated.     

Analysis of a dynamic unilateral contact problem for a cracked viscoelastic body

This work deals with the mathematical analysis of a dynamic unilateral contact problem with friction for a cracked viscoelastic body. We consider here a Kelvin-Voigt viscoelastic material and a nonlocal friction law. To prove the existence of a solution to the unilateral problem with friction, an auxiliary penalized problem is studied. Several estimates on the penalized solutions are given, which enable us to pass to the limit by using compactness results. Received: February 16, 2005     

On crack propagation shapes in elastic bodies

We consider a two dimensional elastic isotropic body with a curvilinear crack. The formula for the derivative of the energy functional with respect to the crack length is discussed. It is proved that this derivative is independent of the crack path provided that we consider quite smooth crack propagation shapes. An estimate for the derivative of the energy functional being uniform with respect to the crack propagation shape is derived.     

Variational formulation of a modified couple stress theory and its application to a simple shear problem

A variational formulation is provided for the modified couple stress theory of Yang et al. by using the principle of minimum total potential energy. This leads to the simultaneous determination of the equilibrium equations and the boundary conditions, thereby complementing the original work of Yang et al. where the boundary conditions were not derived. Also, the displacement form of the modified couple stress theory, which is desired for solving many problems, is obtained to supplement the existing stress-based formulation. All equations/expressions are presented in tensorial forms that are coordinate-invariant. As a direct application of the newly obtained displacement form of the theory, a simple shear problem is analytically solved. The solution contains a material length scale parameter and can capture the boundary layer effect, which differs from that based on classical elasticity. The numerical results reveal that the length scale parameter (related to material microstructures) can have a significant effect on material responses.      

A partitioned strongly coupled fluid-structure interaction method to model heart valve dynamics

A new method for the calculation of fluid-structure interaction (FSI) of highly flexible bodies is presented. This innovative algorithm demonstrates the strong coupling of a commercial computational fluid dynamics code with an in-house coded structural solver. The strong response of the pressure distribution to the displacement can be approximated by a reduced order model for the fluid solver. The Jacobian of this reduced order model is then used in the structural solver to obtain a stable and full implicit iteration scheme. The method is demonstrated on a 2D model of a flexible aortic valve during the cardiac cycle. Furthermore, the model is able to calculate shear stresses on the leaflet.     

Harmonic wave generation in non linear thermoelasticity by variational iteration method and Adomian's method

This paper applies the variational iteration method and Adomian's decomposition method to solve numerically the harmonic wave generation in a nonlinear, one-dimensional elastic half-space model subjected initially to a prescribed harmonic displacement. The results show that the variational iteration method is much easier, more convenient, and more stable and efficient than Adomian decomposition method.     

An evolutionary game based particle swarm optimization algorithm

Particle swarm optimization (PSO) is an evolutionary algorithm used extensively. This paper presented a new particle swarm optimizer based on evolutionary game (EGPSO). We map particles’ finding optimal solution in PSO algorithm to players’ pursuing maximum utility by choosing strategies in evolutionary games, using replicator dynamics to model the behavior of particles. And in order to overcome premature convergence a multi-start technique was introduced. Experimental results show that EGPSO can overcome premature convergence and has great performance of convergence property over traditional PSO.     

Covariance in linearized elasticity

In this paper we covariantly obtain all the governing equations of linearized elasticity. Our motivation is to see if one can make a connection between invariance (covariance) properties of the (global) balance of energy in nonlinear elasticity and those of its counterpart in linear elasticity. We start by proving a Green-Naghdi-Rivilin theorem for linearized elasticity. We do this by first linearizing energy balance about a given reference motion and then by postulating its invariance under isometries of the Euclidean ambient space. We also investigate the possibility of covariantly deriving a linearized elasticity theory, without any reference to the local governing equations, e.g. local balance of linear momentum. In particular, we study the consequences of linearizing covariant energy balance and covariance of linearized energy balance. We show that in both cases, covariance gives all the field equations of linearized elasticity.      

A remark on the junction in a thin multi-domain: the non convex case

Our aim consists of studying, in the spirit of Gamma convergence, a dimension reduction problem for a multi-domain filled of either an hyperelastic material or a non simple grade-two one. We derive asymptotically the limit energy density starting from a sample described trough non convex bulk energy densities, depending either on the first or second order derivative of the displacement.      

Mean field modeling of isotropic random Cauchy elasticity versus microstretch elasticity

We show that the averaged response of random isotropic Cauchy elastic material can be described analytically. It leads to a higher gradient model with explicit expressions for the dependence on the second derivatives of the mean field. A subsequent penalty formulation coincides with a linear elastic micro-stretch model with specific choice of constitutive parameters, depending only on the average cut-off length (the internal characteristic length scale L_c > 0). Thus the microstretch displacement field can be interpreted as an approximated mean field response for these parameter ranges. The mean field free energy in this micro-stretch formulation is not uniformly pointwise positive, nevertheless, the model is well posed.