A mixed finite volume formulation for the solution of gradient elasticity problems

The present works aims at solving the equations of dipolar gradient elasticity via the finite volume method. Initially, a general, mixed, finite volume formulation is stated. As a first approach, the equations of 1D and 2D gradient elasticity are solved via the proposed method. Numerical implementation shows that the suggested model is in excellent agreement with analytic solutions. The approach seems to be promising for extension to 3D problems also.     

A new 3-D finite element for nonlinear elasticity using the theory of a Cosserat point

The theory of a Cosserat point has been used to formulate a new 3-D finite element for the numerical analysis of dynamic problems in nonlinear elasticity. The kinematics of this element are consistent with the standard tri-linear approximation in an eight node brick-element. Specifically, the Cosserat point is characterized by eight director vectors which are determined by balance laws and constitutive equations. For hyperelastic response, the constitutive equations for the director couples are determined by derivatives of a strain energy function. Restrictions are imposed on the strain energy function which ensure that the element satisfies a nonlinear version of the patch test. It is shown that the Cosserat balance laws are in one-to-one correspondence with those obtained using a Bubnov–Galerkin formulation. Nevertheless, there is an essential difference between the two approaches in the procedure for obtaining the strain energy function. Specifically, the Cosserat approach determines the constitutive coefficients for inhomogeneous deformations by comparison with exact solutions or experimental data. In contrast, the Bubnov–Galerkin approach determines these constitutive coefficients by integrating the 3-D strain energy function using the kinematic approximation. It is shown that the resulting Cosserat equations eliminate unphysical locking, and hourglassing in large compression without the need for using assumed enhanced strains or special weighting functions.     

几何精确NURBS有限元中边界条件施加方式对精度影响的三维计算分析

非均匀有理B样条(NURBS)有限元法把计算机辅助几何设计(CAGD)中的NURBS几何构形方法与有限元方法有机结合起来,有效消除了有限元离散模型的几何误差,提高了计算精度。但是由于NURBS基函数不是插值函数,直接在控制节点上施加位移边界条件会引起较大误差。本文详细讨论了NURBS基函数的插值特性,在NURBS有限元分析中采用罚函数法施加位移边界条件,提高了收敛率和计算精度。结合典型三维弹性力学问题,对两种施加位移边界条件的方法进行了对比和分析。计算结果表明,直接施加位移边界条件会导致收敛率和精度的明显降低,而基于罚函数法的NURBS有限元分析则能达到最优收敛率,并具有更高的精度。     

Variational formulation of a simplified strain gradient elasticity theory and its application to a pressurized thick-walled cylinder problem

A detailed variational formulation is provided for a simplified strain gradient elasticity theory by using the principle of minimum total potential energy. This leads to the simultaneous determination of the equilibrium equations and the complete boundary conditions of the theory for the first time. To supplement the stress-based formulation, the coordinate-invariant displacement form of the simplified strain gradient elasticity theory is also derived anew. In view of the lack of a consistent and complete formulation, derivation details are included for the tutorial purpose. It is shown that both the stress and displacement forms of the simplified strain gradient elasticity theory obtained reduce to their counterparts in classical elasticity when the strain gradient effect (a measure of the underlying material microstructure) is not considered. As a direct application of the newly obtained displacement form of the theory, the problem of a pressurized thick-walled cylinder is analytically solved. The solution contains a material length scale parameter and can account for microstructural effects, which is qualitatively different from Lamé’s solution in classical elasticity. In the absence of the strain gradient effect, this strain gradient elasticity solution reduces to Lamé’s solution. The numerical results reveal that microstructural effects can be large and Lamé’s solution may not be accurate for materials exhibiting significant microstructure dependence.     

A formulation of anisotropic continuum elastoplasticity at finite strains. Part I: Modelling

A constitutive model for anisotropic elastoplasticity at finite strains is developed together with its numerical implementation. An anisotropic elastic constitutive law is described in an invariant setting by use of structural tensors and the elastic strain measure C_e. The elastic strain tensor as well as the structural tensors are assumed to be invariant in relation to superimposed rigid body rotations. An anisotropic Hill-type yield criterion, described by a non-symmetric Eshelby-like stress tensor and further structural tensors, is developed, where use is made of representation theorems for functions with non-symmetric arguments. The model also considers non-linear isotropic hardening. Explicit results for the specific case of orthotropic anisotropy are given. The associative flow rule is employed and the features of the inelastic flow rule are discussed in full. It is shown that the classical definition of the plastic material spin is meaningless in conjunction with the present formulation. Instead, the study motivates an alternative definition, which is based on the demand that such a quantity must be dissipation-free, as the plastic material spin is in the case of isotropy. Equivalent spatial formulations are presented too. The full numerical treatment is considered in Part II.     

The conjugate gradient method and block iterative method for penalty finite element of three dimensional navier-stokes equation

A conjugate gradient and block iterative algorithm for element solution of penalty variational form of Navier-Stokes equations are presented. Because the algorithm of solving single variable minimizing problem is simplified, the computing time is greatly saved. In this paper numerical examples are also provided.     

Covariant principal axis formulation of associated coupled thermoplasticity at finite strains and its numerical implementation

In this work we discuss a covariant formulation of the finite strain viscoplasticity in a fully coupled thermomechanical setting. The formulation is presented within the framework of the principal axis methodology, which leads to a very efficient numerical implementation. Several numerical simulations, dealing with fully coupled thermomechanical response at large viscoplastic strains and including both strain localization and cyclic loading cases, are presented in order to illustrate a very satisfying performance of the proposed methodology.     

APPROXIMATION THEORY OF THREE DIMENSIONAL ELASTIC PLATES AND ITS BOUNDARY CONDITIONS WITHOUT USING KIRCHHOFF－LOVE ASSUMPTIONS

