Moderately large forced oblique vibrations of elastic- viscoplastic deteriorating slightly curved beams

Summary An integral equation formulation for the dynamic biaxial response of slightly curved elastic-viscoplastic beams is presented in the context of a multiple field analysis, which takes into account the geometrically nonlinear influence of moderately large deflections. Materials are considered in the regime of rate-dependent plasticity and are subjected to accumulated ductile damage. The latter is modeled by the growth of voids in the plastic zones of an initially porous elastic material. Inelastic defects of the material are considered in the linear elastic background beam by a second imposed strain field (eigenstrains). Geometrically nonlinear effects of large deflections under conditions of immovable supports are approximately taken into account. By inspection, they render another “strain field” to be imposed on the linear background beam. Superposition applies in the linear elastic background in an incremental formulation. Linear methods, as those based on Green's functions and Duhamel's integral, are used to account for the given loads as well as for the resultants of the imposed strain fields. The intensity and the distribution of the imposed strain fields are calculated incrementally in a time-stepping procedure. They are determined by the constitutive law and by application of the nonlinear geometric relations. The numerical procedure resulting from the multiple fields in the elastic background is illustrated for two cases: (1) a preloaded viscoplastic beam of rectangular cross section is subjected to oblique flexural vibrations when forced by a sinusoidal load, and (2) an I-beam with a prescribed initial curvature is severely impacted and thus driven into the plastic regime. Accepted for publication 22 November 1996     

On thermoelastostatics of composites with nonlocal properties of constituents II. Estimation of effective material and field parameters

One considers a linear thermoelastic composite medium, which consists of a homogeneous matrix containing a statistically homogeneous random set of ellipsoidal uncoated or coated heterogeneities. It is assumed that the stress–strain constitutive relations of constituents are described by the nonlocal integral operators, whereas the equilibrium and compatibility equations remain unaltered as in classical local elasticity. The general integral equations connecting the stress and strain fields in the point being considered and the surrounding points are obtained. The method is based on a centering procedure of subtraction from both sides of a known initial integral equation their statistical averages obtained without any auxiliary assumptions such as, e.g., effective field hypothesis implicitly exploited in the known centering methods. In a simplified case of using of the effective field hypothesis for analyzing composites with one sort of heterogeneities, one proves that the effective moduli explicitly depend on both the strain and stress concentrator factor for one heterogeneity inside the infinite matrix and does not directly depend on the elastic properties (local or nonlocal) of heterogeneities. In such a case, the Levin’s (1967) formula in micromechanics of composites with locally elastic constituents is generalized to their nonlocal counterpart. A solution of a volume integral equation for one heterogeneity subjected to inhomogeneous remote loading inside an infinite matrix is proposed by the iteration method. The operator representation of this solution is incorporated into the new general integral equation of micromechanics without exploiting of basic hypotheses of classical micromechanics such as both the effective field hypothesis and “ellipsoidal symmetry” assumption. Quantitative estimations of results obtained by the abandonment of the effective field hypothesis are presented.     

A new class of micro–macro models for elastic–viscoplastic heterogeneous materials

The determination of the effective behavior of heterogeneous materials from the properties of the components and the microstructure constitutes a major task in the design of new materials and the modeling of their mechanical behavior. In real heterogeneous materials, the simultaneous presence of instantaneous mechanisms (elasticity) and time dependent ones (non-linear viscoplasticity) leads to a complex space–time coupling between the mechanical fields, difficult to represent in a simple and efficient way. In this work, a new self-consistent model is proposed, starting from the integral equation for a translated strain rate field. The chosen translated field is the (compatible) viscoplastic strain rate of the (fictitious) viscoplastic heterogeneous medium submitted to a uniform (unknown) boundary condition. The self-consistency condition allows to define these boundary conditions so that a relative simple and compact strain rate concentration equation is obtained. This equation is explained in terms of interactions between an inclusion and a matrix, which lead to interesting conclusions. The model is first applied to the case of two-phase composites with isotropic, linear and incompressible viscoelastic properties. In that case, an exact self-consistent solution using the Laplace–Carson transform is available. The agreement between both approaches appears quite good. Results for elastic–viscoplastic BCC polycrystals are also presented and compared with results obtained from Kröner–Weng's and Paquin et al. (Arch. Appl. Mech. 69 (1999) 14)'s model.     

