Modelling of Shockwave Propagation in Orthotropic Materials

2013 ◽  
Vol 315 ◽  
pp. 557-561 ◽  
Author(s):  
Mohd Khir Mohd Nor ◽  
Rade Vignjevic ◽  
James Campbell

Modelling of shockwave propagation in orthotropic materials requires an appropriate description of material behaviour within elastic and plastic regimes. To deal with this issues, a finite strain constitutive model for orthotropic materials was developed within a consistent thermodynamic framework of irreversible process in this paper. The important features of this material model are the multiplicative decomposition of the deformation gradient and a Mandel stress tensor combined with the new stress tensor decomposition generalised for orthotropic materials. The elastic free energy function and the yield function are defined within an invariant theory by means of the introduction of the structural tensors. The plastic behaviour is characterised within the associative plasticity framework using the Hills yield criterion. The complexity was further extended by coupling the formulation with the equation of state (EOS) to control the response of the material to shock loading. This material model which was developed and integrated in the isoclinic configuration provides a unique treatment for elastic and plastic anisotropy. The effects of elastic anisotropy are taken into account through the stress tensor decomposition and plastic anisotropy through yield surface defined in the generalized deviatoric plane perpendicular to the generalised pressure. To test its ability to describe shockwave propagation, the new material model was implemented into the LLNL-DYNA3D code. The results generated by the proposed material model were compared against the experimental Plate Impact test data of Aluminium Alloy 7010. A good agreement between experimental and simulation was obtained for two principal directions of material orthotropy.

2013 ◽  
Vol 315 ◽  
pp. 635-639 ◽  
Author(s):  
Mohd Khir Mohd Nor ◽  
Rade Vignjevic ◽  
James Campbell

The accuracy and reliability of the new stress tensor decomposition to capture the plasticity behaviour of orthotropic materials under plane-stress conditions was examined in this paper. No experiment was required to perform this work. Therefore, the suitable, published paper which provides a relevant test result and sufficient material properties to characterise the new stress tensor decomposition, was used. This new stress tensor decomposition was used to presents a new yield criterion for orthotropic sheet metals under plane-stress conditions in this work. This was done by assuming the yield surface to be circular in the new deviatoric plane. The predictions of the new effectice stress expression were then compared with the experimental data of 6000 series aluminium alloy sheet (A6XXX-T4) and Al-killed cold-rolled steel sheet SPCE. The predicted new yield surfaces are in good agreement with respect to the experimental data for two materials (A6XXX-T4 and SPCE).


Author(s):  
Ivaylo N. Vladimirov ◽  
Stefanie Reese

In this paper, a finite strain constitutive model for evolving elastic and plastic anisotropy combining nonlinear isotropic and kinematic hardening is presented. The evolution of elastic anisotropy is described by defining the Helmholtz free energy as an isotropic function of a family of evolving structure tensors. In addition, plastic anisotropy is modelled via the dependence of the yield surface on the same family of structure tensors. Exploiting the Clausius-Duhem inequality leads to the interesting result that all tensor-valued internal variables are symmetric. Thus, the integration of the evolution equations can be efficiently performed by means of an implicit algorithm that automatically retains the symmetry of the internal variables in every time step. The material model is used as a user material subroutine UMAT infinite element package ABAQUS/Standard, by means of which the occurrence of earing during cylindrical deep drawing is simulated.


Symmetry ◽  
2021 ◽  
Vol 13 (10) ◽  
pp. 1941
Author(s):  
Sergei Alexandrov ◽  
Yun-Che Wang ◽  
Lihui Lang

Plastic anisotropy significantly influences the behavior of structures subjected to various loading conditions. The extremum principles in the theory of rigid plastic solids are convenient and reliable tools for plastic design. The present paper combines the upper bound theorem and Hill’s quadratic yield criterion for orthotropic materials to evaluate the plastic collapse load of a highly undermatched welded tensile panel with a crack in the weld. The base material is supposed to be rigid. The shape of the crack is quite arbitrary. The orientation of the principal axes of anisotropy varies through the thickness of the weld. The upper bound solution is based on an exact solution for a layer of an anisotropic material. This feature of the upper bound solution is advantageous for increasing its accuracy. A numerical treatment is only necessary to find the solution for the uncracked specimen. This specimen has two axes of symmetry, which simplifies the solution. Simple analytic formulae transform this solution into a solution for the cracked specimens with one axis of symmetry and no symmetry. It is shown that the through-thickness distribution of anisotropic properties significantly affects the limit load.


2013 ◽  
Vol 554-557 ◽  
pp. 2330-2337
Author(s):  
Ivaylo N. Vladimirov ◽  
Stefanie Reese

Sheet metals exhibit anisotropic plastic behavior due to the large plastic deformations that occur during the rolling of the sheet and which induce texture and are responsible for the initial anisotropy. There exist various possibilities to introduce plastic anisotropy into the finite element modelling of sheet metal forming. The initial yield anisotropy can be incorporated either through an anisotropic yield surface or directly by means of a crystallographic texture model. Here, one basically differentiates between empirical and phenomenological anisotropic yield function equations, where the anisotropy coefficients can be obtained from mechanical tests, and texture-based models the coefficients of which are directly determined based on experimentally obtained orientation distributions. Another type of anisotropy that can be usually found in anisotropic materials is the elastic anisotropy. In metal plasticity one often considers the effect of elastic anisotropy significantly smaller than the effect of plastic anisotropy. Consequently, elastic isotropic expressions are often used for elastic stored energy functions with anisotropic yield criteria. However, the influence of elastic anisotropy in the elastoplastic behavior can be very important especially during elastic recovery processes during unloading after forming and springback. This research focuses, therefore, on the study of the influence of elastic anisotropy on the amount of springback in bending processes such as e.g. unconstrained bending. We discuss a finite strain material model for evolving elastic and plastic anisotropy combining nonlinear isotropic and kinematic hardening. The evolution of elastic anisotropy is described by representing the Helmholtz free energy as a function of a family of evolving structure tensors. In addition, plastic anisotropy is modelled via the dependence of the yield surface on the same family of structure tensors. Exploiting the dissipation inequality leads to the interesting result that all tensor-valued internal variables are symmetric. Thus, the integration of the evolution equations can be efficiently performed by means of an algorithm that automatically retains the symmetry of the internal variables in every time step. The material model has been implemented as a user material subroutine UMAT into the commercial finite element software ABAQUS/Standard and has been applied to the simulation of springback of unconstrained bending.


