Spline-based hyperelasticity for transversely isotropic incompressible materials

In this paper we present a spline-based hyperelastic model for incompressible transversely isotropic solids. The formulation is based on the Sussman-Bathe model for isotropic hyperelastic materials. We extend this model to transversely isotropic materials following a similar procedure. Our formulation is able to exactly represent the prescribed behavior for isotropic hyperelastic solids, recovering the Sussman-Bathe model, and to exactly or closely approximate the prescribed behavior for transversely isotropic solids. We have employed our formulation to predict, very accurately, the experimental results of Diani et al. for a transversely isotropic hyperelastic nonlinear material.

Finite elastic deformations of transversely isotropic circular cylindrical tubes

We consider the finite radially symmetric deformation of a circular cylindrical tube of a homogeneous transversely isotropic elastic material subject to axial stretch, radial deformation and torsion, supported by axial load, internal pressure and end moment. Two different directions of transverse isotropy are considered: the radial direction and an arbitrary direction in planes normal locally to the radial direction, the only directions for which the considered deformation is admissible in general. In the absence of body forces, formulas are obtained for the internal pressure, and the resultant axial load and torsional moment on the ends of the tube in respect of a general strain-energy function. For a specific material model of transversely isotropic elasticity, and material and geometrical parameters, numerical results are used to illustrate the dependence of the pressure, (reduced) axial load and moment on the radial stretch and a measure of the torsional deformation for a fixed value of the axial stretch.

Modelización constitutiva y computacional del daño y la fractura de materiales compuestos

En el trabajo se definen modelos constitutivos que permiten reproducir el proceso de fallo de estructuras de materiales compuestos en distintas escalas bajo cargas estáticas. Se define un modelo constitutivo para determinar la respuesta de estructuras de materiales compuestos mediante la teoría de laminados. El modelo es validado mediante un programa de ensayos experimentales con probetas con un agujero central geométricamente similares. Se muestra la capacidad del modelo de detectar el efecto tamaño. Se define un modelo constitutivo para materiales transversalmente isótropos bajo estados tridimensionales de tensión. El modelo se valida analizando numéricamente el proceso de agrietamiento de la matriz. Finalmente se desarrolla un modelo analítico para determinar el agrietamiento de la matriz y la delaminación entre las capas.

Flow and fracture behaviour of high performance alloys

Based on our needs, that is to say, through precise simulation of the impact phenomena that may occur inside a jet engine turbine with an explicit non-linear finite element code, four new material models are postulated. Each one of is calibrated for four high-performance alloys that can be encountered in a modern jet engine. A new uncoupled material model for high strain and ballistic is proposed. Based on a Johnson-Cook type model, the proposed formulation introduces the effect of the third deviatoric invariant by means of three different Lode angle dependent functions. The Lode dependent functions are added to both plasticity and failure models. The postulated model is calibrated for a 6061-T651 aluminium alloy with data taken from the literature. The fracture pattern predictability of the JCX material model is shown performing numerical simulations of various quasi-static and dynamic tests. As an extension of the above-mentioned model, a modification in the thermal softening behaviour due to phase transformation temperatures is developed (JCXt). Additionally, a Lode angle dependent flow stress is defined. Analysing the phase diagram and high temperature tests performed, phase transformation temperatures of the FV535 stainless steel are determined. The postulated material model constants for the FV535 stainless steel are calibrated. A coupled elastoplastic-damage material model for high strain and ballistic applications is presented (JCXd). A Lode angle dependent function is added to the equivalent plastic strain to failure definition of the Johnson-Cook failure criterion. The weakening in the elastic law and in the Johnson-Cook type constitutive relation implicitly introduces the Lode angle dependency in the elastoplastic behaviour. The material model is calibrated for precipitation hardened Inconel 718 nickel-base superalloy. The combination of a Lode angle dependent failure criterion with weakened constitutive equations is proven to predict fracture patterns of the mechanical tests performed and provide reliable results. A transversely isotropic material model for directionally solidified alloys is presented. The proposed yield function is based a single linear transformation of the stress tensor. The linear operator weighs the degree of anisotropy of the yield function. The elastic behaviour, as well as the hardening, are considered isotropic. To model the hardening, a Johnson-Cook type relation is adopted. A material vector is included in the model implementation. The failure is modelled with the Cockroft-Latham failure criterion. The material vector allows orienting the reference orientation in any other that the user may need. The model is calibrated for the MAR-M 247 directionally solidified nickel-base superalloy.

