A finite element model of skeletal muscles

The present paper surveys recent developments in constitutive and computational modelling of skeletal muscles, concerning mainly the generalization to two- and three-dimensional (2D, 3D) continuum deformation analysis of typical one-dimensional (1D) Hill-type muscle models. Extending our previous work in the field and recent contributions by other authors, we describe a constitutive model for skeletal muscles that incorporates all the features of the 3 typical elements (parallel elastic, series elastic and contractile elements) in Hill's muscle model. In particular the proposed incompressible transversely isotropic model incorporates: a multiplicative split of the fibre stretch into contractile and (series) elastic stretches; the possibility of energy storage in the series elastic element; the dependence of the contractile stress on the strain rate; the governing equation of activation dynamics, so that general histories of neural stimulation may be taken as input data. The resulting 2D or 3D constitutive equations are implemented as user subroutines in the large deformation finite element software package ABAQUS. Simple numerical tests are presented and discussed, as well as an example that involves passive or active deformations of a pelvic floor muscle using shell finite elements.     

A finite deformation constitutive model for shape memory polymers based on Hencky strain

In many engineering applications, shape memory polymers (SMPs) usually undergo arbitrary thermomechanical loadings at finite deformation. Thus, development of 3D constitutive models for SMPs within the finite deformation regime has attracted a great deal of interest. In this paper, based on the classical framework of thermodynamics of irreversible processes, employing the logarithmic (or Hencky) strain as a more physical measure of strain, a 3D large-strain macromechanical model is presented. In the constitutive model development, we adopt a multiplicative decomposition of the deformation gradient into elastic and stored parts. In addition, employing the averaging scheme, the logarithmic elastic strain tensor is decomposed into the rubbery and glassy parts. The evolution equations for internal variables are introduced for both cooling and heating processes. The time-discrete form of the proposed model in the implicit form is also presented. Comparing the predicted results with experimental data reported in the literature, the model is validated. Finally, using the finite element method, two boundary value problems e.g., a 3D beam and a medical stent made of SMPs are numerically simulated.     

A finite element model for contact analysis of multiple Cosserat bodies

The objective of this paper is to develop a finite element model for multi-body contact analysis of Cosserat materials. Based on the parametric virtual work principle, a quadratic programming method is developed for finite element analysis of contact problems. The contact problem with friction between two Cosserat bodies is treated in the same way as in plastic analysis. The penalty factors, that are normally introduced into the algorithm for contact analysis, have a direct influence on accuracy of solution. There is no available rule for choosing a reasonable value of these factors for simulation of contact problems of Cosserat materials, and they are therefore cancelled through a special technique so that the numerical results can be of high accuracy. Compared with the conventional work on Cosserat elasticity, the newly developed model is on the contact analysis of the Cosserat materials and is seldom found in the existing literatures. Four examples are computed to illustrate the validity and importance of the model developed.     

Three-dimensional finite element analysis of finite deformation micromorphic linear isotropic elasticity

A finite deformation micromorphic materially linear isotropic elastic model is formulated and implemented for three dimensional finite element analysis. The model is based on the kinematics, balance equations and thermodynamic equations proposed by Eringen and Suhubi (1964). The constitutive equations are calculated in the reference configuration, and the resulting stresses are mapped to the current configuration. The balance of linear momentum and the balance of first moment of momentum are linearized to construct the consistent tangent for three dimensional finite element implementation for solution by the Newton–Raphson method. Three dimensional numerical examples are analyzed to demonstrate preliminarily the implementation.     

Mechanical analysis of 3D braided composites: a finite element model

The analysis of 3D braided composites is more difficult due to its complex microstructure. A new type of finite element method is developed to predict the effective moduli and the local stress within 3D braided composites under the 3D mechanical loading. To verify the present method, the material properties of undamaged 3D braided composites predicted in this paper are compared with the previous work. To demonstrate this method, some examples are analyzed.     

