A modified basis reduction method for limited kinematic hardening shakedown analysis under complex loads

A novel extension of the basis reduction method for kinematic hardening shakedown problem is presented. Firstly, the basis reduction method is implemented based on the modified Newton–Raphson (N-R) method. Then a new technique for the construction of back stress field is introduced, where the simultaneous influence of multiple load corners in shakedown is taken into consideration. Finally, two typical numerical examples are investigated. The results compared with previous works in literatures demonstrated that the proposed method is accurate and the performance in reducing of computation time is significant.     

Thermodynamic based model for the evolution equation of the backstress in cyclic plasticity

A nonlinear kinematic hardening rule is developed here within the framework of thermodynamic principles. The derived kinematic hardening evolution equation has three distinct terms: two strain hardening terms and a dynamic recovery term that operates at all times. The proposed hardening rule, which is referred in this paper as the FAPC (Fredrick and Armstrong–Phillips–Chaboche) kinematic hardening rule, shows a combined form of the Frederick and Armstrong backstress evolution equation, Phillips evolution equation, and Chaboche series rule. A new term is incorporated into the Frederick and Armstrong evolution equation that appears to have agreement with the experimental observations that show the motion of the center of the yield surface in the stress space is directed between the gradient to the surface at the stress point and the stress rate direction at that point. The model is further modified in order to simulate nonproportional cyclic hardening by proposing a measure representing the degree of nonproportionality of loading. This measure represents the topology of the incremental stress path. Numerically, it represents the angle between the current stress increment and the previous stress increment, which is interpreted through the material constants of the kinematic hardening evolution equation. This new kinematic hardening rule is incorporated in a material constitutive model based on the von Mises plasticity type and the Chaboche isotropic hardening type. Numerical integration of the incremental elasto-plastic constitutive equations is based on a simple semi-implicit return-mapping algorithm and the full Newton–Raphson iterative method is used to solve the resulting nonlinear equations. Experimental simulations are conducted for proportional and non-proportional cyclic loadings. The model shows good correlation with the experimental results.     

Asymptotic integration of elasto-plastic constitutive relation for isotropic hardening material and its computational precision

In this paper, an asymptotic expansion solution of the constitutive equation of hardening materials is presented. Its 1st asymptotic integration can give an approximate one with good enough accuracy and the second asymptotic one improves the precision of solutions further. The steps of its algorithms are fairly simple and clear, and its computational workload is considerably reduced. It can be easily incorporated into a general purpose finite element program.The Chinese original of this article was published in the Chinese edition ofActa Mechanica Solida Sinica, No. 1, 1986.     

A combined isotropic-kinematic hardening model for the simulation of warm forming and subsequent loading at room temperature

A coupled isotropic-kinematic hardening material model was developed based on phenomenological observations of performed two stage experiments on a medium carbon steel – SAE 1144, where the first deformation is performed at elevated temperatures and the second deformation at room temperature. Above all, deformations with orthogonal loading at various temperatures were investigated in order to determine the influence of the loading direction as well as of the temperature. Bergström’s theory of work hardening as well as the nonlinear kinematic hardening of an Armstrong–Frederick type were used as a basis for the model development. In the proposed model a relationship between material coefficients of the classical Bergström model and temperature was investigated. The aim of the new material model was to introduce the least possible amount of new parameters as well as to facilitate the mathematical determination of parameters during the fitting of the model with experimental data. The developed model was implemented in an in-house FE-Code in order to simulate the material behavior due to the dynamic strain aging and the hardening behavior after the dynamic strain aging process. Representative simulation results were compared with the experimental data in order to validate the efficiency and the application range of the model.     

循环硬化材料本构模型的隐式应力积分和有限元实现

针对新发展的、能够描述循环硬化行为应变幅值依赖性的粘塑性本构模型,讨论了它的数值实现方法。首先,为了能够对材料的循环棘轮行为（Ratcheting）和循环应力松弛现象进行描述,对已有的本构模型进行了改进;然后,在改进模型的基础上,建立了一个新的、全隐式应力积分算法,进而推导了相应的一致切线刚度（Consistent Tangent Modulus）矩阵的表达式;最后,通过ABAQUS用户材料子程序UMAT将上述本构模型进行了有限元实现,并通过一些算例对一些构件的循环变形行为进行了有限元数值模拟,讨论了该类本构模型有限元实现的必要性和合理性。     

Viscoplastic behaviour of a medium density polyethylene (MDPE): Constitutive equations based on double nonlinear deformation model

The viscoplastic behaviour of a medium density ethylene–butene copolymer (MDPE) is investigated by using samples cut out from thick-walled MDPE pipe. Extensive experimental work has been performed to characterise the nonlinear time-dependent response of such semi-crystalline thermoplastic material. Tests were carried out at 60 °C, on smooth tensile, full axisymmetrically notched creep tensile (FNCT) and double edge notched tensile (DENT) specimens.     

基于精细积分技术的非线性动力学方程的同伦摄动法

