Adaptive FEM with Stabilized Elements for Parameter Identification of Incompressible Hyperelastic Materials

This work is concerned with the identification of material parameters for isotropic, incompressible hyperelastic material models. Besides the principal stretch-based strain-energy function by Ogden an invariant-based strain-energy function by Rivlin/Saunders is considered for which parameter sensitivities are derived. The identification is formulated as a least-squares minimization problem based on the finite element method to account for inhomogeneous states of stresses and strains in the experimental data which is obtained by optical measurements. For the finite element method low-order tetrahedral elements in a mixed displacement-pressure formulation with stabilization are considered. Special attention is payed to an adaptive mesh-refinement based on a goal-oriented a posteriori error indicator to gain reliable material parameters. To approximate error terms an element-wise recovery technique based on enhanced gradients is introduced. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Quantifying the uncertainty in a hyperelastic soft tissue model with stochastic parameters

We present a simple open-source semi-intrusive computational method to propagate uncertainties through hyperelastic models of soft tissues. The proposed method is up to two orders of magnitude faster than the standard Monte Carlo method. The material model of interest can be altered by adjusting few lines of (FEniCS) code. The method is able to (1) provide the user with statistical confidence intervals on quantities of practical interest, such as the displacement of a tumour or target site in an organ; (2) quantify the sensitivity of the response of the organ to the associated parameters of the material model. We exercise the approach on the determination of a confidence interval on the motion of a target in the brain. We also show that for the boundary conditions under consideration five parameters of the Ogden–Holzapfel-like model have negligible influence on the displacement of the target zone compared to the three most influential parameters. The benchmark problems and all associated data are made available as supplementary material.     

A New Unified Approach for the Simulation of a Wide Class of Directional Distributions

The need for effective simulation methods for directional distributions has grown as they have become components in more sophisticated statistical models. A new acceptance–rejection method is proposed and investigated for the Bingham distribution on the sphere using the angular central Gaussian distribution as an envelope. It is shown that the proposed method has high efficiency and is also straightforward to use. Next, the simulation method is extended to the Fisher and Fisher–Bingham distributions on spheres and related manifolds. Together, these results provide a widely applicable and efficient methodology to simulate many of the standard models in directional data analysis. An R package simdd, available in the online supplementary material, implements these simulation methods.     

Computation of energy dissipation in visco-elastic materials at finite deformation

In this contribution, the derivation of the energy dissipation rate in generalized visco-elastic material models with internal stress-type variables and linear evolution equations is outlined. The approximated dissipation rate is computed from a positive quadratic form of the nonlinear non-equilibrium stresses and the inverse of the consistent material tangent tensor. The presented method is used to compute the energy dissipation of visco-elastic rubber material in a large scale application of a steady state rolling tire structure. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Micromechanical modeling of twinning

There exist material models that incorporate mechanical twinning in a homogenized sense, or consider specific aspects, like grain refinement or texture evolution. However, since the RVE-technique became a standard method, it is possible to obtain more detailed predictions based on micromechanical models. In this work, an approach based on a nonconvex elastic potential and the corresponding results of FE calculations are presented. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A non-parametric fluid-equivalent approach for the acoustic characterization of rigid porous materials

The acoustic characterization of porous materials with rigid solid frame plays a key role in the prediction of the acoustic behavior of any dynamic system that incorporates them. In order to obtain an accurate prediction of its frequency-dependent response, a suitable choice of the parametric models for each material is essential. However, such models could be inadequate for a given material or only valid in a specific frequency range. In this work, a novel non-parametric methodology is proposed for the characterization of the acoustic properties of rigid porous materials. Unlike most widespread methodologies, this technique is based on the solution of a sequence of frequency-by-frequency well-posed inverse problems, thus increasing the characterization accuracy. Once a reduced number of experimental measurements is available, the proposed method avoids the a priori choice of a parametric model.     

Instability prevention at the modeling of large elasto-plastic strains under cyclic loading