A new approach for the simulation of skeletal muscles using the tool of statistical mechanics

The structure of a skeletal muscle is dominated by its hierarchical architecture in which thousands of muscle fibres are arranged within a connective tissue network. The single muscle fibres consist of many force‐producing cells, known as sarcomeres. These micro biological engines are part of a motor unit and contribute to the contraction of the whole muscle. There are a lot of questions concerning the optimisation of muscle strength and agility. Standard experimental investigations are not sufficient to answer these questions because they do not supply enough information. Additionally, these methods are limited because not enough material for testing is accessible. To overcome these problems, numerical testing tools can be an adequate alternative. From the mechanical point of view the material behaviour of muscles is highly non‐linear. They undergo large deformations in space, thereby changing their shape significantly, so that geometrical nonlinearity has to be considered. Many authors use continuum‐based approaches in combination with the finite element method to describe such material behaviour. However, models of this kind require realistic constitutive relations between stress and strain which are difficult to determine in an inhomogeneous material. Furthermore, biomechanical information cannot be fully exploited in these models. The present approach is crucially based on the use of the finite element method. The material behaviour of the muscle is split into a so‐called active and a passive part. To describe the passive part special unit cells consisting of one tetrahedral element and six truss elements have been derived. Embedded into these unit cells are further truss elements which represent bundles of muscle fibres. Depending on the discretisation it is possible to simulate the material behaviour of e.g. artery walls characterised by oriented fibres or soft tissue including only non‐oriented fibres. In summary, the present concept has the advantage that a three‐dimensional model is developed which allows us take into account many physiological processes at the micro level.     

A cohesive finite element for quasi-continua

In this paper, a cohesive finite element method (FEM) is proposed for a quasi-continuum (QC), i.e. a continuum model that utilizes the information of underlying atomistic microstructures. Most cohesive laws used in conventional cohesive FEMs are based on either empirical or idealized constitutive models that do not accurately reflect the actual lattice structures. The cohesive quasi-continuum finite element method, or cohesive QC-FEM in short, is a step forward in the sense that: (1) the cohesive relation between interface traction and displacement opening is now obtained based on atomistic potentials along the interface, rather than empirical assumptions; (2) it allows the local QC method to simulate certain inhomogeneous deformation patterns. To this end, we introduce an interface or discontinuous Cauchy–Born rule so the interfacial cohesive laws are consistent with the surface separation kinematics as well as the atomistically enriched hyperelasticity of the solid. Therefore, one can simulate inhomogeneous or discontinuous displacement fields by using a simple local QC model. A numerical example of a screw dislocation propagation has been carried out to demonstrate the validity, efficiency, and versatility of the method. An erratum to this article can be found at     

Discrete cohesive zone model for mixed-mode fracture using finite element analysis

The discrete cohesive zone model (DCZM) is implemented using the finite element (FE) method to simulate fracture initiation and subsequent growth when material non-linear effects are significant. Different from the widely used continuum cohesive zone model (CCZM) where the cohesive zone model is implemented within continuum type elements and the cohesive law is applied at each integral point, DCZM uses rod type elements and applies the cohesive law as the rod internal force vs. nodal separation (or rod elongation). These rod elements have the provision of being represented as spring type elements and this is what is considered in the present paper. A series of 1D interface elements was placed between node pairs along the intended fracture path to simulate fracture initiation and growth. Dummy nodes were introduced within the interface element to extract information regarding the mesh size and the crack path orientation. To illustrate the DCZM, three popular fracture test configurations were examined. For pure mode I, the double cantilever beam configuration, using both uniform and biased meshes were analyzed and the results show that the DCZM is not sensitive to the mesh size. Results also show that DCZM is not sensitive to the loading increment, either. Next, the end notched flexure for pure mode II and, the mixed-mode bending were studied to further investigate the approach. No convergence difficulty was encountered during the crack growth analyses. Therefore, the proposed DCZM approach is a simple but promising tool in analyzing very general two-dimensional crack growth problems. This approach has been implemented in the commercial FEA software ABAQUS^® using a user defined subroutine and should be very useful in performing structural integrity analysis of cracked structures by engineers using ABAQUS^®.     

