Effect of fiber orientation on the structure of shock waves in carbon fiber-epoxy composites

Using an analytical relation between the Hugoniot (anisotropic and isotropic) states and other thermodynamic (anisotropic and isotropic) states at high pressures, the effect of fiber orientation on the structure of shock waves in carbon fiber-epoxy composites of various symmetry is investigated. A correct nonlinear model of propagation of shock waves in anisotropic materials is proposed, which employs the conception of total generalized pressure and the pressure corresponding to the thermodynamic response, i.e., to the equation of state. The equation generalizes the nonlinear Hugoniot equation to anisotropic materials and is reduced to the classical variant in the case of isotropy. Invoking the relations of nonlinear anisotropic solids and the generalized decomposition of stress tensor, the double structure of shock waves, consisting of nonlinear anisotropic and isotropic elastic parts, is examined. The numerical calculations of Hugoniot levels of stress agree well with experimental data for a carbon fiber-epoxy composite selected.     

Anisotropic elasto-plastic damage formulation for finite strains using a damage structural tensor and the effective stress concept

Using the strain equivalence principle and the effective stress concept an anisotropic finite strain damage model is proposed as a direct extension of the classical isotropic LEMAITRE damage model to the anisotropic finite strain case. The damage tensor is included as a structural tensor in the complementary energy potential. This approach allows to consider a wide range of anisotropic damage phenomena on the level of continuum mechanics. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A phenomenological method of calculating the residual stresses and plastic deformations in a hollow surface-hardened cylindrical sample

A phenomenological method is proposed for calculating the residual stress and plastic deformation fields in a hollow surface-hardened cylindrical sample. Versions of the hardening are considered that lead to isotropy and anisotropy in the plastic deformations in the surface layer. A hardening anisotropy parameter is introduced that relates the axial and circumferential components of the residual plastic deformation tensor. The experimentally determined axial and/or circumferential components of the residual plastic stress tensor are used as the initial information. The tensor fields of the residual stresses and deformations are constructed assuming the hypothesis of surface hardening anisotropy and the absence of secondary plastic compression deformations and that the tangential components of the residual stress tensor and the plastic incompressibility of the material are small. A technique is developed for identifying the parameters of the proposed method. The adequacy is checked using experimental data for test pieces of type 45 and 12X18H10T steels hardened by hydro-shot blasting treatment and of type 45 steel hardened by treatment with a roller. Good agreement is observed between the calculated and experimental results. It is noted that the anisotropic hardening procedure leads to a substantial difference between the circumferential and axial components of the residual stresses in the hardened layer, unlike the case of isotropic hardening where they are practically identical.     

Anisotropic stress-softening model for compressible solids

A general constitutive theory for anisotropic stress softening in compressible solids is presented. The constitutive equation describes anisotropic strain induced behaviour of an initially “isotropic” virgin material. Parameters which characterise damage are proposed together with a concept of damage function. In order to develop an anisotropic stress-softening theory for compressible materials in close parallel to a recent incompressible anisotropic theory, the right stretch tensor is decomposed into its isochoric and dilatational parts. The ’free’ energy is expressed as a function of the dilatation, modified principal stretches, a volume change parameter and invariants of the dyadic products of the principal directions of the right stretch tensor and two structural tensors. A class of free energy functions is discussed and a special form of this class which satisfies the Clausius–Duhem inequality is proposed. Results of the theory applied to uniaxial tension, bulk compression and simple shear deformations are given. A sequence of deformations involving shear, hydrostatic-compression and hydrostatic-tension deformations is also investigated. In the case of hydrostatic-tension deformation, the stress softening is due to cavitation damage. The theoretical results obtained are consistent with expected behaviour and compare well with experimental data.     

非局部Kelvin粘弹性直杆受轴向拉力作用的应变分析

在Kelvin粘弹性体模型中引入非局晨应力应变关系,得到了粘弹性体的非局部本构方程,研究了符合该种本构关系的直杆受到轴向拉力作用的应变响应问题．首先通过变换将应变响应的求解问题转化为Volterra积分方程形式,然后采用对称的指数型核函数,利用Neumann级数展开求解了Volterra积分方程,得到了直杆的应变场．数值算例的计算结果显示了直杆受轴向拉力作用后的蠕变过程,当时间趋近无穷大时,计算结果则退化为非局部弹性计算结果．     

关于损伤张量的阶次

本文首先讨论了较为广泛的连续介质材料的应力变形本构关系,得到了通常以泛函表示的应力变形本构关系的张量表达式.以此为基础,研究了各向异性材料各向异性损伤时,无论从连续介质力学模型出发还是从缺陷模型出发,描述损伤的张量都存在最高阶次的限制;指出了在什么条件下,损伤变量可用低于最高阶次的张量来描述.     

On the dependence of the functions of state of a deformable body on the measure of rotation

