Configurational Forces in Crystal Plasticity

The surface morphology of micro machined surfaces depends on the heterogeneous microstructure. A crystal plasticity model is used to describe the plastic deformation in cp-titanium with its hcp crystal structure. Therefore the basal and prismatic slip systems are taken into account. Furthermore, self and latent hardening are considered. The rate dependency is motivated by a visco plastic evolution law. The cutting process of cp-titanium is modeled within the concept of configurational forces for a standard dissipative media. This framework is implemented into the finite element method. An example illustrates the effects of the microstructure on plastic deformation and configurational forces. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Material Forces in elasto-plastic Materials

In this contribution the concept of configurational forces, also called material forces, is applied to rate–independent, elasto–plastic materials. The theory of configurational forces is briefly recast. Zones of plastic deformation can be interpreted as distributed inhomogeneities. With this background the theory of configurational forces can be applied in many situations, including plastic zones at crack tips, elastic inclusions in elasto–plastic materials and localized deformation. The numerical evaluation is done with the Finite Element Method. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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)     

Configurational Forces in Cutting Processes of Microstructured Titanium

The micro cutting process of microstructured material is influenced by the heterogeneity of the microstructure. In case of cp-titanium an hcp crystal structure is present. Therefore the material response is described with elastic anisotropy and crystal plasticity with the prismatic and basal slip systems of cp-titanium. Also self and latent hardening are considered. The rate dependency is taken into account by a visco-plastic evolution law. In order to investigate a micro cutting process the concept of configurational forces for standard dissipative media is used. Within the finite element method an example illustrates the effects of the heterogeneity and the grain boundary on the configurational force. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Computational homogenization of materials with small deformation to determine configurational forces

In this work the mechanical boundary condition for the micro problem in a two-scaled homogenization using a FE² approach is discussed. The strain tensor is often used in the literature for small deformation problem to determine the boundary conditions for the boundary value problem on the micro level. This strain tensor based boundary condition gives consistent homogenized mechanical quantities, e.g. stress tensor and elasticity tensor, but the present work points out that it leads to unphysical homogenized configurational forces. Instead, we propose a displacement gradient based boundary condition for the micro problem. Results show that the displacement gradient based boundary condition can give not only the consistent homogenized mechanical quantities but also the appropriate homogenized configurational forces. The interpretation of the displacement gradient based boundary condition is discussed. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Strategies for the Computation of Configurational Forces in Dissipative Media

Configurational forces can be interpreted as driving forces on material inhomogeneities such as crack tips. In dissipative media the total configurational force on an inhomogeneity consists of an elastic contribution and a contribution due to the dissipative processes in the material. For the computation of discrete configurational forces acting at the nodes of a finite element mesh, the elastic and dissipative contributions must be evaluated at integration point level. While the evaluation of the elastic contribution is straightforward, the evaluation of the dissipative part is faced with certain difficulties. This is because gradients of internal variables are necessary in order to compute the dissipative part of the configurational force. For the sake of efficiency, these internal variables are usually treated as local history data at integration point level in finite element (FE) implementations. Thus, the history data needs to be projected to the nodes of the FE mesh in order to compute the gradients by means of shape function interpolations of nodal data as it is standard practice. However, this is a rather cumbersome method which does not easily integrate into standard finite element frameworks. An alternative approach which facilitates the computation of gradients of local history data is investigated in this work. This approach is based on the definition of subelements within the elements of the FE mesh and allows for a straightforward integration of the configurational force computation into standard finite element software. The suitability and the numerical accuracy of different projection approaches and the subelement technique are discussed and analyzed exemplarily within the context of a crystal plasticity model. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Microcrack Evolution in Functionally Graded Material under Dynamic Loading

The damage process of metal-ceramic functionally graded material (FGM) is investigated. The microcrack evolution in a layered structure is analyzed using a numerical simulation of stresses and configurational forces. The modelling of an FGM of alumina ceramic and a metallic phase with gradually changing volume fraction of alumina is performed. A structure of two different layers bonded to a substrate is simulated. The stiffness and density of the three materials are varying. The evolution of configurational forces is simulated. The influence of the crack length on the crack driving force is studied for the case of a stress wave loading. The stress loading is applied in the horizontal direction as a dead load. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On a numerical scheme for curved crack propagation based on configurational forces and maximum dissipation