有限元简化模型应用于超声换能器模拟分析

依据有限元方法的基本物理思想,在某些不需要计算辐射声场的准确声学参数和波束特性的工程应用方面,对流体模型进行充分简化,提出了简化模型处理的有效方法,利用该方法对超声换能器进行模拟分析,并进行了样品的制作和测试,实测结果与模型简化分析处理的结果基本一致。可以证明,用该方法进行换能器的优化设计是可行和高效的。     

Simulation of the cabling process for Rutherford cables: An advanced finite element model

In all existing large particle accelerators (Tevatron, HERA, RHIC, LHC) the main superconducting magnets are based on Rutherford cables, which are characterized by having: strands fully transposed with respect to the magnetic field, a significant compaction that assures a large engineering critical current density and a geometry that allows efficient winding of the coils. The Nb₃Sn magnets developed in the framework of the HL-LHC project for improving the luminosity of the Large Hadron Collider (LHC) are also based on Rutherford cables. Due to the characteristics of Nb₃Sn wires, the cabling process has become a crucial step in the magnet manufacturing. During cabling the wires experience large plastic deformations that strongly modify the geometrical dimensions of the sub-elements constituting the superconducting strand. These deformations are particularly severe on the cable edges and can result in a significant reduction of the cable critical current as well as of the Residual Resistivity Ratio (RRR) of the stabilizing copper. In order to understand the main parameters that rule the cabling process and their impact on the cable performance, CERN has developed a 3D Finite Element (FE) model based on the LS-Dyna® software that simulates the whole cabling process. In the paper the model is presented together with a comparison between experimental and numerical results for a copper cable produced at CERN.     

A finite element approach to anisotropic damage of ductile materials in large deformations Part II: finite element formulation and applications

A finite element analysis model for material and geometrical non-linearities due to large plastic deformations of ductile materials is presented using the continuum damage mechanics approach. To overcome limitations of the conventional plastic analysis, a fourth-order tensor damage, defined in Part I of this paper to represent the stiffness degradation in the finite strain regime, is incorporated. General forms of an updated Lagrangian (U.L.) finite element procedure are formulated to solve the governing equations of the coupled elastic–plastic-damage analysis, and a computer program is developed for two-dimensional plane stress/strain problems. A numerical algorithm to treat the anisotropic damage is proposed in addition to the non-linear incremental solution algorithm of the U.L. formulation. Selected examples, compared with published results, show the validity of the presented finite element approach. Finally, the necking phenomenon of a plate with a hole is studied to explore plastic damage in large strain deformations. This revised version was published online in August 2006 with corrections to the Cover Date.     

Modelling cohesive crack growth using a two-step finite element-scaled boundary finite element coupled method

A two-step method, coupling the finite element method (FEM) and the scaled boundary finite element method (SBFEM), is developed in this paper for modelling cohesive crack growth in quasi-brittle normal-sized structures such as concrete beams. In the first step, the crack trajectory is fully automatically predicted by a recently-developed simple remeshing procedure using the SBFEM based on the linear elastic fracture mechanics theory. In the second step, interfacial finite elements with tension-softening constitutive laws are inserted into the crack path to model gradual energy dissipation in the fracture process zone, while the elastic bulk material is modelled by the SBFEM. The resultant nonlinear equation system is solved by a local arc-length controlled solver. Two concrete beams subjected to mode-I and mixed-mode fracture respectively are modelled to validate the proposed method. The numerical results demonstrate that this two-step SBFEM-FEM coupled method can predict both satisfactory crack trajectories and accurate load-displacement relations with a small number of degrees of freedom, even for crack growth problems with strong snap-back phenomenon. The effects of the tensile strength, the mode-I and mode-II fracture energies on the predicted load-displacement relations are also discussed.     

