Nonlinear,finite deformation,finite element analysis

The roles of the consistent Jacobian matrix and the material tangent moduli, which are used in nonlinear incremental finite deformation mechanics problems solved using the finite element method, are emphasized in this paper, and demonstrated using the commercial software ABAQUS standard. In doing so, the necessity for correctly employing user material subroutines to solve nonlinear problems involving large deformation and/or large rotation is clarified. Starting with the rate form of the principle of virtual work, the derivations of the material tangent moduli, the consistent Jacobian matrix, the stress/strain measures, and the objective stress rates are discussed and clarified. The difference between the consistent Jacobian matrix (which, in the ABAQUS UMAT user material subroutine is referred to as DDSDDE) and the material tangent moduli (C^e) needed for the stress update is pointed out and emphasized in this paper. While the former is derived based on the Jaumann rate of the Kirchhoff stress, the latter is derived using the Jaumann rate of the Cauchy stress. Understanding the difference between these two objective stress rates is crucial for correctly implementing a constitutive model, especially a rate form constitutive relation, and for ensuring fast convergence. Specifically, the implementation requires the stresses to be updated correctly. For this, the strains must be computed directly from the deformation gradient and corresponding strain measure (for a total form model). Alternatively, the material tangent moduli derived from the corresponding Jaumann rate of the Cauchy stress of the constitutive relation (for a rate form model) should be used. Given that this requirement is satisfied, the consistent Jacobian matrix only influences the rate of convergence. Its derivation should be based on the Jaumann rate of the Kirchhoff stress to ensure fast convergence; however, the use of a different objective stress rate may also be possible. The error associated with energy conservation and work-conjugacy due to the use of the Jaumann objective stress rate in ABAQUS nonlinear incremental analysis is viewed as a consequence of the implementation of a constitutive model that violates these requirements.     

Exploitation of Configurational Forces in Extended Nonlocal Continua with Microstructure

We investigate aspects of the application of configurational forces in extended nonlocal continua with microstructure. Focussing on multifield approaches to gradient–type inelastic solids, the coupled problem is governed by the macroscopic deformation field, while nonlocal inelastic effects on the microstructure are described by a family of order parameter fields. The dual macro– and micro–field equations are derived within an incremental variational framework. Using an incremental principle, due to the variation with respect to the material position, an additional balance in the material space appears with the dual macro–micro–balances in the physical space. In view of the numerical implementation of this coupled problem by a finite element method, the incremental variational framework is recast into a discrete format in terms of discrete macro– and micro–physical nodal forces and configurational nodal forces. Applying a staggered solution scheme, the configurational branch is used as a postprocessing procedure with all the ingredients known from the solution of the coupled macro–micro–problem. The procedure is implemented for a nonlocal, viscous damage model. The consequences with regard to the configurational nodal forces are assessed by means of a numerical example. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling of bainitic phase transformation

In this work, we present a macroscopic material model for simulation of austenite to bainite and of austenite to martensite transformations accompanied by transformation-induced plasticity (TRIP), which is an important phenomenon in metal forming processes. In order to account for the incubation time the model considers nucleation of the bainite phase. When this quantity attains a barrier term, growth of bainite volume fraction is started. The model formulation allows for individual evolutions of upper and lower bainite. Furthermore, the numerical implementation of the constitutive equations into a finite element program is described. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Identification of Finite-Motion Parameters of Serial Chains via Measurement of 6 × 6 Stiffness Matrices

This paper describes a procedure for identifying geometric and stiffness parameters of a mechanical serial chain of know structure by measuring spatial 6×6 stiffness matrices at different positions. The method uses standard optimization routines to determine model parameters such that the model stiffness matrix features in the Frobenius norm the closest distance possible to the measured matrix. From this local identification, a rough model of parameters of finite-motion is created, from which new measuring positions are guessed. By applying this step repeatedly, a model for finite-displacement parameters can be obtained by a sequence of small force-displacement tests. The method is tested with a dummy device consisting of a revolute joint connecting two rigid links dressed with soft material to mimic for example muscle masses of a surrogate mechanism for the elbow joint of a human arm. Two robots grasping the upper and lower arm generate the motion while the force measurement is carried out by a six-axis force sensor. This makes the method potentially suitable for detecting anatomical parameters by in-vivo measurements. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Efficient crack growth analyzes by combining fast methods for BEM

The combination of fast methods for the boundary element method (BEM) for efficient crack growth analyzes is presented. Due to the nonlinearity of fatigue crack growth an incremental procedure has to be applied. Within each increment a stress analysis is needed. Based on the asymptotic stress field the stress intensity factors (SIFs) are calculated by an extrapolation method. Then, a new crack front is determined by a reliable 3D crack growth criterion. Finally, the numerical model has to be updated for the next increment. The time dominant factor in each increment is the computation of the stress field. Due to the stress concentration problem the BEM is utilized. To speed-up the calculation several independent fast methods are exploited. An algebraic technique is the adaptive cross approximation (ACA) method which is acting on the system matrix itself. The application of the substructure technique leads to a blockwise band matrix and therefore to reduced memory requirements. Further savings in memory and computation time are reached by modelling cracks with the dual discontinuity method (DDM) and using the ACA method in each substructure. The efficiency of the combined methods is shown by a complex industrial example. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling the evolution of microstructures in finite plasticity