Nonlinear rheology of concentrated spherical silica suspensions: 3. Concentration dependence

Nonlinear rheology was examined for concentrated suspensions of spherical silica particles (with radius of 40 nm) in viscous media, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture and pure ethylene glycol. The particles were randomly and isotropically dispersed in the media in the quiescent state, and their effective volume fraction φ_eff ranged from 0.36 to 0.59. For small strains, the particles exhibited linear relaxation of the Brownian stress σ_B due to their diffusion. For large step strains γ, the nonlinear relaxation modulus G(t,γ) exhibited strong damping and obeyed the time-strain separability. This damping was related to γ-insensitivity of strain-induced anisotropy in the particle distribution that resulted in decreases of σ_B/γ. The damping became stronger for larger φ_eff. This φ_eff dependence was related to a hard-core volume effect, i.e., strain-induced collision of the particles that is enhanced for larger φ_eff. Under steady/transient shear flow, the particles exhibited thinning and thickening at low and high γ˙, respectively. The thinning behavior was well described by a BKZ constitutive equation using the G(t,γ) data and attributable to decreases of a Brownian contribution, σ_B/γ˙. The thickening behavior, not described by this equation, was related to dynamic clustering of the particles and corresponding enhancement of the hydrodynamic stress at high γ˙. In this thickening regime, the viscosity growth η₊ after start-up of flow was scaled with a strain γ˙t. Specifically, critical strains γ_d and γ_s for the onset of thickening and achievement of the steadily thickened state were independent of γ˙ but decreased with increasing φ_eff. This φ_eff dependence was again related to the hard-core volume effect, flow-induced collision of the particles enhanced for larger φ_eff. Received: 26 June 1998 Accepted: 9 December 1998     

Accounting for the elastic properties of a non-Newtonian material under its viscosimetric flow

This paper considers the deformation and viscoplastic flow of a non-Newtonian material enclosed between coaxial rigid cylindrical surfaces, each of which performs a rotation followed by a stop and a rotation in the opposite direction. The problem is solved using the model of large elastoviscoplastic deformations, in contrast to the classical solutions obtained using the model of a rigid viscoplastic body. The parameters of the viscosimetric process are calculated in both the region of viscoplastic flow developed and the region of elastic deformation. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 2, pp. 143–151, March–April, 2008.     

Finite viscoplasticity of amorphous glassy polymers in the logarithmic strain space

The paper outlines a new constitutive model and experimental results of rate-dependent finite elastic–plastic behavior of amorphous glassy polymers. In contrast to existing kinematical approaches to finite viscoplasticity of glassy polymers, the formulation proposed is constructed in the logarithmic strain space and related to a six-dimensional plastic metric. Therefore, it a priori avoids difficulties concerning with the uniqueness of a plastic rotation. The constitutive framework consists of three major steps: (i) A geometric pre-processing defines a total and a plastic logarithmic strain measures determined from the current and plastic metrics, respectively. (ii) The constitutive model describes the stresses and the consistent moduli work-conjugate to the logarithmic strain measures in an analogous structure to the geometrically linear theory. (iii) A geometric post-processing maps the stresses and the algorithmic tangent moduli computed in the logarithmic strain space to their nominal, material or spatial counterparts in the finite deformation space. The analogy between the formulation of finite plasticity in the logarithmic strain space and the geometrically linear theory of plasticity makes this framework very attractive, in particular regarding the algorithmic implementation. The flow rule for viscoplastic strains in the logarithmic strain space is adopted from the celebrated double-kink theory. The post-yield kinematic hardening is modeled by different network models. Here, we compare the response of the eight chain model with the newly proposed non-affine micro-sphere model. Apart from the constitutive model, experimental results obtained from both the homogeneous compression and inhomogeneous tension tests on polycarbonate are presented. Besides the load–displacement data acquired from inhomogeneous experiments, quantitative three-dimensional optical measurements of the surface strain fields are carried out. With regard to these experimental data, the excellent predictive quality of the theory proposed is demonstrated by means of representative numerical simulations.     