Fully-automatic modelling of cohesive crack growth using a finite element-scaled boundary finite element coupled method

This study develops a method coupling the finite element method (FEM) and the scaled boundary finite element method (SBFEM) for fully-automatic modelling of cohesive crack growth in quasi-brittle materials. The simple linear elastic fracture mechanics (LEFM)-based remeshing procedure developed previously is augmented by inserting nonlinear interface finite elements automatically. The constitutive law of these elements is modelled by the cohesive/fictitious crack model to simulate the fracture process zone, while the elastic bulk material is modelled by the SBFEM. The resultant nonlinear equation system is solved by a local arc-length controlled solver. The crack is assumed to grow when the mode-I stress intensity factor K_I vanishes in the direction determined by LEFM criteria. Other salient algorithms associated with the SBFEM, such as mapping state variables after remeshing and calculating K_I using a “shadow subdomain”, are also described. Two concrete beams subjected to mode-I and mixed-mode fracture respectively are modelled to validate the new method. The results show that this SBFEM-FEM coupled method is capable of fully-automatically predicting both satisfactory crack trajectories and accurate load-displacement relations with a small number of degrees of freedom, even for problems with strong snap-back. Parametric studies were carried out on the crack incremental length, the concrete tensile strength, and the mode-I and mode-II fracture energies. It is found that the K_I ? 0 criterion is objective with respect to the crack incremental length.     

A finite element ductile failure simulation method using stress-modified fracture strain model

This paper proposes a new method to simulate ductile failure using finite element analysis based on the stress-modified fracture strain model. A procedure is given to determine the stress-modified fracture strain as a function of the stress triaxiality from smooth and notched bar tensile tests with FE analyses. For validation, simulated results using the proposed method are compared with experimental data for cracked bar (tensile and bend) tests, extracted from API X65 pipes, and for full-scale burst test of gouged pipes, showing overall good agreements. Advantages in the use of the proposed method for practical structural integrity assessment are discussed.     

Analyses of magneto-electro-elastic plates using a higher order finite element model

In this paper is presented a higher-order model for static and free vibration analyses of magneto-electro-elastic plates, which allows the study of thin and thick plates. The finite element model is a single layer triangular plate/shell element with 24 degrees of freedom for the generalized mechanical displacements. Two degrees of freedom are introduced per each element layer, one corresponding to the electrical potential and the other for magnetic potential. Solutions are obtained for different laminations of the magneto-electro-elastic plate, as well as for the purely elastic plate as a special case. Results are compared with alternative models for static and free vibrations situations.     

An object-oriented finite element framework for multiphysics phase field simulations

The phase field approach is a powerful and popular method for modeling microstructure evolution. In this work, advanced numerical tools are used to create a framework that facilitates rapid model development. This framework, called MARMOT, is based on Idaho National Laboratory’s finite element Multiphysics Object-Oriented Simulation Environment. In MARMOT, the system of phase field partial differential equations (PDEs) are solved simultaneously together with PDEs describing additional physics, such as solid mechanics and heat conduction, using the Jacobian-Free Newton Krylov Method. An object-oriented architecture is created by taking advantage of commonalities in the phase field PDEs to facilitate development of new models with very little effort. In addition, MARMOT provides access to mesh and time step adaptivity, reducing the cost for performing simulations with large disparities in both spatial and temporal scales. In this work, phase separation simulations are used to show the numerical performance of MARMOT. Deformation-induced grain growth and void growth simulations are also included to demonstrate the muliphysics capability.     

A continuum shell finite element model for impact simulation of woven fabrics

