Nonlinear Pseudo-plastic Orthotropic Solids (LAWS28, 50 and 68)

Conventional Nonlinear Pseudo-plastic Orthotropic Solids (LAW28 and LAW50)

These laws are generally used to model honeycomb material structures as crushable foams. The microscopic behavior of this kind of materials can be considered as a system of three independent orthogonal springs. The nonlinear behavior in orthogonal directions can then be determined by experimental tests. The behavior curves are injected directly in the definition of law. Therefore, the physical behavior of the material can be obtained by a simple law. However, the microscopic elasto-plastic behavior of a material point cannot be represented by decoupled unidirectional curves. This is the major drawback of the constitutive laws based on this approach. The cell direction is defined for each element by a local frame in the orthotropic solid property. If no property set is given, the global frame is used.


Figure 1. Local Frame Definition
The Hooke matrix defining the relation between the stress and strain tensors is diagonal, as there is no Poisson's effect:(1)
[ σ 11 σ 22 σ 33 σ 12 σ 23 σ 31 ]=[ E 11 0 0 0 0 0 0 E 12 0 0 0 0 0 0 E 33 0 0 0 0 0 0 G 12 0 0 0 0 0 0 G 23 0 0 0 0 0 0 G 31 ][ ε 11 ε 22 ε 33 ε 12 ε 23 ε 31 ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aadaWadaabaeqabaGaeq4Wdm3aaSbaaSqaaiaaigdacaaIXaaabeaa aOqaaiabeo8aZnaaBaaaleaacaaIYaGaaGOmaaqabaaakeaacqaHdp WCdaWgaaWcbaGaaG4maiaaiodaaeqaaaGcbaGaeq4Wdm3aaSbaaSqa aiaaigdacaaIYaaabeaaaOqaaiabeo8aZnaaBaaaleaacaaIYaGaaG 4maaqabaaakeaacqaHdpWCdaWgaaWcbaGaaG4maiaaigdaaeqaaaaa kiaawUfacaGLDbaacqGH9aqpdaWadaqaauaabeqGgyaaaaaabaGaam yramaaBaaaleaacaaIXaGaaGymaaqabaaakeaacaaIWaaabaGaaGim aaqaaiaaicdaaeaacaaIWaaabaGaaGimaaqaaiaaicdaaeaacaWGfb WaaSbaaSqaaiaaigdacaaIYaaabeaaaOqaaiaaicdaaeaacaaIWaaa baGaaGimaaqaaiaaicdaaeaacaaIWaaabaGaaGimaaqaaiaadweada WgaaWcbaGaaG4maiaaiodaaeqaaaGcbaGaaGimaaqaaiaaicdaaeaa caaIWaaabaGaaGimaaqaaiaaicdaaeaacaaIWaaabaGaam4ramaaBa aaleaacaaIXaGaaGOmaaqabaaakeaacaaIWaaabaGaaGimaaqaaiaa icdaaeaacaaIWaaabaGaaGimaaqaaiaaicdaaeaacaWGhbWaaSbaaS qaaiaaikdacaaIZaaabeaaaOqaaiaaicdaaeaacaaIWaaabaGaaGim aaqaaiaaicdaaeaacaaIWaaabaGaaGimaaqaaiaadEeadaWgaaWcba GaaG4maiaaigdaaeqaaaaaaOGaay5waiaaw2faamaadmaaeaqabeaa cqaH1oqzdaWgaaWcbaGaaGymaiaaigdaaeqaaaGcbaGaeqyTdu2aaS baaSqaaiaaikdacaaIYaaabeaaaOqaaiabew7aLnaaBaaaleaacaaI ZaGaaG4maaqabaaakeaacqaH1oqzdaWgaaWcbaGaaGymaiaaikdaae qaaaGcbaGaeqyTdu2aaSbaaSqaaiaaikdacaaIZaaabeaaaOqaaiab ew7aLnaaBaaaleaacaaIZaGaaGymaaqabaaaaOGaay5waiaaw2faaa aa@8DC6@
An isotropic material may be obtained if:(2)
E 11 = E 22 = E 33 and G 12 = G 23 = G 31 = E 11 2
Plasticity may be defined by a volumic strain or strain dependent yield curve (Figure 2). The input yield stress function is always positive. If the material undergoes plastic deformation, its behavior is always orthotropic, as all curves are independent to each other.


Figure 2. Honeycomb Typical Constitutive Curve

The failure plastic strain may be input for each direction. If the failure plastic strain is reached in one direction, the element is deleted. The material law may include strain rate effects (LAW50) or may not (LAW28).

Cosserat Medium for Nonlinear Pseudo-plastic Orthotropic Solids (LAW68)

Conventional continuum mechanics approaches cannot incorporate any material component length scale. However, a number of important length scales as grains, particles, fibers, and cellular structures must be taken into account in a realistic model of some kinds of materials. To this end, the study of a microstructure material having translational and rotational degrees-of-freedom is underlying. The idea of introducing couple stresses in the continuum modelling of solids is known as Cosserat theory which returns back to the works of brothers Cosserat in the beginning of 20th century. 1 A recent renewal of Cosserat mechanics is presented in several works of Forest 2 3 4 5 A short summary of these publications is presented in this section.

Cosserat effects can arise only if the material is subjected to non-homogeneous straining conditions. A Cosserat medium is a continuous collection of particles that behave like rigid bodies. It is assumed that the transfer of the interaction between two volume elements through surface element dS occurs not only by means of a traction and shear forces, but also by moment vector as shown in Figure 3.


Figure 3. Equilibrium of Cosserat Volume Element
Surface forces and couples are then represented by the generally non-symmetrical force-stress and couple-stress tensors σ i j and μ i j (units MPA and MPa-m):(3)
t i = σ i j n j : m i = μ i j n j

The force and couple stress tensors must satisfy the equilibrium of momentums:

σ i j , j + f i = ρ u ¨ i
(4)
μ i j , j ε i k l σ k l + c i = I φ ¨ i
Where,
f i
Volume forces
c i
Volume couples
ρ
Mass density
I MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamysaaaa@36C4@
Isotropic rotational inertia
ε i k l
Signature of the perturbation (i,k,l)
In the often used couple-stress, the Cosserat micro-rotation is constrained to follow the material rotation given by the skew-symmetric part of the deformation gradient:(5)
φ i = 1 2 ε i j k u j , k

The associated torsion-curvature and couple stress tensors are then traceless. If a Timoshenko beam is regarded as a one-dimensional Cosserat medium, constraint Equation 5 is then the counterpart of the Euler-Bernoulli conditions.

The resolution of the previous boundary value problem requires constitutive relations linking the deformation and torsion-curvature tensors to the force- and couple-stresses. In the case of linear isotropic elasticity, you have:(6)
σ i j = λ e k k δ i j + 2 μ e i j S y m m . + 2 μ c e i j S k e w S y m m . μ i j = α k k k δ i j + 2 β κ i j S y m m . + 2 γ κ i j S k e w S y m m .

Where, e i j S y m m . and e i j S k e w S y m m . are respectively the symmetric and skew-symmetric part of the Cosserat deformation tensor. Four additional elasticity moduli appear in addition to the classical Lamé constants.

Cosserat elastoplasticity theory is also well-established. von Mises classical plasticity can be extended to micropolar continua in a straightforward manner. The yield criterion depends on both force- and couple-stresses:(7)
f ( σ , μ ) = 3 2 ( a 1 s i j s i j + a 2 s i j s j i + b 1 μ i j μ i j + b 2 μ i j μ j i ) R
Where,
s MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4Caaaa@36EE@
Stress deviator
a i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGPbaabeaaaaa@37F6@ and b i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGPbaabeaaaaa@37F6@
Material constants

Cosserat continuum theory can be applied to several classes of materials with microstructures as honeycombs, liquid crystals, rocks and granular media, cellular solids and dislocated crystals.

1 Cosserat E. and Cosserat F., “Theory of Deformable Bodies”, Hermann, Paris, 1909.
2 Forest S. and Sab K., “Cosserat overall modeling of heterogeneous materials”, Mechanics Research Communications, Vol. 25, No. 4, pp. 449-454, 1998.
3 Forest S., Cailletaud G. and Sievert R., “A Cosserat Theory for Elastoviscoplastic Single Crystals at Finite Deformation”, Archives of Mechanics, Vol. 49, pp. 705-736, 1997.
4 Besson J., Bultel F., and Forest S., “Plasticity Cosserat media. Application to particular composites particles”, 4th Symposium Calculation of Structures, CSMA/Teksea, Toulouse, pp. 759-764, 1999.
5 Forest S., “Cosserat Media”, ed. by K.H.J. Buschow, R.W. Cahn, M.C. Flemings, B. Ilschner, E.J. Kramer and S. Mahajan, Encyclopedia of Materials, Science and Technology, Elsevier, 2001.