Identification of the viscoelastic parameters of a polymer model by the aid of a MCMC method

A procedure to identify the viscoelastic material parameters of a solid amorphous polymer and to estimate their values is presented. Stress–strain material data is obtained for the polymer by a compression experiment. The material behavior of the polymer is modeled according to the generalized Maxwell model, which is fitted to the experimental data by the method of least squares to obtain a first approximation for the model parameters. The identification of the model parameters is completed by a Markov chain Monte Carlo (MCMC) method, which generates the probability distributions of the relevant parameters of the material. The utilized MCMC method enables us to determine a suitable complexity (i.e., the number of Maxwell elements) for the generalized Maxwell model, so that the model best fits the data and, simultaneously, leads to an identifiable set of parameters. The numerical results imply that the uniqueness of the solution is lost when the number of model parameters becomes redundant.     

Thermomechanical Behavior of Electrically Conducting Solids Exposed to an External Electromagnetic Field

This paper describes a procedure for the mathematical and numerical modeling of the thermomechanical behavior of electrically conducting solid bodies exposed to an external electromagnetic field. The constitutive equations for the electromagnetic field are the Maxwell equations written for the region of the solid body and the ambient medium. The stress-strain state of the solid is described using the relations for nonisothermal elastoplastic flow. The effects of the electromagnetic field on the heat-transfer and deformation processes are taken into account via heat release and ponderomotive forces, respectively. The relations between the electric and magnetic inductions and the corresponding field strengths are considered nonlinear. All physicomechanical parameters of the body material are temperature dependent.__________Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 5, pp. 14–26, September–October, 2005.     

Phase equilibrium in non-fluids and non-fluid mixtures

A criterion of phase equilibrium for mixtures of materials with arbitrary symmetry (e.g. between solid and fluid or two solid mixture phases) is deduced using a rational thermodynamics approach. This criterion, known also as Maxwell relation, is expressed via the difference of chemical potential tensors (Eshelby tensors) on the singular surface dividing the bulk phases.The thermomechanical balance equations, the entropy inequality and the Maxwell relation for phase equilibrium are given first for the case of pure (one-constituent) materials of arbitrary symmetry and then for the case of mixtures (including chemically reacting ones) of arbitrary symmetry.In the special case of fluids it is shown that the chemical potential tensors reduce to the classical scalar chemical potentials and the Maxwell relations to the classical thermodynamic criterion for the phase equilibrium.     

Consequences of the Maxwell relation for anti-plane shear deformations of an elastic solid

Previous numerical work has predicted chaotic distributions of phases in certain deformations of an elastic solid whose stress-strain relation in simple shear is nonmonotone. The present work provides an interpretation of these results making use of elastic stability considerations. For a specific material, a conformal mapping technique is shown to generate the totality of all deformations satisfying the Maxwell relation. It is shown that some boundary value problems have no solutions which obey the Maxwell relation. Such problems may have associated with them an infinite sequence of deformations (corresponding to a minimizing sequence in variational calculus) which satisfies the Maxwell relation in the sense of a limit. The properties of such sequences are examined in detail, and it is shown that the chaotic numerical solutions may be interpreted as terms in this sequence.     

On a principle of virtual work for thermo-elastic bodies

A principle of virtual work is proposed for thermo-elastic bodies. From it are derived the equations of motion, the Cauchy stress principle and the Gibbs relations. The principle is also used to analyse the response of internally constrained bodies.     

Effect of Maxwell stress on a moving crack with polarization saturation region in ferroelectric solid

The focus of this work is on a generalized two-dimensional problem of a crack moving in a piezoelectric solid subjected to uniform electrical load at infinity. The novel point includes that the electric field inside the crack is taken into account when polarization saturation region exists. Based on the extended Stroh formalism and complex function method, explicit expressions of both the stress fields in the solid and electric fields inside the crack are derived by using semi-permeable crack model, respectively. Effect of Maxwell stress along the crack surface is investigated and the results are illustrated graphically. It is shown that the moving speed of the crack cannot exceed the lowest bulk wave speed. It is also found that the medium properties inside the crack and surrounding the ferroelectric solid at infinity directly affect the Maxwell stress, and as a result the Maxwell stresses are remarkable and cannot be ignored under different electric load.     

