Constitutive behaviour of mixed mode loaded adhesive layer

Mixed mode testing of adhesive layer is performed with the Mixed mode double Cantilever Beam specimen. During the experiments, the specimens are loaded by transversal and/or shear forces; seven different mode mixities are tested. The J-integral is used to evaluate the energy dissipation in the failure process zone. The constitutive behaviour of the adhesive layer is obtained by a so called inverse method and fitting an existing mixed mode cohesive model, which uses a coupled formulation to describe a mode dependent constitutive behaviour. The cohesive parameters are determined by optimizing the parameters of the cohesive model to the experimental data. A comparison is made with the results of two fitting procedures. It is concluded that the constitutive properties are coupled, i.e. the peel and shear stress depend on both the peel and shear deformations. Moreover, the experiments show that the critical deformation in the peel direction is virtually independent of the mode mixity.     

Orthotropic elastic–plastic material model for paper materials

Paper and paperboard generally exhibit anisotropic and non-linear mechanical material behaviour. In this work, the development of an orthotropic elastic–plastic constitutive model, suitable for modelling of the material behaviour of paper is presented. The anisotropic material behaviour is introduced into the model by orthotropic elasticity and an isotropic plasticity equivalent transformation tensor. A parabolic stress–strain relation is adopted to describe the hardening of the material. The experimental and numerical procedures for evaluation of the required material parameters for the model are described. Uniaxial tensile testing in three different inplane material directions provides the calibration of the material parameters under plane stress conditions. The numerical implementation of the material model is presented and the model is shown to perform well in agreement with experimentally observed mechanical behaviour of paper.     

A constitutive model of aluminum foam for crash simulations

A new constitutive model for metallic foams is developed to overcome the deficiency of existing models in commercial finite element codes. The proposed constitutive model accounts for volume changes under hydrostatic compression and combines the hydrostatic pressure and von Mises stress into one yield function. The change of the compressibility of the metallic foam is handled in the constitutive model by allowing for shape changes of the yield surface in the hydrostatic pressure-von Mises stress space. The backward Euler method is adopted to integrate the constitutive equations to achieve numerical accuracy and stability. The model is implemented into LS-DYNA as a user-defined subroutine, verified with existing solutions, and validated with foam testing data. The verified and validated model is then utilized in the crushing simulations of foam-filled columns with square and hexagonal cross-sections. Two constitutive models are studied: the first using an exponential function to describe the relationship of plastic Poisson's ratio with respect to true strain and the second using linear interpolation function as an alternative approximation. The new foam model provides satisfactory prediction of crushing forces and deformed shapes compared to experimental results. Additionally, the new foam model was shown to have better numerical stability and accuracy than existing LS-DYNA built-in material models.     

A REVIEW ON RESEARCH ON DEVELOPMENT OF LIGAMENT CONSTITUTIVE RELATIONS ON MACRO, MESO, AND MICRO LEVELS

Ligaments are densely connective soft tissues capable of maintaining stability and function of knee joint.As an important factor,the constitutive relation of ligament would afect its biomechanics and further play an essential role in the research on ligament injury,healing and treatment.The objective of this paper is to provide an overview of the current research on ligament constitutive relations on the macro,meso,and micro levels as well as the anatomy and histological structure of ligament.Some studies of biomechanical behaviors during ligament injury and healing periods have also been investigated.Based on the research on ligament constitutive relation in the past three decades,a discussion of some research perspectives is also presented,such as a validated accurate measuring method of in situ strain in ligament,a new constitutive relation involving the distribution of ultra-structural properties,and a rational estimation of ligament injury and healing process by the change of its ultra-structural or histological characteristics.     

Characterizing void nucleation in a damage-based constitutive model using notched tensile sheet specimens

An experimental and numerical test programme was conducted to investigate damage-induced ductile fracture in notched tensile sheet specimens of an aluminum–magnesium alloy. An upper bound, damage-based constitutive model was employed to estimate the formability of the material over a range of stress states found in sheet metal forming. Stress- and strain-based nucleation models are evaluated to characterize damage initiation and fracture of the material. The ligament strain-to-failure, elongation-to-failure and load–displacement curves can be captured using either nucleation rule. The advantages of each nucleation model are discussed in relation to quantitative measurements of damage available in the literature. Stress-based nucleation provides a promising approach for characterizing the nucleation behaviour over a range of stress states compared to strain-based nucleation.     

NUMERICAL CALCULATION OF 2D,UNSTEADY FLOW IN CENTRIFUGAL PUMPS: IMPELLER AND VOLUTE INTERACTION

A numerical model is developed for calculating the two-dimensional, unsteady, incompressible and turbulent flow within the rotating impeller and stationary volute of an industrial centrifugal pump. The objective is the investigation and comprehension of the instantaneous behaviour of centrifugal pumps, aiming at the reduction of vibrations, radial forces and hydraulic noise. The computation is performed within a blade-to-blade streamtube for the impeller and a tube normal to the axis of rotation for the volute. The equations to be solved are the unsteady Reynolds-averaged Navier–Stokes equations along with the continuity equation and the unsteady κ–ϵ equations for turbulence modelling. The finite volume method is applied for space discretization and an implicit scheme for time discretization. A multidomain overlapping grid technique is used for matching together the relative flow field calculated within the rotating impeller and the absolute one calculated within the stationary volute. In this way the impeller and volute interaction is directly taken into account. The numerical model is validated for a centrifugal pump of N _q=32 under design flow conditions. Comparisons between calculation and measurements show fairly good agreement.     