2010 ◽  
Vol 638-642 ◽  
pp. 3821-3826 ◽  
Author(s):  
Sergei Alexandrov

The limit load is an essential input parameter of flaw assessment procedures. The present paper deals with an effect of plastic anisotropy on its value. An upper bound solution for three-dimensional deformation of a highly under-matched welded specimen subject to tension is proposed. The base material is assumed to be rigid, and the weld material obeys Hill’s quadratic yield criterion for orthotropic materials. It is demonstrated that it is crucial to account for both plastic anisotropy and three dimensionality of deformation in limit load calculations for flaw assessment procedures.


Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. C85-C97 ◽  
Author(s):  
Nepomuk Boitz ◽  
Anton Reshetnikov ◽  
Serge A. Shapiro

Radiation patterns of earthquakes contain important information on tectonic strain responsible for seismic events. However, elastic anisotropy may significantly impact these patterns. We systematically investigate and visualize the effect of anisotropy on the radiation patterns of microseismic events. For visualization, we use a vertical-transverse-isotropic (VTI) medium. We distinguish between two different effects: the anisotropy in the source and the anisotropy on the propagation path. Source anisotropy mathematically comes from the matrix multiplication of the anisotropic stiffness tensor with the source strain expressed by the potency tensor. We analyze this effect using the corresponding radiation pattern and the moment tensor decomposition. Propagation anisotropy mathematically comes from the deviation between the polarization and the propagation direction of a quasi P-wave in an anisotropic medium. We investigate both effects separately by either assuming the source to be anisotropic and the propagation to be isotropic or vice versa. We find that both effects have a significant impact on the radiation pattern of a pure-slip source. Finally, we develop an alternative visualization of source mechanisms by plotting beach balls proportional to their potency tensors. For this, we multiply the potency tensor with an isotropic elasticity tensor having the equivalent shear modulus [Formula: see text] and [Formula: see text]. In this way, we visualize the tectonic deformation in the source, independently of the rock anisotropy.


2010 ◽  
Vol 77 (6) ◽  
Author(s):  
Min Kyoo Kang ◽  
Rui Huang

A hydrogel consists of a cross-linked polymer network and solvent molecules. Depending on its chemical and mechanical environment, the polymer network may undergo enormous volume change. The present work develops a general formulation based on a variational approach, which leads to a set of governing equations coupling mechanical and chemical equilibrium conditions along with proper boundary conditions. A specific material model is employed in a finite element implementation, for which the nonlinear constitutive behavior is derived from a free energy function, with explicit formula for the true stress and tangent modulus at the current state of deformation and chemical potential. Such implementation enables numerical simulations of hydrogels swelling under various constraints. Several examples are presented, with both homogeneous and inhomogeneous swelling deformation. In particular, the effect of geometric constraint is emphasized for the inhomogeneous swelling of surface-attached hydrogel lines of rectangular cross sections, which depends on the width-to-height aspect ratio of the line. The present numerical simulations show that, beyond a critical aspect ratio, creaselike surface instability occurs upon swelling.


2012 ◽  
Vol 504-506 ◽  
pp. 77-82 ◽  
Author(s):  
Liana Paraianu ◽  
Dan Sorin Comsa ◽  
Ioan Pavel Nicodim ◽  
Ioan Ciobanu ◽  
Dorel Banabic

The accuracy of the forming limit curves predicted by the Marciniak-Kuczynski model depends on the type and flexibility of the constitutive equations used to describe the mechanical response of the sheet metal. From this point of view, the yield criterion has the most significant influence. The paper presents a comparative analysis referring to the quality of the forming limit curves predicted by the Marciniak-Kuczynski model for the case when the plastic anisotropy of a DC04 sheet metal is described by the BBC2005 yield criterion. The coefficients included in the expression of the BBC2005 equivalent stress are evaluated using different identification strategies (with 4, 6, 7, and 8 mechanical parameters). The forming limit curves predicted by the Marciniak-Kuczynski model in each of the cases previously mentioned are compared with experimental data.


2018 ◽  
Vol 85 (12) ◽  
Author(s):  
Kelin Chen ◽  
Stelios Kyriakides ◽  
Martin Scales

The shear stress–strain response of an aluminum alloy is measured to a shear strain of the order of one using a pure torsion experiment on a thin-walled tube. The material exhibits plastic anisotropy that is established through a separate set of biaxial experiments on the same tube stock. The results are used to calibrate Hill's quadratic anisotropic yield function. It is shown that because in simple shear the material axes rotate during deformation, this anisotropy progressively reduces the material tangent modulus. A parametric study demonstrates that the stress–strain response extracted from a simple shear test can be influenced significantly by the anisotropy parameters. It is thus concluded that the material axes rotation inherent to simple shear tests must be included in the analysis of such experiments when the material exhibits anisotropy.


Sign in / Sign up

Export Citation Format

Share Document