On the eigenvalue and eigentensor problem for a tensor of even rank

In more detail as in [1], we study tensor modules of even orders and the eigenvalue and eigentensor problem for tensors of any even rank. The canonical representation of a tensor in the module ?_2p (Ω) is given. We present several statements and theorems about the eigentensors for tensors of even rank, and for adjoint, normal, Hermitian, and unitary tensors in a module of even order. We note that the eigenvalue and eigentensor problem for the tensor of elastic moduli was studied by Ya. Rykhlevskii in 1983–1984. Earlier, it was studied for tensors of any even rank by I. N. Vekua.     

A discussion about scale invariants for tensor functions

It is found that in some cases the complete and irreducible scale invariants given by Ref.[1] are not independent. There are some implicit functional relations among them. The scale invariants for two different cases are calculated. The first case is an arbitrary second order tensor. The second case includes a symmetric tensor, an antisymmetric tensor and a vector. By using the eigentensor notation it is proved that in the first case there are only six independent scale invariants rather than seven as reported in Ref.[1] and in the second case there are only nine independent scale invariants which are less than that obtained in Ref.[1].     

n the fourth order tensor valued function of the stress in return map algorithm

针对各向同性材料,基于一组相互正交的基张量,建立了一套有 效的相关运算方法. 基张量中的两个分别是归一化的二阶单位张量和偏应力张量,另一个则 使用应力的各向同性二阶张量值函数经过归一化构造所得,三者共主轴. 根据张量函数表示 定理,本构方程和返回映射算法中所涉及到的应力的二阶、四阶张量值函数及其逆都由这组 基所表示. 推演结果表明：这些张量之间的运算,表现为对应系数矩阵之间的简单 关系. 其中,四阶张量求逆归结为对应的3\times3系数矩阵求逆,它对二阶张量的变换 则表现为该矩阵对3times 1列阵的变换. 最后,对这些变换关系应用于返回映 射算法的迭代格式进行了相关讨论.     

On general form of Navier-Stokes equations and implicit factored scheme

A general weak conservative form of Navier-Stokes equations expressed with respect to non-orthogonal Curvilinear coordinates and with primitive variables was obtained by using tensor analysis technique, where the contravariant and covariant velocity components were employed. Compared with the current coordinate transformation method, the established equations are concise and forthright, and they are more convenient to be used for solving problems in body-fitted curvilinear coordinate system. An implicit factored scheme for solving the equations is presented with detailed discussions in this paper. For n-dimensional flow the algorithm requires n-steps and for each step only a block tridiagonal matrix equation needs to be solved. It avoids inverting the matrix for large systems of equations and enhances the speed of arithmetic. In this study, the Beam- Warming’s implicit factored schceme is extended and developed in non-orthogonal curvilinear coordinate system.     

The influence of matrix viscoelasticity on the rheology of polymer blends

We examine the effects of matrix phase viscoelasticity on the rheological modeling of polymer blends with a droplet morphology. Two contravariant, second-rank tensor variables are adopted along with the translational momentum density of the fluid to account for viscoelasticity of the matrix phase and the ellipsoidal droplet shapes. The first microstructural variable is a conformation tensor describing the average extension and orientation of the molecules in the matrix phase. The other microstructural variable is a configuration tensor to account for the average shape and orientation of constant-volume droplets. A Hamiltonian framework of non-equilibrium thermodynamics is then adopted to derive a set of continuum equations for the system variables. This set of equations accounts for local conformational changes of the matrix molecules due to droplet deformation and vice versa. The model is intended for dilute blends of both oblate and prolate droplets, and droplet breakup and coalescence are not taken into account. Only the matrix phase is considered as viscoelastic; i.e., the droplets are assumed to be Newtonian. The model equations are solved for various types of homogeneous deformations, and microstructure/rheology relationships are discussed for transient and steady-state conditions. A comparison with other constrained-volume rheological models and experimental data is made as well.     

Thermomechanics of Swelling Biopolymeric Systems

A two-scale theory for the swelling biopolymeric media is developed. At the microscale, the solid polymeric matrix interacts with the solvent through surface contact. The relaxation processes within the polymeric matrix are incorporated by modeling the solid phase as viscoelastic and the solvent phase as viscous at the mesoscale. We obtain novel equations for the total stress tensor, chemical potential of the solid phase, heat flux and the generalized Darcy's law all at the mesoscale. The constitutive relations are more general than those previously developed for the swelling colloids. The generalized Darcy's law could be used for modeling non-Fickian fluid transport over a wide range of liquid contents. The form of the generalized Fick's law is similar to that obtained in earlier works involving colloids. Using two-variable expansions, thermal gradients are coupled with the strain rate tensor for the solid phase and the deformation rate tensor for the liquid phase. This makes the experimental determination of the material coefficients easier and less ambiguous.     