Discontinuous solutions and boundary value problems for non-linearly elastic fibre-reinforced materials

En este trabajo se han analizado varios problemas en el contexto de la elasticidad no lineal basándose en modelos constitutivos representativos. En particular, se han analizado problemas relacionados con el fenómeno de perdida de estabilidad asociada con condiciones de contorno en el caso de material reforzados con fibras. Cada problema se ha formulado y se ha analizado por separado en diferentes capítulos. En primer lugar se ha mostrado el análisis del gradiente de deformación discontinuo para un material transversalmente isótropo, en particular, el modelo del material considerado consiste de una base neo-Hookeana isótropa incrustada con fibras de refuerzo direccional caracterizadas con un solo parámetro. La solución de este problema se vincula con instabilidades que dan lugar al mecanismo de fallo conocido como banda de cortante. La perdida de elipticidad de las ecuaciones diferenciales de equilibrio es una condición necesaria para que aparezca este tipo de soluciones y por tanto las inestabilidades asociadas. En segundo lugar se ha analizado una deformación combinada de extensión, inación y torsión de un tubo cilíndrico grueso donde se ha encontrado que la deformación citada anteriormente puede ser controlada solo para determinadas direcciones de las fibras refuerzo. Para entender el comportamiento elástico del tubo considerado se ha ilustrado numéricamente los resultados obtenidos para las direcciones admisibles de las fibras de refuerzo bajo la deformación considerada. En tercer lugar se ha estudiado el caso de un tubo cilíndrico grueso reforzado con dos familias de fibras sometido a cortante en la dirección azimutal para un modelo de refuerzo especial. En este problema se ha encontrado que las inestabilidades que aparecen en el material considerado están asociadas con lo que se llama soluciones múltiples de la ecuación diferencial de equilibrio. Se ha encontrado que el fenómeno de instabilidad ocurre en un estado de deformación previo al estado de deformación donde se pierde la elipticidad de la ecuación diferencial de equilibrio. También se ha demostrado que la condición de perdida de elipticidad y ^W=2 = 0 (la segunda derivada de la función de energía con respecto a la deformación) son dos condiciones necesarias para la existencia de soluciones múltiples. Finalmente, se ha analizado detalladamente en el contexto de elipticidad un problema de un tubo cilíndrico grueso sometido a una deformación combinada en las direcciones helicoidal, axial y radial para distintas geotermias de las fibras de refuerzo . In the present work four main problems have been addressed within the framework of non-linear elasticity based on representative constitutive models. Namely, problems related to the loss of stability phenomena associated with boundary value problems for fibre-reinforced materials. Each of the considered problems is formulated and analysed separately in different chapters. We first start with the analysis of discontinuous deformation gradients for a transversely isotropic material under plane deformation. In particular, the material model is an augmented neo-Hookean base with a simple unidirectional reinforcement characterised by a single parameter. The solution of this problem is related to material instabilities and it is associated with a shear band-type failure mode. The loss of ellipticity of the governing differential equations is a necessary condition for the existence of these material instabilities. The second problem involves a detailed analysis of the combined non-linear extension, inflation and torsion of a thick-walled circular cylindrical tube where it has been found that the aforementioned deformation is controllable only for certain preferred directions of transverse isotropy. Numerical results have been illustrated to understand the elastic behaviour of the tube for the admissible preferred directions under the considered deformation. The third problem deals with the analysis of a doubly fibre-reinforced thickwalled circular cylindrical tube undergoing pure azimuthal shear for a special class of the reinforcing model where multiple non-smooth solutions emerge. The associated instability phenomena are found to occur prior to the point where the nominal stress tensor changes monotonicity in a particular direction. It has been also shown that the loss of ellipticity condition that arises from the equilibrium equation and ^W=2 = 0 (the second derivative of the strain-energy function with respect to the deformation) are equivalent necessary conditions for the emergence of multiple solutions for the considered material. Finally, a detailed analysis in the basis of the loss of ellipticity of the governing differential equations for a combined helical, axial and radial elastic deformations of a fibre-reinforced circular cylindrical tube is carried out.

Study of an anisotropic polymeric cellular material under compression loading

