Visco-electro-elastic models of fiber-distributed active tissues

We present a constitutive model for stochastically distributed fiber reinforced visco-active tissues, where the behavior of the reinforcement depends on the relative orientation of the electric field. Following our previous works, for the passive behaviors we adopt a second order approximation of the strain energy density associated to the parameters of the fiber distribution. Consistently, we also assume that the active behavior accounts for the stochastic distribution of the fibers. The ensuing mechanical quantities result to be dependent on two average structure tensors. We introduce an extended Helmholtz free energy density characterized by the inclusion of a directional active potential, dependent on a stochastic anisotropic permittivity tensor. The permittivity tensor is expanded in Taylor series up to the second order, allowing to obtain an approximated active potential with the same structure of the passive Helmholtz free energy density. In particular, the explicit expression of active stress and stiffness are dependent on the two average structure tensors that characterize the passive response. Anisotropy follows from the fiber distribution and inherits its stochastic nature through statistics parameters. The active fiber distributed model is extended here to viscous materials by including the contribution of a dual dissipation potential in the variational formulation of the constitutive updates. Additionally, we present a computational example of application of the electro-viscous-mechanical material model by simulating peristaltic contractions on a portion of human intestine.     

Effects of matrix inelasticity on the overall hysteretic behavior of TiNi-SMA fiber composites

The effects of the inelastic deformation of the matrix on the overall hysteretic behavior of a unidirectional titanium–nickel shape-memory alloy (TiNi-SMA) fiber composite and on the local pseudoelastic response of the embedded SMA fibers are studied under the isothermal loading and unloading condition. The multiaxial phase transformation of the SMA fibers is predicted using the phenomenological constitutive equations which can describe the two-step deformation due to the rhombohedral and martensitic transformations, and the inelastic behavior of the matrix material using the standard nonlinear viscoplastic model. The average behavior of the SMA composite is evaluated with the micromechanical method of cells. It is observed that the inelastic deformation of the matrix due to prior tension results in a compressive stress in the matrix after unloading of the SMA composite and this residual stress impedes the complete recovery of the pseudoelastic strain of the SMA fibers. This explains that a closed hysteresis behavior of the SMA composite is no longer observed in contrast with the case that an elastic behavior of matrix is assumed. The predicted local stress–strain behavior indicates that the cyclic response of matrix is crucial to the design of the hysteretic performance of the SMA composite under the repeated loading conditions.     

Residual Strains using Integrated Continuous Fiber Optic Sensing in Thermoplastic Composites and Structural Health Monitoring

The evolution of spatially resolved internal strain/stress during the manufacturing of thermoplastic composites and subsequent relaxation from water intake are evaluated using an in-situ fiber optic sensor corresponding to a coated optical glass fiber with a nominal diameter of 160 μm. Unidirectional carbon fiber-polyamide 6 composites are produced using compression molding with an embedded fiber optic for strain measurement. The distributed fiber optic based strain sensor is placed in an arrangement to capture 0, 45, and 90° strains in the composite to resolve in-plane strain tensor. Strains are monitored in the direction of fiber optic sensor along its length at high resolution during the various stages of compression molding process. Results indicate considerable internal strains leading to residual stress at the end of processing step along the off-axis (45°) and transverse (90°) directions, and small strains in the carbon fiber pre-preg (0°) direction. At the end of compression molding process, an average of 7000 and 10,000 compressive micro-strains are obtained for residual state of strain in the off-axis and transverse direction. Since water/moisture infusion affects the mechanical properties of polyamide-6 matrix resin, these composite panels with embedded sensors targeted for marine applications are monitored in a water bath at 40 °C simulating accelerated testing conditions. Using the same fiber optic sensor based technique, the strain relaxation was observed during water uptake demonstrating in-situ strain monitoring during both manufacturing and subsequent composite implementation/application environment. The technique presented in this paper shows the potential of optimizing time-temperature-pressure protocols typically utilized in thermoplastic manufacturing, and continuous life-cycle monitoring of composite materials using a small diameter and inexpensive distributed fiber optic sensing.     

短纤维增强三元乙丙橡胶横观各向同性黏-超弹性本构模型

短纤维增强三元乙丙橡胶包覆薄膜,是一种应用于固体火箭发动机缠绕包覆装药的新型复合材料.为了描述其在工作过程中受振动、冲击等载荷作用时的力学行为,基于黏弹性理论和纤维增强连续介质力学理论,提出了一种考虑应变率强化效应的横观各向同性黏-超弹本构模型.模型中应变能函数被分解为超弹性应变能和黏性应变能,其中超弹性应变能包括表征各向同性的橡胶基体应变能和表征各向异性的纤维拉伸应变能,黏性应变能采用表征橡胶和纤维黏性响应的宏观唯象模型.选取表征各应变能的函数形式,经过数学变换、替代、叠加,求解确定最终的应力应变形式,明确模型参数获取的具体步骤,将预测结果与实验结果对比分析,准确性较高.研究表明:该模型能有效预测材料在低应变率下纤维方向为0?～45?的非线性率相关力学特性;模型形式易于实现有限元开发,对固体火箭发动机装药结构完整性分析具有参考价值.     

