An anisotropic flow law for incompressible polycrystalline materials

New and explicit anisotropic constitutive equations between the stretching and deviatoric stress tensors for the two- and three-dimensional cases of incompressible polycrystalline materials are presented. The anisotropy is assumed to be driven by an Orientation Distribution Function (ODF). The polycrystal is composed of transversally isotropic crystallites, the lattice orientation of which can be characterized by a single unit vector. The proposed constitutive equations are valid for any frame of reference and for every state of deformation. The basic assumption of this method is that the principle directions of the stretching and of the stress deviator are the same in the isotropic as well as in the anisotropic case. This means that the proposed constitutive laws are able to model the effects of anisotropy only via a change of the fluidity due to a change of the ODF. Such an assumption is justified to guarantee that, besides knowledge of the parameters involved in the isotropic constitutive equation, the anisotropic material response is completely characterized by only one additional parameter, a type of enhancement factor. Explicit comparisons with experimental data are conducted for Ih–ice. Dedicated to Prof. L.W. Morland on the occasion of his 70^th birthday Received: July 6, 2004; revised: November 8, 2004     

The limit equilibrium of an anisotropic medium under the general plasticity condition

A theory of the limit equilibrium of an anisotropic medium under the general plasticity condition in the plane strain state is developed. The proposed yield criterion (the limit equilibrium condition) is obtained by combining the von Mises–Hill yield criterion of an ideally plastic anisotropic material and Prandtl's limit equilibrium condition for a medium under the general plasticity law. It is shown that the problem is statically determinate, i.e., if the boundary conditions are specified in stresses, the stress state in plastic region can only be obtained using equilibrium equations. It is established that the equations describing the stress state are hyperbolic and have two families of characteristic curves that intersect at variable angles. In deriving the equations describing the velocity field, the material is assumed to be rigid plastic, and the associated law of flow is applied. It is shown that the equations for the velocities are also hyperbolic, and their characteristic curves are identical with those of the equations for stresses. However, the directions of the principal values of the stress and strain rate tensors are different due to the anisotropy of the material. The characteristic directions differ from the isotropic case in that the normal and tangential components of the stress tensor do not satisfy the limit conditions. It is established that the equations obtained allow of partial solutions, and in this case, at least one family of characteristic curves consists of straight lines. The conditions along the lines of discontinuity of the velocity are investigated, and it is shown that, as in the isotropic case, these are characteristic curves of the system of governing equations. In the anisotropic formulation, the well-known Rankine problem of the limit state of a ponderable layer is solved. From an analysis of the velocity field it is shown that plastic flow of the entire layer is possible only for a slope angle equal to the angle of internal friction. For slope angles less than the angle of internal friction, the solutions obtained are solutions of problems of the pressure of the medium on the retaining walls. The change in this pressure as a function of the parameters of anisotropy is investigated, and turns out to be significant.     

Polyconvex energy densities with applications to anisotropic thin shells

This work shows the capability to simulate anisotropic thin shells using polyconvex energy densities by analyzing numerical examples. The variational framework is based on the enhanced assumed strain formulation. The iterative enforcement of the zero normal stress condition at the integration points allows consideration of arbitrary three–dimensional constitutive equations. We consider an additive structure of the energy decoupled in an isotropic part for a matrix material and the superposition of transversely isotropic parts for embedded fiber families and focus on polyconvex strain energy functions. In two representative examples we document anisotropy effects. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Robust Schur complement method for strongly anisotropic elliptic equations

We present robust and asymptotically optimal iterative methods for solving 2D anisotropic elliptic equations with strongly jumping coefficients, where the direction of anisotropy may change sharply between adjacent subdomains. The idea of a stable preconditioning for the Schur complement matrix is based on the use of an exotic non‐conformal coarse mesh space and on a special clustering of the edge space components according to the anisotropy behavior. Our method extends the well known BPS interface preconditioner [2] to the case of anisotropic equations. The technique proposed also provides robust solvers for isotropic equations in the presence of degenerate geometries, in particular, in domains composed of thin substructures. Numerical experiments confirm efficiency and robustness of the algorithms for the complicated problems with strongly varying diffusion and anisotropy coefficients as well as for the isotropic diffusion equations in the ‘brick and mortar’ structures involving subdomains with high aspect ratios. Copyright © 1999 John Wiley & Sons, Ltd.     

Assessing the mechanical responses for anisotropic multi-layered medium under harmonic moving load by Spectral Element Method (SEM)