We show that the measures of strain and initial values of vector and tensor state parameters are divided into subjective nonrotational and objective rotational. Representations of the functions of state are divided in a similar way, and only objective ones do not depend explicitly on the measure of rotation of material axes. We have constructed relations for the reduction of rotation-free differentials of the functions of state to expressions in terms of the tensors of infinitesimal strain and rotation. On this basis, we have obtained objective representations for stress tensor in terms of the derivatives of the state potential of an anisotropic material. The results obtained concern the nonlinear mechanics of initially stressed bodies.     

Non-linear relations of anisotropic elasticity and a particular postulate of isotropy

A general approach to the construction of six-dimensional images of strain processes is proposed with the introduction of a vector basis which, in special cases, is identical to the well-known bases of A. A. Il’yushin, V. V. Novozhilov and Ye. I. Shemyakin and S. A. Khristianovich. The analysis of the properties of materials is based on the use of the concept of characteristic elastic states, which was introduced in the papers of J. Rychlewski. In the case of an isotropic material and four types of anisotropic materials belonging to the cubic, hexagonal, trigonal and tetragonal systems, characteristic subspaces, corresponding to the multiple eigenvalues of the elasticity tensor are defined in a six-dimensional space. In accordance with Hooke's law, the components of the stress and strain vectors in these subspaces preserve their axial alignment for any of their orthogonal transformations. The particular postulate of isotropy, formulated by Il’yushin, is therefore satisfied by definition within the framework of isotropic characteristic subspaces for linear elastic materials. An extension of the particular postulate to strain processes in non-linear anisotropic materials is proposed, on the basis of which a general form of constitutive relations is obtained containing a minimum number of experimentally determinable material functions.     

Simple loading of a polymer material with nonlinear creep

Nonlinear tensor relations between strain, stress, and time are examined for a memory-type medium using degenerate kernels. The material parameters are determined from creep tests in a simple state of stress. Expressions for the strain associated with a complex state of stress and simple loading, found on the basis of the local strains theory, are in satisfactory agreement with the experimental data obtained for specimens of high-density polyethylene.Mekhanika Polimerov, Vol. 3, No. 2, pp. 236–242, 1967     

Monte Carlo approximate tensor moment simulations

An algorithm to generate samples with approximate first‐order, second‐order, third‐order, and fourth‐order moments is presented by extending the Cholesky matrix decomposition to a Cholesky tensor decomposition of an arbitrary order. The tensor decomposition of the first‐order, second‐order, third‐order, and fourth‐order objective moments generates a non‐linear system of equations. The algorithm solves these equations by numerical methods. The results show that the optimization algorithm delivers samples with an approximate residual error of less than 10¹⁶ between the components of the objective and the sample moments. The algorithm is extended for a n‐th‐order approximate tensor moment version, and simulations of non‐normal samples replicated from distributions with asymmetries and heavy tails are presented. An application for sensitivity analysis of portfolio risk assessment with Value‐at‐Risk (VaR) is provided. A comparison with previous methods available in the literature suggests that the methodology proposed reduces the error of the objective moments in the generated samples. ? ? JEL Classification: C14, C15, G32.
Copyright © 2016 John Wiley & Sons, Ltd.     

对称张量的分解和它的应用

本文把任一对称张量分解成两个张量的和,其中之一是“应力型”张量,另一个是“应变型”张量.对称张量空间被分解成两个直交子空间的直和.并用几何语言证明了弹性力学的几个基本原理.     

The representation of the displacement gradient of isotropic elastic body in terms of the piola stress tensor : PMM vol. 40, n 6, 1976, pp. 1070–1077

The representation of the displacement gradient of an isotropic elastic body is analyzed. It is shown on the basis of a single controlling inequality and a polar expansion of the Piola tensor that such representation has generally four branches. The mechanical meaning and the nature of that ambiguity is explained. It is established that when the angles of turn of material fibers are not excessively large, only one of the four branches is obtained. Particular cases in which the nature of ambiguity is more complex are investigated. It is noted that in many practical problems the representation of the displacement gradient by the Piola stress tensor is unambiguous.The considered problem is associated with the variational principle of complementary energy in the nonlinear theory of elasticity, where the statistically feasible fields of the asymmetric Piola stress tensor is varied [1], A method was proposed there for expressing the displacement gradient in terms of the Piola stress tensor for an isotropic elastic body. Later the concept of complementary energy and the representation of the strain gradient in terms of the Piola stress tensor were considered in [2, 3]. Examples of the use of the complementary energy concept are given in [2] and the case of an anisotropic body is considered in [3], These investigations disclosed that the considered representation of the strain tensor leads to ambiguity, but the character and nature of the ambiguity were not fully investigated.     

Quasistatic anti-plane motion in the simplest theory of nonlinear viscoelasticity