This paper emphasizes the influence of micro mechanisms of failure of a cellular material on its phenomenological response. Most of the applications of cellular materials comprise a compression loading. Thus, the study focuses on the influence of the anisotropy in the mechanical behavior of cellular material under cyclic compression loadings. For this study, a Digital Image Correlation (DIC) technique (named Correli) was applied, as well as SEM (Scanning Electron Microscopy) images were analyzed. The experimental results are discussed in detail for a closed-cell rigid poly (vinyl chloride) (PVC) foam, showing stress-strain curves in different directions and why the material can be assumed as transversely isotropic. Besides, the present paper shows elastic and plastic Poisson's ratios measured in different planes, explaining why the plastic Poisson's ratios approach to zero. Yield fronts created by the compression loadings in different directions and the influence of spring-back phenomenon on hardening curves are commented, also.

Analytical formulae for electromechanical effective properties of 3-1 longitudinally porous piezoelectric materials

A unidirectional fiber composite is considered here, the fibers of which are empty cylindrical holes periodically distributed in a transversely isotropic piezoelectric matrix, The empty-fiber cross-section is circular and the periodicity is the same in two directions at an angle pi/2 or pi/3. Closed-form formulae for all electromechanical effective properties of these 3-1 longitudinally periodic porous piezoelectric materials are presented. The derivation of such expressions is based on the asymptotic homogenization method as a limit of the effective properties of two-phase transversely isotropic parallel fiber-reinforced composites when the fibers properties tend to zero. The plane effective coefficients satisfy the corresponding Schulgasser-Benveniste-Dvorak universal type of relations, A new relation among the antiplane effective constants from the solutions of two antiplane strains and potential local problems is found. This relation is valid for arbitrary shapes of the empty-fiber cross-sections. Based on such a relation, and using recent numerical results for isotropic conductive composites, the antiplane effective properties are computed for different geometrical shapes of the empty-fiber cross-section. Comparisons with other analytical and numerical theories are presented. (c) 2008 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

Orthogonal packing of identical rectangles within isotropic convex regions

A mixed integer continuous nonlinear model and a solution method for the problem of orthogonally packing identical rectangles within an arbitrary convex region are introduced in the present work. The convex region is assumed to be made of an isotropic material in such a way that arbitrary rotations of the items, preserving the orthogonality constraint, are allowed. The solution method is based on a combination of branch and bound and active-set strategies for bound-constrained minimization of smooth functions. Numerical results show the reliability of the presented approach. (C) 2010 Elsevier Ltd. All rights reserved.

Effective Properties Evaluation for Smart Composite Materials

The purpose of this article is to present a method which consists in the development of unit cell numerical models for smart composite materials with piezoelectric fibers made of PZT embedded in a non-piezoelectric matrix (epoxy resin). This method evaluates a globally homogeneous medium equivalent to the original composite, using a representative volume element (RVE). The suitable boundary conditions allow the simulation of all modes of the overall deformation arising from any arbitrary combination of mechanical and electrical loading. In the first instance, the unit cell is applied to predict the effective material coefficients of the transversely isotropic piezoelectric composite with circular cross section fibers. The numerical results are compared to other methods reported in the literature and also to results previously published, in order to evaluate the method proposal. In the second step, the method is applied to calculate the equivalent properties for smart composite materials with square cross section fibers. Results of comparison between different combinations of circular and square fiber geometries, observing the influence of the boundary conditions and arrangements are presented.

Designing piezoresistive plate-based sensors with distribution of piezoresistive material using topology optimization