A new computational approach is developed to predict the impact behaviour of fabric panels based on the detailed response of the smallest repeating unit (unit cell) in the fabric. The unit cell is constructed and calibrated using measured geometrical (weave architecture, crimp, voids, etc.) and mechanical properties of the fabric. A pre-processor is developed to create a 3D finite element mesh of the unit cell using the measured fabric cross-sectional micro-images. To render an efficient method for simulation of multi-layer packs, these unit cells are replaced with orthotropic shell elements that have similar macroscopic (smeared) mechanical properties as the unit cell. The aim is to capture the essence of the response of a unit cell in a single representative shell element, which would replace the more complicated and numerically costly 3D solid model of the yarns in a crossover. The 3D finite element analysis of the unit cell is used to provide a baseline mechanical response for calibrating the constitutive model in the equivalent shell representation. This shell element takes advantage of a simple physics-based analytical relationship to predict the behaviour of the fabric's warp and weft yarns under general applied displacements in these directions. The analytical model is implemented in the commercial explicit finite element code, LS-DYNA, as a user material routine (UMAT) for shell elements. Layers of fabric constructed from these specialized elements are stacked together to create fabric targets that are then analysed under projectile impact. This approach provides an efficient numerical model for the dynamic analysis of multi-layer fabric structures while taking into account several geometrical and material attributes of the yarns and the fabric.     

A finite element model for liquid phase electroepitaxy

A finite element numerical simulation model for the liquid phase electroepitaxial growth process of gallium arsenide is presented. The basic equations obtained from the fundamental principles of electrodynamics of continua, the constitutive equations for the liquid and solid phases derived from a rational thermodynamic theory, and the associated interface and boundary conditions are presented for a two-dimensional axisymmetric growth cell configuration. The field equations are solved numerically by an adaptive finite element procedure. The effect of moving interfaces is taken into account. Numerical simulations are carried out for different convection levels by changing the value of the gravitational constant. Results show that convection has significant effect on the growth process under normal gravity conditions and results in thickness non-uniformity of the grown layers. The thickness non-uniformity leads to curved interfaces of growth and dissolution, which enhance convection.     

A scaled boundary finite element formulation for dynamic elastoplastic analysis

This study presents the development of the scaled boundary finite element method (SBFEM) to simulate elastoplastic stress wave propagation problems subjected to transient dynamic loadings. Material nonlinearity is considered by first reformulating the SBFEM to obtain an explicit form of shape functions for polygons with an arbitrary number of sides. The material constitutive matrix and the residual stress fields are then determined as analytical polynomial functions in the scaled boundary coordinates through a local least squares fit to evaluate the elastoplastic stiffness matrix and the residual load vector semianalytically. The treatment of the inertial force within the solution of the nonlinear system of equations is also presented within the SBFEM framework. The nonlinear equation system is solved using the unconditionally stable Newmark time integration algorithm. The proposed formulation is validated using several benchmark numerical examples.     

A mixed finite element method for beam and frame problems

 In this work we consider solutions for the Euler-Bernoulli and Timoshenko theories of beams in which material behavior may be elastic or inelastic. The formulation relies on the integration of the local constitutive equation over the beam cross section to develop the relations for beam resultants. For this case we include axial, bending and shear effects. This permits consideration in a direct manner of elastic and inelastic behavior with or without shear deformation. A finite element solution method is presented from a three-field variational form based on an extension of the Hu–Washizu principle to permit inelastic material behavior. The approximation for beams uses equilibrium satisfying axial force and bending moments in each element combined with discontinuous strain approximations. Shear forces are computed as derivative of bending moment and, thus, also satisfy equilibrium. For quasi-static applications no interpolation is needed for the displacement fields, these are merely expressed in terms of nodal values. The development results in a straight forward, variationally consistent formulation which shares all the properties of so-called flexibility methods. Moreover, the approach leads to a shear deformable formulation which is free of locking effects – identical to the behavior of flexibility based elements. The advantages of the approach are illustrated with a few numerical examples. Dedicated to the memory of Prof. Mike Crisfield, for his cheerfulness and cooperation as a colleague and friend over many years.