The article is concerned with the mechanical responses of anisotropic multi-layered medium under harmonic moving load. An analytical solution for two-dimensional anisotropic multi-layered medium subjected to harmonic moving load is devoted via Spectral Element Method (SEM), while the anisotropic property is approximated as transverse isotropy. Starting with the constitutive equations of transversely isotropic body and the governing equations of motion based on the loading properties. The analytical spectral elements in the wavenumber domain are obtained according to the principle of wave superposition and Fourier transformation. Then, the spectral global stiffness matrix of the multi-layered medium is derived by assembling the nodded stiffness matrices of all layers depended on the different interlayer conditions between the adjacent layers, i.e. sliding and bonded. The corresponding analytical solutions are achieved by taking the Fourier series and Inverse Fast Fourier Transform (IFFT) algorithm. Finally, some examples are given to validate the accuracy of the proposed analytical solution, and to demonstrate the impact of both anisotropy, top layer thickness, interlayer conditions, and loading properties (velocity and natural frequency) on the mechanical response of the multi-layered medium.     

Modeling of induced anisotropy at large deformations for polymers

In this work we develop a model to describe the induced plasticity of polymers at large deformations. Polymers such as stretch films exhibit a pronounced strength in the loading direction. The undeformed state of the films is isotropic, whereas after the uni-axial loading the material becomes anisotropic. In order to consider this induced anisotropy during the stretch process, a spectral decomposition of the inelastic right CAUCHY-GREEN tensor is done. Therefore, the yield function can be formulated as a function of the anisotropic tensor, where again the anisotropic tensor is a function of the maximum eigenvalue. A backward EULER scheme is used for updating the evolution equations, and the algorithmic tangent operator is derived. The numerical implementation of the resulting set of constitutive equations is used in a finite element program and for parameter identification. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

页岩储层各向异性对水平井坍塌压力的影响

为揭示层理页岩强度各向异性对水平井井壁坍塌压力的影响规律,对取自四川盆地龙马溪组页岩钻井不同层理角度的岩芯,通过岩芯观测、偏光显微镜、扫描电镜试验分别从宏观和微观角度研究了页岩层理特征,并基于单轴抗压强度研究了页岩强度各向异性特征.层理结构发育导致页岩呈现出显著的各向异性特征,粗略地将页岩视为各向同性体及简单采用Mohr-Coulomb准则作为破坏判据,使得预测的维持井壁稳定的坍塌压力不能满足安全钻井的需要.该文在研究层理地层岩石力学特性的基础上,采用横观各向同性地层井周应力分析模型,研究了弹性模量及地应力各向异性比对井周应力的影响;并结合考虑中间主应力的Mogi-Coulomb判据,分析了井壁坍塌压力对弹性模量E和Poisson(泊松)比ν各向异性比值变化的敏感性.实例分析结果表明:弹性模量和水平地应力各向异性比值的变化会对井周应力产生较大影响;Poisson比各向异性对井壁坍塌压力没有明显影响,而弹性模量比值的变化对井壁坍塌压力影响较大;横观各向同性地层中井壁坍塌压力与储层水平和垂直方向的硬脆性强弱关系有关,弹性力学参数各向异性既有可能利于井壁稳定,也可能更易导致井壁的剪切破坏;在实际工程中应根据储层岩石水平和垂直方向的硬脆性强弱关系,确定安全钻井液密度窗口.该结果对现场施工具有很好的指导意义.     

Non-linear radial oscillations of a transversely isotropic hyperelastic incompressible tube

The constitutive equation for a transversely isotropic incompressible hyperelastic material is written in a covariant form for arbitrary orientation of the anisotropic director. Three non-linear differential equations are derived for radial oscillations in radial, tangential and longitudinal transversely isotropic thin-walled cylindrical tubes of generalised Mooney-Rivlin material. A Lie point symmetry analysis is performed. The conditions on the strain-energy function and on the net applied surface pressure for Lie point symmetries to exist are determined. For radial and tangential transversely isotropic tubes the differential equations are reduced to Abel equations of the second kind. Radial oscillations in a longitudinal transversely isotropic tube and in an isotropic tube are described by the Ermakov-Pinney equation.     

准各向同性复合材料弹性和强度的各向异性效应(Ⅰ)──模型*