Incremental Magnetoelastic Deformations,with Application to Surface Instability

In this paper the equations governing the deformations of infinitesimal (incremental) disturbances superimposed on finite static deformation fields involving magnetic and elastic interactions are presented. The coupling between the equations of mechanical equilibrium and Maxwell’s equations complicates the incremental formulation and particular attention is therefore paid to the derivation of the incremental equations, of the tensors of magnetoelastic moduli and of the incremental boundary conditions at a magnetoelastic/vacuum interface. The problem of surface stability for a solid half-space under plane strain with a magnetic field normal to its surface is used to illustrate the general results. The analysis involved leads to the simultaneous resolution of a bicubic and vanishing of a 7×7 determinant. In order to provide specific demonstration of the effect of the magnetic field, the material model is specialized to that of a “magnetoelastic Mooney–Rivlin solid”. Depending on the magnitudes of the magnetic field and the magnetoelastic coupling parameters, this shows that the half-space may become either more stable or less stable than in the absence of a magnetic field.      

Iteration method in linear elasticity of random structure composites containing heterogeneities of noncanonical shape

We consider a linearly elastic composite medium, which consists of a homogeneous matrix containing a statistically inhomogeneous random set of heterogeneities of arbitrary shape. The general integral equations connecting the stress and strain fields in the point being considered with the stress and strain fields in the surrounding points are obtained for the random fields of heterogeneities. The method is based on a recently developed centering procedure where the notion of a perturbator is introduced and statistical averages are obtained without any auxiliary assumptions such as, e.g., effective field hypothesis implicitly exploited in the known centering methods. Effective elastic moduli and the first statistical moments of stresses in the heterogeneities are estimated for statistically homogeneous composites with the general case of both the shape and inhomogeneity of the heterogeneities moduli. The explicit new representations of the effective moduli and stress concentration factors are built by the iteration method in the framework of the quasicristallite approximation but without basic hypotheses of classical micromechanics such as both the EFH and “ellipsoidal symmetry” assumption. Numerical results are obtained for some model statistically homogeneous composites reinforced by aligned identical homogeneous heterogeneities of noncanonical shape. Some new effects are detected that are impossible in the framework of a classical background of micromechanics.     

Electroelastic field for an impermeable anti-plane shear crack in a piezoelectric ceramics plate

Electroelastic behavior of a cracked piezoelectric ceramics plate subjected to four cases of combined mechanical-electrical loads is analyzed. The integral transform method is applied to convert the problem involving an impermeable anti-plane crack to dual integral equations. Solving the resulting equations, the explicit analytic expressions for electroelastic field along the crack line and the intensity factors of relevant quantities near the crack tip and the mechanical strain energy release rate are obtained. The known results for an infinite piezoelectric ceramics plane containing an impermeable anti-plane crack are recovered from the present results only if the thickness of the plate h → ∞. Biography: LI Xian-fang (1964-)     

ASYMPTOTIC ANALYSIS OF MODE Ⅱ STATIONARY GROWTHCRACK ON ELASTIC-ELASTIC POWER LAWCREEPING BIMATERIAL INTERFACE

A mechanical model was established for mode Ⅱ interfacial crack static growingalong an elastic-elastic power law creeping bimaterial interface. For two kinds of boundaryconditions on crack faces, traction free and frictional contact, asymptotic solutions of thestress and strain near tip-crack were given. Results derived indicate that the stress andstrain have the same singularity, there is not the oscillatory singularity in the field; thecreep power-hardening index n and the ratio of Young‘s module notably influence the crack-tip field in region of elastic power law creeping material and n only influences distribution ofstresses and strains in region of elastic material. When n is bigger, the creepingdeformation is dominant and stress fields become steady, which does not change with n.Poisson‘s ratio does not affect the distributing of the crack-tip field.     