A numerical scheme is presented to predict crack trajectories in two dimensional components. First a relation between the curvature in mixed–mode crack propagation and the corresponding configurational forces is derived, based on the principle of maximum dissipation. With the help of this, a numerical scheme is presented which is based on a predictor–corrector method using the configurational forces acting on the crack together with their derivatives along real and test paths. With the help of this scheme it is possible to take bigger than usual propagation steps, represented by splines. Essential for this approach is the correct numerical determination of the configurational forces acting on the crack tip. The methods used by other authors are shortly reviewed and an approach valid for arbitrary non–homogenous and non–linear materials with mixed–mode cracks is presented. Numerical examples show, that the method is a able to predict the crack paths in components with holes, stiffeners etc. with good accuracy. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation of domain structures in cracked ferroelectric materials

The domain structure around a crack tip plays a significant role in the fracture behavior of ferroelectrics. A continuum phase field model is used to investigate the microstructure at the crack front. The concept of the Eshelby momentum tensor and configurational forces is then generalized to account for the contributions of the polarization term. Implementation of the generalized configurational force in the Finite Element code enables us to numerically obtain the driving force at the crack tip, which corresponds to the crack-tip energy release rate. Calculations show that additional positive electric fields tend to prohibit crack growth, whereas additional negative electric fields tend to promote crack growth. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On Configurational-Force-Driven Crack Propagation and Adaptive Mesh Refinement in Finite Inelasticity

We introduce a consistent variational framework for inelasticity at finite strains, yielding dual balances in physical and material space as the Euler equations. The formulation is employed for the simultaneous usage of configurational forces as both driving forces for crack propagation as well as h-adaptive mesh refinement. The theoretical basis builds upon a global balance of internal and external power, where the mechanical response is exclusively governed by two scalar functions, the free energy function and a dissipation potential. The resulting variational structure is exploited in the context of fracture mechanics and yields evolution equations for internal variables. In the discrete setting, we present a geometry model fully separated from the finite element mesh structure that represents structural changes of the material configuration due to crack propagation. Advanced meshing algorithms provide an optimal discretization at the crack tip. Local and global criteria are obtained via error estimators based on configurational forces being interpreted as indicators of an energetic misfit due to an insufficient discretization. The numerical handling is decomposed into a staggered algorithm scheme for the dual set of equilibrium equations in material and physical space and efficient mesh generation tools. Exemplary numerical examples are considered to illustrate the method and to underline the effects of inelastic material behaviour in the presented context. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An Approach for Micromechanical Modelling of Damage Mechanisms

In the framework of numerical simulation of damaged materials, softening behaviour represents an important topic. Thereby, the decrease of stiffness is mainly caused by the evolution of microvoids. In contrast to the established phenomenological damage approaches, the explicit consideration of effects on the micro scale can lead to an improved approximation quality. In this work, we discuss an approach to describe microstructural evolution. Based on a two phase micro model representing the macroscopical material behaviour, the structural evolution on the micro scale will be modelled based on configurational forces. Besides some theoretical basics on configurational forces at two phase systems and the definition of suitable evolution laws, we present an application of this approach on void growth process in rubberlike material. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Fracture Mechanics of Rubber Material

This contribution presents an approach based on the material force method which allows to evaluate the fracture sensitivity of rubber material. In order to account for the rate–independent dissipation of filled elastomers subjected to low rates, an endochronic plasticity formulation is introduced. It is shown how configurational forces under consideration of material body forces correspond with the elastic and dissipative energy and how they can be related to energy release rates measured in test configurations. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation of micro-cutting considering finite deformation crystal plasticity

Micro-machining processes on metalic microstructures are influenced by the crystal structure, i. e. the grain orientation. Furthermore, the chip formation underlies large deformations. To perform finite element simulations of micro-cutting processes, a large deformation material model is necessary in order to model the hyperelastic and finite plastic material behaviour. In the case of cp-titanium material with hcp-crystal structure the anisotropic behaviour must be considered by an appropriate set of slip planes and slip directions. In the present work the impact of the grain orientation on the plastic deformation is demonstrated by means of finite element simulations of a finite deformation single slip crystal plasticity model. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Application of Material Forces to Hyperelastodynamics