实验观测表明,准各向同性材料,如N轴纤维增强复合材料层合板和编织材料,其面内刚度和强度具有不同程度的方向性,且强度的各向异性程度往往明显高于弹性性质的各向异性程度。本文根据张量函数表示理论所提出的本构方程和强度准则的一般模型,结合有关实验数据,分析了材料弹性性能和强度的非各向同性效应。具体给出了几类本构模型和强度准则的特殊形式并讨论了本文所得到结果的若干力学性质。本文第Ⅱ部分具体讨论了含单个椭圆孔或裂纹的无限大板的有关强度的各向异性效应,并用细观力学方法检验了本文的模型。     

Nonlinear analysis and elastic-plastic behavior of anisotropic structures

An elastic-plastic analysis of anisotropic work-hardening materials based on a quadratic approximation of the Tsai-Wu criterion is presented. General expressions for the anisotropic parameters in the yield condition are derived for initial and subsequent yielding. Particularly, the plastic constitutive relations are expressed by means of both the flow theory as well as the deformation theory extended to anisotropic plasticity. The numerical algorithms are based upon the notion of a return mapping procedure and a consistent tangent operator valid for anisotropic elastic-plastic materials including work-hardening effects is developed. The solution equations are evaluated by consistent linearization of a nonlinear variational principle and a Newton-Raphson scheme is adopted for the iterative solution of the nonlinear problems. Numerical examples exhibit the reliable performance of the proposed algorithm in some practical calculations. The effects of anisotropy and the differences between flow and deformation theories in the obtained solutions are discussed and compared with available numerical results.     

Instability of solitary waves in non-linear composite media

A problem on the dynamic instability of soliton solutions (solitary waves) of Hamilton's equations, describing plane waves in non-linear elastic composite media with or without anisotropy, is considered. In the anisotropic case, there are two two-parameter families of solitary waves: fast and slow and, when there is no anisotropy, there is one three-parameter family. A classification of the instability of solitary waves of the fast family in the anisotropic case and of representatives of families of solitary waves, the velocities of which lie outside of the range of stability when there is anisotropy and when there is no anisotropy, is given on the basis of a numerical solution of a Cauchy problem for the model equations. In this paper, the fundamental equations describing plane waves in non-linear, anisotropic, elastic composites are derived, the Hamilton form of the basic equations is presented, the symmetries in the anisotropic and isotropic cases are considered, the conserved quantities and the soliton solutions, that is, the solitary waves are presented, the nature of the instability of representatives of all three families is investigated, brief formulation of the results is presented and problems of the instability of the fast family in the anisotropic case and of representatives of the families, the velocities of which lie outside of the range of stability in the presence and absence of anisotropy (explosive instability), are discussed.     

Water reservoir loading of long anisotropic valleys with irregular topographies

The nature of the stress field induced by water reservoir loading of long anisotropic valleys is considered using an analytical method proposed earlier by the authors. The rock mass is modeled as a linearly elastic, transversely isotropic or nearly isotropic and homogeneous continuum that deforms under a condition of plane strain. The method is used to determine stresses in symmetric and asymmetric valleys induced by water reservoir loading with or without gravity. Numerical examples are presented to analyze the disturbance associated with a water reservoir on the existing gravitational stress field. They show that the presence of a water reservoir does not much affect the horizontal and vertical gravity-induced normal stresses except near the valley bottom. There the effect of the water reservoir is to reduce the magnitude of the gravity-induced horizontal tensile stress and to increase the magnitude of the gravity-induced vertical compressive stress. On the other hand water reservoir loading affects the gravity-induced shear stress in a relatively large region below the valley. Finally a parametric study is presented on the effect of

1. (1) water reservoir height,

2. (2) degree of anisotropy, and

3. (3) orientation of anisotropy on the magnitude and distribution of stresses induced by water reservoir loading only.

     

各向异性液体-多孔饱和介质在机械荷载作用下的弹性动力分析

将一个各向异性液体-多孔饱和介质的弹性动力分析,归结为一个横观各向同性液体-多孔饱和介质在机械荷载作用下的变形问题.自然界中有些物理问题,仅在一个方向发生变形,例如,与变形结构和变形柱有关的问题.土力学中,通常假设只有竖向沉降,从而归结为一维多孔弹性模型.采用各向异性液体-多孔饱和介质的一维变形模型,研究了在不同时间和距离下扰动的变化.给出了在不同类型荷载作用下,介质的各向异性对位移分布和应力分布的影响.     

Non-linear laws of fluid flow through anisotropic porous media