Micro–macro analysis of shape memory alloy composites

The aim of the paper is to develop a micro–macro approach for the analysis of the mechanical behavior of composites obtained embedding long fibers of Shape Memory Alloys (SMA) into an elastic matrix. In order to determine the overall constitutive response of the SMA composites, two homogenization techniques are proposed: one is based on the self-consistent method while the other on the analysis of a periodic composite. The overall response of the SMA composites is strongly influenced by the pseudo-elastic and shape memory effects occurring in the SMA material. In particular, it is assumed that the phase transformations in the SMA are governed by the wire temperature and by the average stress tensor acting in the fiber. A possible prestrain of the fibers is taken into account in the model.Numerical applications are developed in order to analyze the thermo-mechanical behavior of the SMA composite. The results obtained by the proposed procedures are compared with the ones determined through a micromechanical analysis of a periodic composite performed using suitable finite elements.Then, in order to study the macromechanical response of structural elements made of SMA composites, a three-dimensional finite element is developed implementing at each Gauss point the overall constitutive laws of the SMA composite obtained by the proposed homogenization procedures. Some numerical applications are developed in order to assess the efficiency of the proposed micro–macro model.     

An anisotropic rotary diffusion model for fiber orientation in short- and long-fiber thermoplastics

The Folgar–Tucker model, which is widely-used to predict fiber orientation in injection-molded composites, accounts for fiber–fiber interactions using isotropic rotary diffusion. However, this model does not match all aspects of experimental fiber orientation data, especially for composites with long discontinuous fibers. This paper develops a fiber orientation model that incorporates anisotropic rotary diffusion. From kinetic theory we derive the evolution equation for the second-order orientation tensor, correcting some errors in earlier treatments. The diffusivity is assumed to depend on a second-order space tensor, which is taken to be a function of the orientation state and the rate of deformation. Model parameters are selected by matching the experimental steady-state orientation in simple shear flow, and by requiring stable steady states and physically realizable solutions. Also, concentrated fiber suspensions align more slowly with respect to strain than models based on Jeffery's equation, and we incorporate this behavior in an objective way. The final model is suitable for use in mold filling and other flow simulations, and it gives improved predictions of fiber orientation for injection molded long-fiber composites.     

An improved Jones-Nelson nonlinear material model and rubber composite lamina analysis

A Jones-Nelson model has been applied to depict nonlinear stress-strain relations of composite laminae, where mechanical properties were expressed by strain energy density. The nonlinear material matrix is only a function of the strain energy density. Then a material model could be conveniently applied under complex stress condition. In this paper, by introducing large displacement stress-strain measurement and varying-Poisson's ratio idea, an improved Jones-Nelson material model is presented, where the expanding problem of material properties and convergence problems are overcome. Meanwhile a discuss of the reorientation of fiber and a material nonlinear analysis of rubber composite lamina under super large deformation conditions are made. The prediction results of improved material model are in fairly good agreement with those of the experiments.     

Determination of the effective elastic–plastic response of metal–ceramic composites

The mechanical response of metal–ceramic composites is analysed through a homogenization model accounting for the mechanical behaviour of the constituent materials. In order to achieve this purpose a nonlinear homogenization method based on the phase field approach has been suitably implemented into a numerical code. A prescribed homogenized strain state is applied to a unit volume element of a metal–ceramic composite with proportional loading in which all components of the strain tensor are proportional to one scalar parameter. The mechanical response of the material has been modeled by considering a von Mises plasticity model for the metal phase and a Drucker–Prager associative elastic–plastic material model for the ceramic phase. A two stages plasticity has been obtained in which inelastic strain develops in the metal phase followed by a fully plastic response. A comparison with a finite element model of the stress–strain response of an axisymmetric unit cell has been carried out with the purpose to validate the homogenization based modeling presented in the paper. Plastic parameters of a Drucker–Prager yield surface for the homogenized composite have been calculated at different materials compositions. Associative Drucker–Prager plasticity has been found to be accurate for high ceramic content.     

A structural multi-mechanism constitutive equation for cerebral arterial tissue