The cyclic instability phenomenon is investigated at the modeling of large elasto-plastic strains. The possible causes of the cyclic instability and conditions ensuring cyclical stability of elasto-plastic models are analyzed for the case of large strains. Among the possible causes of the cyclic instability the following are considered: the method of strain decomposition on elastic and plastic parts; the constitutive law for the elastic deformation (hypo- and hyper-elasticity); the constitutive equation for the plastic deformation; the constitutive relation for the plastic spin; kinematic hardening law, in particular, the type of the objective rate in the generalized Prager's law. Predictions of 50 various models of the elasto-plastic material have been compared in order to find the causes of the cyclic instability. Two test problems are considered: cyclic simple shear, combined cyclic simple shear and tension-compression. Results of numerical experiments for the various material models are presented and discussed. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Micromechanical strength analysis of fiber reinforced plastics

Strength of an unidirectional lamina is computed with a representative volume element. An approximative “voxel” meshing method is used in conjunction with continuum damage mechanics to simulate crack growth in the RVE. Two material models for the nonlinear material behaviour of the epoxy resin are compared. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Parallel solution of volume-coupled multi-field problems using an Abaqus-PANDAS software interface

Numerical simulations have proven to be a powerful tool in several engineering disciplines, such as mechanical, civil and biomechanical engineering, and are thus widely used. However, the reliability of the simulations strongly relies on the governing material model. These models are usually developed in academic or industrial research projects and are implemented into dedicated software packages to proof their concepts. A transfer of these models from the research into a production-related environment is often time consuming and prone to failures, and therefore a costly task. The present work introduces a general interface between the research code PANDAS, which is a dedicated multi-field finite-element solver based on a monolithic solution strategy, and the commercial finite-element package Abaqus. The coupling is based on the user-defined element subroutine (UEL) of Abaqus. This procedure, on the one hand, allows for a straight-forward embedding of the PANDAS material models into Abaqus. On the other hand, it provides, in comparison to the native UEL subroutine of Abaqus, a user-friendly programming environment for user-defined material models with an extended number of degrees of freedom. Furthermore, the coupling also supports the parallel-analysis capabilities for large-scale problems on high-performance computing clusters. The Abaqus-PANDAS linkage can be applied to various coupled multi-field problems. However, the present contribution addresses, in particular, volume-coupled multi-field problems as they arise when proceeding from the Theory of Porous Media (TPM) as a modelling framework. For instance, it can be used to model partially or fully saturated soils, or chemically or electro-chemically driven swelling phenomena as they appear, for example, within hydrogels. Additionally, discontinuities, such as cracks, can be described for instance via phase-field models or by the extended finite-element method (XFEM). (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Acoustical Performance of Concreted Wood Fiber Materials

In this paper, an inverse acoustical characterization for determination of porosity, flow resistivity, tortuosity and other parameters of a concreted wood material was proposed. Sound absorption measurements for a concreted wood fiber material were made in an impedance tube. Four acoustical models from literature were briefly described. Each model was fitted to the experimental data in order to compute their unknown parameters. Using a least-square method, the experimental data (sound absorption coefficient) was used to optimize the acoustical models by best-fitting the unknown parameters. The results of this paper show that the presented technique leads to reliable estimates of the physical parameters of the material. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation du drapé des tissus par maillages adaptatifsAdaptive meshing for cloth draping

Numerical simulation methods regarding fabric and cloth draping are generally based on mechanical models. These models are usually applied to uniform grids representing the true geometry of the fabric. Fabrics being a very flexible material, wrinkles appear on its surface when submitted to free or constrained motion. The main problem of the simulation is to represent realistically cloth surface motion. This is strongly dependent on the surface discretization. We present a new cloth animation scheme based on adaptive surface discretization. It can be seen as a multi-grid method which allows us to obtain realistic simulations. We propose also a new mechanical model well suited to our adaptive meshing strategy. A numerical example is given to show the efficiency of the method. To cite this article: J. Villard et al., C. R. Acad. Sci. Paris, Ser. I 335 (2002) 561–566.     

Gradient‐based optimization of a viscoelastic Ogden type model for cellular polymers

The parameter identification is the interface between a theoretical material model and its application in numerical computations. Only by an accurate identification of the theoretically introduced material parameters, an applicable simulation of the material is achieved. An increasing standard of the parameter identification is set by the requirements of complex material models used in computer‐aided engineering. A common identification strategy is a gradient‐based optimization of a least‐squares functional, e. g. the sequential quadratic programming (SQP) technique. In this paper, the SQP method is used to optimize material models of cellular polymers. In particular, the optimization is shown for a viscoelastic polyurethane (PU) foam. Due to the high‐grade nonlinear material behaviour, the foam is modelled by a finite viscoelastic Ogden type law in the framework of the Theory of Porous Media (TPM).     

A cellular-automata method for studying porous media properties

A cellular-automata (CA) approach for investigating properties of porous media with tortuous channels and different smoothness of pore walls is proposed. This approach is aimed at combining two different CA models: the first one is intended for constructing the morphology of a porous material; the second, for simulating a fluid flow through it. The porous media morphology is obtained as a result of evolution of a cellular automaton, forming a “steady pattern.” The result is then used for simulating a fluid flow through a porous medium by applying the Lattice Gas CA model. The method has been tested on a small fragment of a porous material and implemented for investigating a carbon electrode of a hydrogen fuel cell on a multiprocessor cluster.     

Estimating the cycle time of three-stage material handling systems

Because of high investment costs, the productivity of material handling systems must be accurately estimated before various logistical, industrial, or transportation systems can be implemented. This paper proposes analytical models for three-stage material handling systems. Two possible approaches to the estimation of the productivity of three-stage material handling systems are considered: one, the continuous Markov chain model, and two, approximated mathematical models. The approximated models are based on the probability theory and permit very accurate calculations of the compound cycle time in cases when the probability distribution of the “technical” cycle times of each stage is known. Finally, some numerical results obtained by the proposed models are compared with those results obtained by a simulation study.     

颜色读数与物质浓度关系的数学建模

比色法作为一种检测物质浓度的常用方法,是通过待测物质溶液颜色读数来推测所含物质浓度,根据所提供的物质浓度与不同色卡颜色读数的试验数据,对物质浓度和颜色读数进行相关性分析,拟合颜色读数与物质浓度的指数模型和建立灰色预测模型,对这两种模型的拟合值与原始值的误差进行比较分析,给出了颜色读数与物质浓度关系的优化模型.     

A finite element framework for the numerical implementation of boundary potentials

This contribution deals with the implications of boundary potential energies on deformational mechanics in the framework of the finite element method at finite strains. The common material models in continuum mechanics are taking the bulk into account, nevertheless, neglecting the boundary. However, boundary effects sometimes play a dominant role in the material behavior, e.g. surface tension in fluids. The boundary potentials, in general, are allowed to depend not only on the boundary deformation gradient but also on the spatial surface–normal / curve–tangent, as well. For the finite element implementation, a suitable curvilinear coordinate system attached to the boundary is defined and corresponding geometrical and kinematical derivations are carried out. Afterwards, the discretization of the generalized weak formulation, including boundary potentials, is carried out and finally numerical examples are presented to demonstrate the boundary effects due to the different proposed material behavior. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The singular sources method for cracks

The singular sources method is given to detect the shape of a thin infinitely cylindrical obstacle from a knowledge of the TM‐polarized scattered electromagnetic field in large distance. The basic idea is based on the singular behaviour of the scattered field of the incident point source on the cross‐section of the cylinder. We assume that the scatterer is a perfect conductor which is possibly coated by a material and investigate two models with different boundary conditions. Also we give a uniqueness proof for the shape reconstruction. Copyright © 2006 John Wiley & Sons, Ltd.     

Noise Problem in Semiconductors

This paper deals with the problem of calculating the variancesand covariances of the fluctuating numbers of conduction electronsin semiconductors, since the noise power associated with thecurrents passing through the material can be obtained in termsof these quantities. The usual finite difference equation isobtained for the probability distribution of these electronnumbers, and it is shown how present methods of calculationare too cumbersome to apply to realistic semiconductor models.In a new calculation it is shown that it may be possible toobtain a better signal-to-noise ratio in the material by introducinga critical number of impurities. To study this principle ingeneral a more rapid and straightforward method must be foundof obtaining the required quantities. To achieve this aim anew method is tried—the difference equation is transformedinto a partial differential equation for the probability generatingfunction and is solved exactly for a case in which there isone type of impurity present. The equation has not yet beensolved when the material contains more than one type of impurity.