Modeling of textural anisotropy in granular materials with stochastic micro-polar hypoplasticity

The paper deals with a numerical analysis of the effect of textural anisotropy on the behaviour of cohesionless granular materials with consideration of shear localization. For a simulation of the mechanical behaviour of a granular material during a monotonic deformation path, an isotropic micro-polar hypoplastic constitutive model was used. To describe textural effects, spatially correlated random fields of the initial void ratio were subject to rotation against the horizontal axis. The 2D random fields were generated using a conditional rejection method. The results were compared with those obtained with an anisotropic micro-polar constitutive model for a uniform distribution of the initial void ratio. The calculations were carried out with an initially dense granular specimen during plane strain compression under constant lateral pressure.     

Characterization of pressure-wave propagation during the condensation of R404A and R134a refrigerants in pipe mini-channels that undergo periodic hydrodynamic disturbances

In the present paper, an attempt was undertaken to model the propagation of a pressure wave triggered by periodic hydrodynamic instabilities in the condensation of the R404A and R134a refrigerants in pipe mini-channels. A homogenous transient two-fluid model was used based on balance equations. The model presents the complexity of multi-phase flows. The heat exchange between the phases in the condensation process was calculated using the one-dimensional form of Fourier’s equation. The dependence which defines the interfacial the drag forces for the flow in mini-channels was also taken into consideration. The results of the numerical calculations were verified with experimental investigations and satisfactory compliance was obtained.     

Application of MLPG in Large Deformation Analysis

Two-dimensional large deformation analysis of hyperelastic and elasto-plastic solids based on the Meshless Local Petrov–Galerkin method (MLPG) is presented. A material configuration based the nonlinear MLPG formulation is introduced for the large deformation analysis of both path-dependent and path-independent materials. The supports of the MLS approximation functions cover the same sets of nodes during material deformation, thus the shape function needs to be computed only in the initial stage. The multiplicative hyperelasto-plastic constitutive model is adopted to avoid objective time integration for stress update in large rotation. With this constitutive model, the computational formulations for path-dependent and path-independent materials become identical. Computational efficiency of the nonlinear MLPG method is discussed and optimized in several aspects to make the MLPG an O(N) algorithm. The numerical examples indicate that the MLPG method can solve large deformation problems accurately. Moreover, the MLPG computations enjoy better convergence rate than the FEM under very large particle distortion.The project supported by the National Natural Science Foundation of China (10472051). The English text was polished by Keren Wang.     

DETERMINATION OF QUADRICEPS FORCES IN SQUAT AND ITS APPLICATION IN CONTACT PRESSURE ANALYSIS OF KNEE JOINT

While the quadriceps muscles of human body are quite important to the daily activities of knee joints,the determination of quadriceps forces poses significant challenges since it cannot be measured in ...     

Numerical scheme for accurately capturing gas migration described by 1D multiphase drift flux model

A numerical scheme is designed and implemented to solve a simplified set of equations modeling 1-D multi-phase flow based on drift flux model in an isothermal setup with phase dissolution. The difficulty in obtaining the analytical Jacobian of the fluxes leads to the difficulty in obtaining an efficient linearized Riemann solver which in turn affects the accuracy in capturing the contact wavefront/gas migration. To address this issue a fully explicit second order finite volume solver based on flux corrected transport (FCT) is implemented. The choice of variables used for limiting the fluxes affects the amount of numerical diffusion and an appropriate choice of the gradient in volume fraction is used. Practical test cases while drilling in the oil and gas industry, of gas injection inside a well annulus and shut-in of the vertical well are presented. The results conclude that the FCT solver is better and efficient for accurately capturing gas migration for multi-phase models with phase behaviour involving slip velocities given by algebraic relations.     

An application of a damage constitutive model to concrete at high temperature and prediction of spalling

A characteristic feature of concrete under uniaxial compression is the development of cracks parallel to the loading direction. A damage constitutive model proposed by Ortiz [Ortiz, M., 1985. A constitutive theory for the inelastic behaviour of concrete. Mech. Mater. 4, 67–93] can predict the transverse tensile stress responsible for these cracks by considering the interaction between the aggregate and the mortar and the development of damage in the latter. When concrete is exposed to high temperature, as is the case during fire, the failure mode is thermal spalling. In order to improve the prediction of the stresses involved in this failure Ortiz’s model is extended to account for the effects of high temperature. Published experimental results for uniaxial and biaxial compression at high temperatures are used to calibrate the temperature dependence of some of the material properties. The transient creep strain is accounted for by modifying the constrained thermal strain. The stress analysis is coupled with hygro-thermal analysis of heat, mass transfer and pore pressure build-up. The effect of pore pressure on the damage evolution is modeled by applying a body force in the stress analysis module proportional to the pressure gradient. A numerical example of concrete under fire is solved and the computed results are discussed. Spalling is predicted when the damage variable reaches its maximum value of unity. The predicted depth and time of spalling for a range of variation of permeability and initial liquid water content are presented. They are in good agreement with published experimental results.     