In the framework of numerical simulation, the analysis of fracture and damage mechanisms play an increasing role. In this context, the configurational force method is an useful tool for investigating structures with the finite element method [1] to predict the material behaviour more precisely. In this paper, we discuss the configurational forces in the context of high strain rate loaded hyperelastic structures [6]. Balance of momentum in the material space and the weak formulation are shown with focus on the explicite finite element method. Finally, an application will be presented based on a detailed treatment of structures subjected to short time loading. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Configurational forces in a multiscale approach to piezoelectric materials

Due to the growing interest in determining the macroscopic material response of inhomogeneous materials, computational methods are becoming increasingly concerned with the application of homogenization techniques. In this work, a two-scale classical homogenization of an electro-mechanically coupled material using a FE²-approach is discussed. We explicitly formulated the homogenized coefficients of the elastic, piezoelectric and dielectric tensors for small strain as well as the homogenized remanent strain and remanent polarization. In the homogenization different representative volume elements (RVEs), which capture the micro-structure of the inhomogeneous material, are used to represent the macroscopic material response. Two different schemes are considered. In the first case, domain wall movement is not allowed, but in the second case the movement of the domain walls is taken into account using thermodynamic considerations. Later this technique is used to determine the macroscopic and microscopic configurational forces on defects [2]. These defect situations include the driving force on a crack tip. The effect of the applied electric field on configurational forces at the crack tip is investigated. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Configurational forces on material defects

In micromechanical applications, the behaviour of material defects has been a great subject of interest. In this paper, the idea of Eshelby's energy momentum tensor is briefly reconsidered with respect to problems in solid mechanics. The derived configurational force balance is used to obtain the thermodynamic driving forces acting on centers of dilatation, dislocations, crack tips and interfaces of misfitting precipitates.     

A phase field model for martensitic transformations with a temperature-dependent separation potential

Metallic materials often exhibit a complex microstructure with varying material properties in the different phases. Of major importance in mechanical engineering is the evolution of the austenitic and martensitic phases in steel. The martensitic transformation can be induced by heat treatment or by plastic surface deformation at low temperatures. A two dimensional elastic phase field model for martensitic transformations considering several martensitic orientation variants to simulate the phase change at the surface is introduced in [1]. However here, only one martensitic orientation variant is considered for the sake of simplicity. The separation potential is temperature dependent. Therefore, the coefficients of the Landau polynomial are identified by results of molecular dynamics (MD) simulations for pure iron [1]. The resulting separation potential is applied to analyse the mean interface velocity with respect to temperature and load. The interface velocity is computed by use of the dissipative part to the configurational forces balance as suggested in [3]. The model is implemented in the finite element code FEAP using standard 4-node elements with bi-linear shape functions. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Configurational forces in a micro-cracked body

A framework to model micro-cracked bodies based on the theory of continua with microstructure is presented: assuming a crack embedded in each material element, a field is considered characterising size and orientation of the crack, in addition to the usual placement in the Euclidian space. To deal with crack evolutions, we call upon configurational forces which, in the present case, occur at the level of the microstructure. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

基于非约束模态的中心刚体-Timoshenko梁动力学建模与分析北大核心CSCD

梁的横向变形会导致梁纵向缩短,建模过程中考虑梁横纵变形二次耦合项则存在动力刚化现象,这说明梁的纵向变形会对模型的广义刚度造成影响.对于做旋转运动的梁结构,旋转运动时还会受到离心力的作用而产生轴向拉力,轴向拉力同样也会引起梁的轴向变形,这种影响对粗短梁更加明显.以大范围运动中心刚体-Timoshenko梁模型为研究对象:首先,运用Timoshenko梁理论以及Hamilton原理建立含离心力的动力学模型;其次,引入非约束模态概念,采用Frobenius方法求解非约束模态振型函数以及固有频率;最后,通过数值仿真探究不同恒定转速时非约束模态与约束模态广义刚度的差异和非约束模态条件下离心力对模型的影响.