ＡＰＰＲＯＸＩＭＡＴＩＯＮＴＨＥＯＲＹＯＦＴＨＲＥＥＤＩＭＥＮＳＩＯＮＡＬＥＬＡＳＴＩＣＰＬＡＴＥＳＡＮＤＩＴＳＢＯＵＮＤＡＲＹＣＯＮＤＩＴＩＯＮＳＷＩＴＨＯＵＴＵＳＩＮＧＫＩＲＣＨＨＯＦＦ－ＬＯＶＥＡＳＳＵＭＰＴＩＯＮＳＣｈｉｅｎＷｅｉ－ｚａｎｇ（...     

A homogenization theory of strain gradient single crystal plasticity and its finite element discretization

In this study, a homogenization theory based on the Gurtin strain gradient formulation and its finite element discretization are developed for investigating the size effects on macroscopic responses of periodic materials. To derive the homogenization equations consisting of the relation of macroscopic stress, the weak form of stress balance, and the weak form of microforce balance, the Y-periodicity is used as additional, as well as standard, boundary conditions at the boundary of a unit cell. Then, by applying a tangent modulus method, a set of finite element equations is obtained from the homogenization equations. The computational stability and efficiency of this finite element discretization are verified by analyzing a model composite. Furthermore, a model polycrystal is analyzed for investigating the grain size dependence of polycrystal plasticity. In this analysis, the micro-clamped, micro-free, and defect-free conditions are considered as the additional boundary conditions at grain boundaries, and their effects are discussed.     

Equations of the linear theory of elasticity of anisotropic materials,reduced to three independent wave equations

Novosibirsk. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, No. 6, pp. 143–150, November–December, 1994.     

用差分法与超松弛迭代法求高维平稳FPK方程的解

用不同精度的差分格式将高维平稳FPK方程离散化为线性代数方程组,然后用超松弛迭代法求解该线性代数方程组得到平稳FPK方程的近似解。讨论了不同的差分格式、网格密度及超松弛因子对解精度及收敛速度的影响,并与其他方法的计算精度进行比较,提出用多重网格算法提高计算效率。研究了典型的二维及四维随机系统的稳态响应,算例表明,该算法具有简洁、节省存储量且精度高的特点,是求解高维平稳FPK方程解的有效算法。     

A viscoplasticity theory applied to proportional and non-proportional cyclic loading at small strains

The small strain, isotropic, Viscoplasticity theory Based on Overstress (VBO) is modified so that the experimentally observed, complex cyclic hardening behavior under proportional and non-proportional loading with fixed as well as variable strain amplitude can be reproduced. This is accomplished by formulating a growth law for the isotropic, rate independent stress, a scalar valued state variable of VBO. The non-proportionality measure employed for modeling the effects of loading path, amplitude and prior history is Tanaka's fourth order tensor [Eur. J. Mech. Solids 13 (1994) 155]. Numerical experiments show the responses in step-up and down two amplitude tests and include further hardening after an increase in the strain amplitude. The differences in the responses to proportional and non-proportional loading including circular, square, one step and two steps are demonstrated. History dependence of the hardening is depicted. The cross hardening behavior, sudden increase in the stress level which is followed by softening after the path change are modeled well.     

The second order approximation theory of three dimensional elastic plates and its boundary conditions without using kirchhoff-love assumptions

The first order approximation theory of three dimensional elastic plates and its boundary conditions presented in the previous paper^[1] establishes six differential equations for the solutions of six undetermined functions u_o, u_a, A_(o) and S_(2)a defined in the x, y plane. They can be divided into two groups, each constitutes three equations to calculate u_o, S_(2)a and u_a, A_(o) respectively. Their boundary conditions as well as these equations are derived from the stationary conditions of variations of a functional for this problem based on the generalized variational principle. The solutions given by this theory are close to those given by the classical theory of thin plates as the ratio of thickness h to width a is small. For large ratio, say h/a=0.3 a considerable difference arises between the two theories. It has not been made clear that in what range of the ratio such difference is reasonable to give more precise solutions. In order, to solve this problem, we must study the second order approximation theory. In this paper following the previous one, we shall establish the second order approximation theory by applying the, stationary condition of variations of a functional for this problem based on the generalized variational principle, to derive nine differential equations and the relate boundary conditions, which are used to calculate nine undetermined functions u_o u_a, A_(o), A₍₁₎, S_(2)a and S_(3)a. And the range of the validity of the first order approximation theory can be found out by comparing the second order theory with the first order theory and the classical theory. It should be pointed out here that the equations of, the second order theory can also be divided into two groups to be solved separately, and the procedure of solution is not too complicate to perform as well. Here, we will use the same notations adopted in the previous paper, and not repeat their definitions.     

General solutions and reduction of a system of equations of the linear theory of elasticity to diagonal form

Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, No. 5, pp. 112–122, September–October, 1993.     

Thermodynamic and rate variational formulation of models for inhomogeneous gradient materials with microstructure and application to phase field modeling

In this work,thermodynamic models for the energetics and kinetics of inhomogeneous gradient materials with microstructure are formulated in the context of continuum thermodynamics and material theory.For simplicity,attention is restricted to isothermal conditions.The materials of interest here are characterized by(1) first- and secondorder gradients of the deformation field and(2) a kinematic microstructure field and its gradient(e.g.,in the sense of director,micromorphic or Cosserat microstructure).Material inhomogeneity takes the form of multiple phases and chemical constituents,modeled here with the help of corresponding phase fields.Invariance requirements together with the dissipation principle result in the reduced model field and constitutive relations.Special cases of these include the wellknown Cahn-Hilliard and Ginzburg-Landau relations.In the last part of the work,initial boundary value problems for this class of materials are formulated with the help of rate variational methods.