We consider a very simple model in the framework of differential viscoelastic materials which are isotropic and incompressible. In this model the Cauchy stress tensor is split in an elastic part and a dissipative part. The elastic part is derived from a strain-energy density function only of the first invariant of the Cauchy–Green strain tensor. The dissipative part is like the Navier–Stokes equations: linear in the stretching tensor with a constant viscosity parameter. For this model we provide some time and spatial estimates in the quasistatic approximations for the equations governing anti-plane shear motions. Several explicit examples for specific form of the strain energy are produced. Our results impose analytical restrictions on the mathematical properties of the strain energy to ensure a physical behavior in the creep and recovery experiments. Moreover, we show polynomial decay for the spatial behavior in the class of stress-hardening (or strain-stiffening) materials. For stress-softening materials a Phragmen–Lindelof alternative is proved.     

Topological sensitivity analysis of a multi‐scale constitutive model considering a cracked microstructure

This paper deals with the sensitivity analysis of the macroscopic elasticity tensor to topological microstructural changes of the underlying material. In particular, the microstucture is topologicaly perturbed by the nucleation of a small circular inclusion. The derivation of the proposed sensitivity relies on the concept of topological derivative, applied within a variational multi‐scale constitutive framework where the macroscopic strain and stress at each point of the macroscopic continuum are defined as volume averages of their microscopic counterparts over a representative volume element (RVE) of material associated with that point. We consider that the RVE can contain a number of voids, inclusions and/or cracks. It is assumed that non‐penetration conditions are imposed at the crack faces, which do not allow the opposite crack faces to penetrate each other. The derived sensitivity leads to a symmetric fourth‐order tensor field over the unperturbed RVE domain, which measures how the macroscopic elasticity parameters estimated within the multi‐scale framework changes when a small circular inclusion is introduced at the micro‐scale level. Copyright © 2009 John Wiley & Sons, Ltd.     

A Method to Estimate the Inelastic Response in Notched Structures of Anisotropic Materials with Elastic ‐ Perfectly Plastic Behavior

The fatigue strength of engineering structures is often limited by notches which arise from constructive features or manufacturing defects. The constitutive behavior in notch regions is then characterized by small plastic zones which are contained in an elastic region. To avoid costly plastic calculations, approximate methods have been developed to estimate the inelastic stress‐strain response at the notch tip. One of the first and best known approximate models is that of Neuber which, over the last 40 years, has received considerable attention particularly in connection with fatigue life prediction. Numerous studies have been conducted to the verification and the generalization of the Neuber approach which respect to multiaxiality, cyclic loading and creep conditions. Recently an extension of the Neuber method to anisotropic materials has been proposed in [1] and applied to directionally solidified and single crystal Nickel based superalloys as they are used in high temperature material applications. In this short notice we modify the approach in [1] for the special case of an elastic ‐ perfectly plastic anisotropic material.     

Application of a Viscoelastic Material Model in Electro-Mechanics

In this contribution, the numerical modeling of electro-viscoelastic material is considered. The electro-mechanical problem formulated in terms of a symmetrized stress tensor is extended to a viscoelastic material model. For the incorporation of the viscosity model, the logarithmic strain space setting is utilized which mimics the small strain setting. Therefore a rheological model for viscosity from the geometrically linear theory can be used. Numerical examples for a typical uniaxial tensile test show the capability of the method to demonstrate typical relaxation and creep behavior. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Construction of polyconvex,anisotropic free‐energy functions

The existence of minimizers of some variational principles in finite elasticity is based on the concept of quasiconvexity, introduced by Morrey [6]. This integral inequality is rather complicated to handle. Thus, the sufficient condition for quasiconvexity, the polyconvexity condition in the sense of Ball [1], is a more important concept for practical applications, see also Ciarlet [4] and Dacorogna [5]. In the case of isotropy there exist some models which satisfy this condition. Furthermore, there does not exist a systematic treatment of anisotropic, polyconvex free‐energies in the literature. In the present work we discuss some aspects of the formulation of polyconvex, anisotropic free‐energy functions in the framework of the invariant formulation of anisotropic constitutive equations and focus on transverse isotropy.     

Optimal convergence analysis of mixed finite element methods for fourth‐order elliptic and parabolic problems

By using a special interpolation operator developed by Girault and Raviart (finite element methods for Navier‐Stokes Equations, Springer‐Verlag, Berlin, 1986), we prove that optimal error bounds can be obtained for a fourth‐order elliptic problem and a fourth‐order parabolic problem solved by mixed finite element methods on quasi‐uniform rectangular meshes. Optimal convergence is proved for all continuous tensor product elements of order k ≥ 1. A numerical example is provided for solving the fourth‐order elliptic problem using the bilinear element. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Une justification du modèle de coques de Naghdi

We consider a thin linearly elastic loaded shell allowing non-zero inextensional displacements. Under some assumptions on the loads, we prove that the tangential and normal parts of the stress tensor are small compared with the transverse pan, when the thickness of the shell goes to zero. Besides, the displacement vector and the transverse pan of the stress tensor are of the same order of magnitude with respect to the thickness when the material constituting the shell is Isotropic and homogeneous. The limit model, which is a flexural model, can also be obtained from Naghdi's model but not from Koiter's model. In some cases of anisotropic materials, the displacement vector is of a larger order of magnitude than the stress tensor, when the thickness goes to zero.