Internal energy variational principles and governing equations in electroelastic analysis

The first thermodynamic law contains a universal thermodynamic variational principle. The complete internal energy variational principle in the electroelastic analysis is not discussed in previous papers. In this paper this principle will be discussed. From this principle the simple complete governing equations can be deduced, and the Maxwell stress can be naturally derived from this variational principle. It is shown that the Maxwell stress has slightly different forms determined by using internal energy or electric Gibbs free energy variational principle, but substantially they are the same. In the second-order precision the Maxwell stress is uniquely determined, and its expression has the same form for all deformable and rigid dielectrics. The electroelastic analyses in the dielectric should be studied together with its environment, because the electric field exists in all materials except the ideal conductor. The complete governing equations under finite deformation in the initial configuration are also discussed.     

The extended Padé-Laplace method for efficient discretization of linear viscoelastic spectra

The discretization of linear viscoelastic spectra is valuable as a starting point for non-linear viscoelastic modeling. However, obtaining the parameters of the generalized Maxwell model from linear viscoelastic experiments with naive least squares procedures is known to be an ill-posed problem. A novel technique, the Padé-Laplace method was recently elucidated (Fulchiron et al., 1993) for robustly extracting the parameters of the generalized Maxwell model from stress relaxation experiments, without any a priori assumption about the number of Maxwellian modes. We extend this method for obtaining the Maxwellian modes from dynamic data and discuss the relationship between continuous viscoelastic spectra and the Maxwellian modes obtained by this procedure. Furthermore, the applicability of this method with experimental data in limited time/frequency windows is clarified. Finally, a procedure for assembling the discretized spectrum with the Padé-Laplace method applied to both stress relaxation and dynamic data with typical experimental time/frequency cutoffs is developed.Dedicated to Professor H. Janeschitz-Kriegl on the occasion of his 70th birthday.     

Propagation of electromagnetic wave in the three phases soil media

IntroductionItisofparamountimportancetogetgeotechnicalparameterssuchasmoisture ,compactness,andshearstrengthforgeotechnicalengineeringtransaction .Butthemethodtofetchthesoilsamplesforlaboratorytestcannotavoidbringingthedisturbance ,andmakethetestresultsinconsistenttotheactualstatusofsoil.Inrecentyears,theapplicationofEMwaveinthedetectinghiddentroublesindamsandundergroundfoundationsdevelopedrapidly .Buttherearestillsomefactorsthatrestricttheadvanceofapplication .Firstly ,thesoilconstitutesthr…     

Gravitational motion of a two-particle cluster between two parallel plane solid wallsSédimentation de deux particules solides entre deux plans parallèles

The slow viscous settling migration of a 2-particule cluster between two solid and parallel plane walls is investigated by resorting to a Boundary Element Method. The procedure, valid for arbitrarily-shaped bodies, is presented and preliminary numerical results for both identical spheres and a spheroid-sphere cluster are discussed. To cite this article: L. Pasol, A. Sellier, C. R. Mecanique 334 (2006).     

Viscoelastic solutions for conical indentation