A structural multi-mechanism constitutive equation is developed to describe the nonlinear, anisotropic, inelastic mechanical behavior of cerebral arterial tissue. Elastin and collagen fibers are treated as separate components (mechanisms) of the artery. Elastin is responsible for load bearing at low strain levels while the collagen mechanism is recruited for load bearing at higher strain levels. This work builds on an earlier model in which both the elastin and collagen mechanisms are represented by isotropic response functions [Wulandana, R., Robertson, A.M., 2005. An inelastic multi-mechanism constitutive equation for cerebral arterial tissue. Biomech. Model. Mechan. 4 (4), 235–248]. Here, the anisotropic material response of the wall is introduced through the collagen mechanism which is composed of helically distributed families of fibers. The orientation of these families is described using either a finite number of fiber orientations or a fiber distribution function. The fiber orientation or dispersion function can be prescribed directly from arterial histology data, or, taking a phenomenological approach, based on data fitting from bi-axial measurements. The activation of the collagen mechanism is specified using a new fiber strain based activation criterion. The multi-mechanism constitutive equation is applied to the simple case of cylindrical inflation and material constants are determined based on available inelastic experimental data for cerebral arteries. While the proposed model captures all features of this inelastic data, there is a pressing need for further experiments to refine the model.     

改进的Jones—Nelson物理非线性材料模型与橡胶复合材料层板分析

Ｊｏｎｅｓ－Ｎｅｌｓｏｎ模型是复合材料物理非线性应力应变关系的一种描述方法，它通过建立材料刚度与应变能密度的关系以及考虑纤维增强复合材料的物理非线性问题，非线性的材料矩阵只为应变能的子数，使得材料模型可以方便地应用地复杂应力状态下，通过 在大变形的应变应力变量和变泊松比概念，对Ｊｏｎｅｓ－Ｎｅｌｓｏｎ模型进行了改进，解决了材料特性的扩展问题和收敛问题，同时考虑了大变形中纤维铺设角度的重新取向，使其     

含正交排列夹杂和缺陷材料的等效弹性模量和损伤

研究含正交排列夹杂和缺陷材料的等效弹性模量和损伤,推导了以Eshelby－Mori－Tanaka方法求解多相各向异性复合材料等效弹性模量的简便计算公式,针对含三相正交椭球状夹杂的正交各向异性材料,得到了由细观参量（夹杂的形状、方位和体积分数）表示的等效弹性模量的解析表达式．在此基础上,提出了一个宏细观结合的正交各向异性损伤模型,从而建立了以细观量为参量的含损伤材料的应力应变关系．最后,对影响材料损伤的细观结构参数进行了分析．     

Non-linear strain rate dependent micro-mechanical composite material model for finite element impact and crashworthiness simulation

The present study aims at implementation of a strain rate dependent, non-linear, micro-mechanics material model for laminated, unidirectional polymer matrix composites into the explicit finite element code LSDYNA. The objective is to develop an accurate and simple micro-mechanical, rate dependent material model, which is computationally efficient. Within the model a representative volume cell is assumed. The stress-strain relation including rate dependent effects for the micro-model is derived for both shell elements and solid elements. Micro-failure criterion is presented for each material constituent and failure mode. The implemented model can deal with problems such as impact, crashworthiness, and failure analysis under quasi-static loads. The developed material model has a wide range of applications such as jet engine jackets, armor plates, and structural crashworthiness simulation. The deformation response of two representative composite materials with varying fiber orientation is presented using the described technique. The predicted results compare favorably to experimental values.     

Continuous versus discrete (invariant) representations of fibrous structure for modeling non-linear anisotropic soft tissue behavior

The non-linear anisotropic mechanical response of soft tissue is largely dependent on the structure of the underlying collagen network. Collagen structure has been successfully quantified for various tissue types in terms of a locally defined fiber orientation distribution function. The continuous distribution function derived from structural data can be directly incorporated into an integral representation of the strain energy function for modeling tissue behavior. Alternatively, non-integral (often invariant-based) strain energy functions have been developed in which the collagen network structure is approximated using a discrete set of fiber classes. The advantage of such an approach is increased computational efficiency since the values of the strain energy and its derivatives (e.g. stress) can be evaluated without numerical integration. However, because of the structural simplifications such models are presumably unable to predict mechanical data as accurately as the models which incorporate a continuous orientation distribution function. In this work the ability of discrete versus continuous fiber models to capture the non-linear anisotropic response of soft tissue is critically analyzed. Both unimodal and bimodal fiber distributions are considered. A general formulation has been developed in terms of an arbitrary fiber strain energy function, such that the analysis can be performed for any suitable fiber material model. For tissue structures in which a discrete representation is suitable, techniques are presented for establishing the range of loading conditions in which model accuracy is not significantly compromised, thus justifying the use of an invariant-based modeling approach.     

Folding of fiber composites with a hyperelastic matrix