Kochmann and Hackl introduced in [1] a micromechanical model for finite single crystal plasticity. Based on thermodynamic variational principles this model leads to a non-convex variational problem. Employing the Lagrange functional, an incremental strategy was outlined to model the time-continuous evolution of a first order laminate microstructure. Although this model provides interesting results on the material point level, due to the global minimization in the evolution equations, the calculation time and numerical instabilities may cause problems when applying this model to macroscopic specimens. In order to avoid these problems, a smooth transition zone between the laminates is introduced to avoid global minimization, which makes the numerical calculations cumbersome compared to the model in [1]. By introducing a smooth viscous transition zone, the dissipation potential and its numerical treatment have to be adapted. We obtain rate-dependent time-evolution equations for the internal variables based on variational techniques and show as an example single slip shear. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Optimization the shape of finite element in shot peening process

This paper contains modelling and numerical simulations of shot peening process. The application in Ansys/LS – Dyna programme were elaborated. The phenomena of shot peening process on a typical incremental step were described using a step-by-step incremental procedure, with an updated Lagrangian formulation.Finite elements methods (FEM) and the dynamic explicit method (DEM) were used to obtain the solution. The main purpose of this article is to determinate optimal model of shot peening process (convergence resulting for minimal number of finite elements and their optimal shapes). Examples of calculations were presented. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modelling the mechanical properties of biaxial weft‐knitted fabric reinforced

The subjects of the analysis are glass/epoxy‐composites with biaxial weft knitted fabric reinforcement. This fabric combines the advantages of high in‐plane stiffness and delamination resistance provided by the biaxial and knitted fabric, respectively. The asymptotic homogenisation procedure based on the displacement method in conjunction with a finite element model provides the possibility to obtain the effective material stiffness. The effort to generate a pure 3D‐model of the unit cell would be extremely high. A solution of this problem gives the Binary Model where a simplified geometry of the interior structure is used. Additionally, the Voigt and Reuss bounds for the Young's and shear moduli are found by using the 3D orientation averaging procedure. The results obtained with the Binary Model will be compared to these bounds and the experimental results. The work done for this paper is sponsored by the German Research Community (DFG) in the context of the priority program (SPP) no. 1123. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

识别剪切位点的分组权重矩阵方法

提出识别剪切位点的分组权重矩阵方法．分组权重矩阵是一个子总体数未知的混合权重矩阵模型．为得到未知参数的极大似然估计,采用了EM方法,而子总体数是利用隙统计量确定的．作者给出了该方法的具体计算步骤．将该方法应用于人类和水稻基因组的供子和受子识别,并与现行的主要方法作了比较．结果表明,分组权重矩阵法的识别效果更为精确．     

Identification of optimal material models and parameters in finite strain plasticity

The numerical simulation of the behaviour of a workpiece during manufacturing depends to a large extent on the quality of the applied material model. In this work, a method for the identification of constitutive models and material parameters in engineering applications is proposed. The presented method is used in the setting of optimal experimental design and is based on successive optimization of a set of finite strain plasticity models with kinematic and/or isotropic hardening. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mixed Variational Principles and Robust Finite Element Design of Gradient Plasticity at Finite Strains

The modeling of size effects in elastic-plastic solids, such as the width of shear bands or the grain size dependence in polycrystals, must be based on non-standard theories which incorporate length-scales. This is achieved by models of strain gradient plasticity, incorporating spatial gradients of selected micro-structural fields which describe the evolving dissipative mechanisms. The key aspect of this work is to provide a rigorous incremental variational formulation and mixed finite element design of additive finite gradient plasticity in the logarithmic strain space. We start from a mixed saddle point principle for metric-type plasticity, which is specified for the important model problem of isochoric plasticity with gradient-extended hardening/softening response. To this end, we propose a novel finite element design of the coupled problem incorporating a local-global solution strategy of short- and long-range fields. This includes several new aspects, such as extended Q1P0-type and MINI-type finite elements for gradient plasticity [4]. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A rough set-based incremental approach for learning knowledge in dynamic incomplete information systems

With the rapid growth of data sets nowadays, the object sets in an information system may evolve in time when new information arrives. In order to deal with the missing data and incomplete information in real decision problems, this paper presents a matrix based incremental approach in dynamic incomplete information systems. Three matrices (support matrix, accuracy matrix and coverage matrix) under four different extended relations (tolerance relation, similarity relation, limited tolerance relation and characteristic relation), are introduced to incomplete information systems for inducing knowledge dynamically. An illustration shows the procedure of the proposed method for knowledge updating. Extensive experimental evaluations on nine UCI datasets and a big dataset with millions of records validate the feasibility of our proposed approach.     