将精细积分技术（PIM）和同伦摄动方法（HPM）相结合,给出了一种求解非线性动力学方程的新的渐近数值方法。采用精细积分法求解非线性问题时,需要将非线性项对时间参数按Taylor级数展开,在展开项少时,计算精度对时间步长敏感;随着展开项的增加,计算格式会变得越来越复杂。采用同伦摄动法,则具有相对筒单的计算格式,但计算精度较差,应用范围也限于低维非线性微分方程。将这两种方法相结合得到的新的渐近数值方法则同时具备了两者的优点,既使同伦摄动方法的应用范围推广到高维非线性动力学方程的求解,又使精细积分方法在求解非线性问题时具有较简单的计算格式。数值算例表明,该方法具有较高的数值精度和计算效率。     

A finite strain kinematic hardening constitutive model based on Hencky strain: General framework, solution algorithm and application to shape memory alloys

The logarithmic or Hencky strain measure is a favored measure of strain due to its remarkable properties in large deformation problems. Compared with other strain measures, e.g., the commonly used Green-Lagrange measure, logarithmic strain is a more physical measure of strain. In this paper, we present a Hencky-based phenomenological finite strain kinematic hardening, non-associated constitutive model, developed within the framework of irreversible thermodynamics with internal variables. The derivation is based on the multiplicative decomposition of the deformation gradient into elastic and inelastic parts, and on the use of the isotropic property of the Helmholtz strain energy function. We also use the fact that the corotational rate of the Eulerian Hencky strain associated with the so-called logarithmic spin is equal to the strain rate tensor (symmetric part of the velocity gradient tensor). Satisfying the second law of thermodynamics in the Clausius-Duhem inequality form, we derive a thermodynamically-consistent constitutive model in a Lagrangian form. In comparison with the available finite strain models in which the unsymmetric Mandel stress appears in the equations, the proposed constitutive model includes only symmetric variables. Introducing a logarithmic mapping, we also present an appropriate form of the proposed constitutive equations in the time-discrete frame. We then apply the developed constitutive model to shape memory alloys and propose a well-defined, non-singular definition for model variables. In addition, we present a nucleation-completion condition in constructing the solution algorithm. We finally solve several boundary value problems to demonstrate the proposed model features as well as the numerical counterpart capabilities.     

A combined parametric quadratic programming and precise integration method based dynamic analysis of elastic-plastic hardening/softening problems

The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems. The project supported by the National Key Basic Research Special Foundation (G1999032805), the National Natural Science Foundation of China (19872016, 50178016, 19832010) and the Foundation for University Key Teacher by the Ministry of Education of China     

Analysis of the electromechanical behavior of ferroelectric ceramics based on a nonlinear finite element model

A nonlinear finite element (FE) model based on domain switching was proposed to study the electromechanical behavior of ferroelectric ceramics. The incremental FE formulation was improved to avoid any calculation instability. The problems of mesh sensitivity and convergence, and the efficiency of the proposed nonlinear FE technique have been assessed to illustrate the versatility and potential accuracy of the said technique. The nonlinear electromechanical behavior, such as the hysteresis loops and butterfly curves, of ferroelectric ceramics subjected to both a uniform electric field and a point electric potential has been studied numerically. The results obtained are in good agreement with those of the corresponding theoretical and experimental analyses. Furthermore, the electromechanical coupling fields near (a) the boundary of a circular hole, (b) the boundary of an elliptic hole and (c) the tip of a crack, have been analyzed using the proposed nonlinear finite element method (FEM). The proposed nonlinear electromechanically coupled FEM is useful for the analysis of domain switching, deformation and fracture of ferroelectric ceramics.The project supported by the National Natural Science Foundation of China (10025209, 10132010 and 90208002), the Research Grants of the Council of the Hong Kong Special Administrative Region, China (HKU7086/02E) and the Key Grant Project of the Chinese Ministry of Education (0306)     

On a new extended finite element method for dislocations: Core enrichment and nonlinear formulation

A recently developed finite element method for the modeling of dislocations is improved by adding enrichments in the neighborhood of the dislocation core. In this method, the dislocation is modeled by a line or surface of discontinuity in two or three dimensions. The method is applicable to nonlinear and anisotropic materials, large deformations, and complicated geometries. Two separate enrichments are considered: a discontinuous jump enrichment and a singular enrichment based on the closed-form, infinite-domain solutions for the dislocation core. Several examples are presented for dislocations constrained in layered materials in 2D and 3D to illustrate the applicability of the method to interface problems.     

基于平均值定理和点积分方案的自然单元法及其程序实现

自然单元法是一种基于自然邻接点插值求解偏微分方程的无网格数值方法.它使用Voronoi图或Delaunay三角形作为背景积分网格,使用几何测度构造插值点形函数并形成刚度矩阵.平均值定理定义在未知函数定义域内任何球心(或圆心)的值等于球面(或圆周)上值的平均或加权平均,对于未知函数所满足的平衡方程是充分必要的.因此用平均值定理和点积分方案将求解域内平均应变值由散度定理转化为区域周界上的环路积分,改进传统的积分格式.算例表明,这一积分方案能进一步精简计算量和提高计算效率,是一种自适应的数值计算方法.     