Rheology of lyotropic and thermotropic lamellar phases

We investigate the steady-state rheological behaviour of the lamellar phase of a lyotropic system (CpCl, hexanol, brine) and of a thermotropic system (8CB). Power laws characterize the behaviour of the imposed stress as a function of the measured shear rate and similarities are observed for both systems; the same regime γ˙∼σ ^m with m≈1.7 is obtained at low shear stresses corresponding to a texture of oily streaks oriented in the direction of the flow, as shown by microscopic observations. The “onion state” only exists in the case of dilute samples of the lyotropic lamellar phase; the stress then varies as γ˙∼σ ^m with m≈4.8, as already observed by other groups with different systems. Rheological measurements at different temperatures allow determination of different activation energies relating to the still badly understood processes involved in the different rheological regimes. We propose a model which reproduces the experimental power laws and which is based on an analogy with the theory of high-temperature creep in metals and alloys. Received: 19 October 1999/Accepted: 1 November 1999     

Mechanical Performance of a Recycled Carbon Fibre/PP Composite

A composite made of recycled carbon fibres in recycled polypropylene matrix is studied experimentally to describe the features of the elastic and time dependent nonlinear mechanical behaviour. The properties of the developed material have a large variability to be addressed and understood. It was found that the stress-strain curves in tension are rather nonlinear at low strain rate and the strength is sensitive to strain rate. The elastic properties’ reduction for this composite after loading to high strains is rather limited. More important is that even in the “elastic region” due to viscoelastic effects the slope of loading–unloading curve is not the same and that at higher stress large viscoplastic strains develop and creep rupture is typical. The time and stress dependence of viscoplastic strains was analysed and described theoretically. The viscoelastic response of the composite was analysed using creep compliance, which was found to be slightly nonlinear.     

New integral formulation and self-consistent modeling of elastic–viscoplastic heterogeneous materials

Predicting the overall behavior of heterogeneous materials, from their local properties at the scale of heterogeneities, represents a critical step in the design and modeling of new materials. Within this framework, an internal variables approach for scale transition problem in elastic–viscoplastic case is introduced. The proposed micromechanical model is based on establishing a new system of field equations from which two Navier’s equations are obtained. Combining these equations leads to a single integral equation which contains, on the one hand, modified Green operators associated with elastic and viscoplastic reference homogeneous media, and secondly, elastic and viscoplastic fluctuations. This new integral equation is thus adapted to self-consistent scale transition methods. By using the self-consistent approximation we obtain the concentration law and the overall elastic–viscoplastic behavior of the material. The model is first applied to the case of two-phase materials with isotropic, linear and compressible viscoelastic properties. Results for elastic–viscoplastic two-phase materials are also presented and compared with exact results and variational methods.     

Semi-weight function method on computation of stress intensity factors in dissimilar materials

Semi-weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of continuity across interface and free crack surface, the stress and displacement fields were obtained. The eigenvalue of these fields is lambda. Semi-weight functions were obtained as virtual displacement and stress fields with eigenvalue-lambda. Integral expression of fracture parameters, KⅠ and KⅡ, were obtained from reciprocal work theorem with semi-weight functions and approximate displacement and stress values on any integral path around crack tip. The calculation results of applications show that the semi-weight function method is a simple, convenient and high precision calculation method.     

A non-linear uniaxial integral constitutive equation incorporating rate effects,creep and relaxation