Piezoresistive sensors are commonly made of a piezoresistive membrane attached to a flexible substrate, a plate. They have been widely studied and used in several applications. It has been found that the size, position and geometry of the piezoresistive membrane may affect the performance of the sensors. Based on this remark, in this work, a topology optimization methodology for the design of piezoresistive plate-based sensors, for which both the piezoresistive membrane and the flexible substrate disposition can be optimized, is evaluated. Perfect coupling conditions between the substrate and the membrane based on the `layerwise' theory for laminated plates, and a material model for the piezoresistive membrane based on the solid isotropic material with penalization model, are employed. The design goal is to obtain the configuration of material that maximizes the sensor sensitivity to external loading, as well as the stiffness of the sensor to particular loads, which depend on the case (application) studied. The proposed approach is evaluated by studying two distinct examples: the optimization of an atomic force microscope probe and a pressure sensor. The results suggest that the performance of the sensors can be improved by using the proposed approach.

Experimental analyses of the poly(vinyl chloride) foams' mechanical anisotropic behavior

Assessing a full set of mechanical properties is a rather complicate task in the case of foams, especially if material models must be calibrated with these results. Many issues, for example anisotropy and heterogeneity, influence the mechanical behavior. This article shows through experimental analyses how the microstructure affects different experimental setups and it also quantifies the degree of anisotropy of a poly(vinyl chloride) foam. Monotonic and cyclic experimental tests were carried out using standard compression specimens and non-standard tensile specimens. Results are complemented and compared with the aid of a digital image correlation technique and scanning electron microscopy analyses. Mechanical properties (e.g., elastic and plastic Poisson's ratios) are evaluated for compression and tensile tests, for two different material directions (normal and in-plane). The material is found to be transversely isotropic. Differences in the results of the mechanical properties can be as high as 100%, or even more depending on the technique used and the loading direction. Also, the experimental analyses show how the material's microstructure behavior, like the evolution of the herein identified yield fronts and a spring back phenomenon, can influence the phenomenological response and the failure mechanisms as well as the hardening curves. POLYM. ENG. SCI., 52:2654-2663, 2012. (C) 2012 Society of Plastics Engineers

Effective properties evaluation for smart composite materials

The purpose of this article is to present a method which consists in the development of unit cell numerical models for smart composite materials with piezoelectric fibers made of PZT embedded in a non-piezoelectric matrix (epoxy resin). This method evaluates a globally homogeneous medium equivalent to the original composite, using a representative volume element (RVE). The suitable boundary conditions allow the simulation of all modes of the overall deformation arising from any arbitrary combination of mechanical and electrical loading. In the first instance, the unit cell is applied to predict the effective material coefficients of the transversely isotropic piezoelectric composite with circular cross section fibers. The numerical results are compared to other methods reported in the literature and also to results previously published, in order to evaluate the method proposal. In the second step, the method is applied to calculate the equivalent properties for smart composite materials with square cross section fibers. Results of comparison between different combinations of circular and square fiber geometries, observing the influence of the boundary conditions and arrangements are presented.

Modelo de comportamiento plástico y rotura en aleaciones de solidificación direccional

Esta investigación presenta un modelo de material para aleaciones de solidificación direccional que poseen un comportamiento mecánico transversalmente isótropo. Se han realizado una serie de ensayos de tracción sobre probetas cilíndricas a varias velocidades de deformación y a varias temperaturas sobre la superaleación de base níquel de solidificación direccional MAR-M 247 con objeto de conocer su comportamiento mecánico. Los ensayos se realizaron sobre probetas cilíndricas cuya dirección longitudinal forma 0º y 90º con la de la orientación de crecimiento de los granos. Para representar el comportamiento plástico anisótropo se ha formulado una función de plastificación de forma no cuadrática basada en la transformación lineal de tensores. Con el propósito de simplificar en todo lo posible el modelo se ha considerado un endurecimiento isótropo. Para probar la validez del modelo propuesto se ha implementado el mismo como modelo de material definido por el usuario en el código no lineal de elementos finitos LS-DYNA. In this research a material model for directionally solidified alloys with transversely isotropic mechanic behavior is presented. In order to characterize the mechanical behavior of the Mar-M 247 directionally solidified nickel based superalloy, tensile tests of axisymmetric smooth specimens were performed at various strain rates and temperatures. The specimens were machined making sure that the longitudinal axis of them was forming 0º and 90º with the grain growth orientation. To represent the plastic flow, a non-quadratic anisotropic function based on linear transformation of tensors has been formulated. For the sake of simplicity isotropic strain hardening of the material has been considered. To prove the validity of the model, a material subroutine has been implemented in LS-DYNA non-linear finite element code as a user defined material model.

WYPIWYG Hyperelasticity

In this paper we describe a new promising procedure to model hyperelastic materials from given stress-strain data. The main advantage of the proposed method is that the user does not need to have a relevant knowledge of hyperelasticity, large strains or hyperelastic constitutive modelling. The engineer simply has to prescribe some stress strain experimental data (whether isotropic or anisotropic) in also user prescribed stress and strain measures and the model almost exactly replicates the experimental data. The procedure is based on the piece-wise splines model by Sussman and Bathe and may be easily generalized to transversely isotropic and orthotropic materials. The model is also amenable of efficient finite element implementation. In this paper we briefly describe the general procedure, addressing the advantages and limitations. We give predictions for arbitrary ?experimental data? and also give predictions for actual experiments of the behaviour of living soft tissues. The model may be also implemented in a general purpose finite element program. Since the obtained strain energy functions are analytic piece-wise functions, the constitutive tangent may be readily derived in order to be used for implicit static problems, where the equilibrium iterations must be performed and the material tangent is needed in order to preserve the quadratic rate of convergence of Newton procedures.