The aim of indentation analysis is to link indentation data, typically an indentation force vs. indentation depth curve, P–h, to meaningful mechanical properties of the indented material. While well established for time independent behavior, the presence of a time dependent behavior can strongly affect both the loading and the unloading responses. The paper presents a framework of viscoelastic indentation analysis based on the method of functional equations, developed by Lee and Radok [1960, The contact problem for viscoelastic bodies, J. Appl. Mech. 27, 438–444]. While the method is restricted to monotonically increasing contact areas, we show that it remains valid at the very beginning of the unloading phase as well. Based on this result, it is possible to derive closed form solutions following the classical procedure of functional formulations of viscoelasticity: (1) the identification of the indentation creep function, which is the indentation response to a Heaviside load; and (2) a convolution integral of the load history over the indentation creep function. This is shown here for a trapezoidal loading by a conical indenter on three linear isotropic viscoelastic materials with deviator creep: the 3-parameter Maxwell model, the 4-parameter Kelvin–Voigt model and the 5-parameter combined Kelvin–Voigt–Maxwell model. For these models, we derive closed form solutions that can be employed for the back-analysis of indentation results from the loading and holding period and for the definition of unloading time criteria that ensure that viscous effects are negligible in the unloading response.     

A novel time dependent cohesive zone model for the debonding interface between solid propellant and insulation

A novel time dependent cohesive zone model (CZM) is proposed in this paper based on two main assumptions. Firstly, ultimate cohesive parameters are inherent and fixed for a given non-aging bond interface. The apparent cohesive parameters are time related variables. Secondly, relaxation response of the interface is the main reason for the time dependent traction. Numerical simulation shows that the traction, critical displacement as well as damage initiation displacement will increase with imposed loading rate and parameter λ for single Maxwell box based model. N single Maxwell box connected in parallel construct the N Maxwell box based model, and each Maxwell box bears 1/N traction of the interface. Double cantilever beam (DCB) is utilized to investigate the structure response with the single Maxwell box based model including constant stretch and relaxation test. Quite good agreement between the numerical and experimental reaction force–displacement curves is obtained from stretch test of Double cantilever sandwich beam (DCSB) specimen with four different N Maxwell box based model, especially when the number of the Maxwell box is 7. It is a fact that the model will be more adaptive with more Maxwell box connected in parallel which can be revealed by the verification test.     

Numerical study around the corotational Maxwell model for the viscoelastic fluid flows

We present numerical results for the FEM (finite element method) presented in [Comput. Methods Appl. Mech. Engrg. 191 (2002) 5045–5065]. This method is devoted to the approximation of fluid flows obeying the Oldroyd model. A particularity of this method, is to take into account the purely viscoelastic case, the so-called Maxwell model, important in practice. Numerical results are given for a fluid flowing in an abrupt plane 4 to 1 contraction. We use the corotational Maxwell model as benchmark in the choice of our computations. Results are also given for the upper convected Maxwell model. Interesting effects appear on the velocity profile: a phenomenon of quasi slip at the downstream wall.     

Rheological representation of fractional order viscoelastic material models

We develop rheological representations, i.e., discrete spectrum models, for the fractional derivative viscoelastic element (fractional dashpot or springpot). Our representations are generalized Maxwell models or series of Kelvin-Voigt units, which, however, maintain the number of parameters of the corresponding fractional order model. Accordingly, the number of parameters of the rheological representation is independent of the number of rheological units. We prove that the representations converge to the corresponding fractional model in the limit as the number of units tends to infinity. The representations extend to compound fractional derivative models such as the fractional Maxwell model, fractional Kelvin-Voigt model, and fractional standard linear solid. Computational experiments show that the rheological representations are accurate approximations of the fractional order models even for a small number of units.     

Stabilization of the Rotation Axis of a Solid by Coupled Perfectly Rigid Bodies

The problem of stabilizing the axis of a solid by coupled perfectly rigid bodies (PRBs) is solved. The solid executes a plane-parallel motion. The PRBs can rotate as a single rigid body about the centroidal axis of the solid and counterrotate about its transverse axes through equal angles. There is a particle inside the solid which causes its imbalance. It is established that the principal state (if any) of the system—rotation about the centroidal axis—is stable, whereas the rest (unwanted) states are unstable __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 8, pp. 122–129, August 2005.     

The thermostressed state of electrically conductive bodies subjected to an external electromagnetic field