A Recasting of Anisotropic Poroelasticity in Matrices of Tensor Components

The equations associated with the theory of anisotropic poroelastic materials undergoing small deformations are recast in a matrix notation where the matrices are composed of proper tensor components. Using this notation the compatibility conditions on the components of the strain tensor are expressed in terms of the stress tensor, the pore pressure and the anisotropic elastic coefficients of the medium.     

Positive‐definite q‐families of continuous subcell Darcy‐flux CVD(MPFA) finite‐volume schemes and the mixed finite element method

A new family of locally conservative cell‐centred flux‐continuous schemes is presented for solving the porous media general‐tensor pressure equation. A general geometry‐permeability tensor approximation is introduced that is piecewise constant over the subcells of the control volumes and ensures that the local discrete general tensor is elliptic. A family of control‐volume distributed subcell flux‐continuous schemes are defined in terms of the quadrature parametrization q (Multigrid Methods. Birkhauser: Basel, 1993; Proceedings of the 4th European Conference on the Mathematics of Oil Recovery, Norway, June 1994; Comput. Geosci. 1998; 2 :259–290), where the local position of flux continuity defines the quadrature point and each particular scheme. The subcell tensor approximation ensures that a symmetric positive‐definite (SPD) discretization matrix is obtained for the base member (q=1) of the formulation. The physical‐space schemes are shown to be non‐symmetric for general quadrilateral cells. Conditions for discrete ellipticity of the non‐symmetric schemes are derived with respect to the local symmetric part of the tensor. The relationship with the mixed finite element method is given for both the physical‐space and subcell‐space q‐families of schemes. M‐matrix monotonicity conditions for these schemes are summarized. A numerical convergence study of the schemes shows that while the physical‐space schemes are the most accurate, the subcell tensor approximation reduces solution errors when compared with earlier cell‐wise constant tensor schemes and that subcell tensor approximation using the control‐volume face geometry yields the best SPD scheme results. A particular quadrature point is found to improve numerical convergence of the subcell schemes for the cases tested. Copyright © 2007 John Wiley & Sons, Ltd.     

Effective thermal conductivity of transversely isotropic media with arbitrary oriented ellipsoïdal inhomogeneities

The objective of this work is to investigate the thermal conduction phenomena in transversely isotropic geomaterials or rock-like composites with arbitrary oriented ellipsoïdal inhomogeneities of low aspect ratio. Based on the evaluation of the Green function, we provide here new expressions for the interaction tensor whose knowledge permits to obtain the concentration tensor of the polarization field used itself to evaluate the effective thermal conductivity tensor by homogenization. Some particular cases of the obtained general solution are equally presented, in order to validate the developed formalism. The obtained results are next used to study the effect of matrix anisotropy, pores systems and microstructure-related parameters on the overall effective thermal conductivity in transversely isotropic rocks. A two-step homogenization scheme is developed for the prediction of the initial anisotropy effects and to test the ability of the proposed model in the evaluation of effective thermal conductivity. With the help of an Orientation Distribution Function (ODF) the anisotropy due to the pore systems is also accounted. Numerical applications and comparisons with available experimental data are finally carried out for a partially saturated Opalinus clay and an argillite which are both composed of an argillaceous matrix and multiple solid minerals constituents.     

Overall damage and elastoplastic deformation in fibrous metal matrix composites

A phenomenological constitutive model for fibrous composite materials with a ductile matrix is postulated incorporating damage mechanics with micromechanical behavior. The model is first formulated in an undamaged composite system and then transformed consistently achieved in terms of an overall damage tensor M for the whole composite. In the process of formulating this model, interesting results are obtained demonstrating the necessity of using a non-associated flow rule for plasticity in the damaged composite system together with a Hill's type yield criterion. It is also shown that using a Ziegler-Prager kinematic hardening rule for the ductile matrix leads to a general kinematic hardening rule for the composite that is a combination of a generalized Ziegler-Prager model and a Phillips-type model. Finally, an explicit expression for the elastoplastic stiffness tensor for the damaged composite is obtained.     

Constitutive modeling of complex interfaces based on a differential interfacial energy function

An important theory on the dynamics of complex interfaces is the Doi and Ohta theory where the interfacial contribution to the Cauchy stress tensor is determined from an interfacial conformation tensor. For a uniform deformation field in the Eulerian framework, Doi and Ohta adopted a decoupling approximation to reduce a fourth-order tensor into two second-order tensors and derived a differential equation governing the evolution of the interfacial conformation tensor. In this paper, a different formulation is presented for establishing the Cauchy stress tensor based on a path-independent interfacial energy function. By differentiating this interfacial energy function against a Lagrangian strain tensor, a nearly closed-form solution for the stress tensor was determined, involving only elementary algebraic and matrix operations. From this process, the stress-conformation relation proposed by Doi and Ohta is also confirmed from a thermodynamic perspective. The testing cases with uniaxial elongation and simple shear further showed improved fitting to the analytical or exact solutions.     