<Emphasis Type="Italic">Ex Vivo</Emphasis> Intervertebral Disc Bulging Measurement Using a Fibre Bragg Grating Sensor

Advances using optical fibres as sensors may represent an important contribution for development of minimally invasive techniques in biomedical and biomechanical applications. Concerning spine injuries, intervertebral disc (IVD) degeneration is a major clinical issue since it represents gross structural disruption and it is irreversible. Measuring biomechanical parameters of the IVD should contribute for better understanding on its mechanical response to external applied forces. The purpose of this study was to explore the potential of a Fibre Bragg Grating (FBG) sensor to measure strain caused by bulging of the intervertebral disc under axial compression. Disc bulging is a consequence of IVD compression and a technique to register this behaviour is addressed in this study. Needle-mounted sensors were already used to measure IVD pressure in cadaveric material. In this study we also explored the possibility of using needles only for sensor guiding and positioning leaving sensor directly in contact with the IVD material. An ex vivo porcine dorsal functional spinal unit was instrumented with a FBG sensor and submitted to axial compression. Results suggest the sensor’s ability to measure strain response to load. Bulging of the annulus fibrosus as a consequence of axial compression was confirmed using the FBG sensor. Hysteresis and viscoelastic behaviour were observable suggesting that energy is dissipated by the deformation of the annulus and that unloading time was insufficient for disc recovery. Nevertheless the relatively low strain sensitivity of the sensor as well as signal artefacts caused by transverse loading may constitute a problem in the analysis and interpretation of strain data. The technique may not be suitable for measurement of physiologic bulging being more indicative of the radial force exerted by the annulus.     

Peridynamic States and Constitutive Modeling

A generalization of the original peridynamic framework for solid mechanics is proposed. This generalization permits the response of a material at a point to depend collectively on the deformation of all bonds connected to the point. This extends the types of material response that can be reproduced by peridynamic theory to include an explicit dependence on such collectively determined quantities as volume change or shear angle. To accomplish this generalization, a mathematical object called a deformation state is defined, a function that maps any bond onto its image under the deformation. A similar object called a force state is defined, which contains the forces within bonds of all lengths and orientation. The relation between the deformation state and force state is the constitutive model for the material. In addition to providing a more general capability for reproducing material response, the new framework provides a means to incorporate a constitutive model from the conventional theory of solid mechanics directly into a peridynamic model. It also allows the condition of plastic incompressibility to be enforced in a peridynamic material model for permanent deformation analogous to conventional plasticity theory.      

复杂应力状态下木材力学性能的数值模拟

针对木材复杂的各向异性材料特点,建立了能反映木材正交各向异性弹性、抗拉和抗压强度不等、抗拉或抗剪时发生脆性破坏而受压时发生塑性变形等特性的本构模型。将木材弹性应力—应变关系简化为正交各向异性;选用Yamada-Sun强度准则来判断木材抗压时是否屈服,抗拉或抗剪时是否发生应变软化;通过引入损伤因子和弹性应变能,建立了木材...     

材料非线性微-宏观分析的多尺度方法研究

介绍并比较了近年来在材料非线性微-宏观分析多级数值方法方面的研究工作. 针 对考虑材料内摩擦接触的颗粒材料多尺度计算问题，建立一种基于数值技术的多级分析方 法. 方法的特点是在对材料进行微观分析的基础上建立宏观材料的多尺度非线性数值本构模 型. 而对材料弹塑性多级分析问题，建立了基于转换场技术的算法，采用近似技术建立非线 性分析的本征应变矩阵，使方法具有表达简单与实现方便的特点. 给出了数值算例， 通过比较说明了方法的正确性与有效性.     

Cohesive–frictional interface constitutive model

In the framework of numerical analysis of joined bodies, the present paper is devoted to the constitutive modeling, via an interface kinematic formulation, of mechanical behaviour of internal adhesive layers. The proposed interface constitutive model couples a cohesive behaviour, based on the damage mechanics theory, with a frictional one, defined in a non-associative plasticity framework. Namely, the interface formulation follows the transition of the adhesive material from the sound elastic condition to the fully cracked one. This formulation is able to model, by means of a specific interpretation of the damage variable and in a relevant mathematical setting, the interface intermediate mechanical properties, during the microcracks spreading process up to the discontinuity surface formation (macrocrack). The constitutive modeling is performed in fully compliance with the thermodynamic principles, in order to ensure the thermodynamic consistency requirement. In the present work, various monotonic and cyclic loading conditions are examined in order to show the main features of the constitutive formulation as well as several significant differences with respect to other existing models. Computational efficiency of the interface constitutive model is tested in a numerical application by FEM resolution strategy approach.