A technique is proposed for mathematical and numerical modeling of thermomechanical processes in electrically conductive bodies subjected to an external electromagnetic field. The initial relations for the determination of the electromagnetic field are the Maxwell equations. The stress and strain states of the body are described using the equations of nonisothermal elastoplasticity. The model takes into account the coupling of the electromagnetic and thermal fields. All physical and mechanical parameters of the material depend on temperature. The process of high-temperature induction treatment of a ferromagnetic cylinder is considered as an example __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 12, pp. 13–25, December 2005.     

Linear gravity waves on Maxwell fluids of finite depth

Linear surface gravity waves on Maxwell viscoelastic fluids with finite depth are studied in this paper. A dispersion equation describing the spatial decay of the gravity wave in finite depth is derived. A dimensionless memory (time) number 0 is introduced. The dispersion equation for the pure viscous fluid will be a specific case of the dispersion equation for the viscoelastic fluid as θ=0. The complex dispersion equation is numerically solved to investigate the dispersion relation. The influences of θ and water depth on the dispersion characteristics and wave decay are discussed. It is found that the role of elasticity for the Maxwell fluid is to make the surface gravity wave on the Maxwell fluid behave more like the surface gravity wave on the inviscid fluid.     

A modified SPH method for simulating motion of rigid bodies in Newtonian fluid flows

A weakly compressible smoothed particle hydrodynamics (WCSPH) method is used along with a new no-slip boundary condition to simulate movement of rigid bodies in incompressible Newtonian fluid flows. It is shown that the new boundary treatment method helps to efficiently calculate the hydrodynamic interaction forces acting on moving bodies. To compensate the effect of truncated compact support near solid boundaries, the method needs specific consistent renormalized schemes for the first and second-order spatial derivatives. In order to resolve the problem of spurious pressure oscillations in the WCSPH method, a modification to the continuity equation is used which improves the stability of the numerical method. The performance of the proposed method is assessed by solving a number of two-dimensional low-Reynolds fluid flow problems containing circular solid bodies. Wherever possible, the results are compared with the available numerical data.     

A simple and efficient BEM implementation of quasistatic linear visco-elasticity

A simple yet efficient procedure to solve quasistatic problems of special linear visco-elastic solids at small strains with equal rheological response in all tensorial components, utilizing boundary element method (BEM), is introduced. This procedure is based on the implicit discretisation in time (the so-called Rothe method) combined with a simple “algebraic” transformation of variables, leading to a numerically stable procedure (proved explicitly by discrete energy estimates), which can be easily implemented in a BEM code to solve initial-boundary value visco-elastic problems by using the Kelvin elastostatic fundamental solution only. It is worth mentioning that no inverse Laplace transform is required here. The formulation is straightforward for both 2D and 3D problems involving unilateral frictionless contact. Although the focus is to the simplest Kelvin–Voigt rheology, a generalization to Maxwell, Boltzmann, Jeffreys, and Burgers rheologies is proposed, discussed, and implemented in the BEM code too. A few 2D and 3D initial-boundary value problems, one of them with unilateral frictionless contact, are solved numerically.     

Theoretical and experimental study of the isothermal mechanical behaviour of alloys in the semi-solid state

An isothermal constitutive model for semi-solid alloys based on the concepts of mechanics of continuous media and the theory of mixtures is presented. The model is applicable to semi-solid states obtained either by solidification from liquid state or partial remelting from solid state in which each of the solid and the liquid phases is contiguous. During deformation their behaviours are coupled: the densification of the solid matrix considered as a porous viscoplastic medium saturated with a liquid drives the fluid flow behaviour, and the resulting pressure distribution in the liquid affects in turn the stresses and the densification of the solid. The identification procedure of the model uses two types of mechanical tests: uniaxial compression and drained die pressing (filtration) carried out with A356 alloy. The identification results are then validated using drained triaxial compression.