Mindlin second-gradient elastic properties from dilute two-phase Cauchy-elastic composites Part II: Higher-order constitutive properties and application cases

Starting from a Cauchy elastic composite with a dilute suspension of randomly distributed inclusions and characterized at first-order by a certain discrepancy tensor (see part I of the present article), it is shown that the equivalent second-gradient Mindlin elastic solid: (i) is positive definite only when the discrepancy tensor is negative defined; (ii) the non-local material symmetries are the same of the discrepancy tensor, and (iii) the non-local effective behaviour is affected by the shape of the RVE, which does not influence the first-order homogenized response. Furthermore, explicit derivations of non-local parameters from heterogeneous Cauchy elastic composites are obtained in the particular cases of: (a) circular cylindrical and spherical isotropic inclusions embedded in an isotropic matrix, (b) n-polygonal cylindrical voids in an isotropic matrix, and (c) circular cylindrical voids in an orthotropic matrix.     

General invariant representations of the constitutive equations for isotropic nonlinearly elastic materials

This paper develops general invariant representations of the constitutive equations for isotropic nonlinearly elastic materials. Different sets of mutually orthogonal unit tensor bases are constructed from the strain argument tensor by using the representation theorem and corresponding irreducible invariants are defined. Their relations and geometrical interpretations are established in three dimensional principal space. It is shown that the constitutive law linking the stress and strain tensors is revealed to be a simple relationship between two vectors in the principal space. Relative to two different sets of the basis tensors, the constitutive equations are transformed according to the transformation rule of vectors. When a potential function is assumed to exist, the vector associated with the stress tensor is expressed in terms of its gradient with respect to the vector associated with the strain tensor. The Hill’s stability condition is shown to be that the scalar product of the increment of those two vectors must be positive. When potential function exists, it becomes to be that the 3 × 3 constitutive matrix derived from its second order derivative with respect to the vector associated with the strain must be positive definite. By decomposing the second order symmetric tensor space into the direct sum of a coaxial tensor subspace and another one orthogonal to it, the closed form representations for the fourth order tangent operator and its inversion are derived in an extremely simple way.     

Rheology of polymer blends with matrix-phase viscoelasticity and a narrow droplet size distribution

A Hamiltonian framework of non-equilibrium thermodynamics is adopted to construct a set of dynamical continuum equations for a polymer blend with matrix viscoelasticity and a narrow droplet size distribution that is assumed to obey a Weibull distribution function. The microstructure of the matrix is described in terms of a conformation tensor. The variable droplet distribution is described in terms of two thermodynamic variables: the droplet shape tensor and the number density of representative droplets. A Hamiltonian functional in terms of the thermodynamic variables is introduced and a set of time evolution equations for the system variables is derived. Sample calculations for homogenous flows and constant droplet distribution are compared with data of a PIB/PDMS blend and a HPC/PDMS blend with high viscoelastic contrast. For the PIB/PDMS blend, satisfactory predictions of the flow curves are obtained. Sample calculations for a blend with variable droplet distribution are performed and the effect of flow on the rheology, droplet morphology, and on the droplet distribution are discussed. It is found that deformation can increase or decrease the dispersity of the droplet morphology for the flows investigated herein.     

Finite Bending of a Rectangular Block of an Elastic Hencky Material

Hencky's elasticity model is an isotropic, finite hyperelastic equation obtained by simply replacing the Cauchy stress tensor and the infinitesimal strain tensor in the classical Hooke's law for isotropic infinitesimal elasticity with the Kirchhoff stress tensor and Hencky's logarithmic strain tensor. A study by Anand in 1979 and 1986 indicates that it is a realistic finite elasticity model that is in good accord with experimental data for a variety of engineering materials for moderate deformations. Most recently, by virtue of well-founded physical grounds and rigorous mathematical procedures it has been demonstrated by these authors that this model may be essential to achieving self-consistent Eulerian rate type theories of finite inelasticity, e.g., the J ₂-flow theory for metal plasticity, etc. Its predictions have been studied for some typical deformation modes, including extension, simple shear and torsion, etc. Here we are concerned with finite bending of a rectangular block. We show that a closed-form solution may be obtained. We present explicit expressions for the bending angle and the bending moment in terms of the maximum or minimum circumferential stretch in a general case of compressible deformations for any assigned stretch normal to the bending plane. In particular, simplified results are derived for the plane strain case and for the case of incompressibility. This revised version was published online in June 2006 with corrections to the Cover Date.     