一种求解大型结构动力响应的快速计算方法

针对大型工程结构动力响应求解效率较低的问题,提出了一种基于减基法的快速求解方法。该方法对动力学方程进行时间域积分构建减基空间,利用Galerkin映射向减基空间进行投影得到减缩方程,利用减缩方程快速求得原系统的逼近解,从而大大提高了动力学方程的求解效率。该方法还考虑了Galerkin映射下的奇异性计算,文中算例验证了该...     

Size-dependent effect on biaxial and shear nonlinear buckling analysis of nonlocal isotropic and orthotropic micro-plate based on surface stress and modified couple stress theories using differential quadrature method

The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress,the modified couple stress theory(MCST),and the nonlocal elasticity theories using the differential quadrature method(DQM)is presented.Main advantages of the MCST over the classical theory(CT)are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter.Based on the nonlinear von K′arm′an assumption,the governing equations of equilibrium for the micro-classical plate considering midplane displacements are derived based on the minimum principle of potential energy.Using the DQM,the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained.Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature.A parametric study is conducted to show the effects of the aspect ratio,the side-to-thickness ratio,Eringen’s nonlocal parameter,the material length scale parameter,Young’s modulus of the surface layer,the surface residual stress,the polymer matrix coefficients,and various boundary conditions on the dimensionless uniaxial,biaxial,and shear critical buckling loads.The results indicate that the critical buckling loads are strongly sensitive to Eringen’s nonlocal parameter,the material length scale parameter,and the surface residual stress effects,while the effect of Young’s modulus of the surface layer on the critical buckling load is negligible.Also,considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate.The results show that the critical biaxial buckling load increases with an increase in G12/E2and vice versa for E1/E2.It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude.Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios,it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.     

An unstructured finite element model for incompressible two-phase flow based on a monolithic conservative level set method

We present a robust numerical method for solving incompressible, immiscible two-phase flows. The method extends both a monolithic phase conservative level set method with embedded redistancing and a semi-implicit high-order projection scheme for variable-density flows. The level set method can be initialized conveniently via a simple phase indicator field instead of a signed distance function (SDF). To process the indicator field into a SDF, we propose a new partial differential equation-based redistancing method. We also improve the monolithic level set scheme to provide more accuracy and robustness in full two-phase flow simulations. Specifically, we perform an extra step to ensure convergence to the signed distance level set function and simplify other aspects of the original scheme. Lastly, we introduce consistent artificial viscosity to stabilize the momentum equations in the context of the projection scheme. This stabilization is algebraic, has no tunable parameters and is suitable for unstructured meshes and arbitrary refinement levels. The overall methodology includes few numerical tuning parameters; however, for the wide range of problems that we solve, we identify only one parameter that strongly affects performance of the computational model and provide a value that provides accurate results across all the benchmarks presented. This methodology results in a robust, accurate, and efficient two-phase flow model, which is mass- and volume-conserving on unstructured meshes and has low user input requirements, making it attractive for real-world applications.     

Efficient solution algorithms for thermomechanical contact problems based on a domain/boundary decomposition method

An efficient domain/boundary decomposition method is presented for fully coupled thermomechanical problems with contact boundaries. The whole domain is regarded as a union of subdomains, an interface, and contact interfaces. Penalized variational formulations are performed to connect the interface or contact interfaces with the neighboring subdomains that satisfy continuity constraints on the displacement and temperature fields. As a result, non-linear finite element computations due to the contact boundaries can be localized within a few subdomains or contact interfaces. Therefore, the computational efficiency can be enhanced considerably by devising suitable solution algorithms. A variety of numerical examples were tested to confirm the important features of the new algorithms presented.     

Development of a three-dimensional solver based on the local domain-free discretization and immersed boundary method and its application for incompressible flow problems

Recently, the author and two other coauthors have proposed a two-dimensional hybrid local domain-free discretization and immersed boundary method (LDFD-IBM), which can be used to solve the flow problem with complex geometries. In this paper, the LDFD-IBM is extended to solve a three-dimensional unsteady incompressible flow with the complex computational domain. The technical issues related to the implementation of the LDFD-IBM in three-dimensional problems are discussed in detail, particularly for the discretization of Navier-Stokes equations, mesh strategies for a three-dimensional flow, and the fast algorithm on the identification of the status of mesh nodes (ie, to identify if the mesh node is located in the solid domain, in the fluid domain, or near the immersed boundary). Numerical tests show that the LDFD-IBM can accurately solve three-dimensional incompressible problems with ease.     

Diagonal form fast multipole boundary element method for 2D acoustic problems based on Burton-Miller boundary integral equation formulation and its applications

This paper describes formulation and implementation of the fast multipole boundary element method (FMBEM) for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.