Non-linear laws of fluid flow through anisotropic porous media are written out in invariant tensor form for all crystallographic point symmetry groups. The equations, as is customary in seepage theory [1, 2], are represented by expressions containing the seepage velocity up to and including the third degree. Expressions defining non-linear flow resistances are given and it is shown that, when one transfers from linear to non-linear seepage laws, the symmetry group of the flow properties may change. For example, the isotropic flow properties manifested in Darcy's law may become essentially anisotropic in a non-linear law and display asymmetry, that is, they may be different along one straight line in the positive and negative directions. It is shown that, compared with linear seepage laws for anisotropic media, when flow properties may be defined by just four essentially different types of equation, in non-linear laws the appearance of anisotropy is highly diversified and the number of distinct types of equation increases considerably.     

非线性各向同性弹性材料热应力本构方程

在非线性各向同性弹性体张量形式的本构方程基础上,仅考虑温度初值和增值,按照表示定理,补充考虑温度影响的完备项,建立了非线性各向同性弹性材料完备的多项式形式的热应力本构方程和应变能函数．作为应用举例,利用MATLAB软件,将本构方程与现有文献中高温金属材料单向拉伸和压缩情况下弹性阶段的实验数据进行了拟合,结果表明实验值与所提出的理论模型的结果显示了良好的一致性．     

Pseudo-Rigid Bodies: A Geometric Lagrangian Approach

The pseudo-rigid body model is viewed within the context of continuum mechanics and elasticity theory. A Lagrangian reduction, based on variational principles, is developed for both anisotropic and isotropic pseudo-rigid bodies. For isotropic Lagrangians, the reduced equations of motion for the pseudo-rigid body are a system of two (coupled) Lax equations on so(3)×so(3) and a second-order differential equation on the set of diagonal matrices with a positive determinant. Several examples of pseudo-rigid bodies such as stretching bodies, spinning gas cloud and Riemann ellipsoids are presented.     

正交各向异性材料的混合硬化弹塑性本构方程

建立了混合硬化正交各向异性材料的屈服准则,进而推导了与之相关的塑性流动法则.根据简单应力状态的实验曲线,可得到广义等效应力-应变关系.初始屈服曲面与材料的弹性常数有关,材料退化为各向同性且只考虑各向同性硬化时,屈服函数退化为Huber-Mises屈服函数,相关的本构方程退化为Prandtl-Reuss方程.     

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.     

Robust a posteriori error estimates for conforming discretizations of a singularly perturbed reaction-diffusion problem on anisotropic meshes

In this paper, we derive robust a posteriori error estimates for conforming approximations to a singularly perturbed reaction-diffusion problem on anisotropic meshes, since the solution in general exhibits anisotropic features, e.g., strong boundary or interior layers. Based on the anisotropy of the mesh elements, we improve the a posteriori error estimates developed by Cheddadi et al., which are reliable and efficient on isotropic meshes but fail on anisotropic ones. Without the assumption that the mesh is shape-regular, the resulting mesh-dependent error estimator is shown to be reliable, efficient and robust with respect to the reaction coefficient, as long as the anisotropic mesh sufficiently reflects the anisotropy of the solution. We present our results in the framework of the vertex-centered finite volume method but their nature is general for any conforming one, like the piecewise linear finite element one. Our estimates are based on the usual H(div)-conforming, locally conservative flux reconstruction in the lowest-order Raviart-Thomas space on a dual mesh associated with the original anisotropic simplex one. Numerical experiments in 2D confirm that our estimates are reliable, efficient and robust on anisotropic meshes.     

Flow-induced anisotropic viscosity in short fiber reinforced polymers

The commonly used flow models for fiber reinforced polymers often neglect the flow induced mechanical anisotropy of the suspension. With an increasing fiber volume fraction, this plays, however, an important role. There are some models which count on this effect, they are, however, phenomenological and require a fitted model parameter. In this paper, a micromechanically based constitutive law is proposed which considers the flow induced anisotropic viscosity of the fiber suspension. The introduced viscosity tensor can handle arbitrary anisotropy of the fluid-fiber mixture depending on the actual fiber orientation distribution. A homogenization method for unidirectional structures in contribution with orientation averaging is used to determine the effective viscosity tensor. The motion of rigid ellipsoidal fibers induced by the flow of the matrix material is described by Jeffery's equation. A numerical implementation of the introduced model is applied to representative flow modes. The calculated stress values are analyzed in transient and stationary flow cases. They show a less pronounced anisotropic viscous behaviour in every investigated case compared to the results obtained by the use of the Dinh-Armstrong constitutive law. The reason for the qualitative difference is that the presented model depends on the complete orientation information, while the other one is linear in the fourth-order fiber orientation tensor. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)