A previously proposed first order non-linear differential equation for uniaxial viscoplasticity, which is non-linear in stress and strain but linear in stress and strain rates, is transformed into an equivalent integral equation. The proposed equation employs total strain only and is symmetric with respect to the origin and applies for tension and compression. The limiting behavior for large strains and large times for monotonic, creep and relaxation loading is investigated and appropriate limits are obtained. When the equation is specialized to an overstress model it is qualitatively shown to reproduce key features of viscoplastic behavior. These include: initial linear elastic or linear viscoelastic response: immediate elastic slope for a large instantaneous change in strain rate normal strain rate sensitivity and non-linear spacing of the stress-strain curves obtained at various strain rates; and primary and secondary creep and relaxation such that the creep (relaxation) curves do not cross. Isochronous creep curves are also considered. Other specializations yield wavy stress-strain curves and inverse strain rate sensitivity. For cyclic loading the model must be modified to account for history dependence in the sense of plasticity.     

The second-law analysis of mixed convection in rectangular ducts

 The rate of entropy generation, S˙ ^′ _G[W/mK], is examined both theoretically and numerically for forced and mixed convection in a rectangular duct heated at the bottom. Under fully-developed flow conditions S˙ ^′ _G is expressed in terms of relevant non-dimensional hydrodynamic and thermal parameters. Numerically, it is demonstrated that S˙ ^′ _G is a single, effective parameter to examine both thermal and hydrodynamic fields and their variations. Received on 22 November 1999     

An admissible form of an inelastic constitutive equation—I. General theory

From the viewpoint of irreversible thermodynamics an admissible form of rate-type constitutive equation of inelastic materials is given. The displacement gradient tensor F referred to the temporarily fixed reference frame which coincides with the Euler frame at the instant of the reference time is decomposed linearly into elastic and inelastic parts so that the procedure of formulation is simplified and clarified. The inelastic deformation rate is directly related to the internal production rate of entropy. The existence of an inelastic potential of the usual sense is not assumed, though the result can be understood to include the conventional flow theory based on an inelastic potential. An example of an elastoviscoplastic constitutive equation is given and some properties of yield surfaces are discussed.     

Exact Connections of Warping Fields and Torsional Rigidities of Symmetric Composite Bars

Saint-Venant's torsion of symmetric cylindrical bars consisting of two or four homogeneous phases is studied. A symmetric section is meant that the cross section of the cylindrical bar possesses reflectional symmetry with respect to one or more axes. Each constituent region may have different shear modulus. The idea of the analysis is to superimpose suitably reflected potentials to obtain the torsion solution of the same composite section but with different moduli. For two-phase sections, we show that, if the warping fields for a given symmetric section with phase shear moduli μ₁ and μ₂ are known a priori, then the warping fields for the same configuration but with a different set of constituent moduli μ₁ ^′ and μ₂ ^′ are readily found through simple linear superpositions. Further, suppose that the torsional rigidities T(μ₁,μ₂) and T(μ₁ ^′,μ₂ ^′) for any two sets of phase moduli can be measured by some experimental tests or evaluated through numerical procedures, then the torsional rigidity for any other combinations of constituent moduli T(μ₁ ^′′,μ₂ ^′′) can be exactly determined without any recourse to the field solutions of governing differential equations. Similar procedures can be applied to a 4-phase symmetric section. But the coefficients of superposition are only found for a few branches. Specifically, we find that depending on the conditions of μ and μ^′, admissible solutions can be divided into three categories. When the correspondence between the warping field is known to exist, a link between the torsional rigidities can be established as well. This revised version was published online in June 2006 with corrections to the Cover Date.     

A new constitutive equation for elastoviscoplastic fluid flows

From thermodynamic theory, a new three-dimensional model for elastoviscoplastic fluid flows is presented. It extends both the Bingham viscoplastic and the Oldroyd viscoelastic models. Fundamental flows are studied: simple shear flow, uniaxial elongation and large amplitude oscillatory shear. The complex moduli (G^′,G^″)$(G^{\prime }\text{,}{G}^{″})$ are found to be in qualitative agreement with experimental data for materials that present microscopic network structures and large scale rearrangements. Various fluids of practical interest, such as liquid foams, droplet emulsions or blood, present such elastoviscoplastic behavior: at low stress, the material behaves as a viscoelastic solid, whereas at stresses above a yield stress, the material behaves as a fluid.