This paper presents an experimental and numerical study of the folding behavior of thin composite materials consisting of carbon fibers embedded in a silicone matrix. The soft matrix allows the fibers to microbuckle without breaking and this acts as a stress relief mechanism during folding, which allows the material to reach very high curvatures. The experiments show a highly non-linear moment vs. curvature relationship, as well as strain softening under cyclic loading. A finite element model has been created to study the micromechanics of the problem. The fibers are modeled as linear-elastic solid elements distributed in a hyperelastic matrix according to a random arrangement based on experimental observations. The simulations obtained from this model capture the detailed micromechanics of the problem and the experimentally observed non-linear response. The proposed model is in good quantitative agreement with the experimental results for the case of lower fiber volume fractions but in the case of higher volume fractions the predicted response is overly stiff.     

Strain-space plasticity analysis of fibrous composites

A strain-space formulation of plasticity theory for metal matrix fibrous composites is discussed. Specific results are obtained for a composite system in which the fibers are elastic until failure, while the matrix is of the Mises type, with kinematic hardening. The material model of the composite is based on the assumption that the fiber diameter is vanishingly small, and the fiber volume fraction is finite. Explicit expressions are obtained for instantaneous strain concentration factors in the phases, for the instantaneous overall stiffness, and for the overall loading surface and hardening rule during mechanical straining along a prescribed path.     

基于迟滞行为的2D-SiC/SiC复合材料组份力学性能分析

基于剪滞理论, 建立了单向纤维增强陶瓷基复合材料的加卸载理论模型, 分析了基体长碎块和短碎块对材料迟滞力学行为的不同影响. 通过拉伸循环加卸载试验, 获得了2D-SiC/SiC 复合材料的迟滞应力—应变行为.依据材料基体损伤特点, 将试验结果代入长碎块对应理论推导结果, 计算得到了4 个表征材料组份性能的参数:基体开裂应力为90 MPa, 热残余应力为19 MPa, 界面脱粘能为3.1 Jm2, 界面滑移力为74 MPa. 最后结合少量短碎块的存在对试验结果的影响, 定性分析了计算结果的偏差. 结果表明, 获得的材料组分性能参数具有较小的分散性, 并能够准确表征材料整体的力学行为.      

基于神经网络与光纤传感阵列的结构状态监测方法

以探测机敏复合材料与结构内应变、应力及损伤等物理状态的位置信息为目的．提出了一种新颖的可内埋于材料与结构内作为结构状态监测用的强度调制型光纤传感阵列网络，并采用人工神经网络来处理传感阵列输出的并行分布式传感信号，阐述了适用的Ｋｏｈｏｎｅｎ网模型及其变化形式，给出了仿真实验结果     

A composites-based hyperelastic constitutive model for soft tissue with application to the human annulus fibrosus

This paper presents a composites-based hyperelastic constitutive model for soft tissue. Well organized soft tissue is treated as a composite in which the matrix material is embedded with a single family of aligned fibers. The fiber is modeled as a generalized neo-Hookean material in which the stiffness depends on fiber stretch. The deformation gradient is decomposed multiplicatively into two parts: a uniaxial deformation along the fiber direction and a subsequent shear deformation. This permits the fiber-matrix interaction caused by inhomogeneous deformation to be estimated by using effective properties from conventional composites theory based on small strain linear elasticity and suitably generalized to the present large deformation case. A transversely isotropic hyperelastic model is proposed to describe the mechanical behavior of fiber-reinforced soft tissue. This model is then applied to the human annulus fibrosus. Because of the layered anatomical structure of the annulus fibrosus, an orthotropic hyperelastic model of the annulus fibrosus is developed. Simulations show that the model reproduces the stress-strain response of the human annulus fibrosus accurately. We also show that the expression for the fiber-matrix shear interaction energy used in a previous phenomenological model is compatible with that derived in the present paper.     

A constitutive model for shape memory alloys accounting for thermomechanical coupling

This paper presents a generalized Zaki-Moumni (ZM) model for shape memory alloys (SMAs) [cf. Zaki, W., Moumni, Z., 2007a. A three-dimensional model of the thermomechanical behavior of shape memory alloys. J. Mech. Phys. Solids 55, 2455-2490 accounting for thermomechanical coupling. To this end, the expression of the Helmholtz free energy is modified in order to derive the heat equation in accordance with the principles of thermodynamics. An algorithm is proposed to implement the coupled ZM model into a finite element code, which is then used to solve a thermomechanical boundary value problem involving a superelastic SMA structure. The model is validated against experimental data available in the literature. Strain rate dependence of the mechanical pseudoelastic response is taken into account with good qualitative as well as quantitative accuracy in the case of moderate strain rates and for mechanical results in the case of high strain rates. However, only qualitative agreement is achieved for thermal results at high strain rates. It is shown that this discrepancy is mainly due to localization effects which are note taken into account in our model. Analyzing the influence of the heat sources on the material response shows that the mechanical hysteresis is mainly due to intrinsic dissipation, whereas the thermal response is governed by latent heat. In addition, the variation of the area of the hysteresis loop with respect to the strain rate is discussed. It is found that this variation is not monotonic and reaches a maximum value for a certain value of strain rate.