Incremental self-consistent approach for the estimation of nonlinear material behavior of metal matrix composites

The application of homogenization methods to compute the macroscopic material response of metal matrix composites is a possibility to save memory and computation time in comparison to full field simulations. This paper deals with a method to extend the self-consistent scheme from linear elasticity theory to nonlinear problems. The idea is to approximate the nonlinear problem by an incrementally linear one. Since time discretization of the deformation process implies a certain linearization, we use the algorithmic consistent tangent operator of the composite for defining the linear comparison material in each time step. This is in contrast to classical incremental self-consistent approaches which use continuum tangent or secant operators. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On physics-based crystal plasticity models: Application to virtual material testing of sheet metals

It is possible to pursue a multi-scale modeling approach for sheet forming simulations by applying the concept of virtual material testing to determine the yield surface from the microstructure of a given material. Full-field simulations with phenomenological crystal plasticity models are widely used for this kind of investigations. However, recent developments focus on incorporating physical quantities like dislocation density into these models. In this work, a dislocation density based crystal plasticity model is used to investigate the plastic anisotropy of the deep drawing steel DC04. In particular, we focus on the prediction of R-values, which can be used to calibrate macroscopic plasticity models. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Bifurcation and chaos in escape equation model by incremental harmonic balancing

Periodic motions of the nonlinear system representing the escape equation with cosine and sine parametric excitations and external harmonic excitations are obtained by the incremental harmonic balance (IHB) method. The system contains quadratic stiffness terms. The Jacobian matrix and the residue vector for the type of nonlinearity with parametric excitation are explicitly derived. An arc length path following procedure is used in combination with Floquet theory to trace the response diagram and to investigate the stability of the periodic solutions. The system undergoes chaotic motion for increase in the amplitude of the harmonic excitation which is investigated by numerical integration and represented in terms of phase planes, Poincaré sections and Lyapunov exponents. The interpolated cell mapping (ICM) method is used to obtain the initial condition map corresponding to two coexisting period 1 motions. The periodic motions and bifurcation points obtained by the IHB method compare very well with results of numerical integration.     

Dynamic demand lot-sizing rules for incremental quantity discounts

Numerous studies have developed and compared lot-sizing procedures for finite-horizon dynamic demand, material requirements planning (MRP), environments when either no purchase discounts exist or for the case of all-units quantity discounts. This paper examines lot-sizing rules when product price schedules follow incremental quantity discounts. The optimal (non-discount) procedure and some traditional heuristic procedures are modified to incorporate incremental quantity discounts. We further modify two heuristics with a ‘look-ahead enhancement’ that performs very well under experimentation. Numerical tests revealed the overall best-performing heuristic in this study to be a modified ‘least-unit cost’ method with a look-ahead enhancement. That procedure produced an average cost penalty vs optimal of 0.26%.     

Incremental unknowns in three-dimensional stationary problem

The main purpose of this work is to set up the explicit matrix framework appropriate to three-dimensional partial differential equations by means of the incremental unknowns method. Multilevel schemes of the incremental unknowns are presented in the three space dimensions, and through numerical experiments, we confirm that the incremental unknowns method is efficient and the hierarchical preconditioning based on the incremental unknowns can be applied in a more general form.      

基于剪切效应纤维梁单元的结构非线性有限元数值模拟

基于Euler-Bernoulli梁理论的经典纤维模型忽略了剪切变形给截面带来的影响,为了得到更加精确的梁单元模型,该文基于考虑剪切效应的纤维梁单元,根据Timoshenko梁理论,推导了该纤维梁单元的刚度矩阵,并结合弹塑性增量理论,同时考虑了几何非线性和材料非线性的双重影响,建立了压弯剪复杂应力状态下结构非线性有限元...     

On convex quadratic programs with linear complementarity constraints

The paper shows that the global resolution of a general convex quadratic program with complementarity constraints (QPCC), possibly infeasible or unbounded, can be accomplished in finite time. The method constructs a minmax mixed integer formulation by introducing finitely many binary variables, one for each complementarity constraint. Based on the primal-dual relationship of a pair of convex quadratic programs and on a logical Benders scheme, an extreme ray/point generation procedure is developed, which relies on valid satisfiability constraints for the integer program. To improve this scheme, we propose a two-stage approach wherein the first stage solves the mixed integer quadratic program with pre-set upper bounds on the complementarity variables, and the second stage solves the program outside this bounded region by the Benders scheme. We report computational results with our method. We also investigate the addition of a penalty term y ^T Dw to the objective function, where y and w are the complementary variables and D is a nonnegative diagonal matrix. The matrix D can be chosen effectively by solving a semidefinite program, ensuring that the objective function remains convex. The addition of the penalty term can often reduce the overall runtime by at least 50 %. We report preliminary computational testing on a QP relaxation method which can be used to obtain better lower bounds from infeasible points; this method could be incorporated into a branching scheme. By combining the penalty method and the QP relaxation method, more than 90 % of the gap can be closed for some QPCC problems.