Surface Impedance Tensor and Green’s Function for Weakly Anisotropic Elastic Materials

Recently Fu and Mielke uncovered a new identity that the surface impedance tensor of any anisotropic elastic material has to satisfy. By solving algebraically a matrix equation that follows from the new identity, we derive an explicit expression for the surface impedance tensor, which is correct up to terms linear in the components of the anisotropic part of the elasticity tensor of the material in question. From the well-known relationship between the surface impedance tensor and the Green’s function for infinite space, we obtain an explicit expression for the Green’s function, which is correct up to terms linear in the components of the anisotropic part of the elasticity tensor.      

求解钢压杆极限承载力的有限条塑性系数增量初应力法

本文提出了有限条塑性系数增量初应力法，用于分析钢压杆的弹塑性稳定极限承载力，该法采用分级加载，用有限条法建立结构的增量平衡方程；在塑性范围，引入截面的塑性系数对弹性刚度进行折减得到结构的弹性刚度矩阵；用修正的Ｎｅｗｔｏ－Ｒａｐｈｓｏｎ方法迭代求解，数值结果表明，本法效率较高，与钢压杆试验结果吻合良好，可以考虑残余应力和载荷偏心的影响，可望实现大型超静定结构的极限载力分析。     

Strain gradient solution for a finite-domain Eshelby-type anti-plane strain inclusion problem

A solution for the finite-domain Eshelby-type inclusion problem of a finite elastic body containing an anti-plane strain inclusion of arbitrary cross-sectional shape prescribed with a uniform eigenstrain and a uniform eigenstrain gradient is derived in a general form using a simplified strain gradient elasticity theory (SSGET). The formulation is facilitated by an extended Betti’s reciprocal theorem and an extended Somigliana’s identity based on the SSGET and suitable for anti-plane strain problems. The disturbed displacement field is obtained in terms of the SSGET-based Green’s function for an infinite anti-plane strain elastic body. The solution reduces to that of the infinite-domain anti-plane strain inclusion problem when the boundary effect is not considered. The problem of a circular cylindrical inclusion embedded concentrically in a finite cylindrical elastic matrix undergoing anti-plane strain deformations is analytically solved by applying the general solution, with the Eshelby tensor and its average over the circular cross section of the inclusion obtained in closed forms. This Eshelby tensor, being dependent on the position, inclusion size, matrix size, and a material length scale parameter, captures the inclusion size and boundary effects, unlike existing ones. It reduces to the classical linear elasticity-based Eshelby tensor for the circular cylindrical inclusion in an infinite matrix if both the strain gradient and boundary effects are suppressed. Numerical results quantitatively show that the inclusion size effect can be quite large when the inclusion is small and that the boundary effect can dominate when the inclusion volume fraction is high. However, the inclusion size effect is diminishing with the increase of the inclusion size, and the boundary effect is vanishing as the inclusion volume fraction becomes sufficiently low.     

The existence of real distinct slownesses in photoelastic media

The continuum mathematical models of photoelasticity rely upon the ability of an optical medium, rendered anisotropic by the presence of a stress field, to sustain two transmitted light waves, with distinct slownesses and polarisations. The purpose of this paper is to examine the mathematical models with a view to determining whether they validly support birefringence over a wide range of material properties and impressed stress fields, and to determine conditions under which birefringence is guaranteed by the mathematical models. The classical models of Maxwell and Neumann and the recent finite deformation model of Smith and Rivlin are treated. In each case the secular equation is given and conditions are derived in the form of inequalities for the existence of distinct real slowness solutions. Methods of direct inspection and function extremisation are employed to verify these inequalities, with the aid of a geometric analogy. It is found that for magnetically isotropic media the distinctness of the three eigenvalues of the dielectric tensor is a necessary and sufficient condition for birefringence for all directions of propagation. For magnetically anisotropic media, birefringence is guaranteed by the distinctness of at least two of the eigenvalues of the matrix product of the dielectric tensor and the specific magnetic reluctance tensor.     

Areal moments of inertia revisited: on the distinction between the principal directions

Three commonly used methods to determine the principal moments of inertia of a plane area and their directions are based on: (i) the stationarity condition for the axial moment of inertia, (ii) the eigenvalue analysis, and (iii) Mohr’s circle. In this paper we provide two new derivations, which are based on: (a) the matrix diagonalization and the invariant tensor properties, and (b) the conjugacy property of the moment of inertia vectors. A new general expression is derived which specifies the principal directions of inertia, as well as the directions of the maximum and minimum product of inertia. A comparative study of the five presented approaches is given, which is of interest from both conceptual and methodological points of view. The connection between the deviatoric part of the moment of inertia tensor and Land’s circle of inertia is also given. The presented analysis applies to any two-by-two symmetric second order tensor.