A concurrent micromechanical model for predicting nonlinear viscoelastic responses of composites reinforced with solid spherical particles

A concurrent micromechanical model for predicting nonlinear viscoelastic responses of particle reinforced polymers is developed. Particles are in the form of solid spheres having micro-scale diameters. The composite microstructures are idealized by periodically distributed cubic particles in a matrix medium. Each particle is assumed to be fully surrounded by polymeric matrix such that contact between particles can be avoided. A representative volume element (RVE) is then defined by a single particle embedded in the cubic matrix. A spatial periodicity boundary condition is imposed to the RVE. One eighth unit-cell model with four particle and polymer subcells is generated due to the three-plane symmetry of the RVE. The solid spherical particle is modeled as a linear elastic material. The polymeric matrix follows nonlinear viscoelastic behaviors of thermorheologically simple materials. The homogenized micromechanical relation is developed in terms of the average strains and stresses in the subcells and traction continuity and displacement compatibility at the subcells’ interfaces are imposed. A stress–strain correction scheme is also formulated to satisfy the linearized micromechanical and the nonlinear constitutive relations. The micromechanical model provides three-dimensional (3D) effective properties of homogeneous composite responses, while recognizing microstructural geometries and in situ material properties of the heterogeneous medium. The micromechanical formulation is designed to be compatible with general displacement based finite element (FE) analyses. Experimental data and analytical micromechanical models available in the literature are used to verify the capability of the above micromechanical model for predicting the overall composite behaviors. The proposed micromodel is also examined in terms of computational efficiency and accuracy.     

Materially and geometrically non-linear woven composite micro-mechanical model with failure for finite element simulations

A computational micro-mechanical material model of woven fabric composite material is developed to simulate failure. The material model is based on repeated unit cell approach. The fiber reorientation is accounted for in the effective stiffness calculation. Material non-linearity due to the shear stresses in the impregnated yarns and the matrix material is included in the model. Micro-mechanical failure criteria determine the stiffness degradation for the constituent materials. The developed material model with failure is programmed as user-defined sub-routine in the LS-DYNA finite element code with explicit time integration. The code is used to simulate the failure behavior of woven composite structures. The results of finite element simulations are compared with available test results. The model shows good agreement with the experimental results and good computational efficiency required for finite element simulations of woven composite structures.     

形状记忆合金纤维复合材料的等效力学行为

在Aboudi提出的胞元模型以及Liu等建立的形状记忆合金的本构模型的基础上，由Legendre多项式，假设每个子胞元的位移场、应变场和应力场，再由子胞元间交界面的应力连续条件和外荷载边界条件推导出基体为弹塑性材料的形状记忆合金纤维复合材料的胞元模型；模拟了呈周期对称的形状记忆合金纤维复合材料受轴向单向拉伸、横向拉伸和横向剪切荷载作用下的等效力学行为，与有限元解进行了比较，结果基本一致。与有限元法比较起来，本文推导出的形状记忆合金纤维复合材料的胞元模型更具高效性。     

Nonlinear finite element analysis of inflatable beams made from orthotropic woven fabric

This paper was devoted to the three-dimensional nonlinear finite element analysis of inflatable beams. The beams under consideration are made of modern textile materials and can be used as a load-bearing beams or arches when inflated. A 3D Timoshenko beam with a homogeneous orthotropic woven fabric (OWF) was proposed. The model took into account the geometric nonlinearities and the follower force resulting from the inflation pressure. The use was made of the usual total Lagrangian form of the virtual work principle to perform the nonlinear equilibrium equations which were discretized by the finite element method. Two kinds of solutions were then investigated: finite elements solutions for linearized problems which were obtained by the means of the linearization around the prestressed reference configuration of the nonlinear equations and nonlinear finite element solutions which were performed by the use of an optimization algorithm based on the Quasi-Newton method. As an example, the bending problem of a cantilever inflated beam under concentrated load was considered and the deflection results improve the existing theoretical models. As these beams are made from fabric, the beam models were validated through their comparison with a 3D thin-shell finite element model. The influence of the material effective properties and the inflation pressure on the beam response was also investigated through a parametric study. The finite elements solutions for linearized problems were found to be close to the theoretical results existing in the literature. On the other hand, the results for the nonlinear finite element model were shown to be close to the results for the linearized finite elements model in the case of high mechanical properties and the nonlinear finite element model was used to improve the linearized model when the mechanical properties of the fabric are low.     

Calibrating a J2 plasticity material model using a 2D inverse finite element procedure

Material models are the key ingredients to accurately capture the global mechanical response of structural systems. The use of finite element analysis has proven to be effective in simulating nonlinear engineering applications. However, the choice of the appropriate material model plays a big role in the value of the numerical predictions. Such models are not expected to exactly reproduce global experimental response in all cases. Alternatively, the measured global response at specific domain or surface points can be used to guide the nonlinear analysis to successively extract a representative material model. By selecting an initial set of stress–strain data points, the load–displacement response at the monitoring points is computed in a forward incremental analysis without iterations. This analysis retains the stresses at the integration points. The corresponding strains are not accurate since the computed displacements are not anticipated to match the measured displacements at the monitoring points. Therefore, a corrective incremental displacement analysis is performed at the same load steps to adjust for displacements and strains everywhere by matching the measured displacements at the monitoring points. The stress–strain vectors at the most highly stressed integration point are found to establish an improved material model. This model is used within a multi-pass incremental nonlinear finite element analysis until the discrepancy between the measured and the predicted structural response at the monitoring points vanishes. The J2 flow theory of plasticity is used as a constitutive framework to build the tangent elastic–plastic matrices. The applicability of the proposed approach is demonstrated by solving 2D inverse continuum problems. The comparisons presented support the effectiveness of the proposed approach in accurately calibrating the J2 plasticity material model for such problems.     

NONLINEAR MICRO－MECHANICAL MODEL FOR PLAIN WOVEN FABRIC

The warp yarns and weft yarns of plain woven fabric which, being the principal axes of material of fabric, are orthogonal in the original configuration, but are obliquely crossed in the deformed configuration in general. The orthotropic constitutive model is unsuitable for fabric. In the oblique principal axes system the relations between loaded stress vectors and stress tensor are investigated, the stress fields of micro-weaving structures of fabric due to pure shear are carefully studied and, finally, a nonlinear micro-mechanical model for plain woven fabric is proposed. This model can accurately describe the nonlinear mechanical behavior of fabric observed in experiments. Under the assumption of small deformation and linearity of mechanical properties of fabric the model will degenerate into the existing linear model.     

A cell-based constitutive model for embryonic epithelia and other planar aggregates of biological cells

Embryonic epithelia are shown to have in common with plastic materials a number of key characteristics, including fabric evolution, “yielding”, “particle” (cell) rearrangement, energy dissipation, dependence of stress on fabric and irreversibility of deformation. The strains apparent at the tissue level can be large (several hundred percent over the course of 5–10 h), and are possible because of in-plane cell rearrangement. We propose a cell-based constitutive model, the first of its kind, to relate in-plane stresses, tissue deformations, evolution of cellular fabric (cell size, shape and orientation), mitosis and cell rearrangement. The governing equations are based on results from finite element models, statistical mechanics analyses and experiments. The constitutive model overcomes drawbacks of existing finite element models where cells are modeled using multiple elements, and it confirms that tissue fabric is a primary determinant of stress and deformation. Fabric predictions made using the model are as good as the available data, even when strain histories are complex or multiple biological processes are active simultaneously. The model provides insights into the mechanics of embryonic epithelia and other labile biological tissues, and it sets the stage for future computational studies of whole embryos.     

Local elastodynamic stresses in the unit cell of a woven fabric composite

Summary A theoretical study of the local elastodynamic stresses of woven fabric composites under dynamic loadings is presented in this article. The analysis focuses on the unit cell of an orthogonal woven fabric composite, which is composed of two sets of mutually orthogonal yarns of either the same fiber (nonhybrid fabric) or different fibers (hybrid fabric) in a matrix material. Using the mosaic model for simplifying woven fabric composites and a shear lag approach to account for the inter-yarn deformation, a one-dimensional analysis has been developed to predict the local elastodynamic and elastostatic behavior. The initial and boundary value problems are formulated and then solved using Laplace transforms. Closed form solutions of the dynamic displacements and stresses in each yarn and the bond shearing stresses at the interfaces between adjacent yarns are obtained in the time domain for any type of in-plane impact loadings. When time tends to infinity, the dynamic solutions approach to their corresponding static solutions, which are also developed in this article. Solutions of certain special cases are identical to those reported in the literature. Lastly, the dynamic stresses and bond shearing stresses of plain weave composites subjected to step uniform impacts are presented and discussed as an example of the general analytical model. Received 3 May 1999; accepted for publication 22 September 1999     

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.     

Impact of woven fabric: Experiments and mesostructure-based continuum-level simulations

Woven fabric is an increasingly important component of many defense and commercial systems, including deployable structures, restraint systems, numerous forms of protective armor, and a variety of structural applications where it serves as the reinforcement phase of composite materials. With the prevalence of these systems and the desire to explore new applications, a comprehensive, computationally efficient model for the deformation of woven fabrics is needed. However, modeling woven fabrics is difficult due, in particular, to the need to simulate the response both at the scale of the entire fabric and at the meso-level, the scale of the yarns that compose the weave. Here, we present finite elements for the simulation of the three-dimensional, high-rate deformation of woven fabric. We employ a continuum-level modeling technique that, through the use of an appropriate unit cell, captures the evolution of the mesostructure of the fabric without explicitly modeling every yarn. Displacement degrees of freedom and degrees of freedom representing the change in crimp amplitude of each yarn family fully determine the deformed geometry of the mesostructure of the fabric, which in turn provides, through the constitutive relations, the internal nodal forces. In order to verify the accuracy of the elements, instrumented ballistic impact experiments with projectile velocities of 22-550 m/s were conducted on single layers of Kevlar^® fabric. Simulations of the experiments demonstrate that the finite elements are capable of efficiently simulating large, complex structures.     

A continuum constitutive model for the mechanical behavior of woven fabrics

We propose a new approach for developing continuum models for the mechanical behavior of woven fabrics in planar deformation. We generate a physically motivated continuum model that can both simulate existing fabrics and predict the behavior of novel fabrics based on the properties of the yarns and the weave. The approach relies on the selection of a geometric model for the fabric weave, coupled with constitutive models for the yarn behaviors. The fabric structural configuration is related to the macroscopic deformation through an energy minimization method, and is used to calculate the internal forces carried by the yarn families. The macroscopic stresses are determined from the internal forces using equilibrium arguments. Using this approach, we develop a model for plain weave ballistic fabrics, such as Kevlar®, based on a pin-joined beam geometry. We implement this model into the finite element code ABAQUS and simulate fabrics under different modes of deformation. We present comparisons between model predictions and experimental findings for quasi-static modes of in-plane loading.     

Positive‐definite q‐families of continuous subcell Darcy‐flux CVD(MPFA) finite‐volume schemes and the mixed finite element method

A new family of locally conservative cell‐centred flux‐continuous schemes is presented for solving the porous media general‐tensor pressure equation. A general geometry‐permeability tensor approximation is introduced that is piecewise constant over the subcells of the control volumes and ensures that the local discrete general tensor is elliptic. A family of control‐volume distributed subcell flux‐continuous schemes are defined in terms of the quadrature parametrization q (Multigrid Methods. Birkhauser: Basel, 1993; Proceedings of the 4th European Conference on the Mathematics of Oil Recovery, Norway, June 1994; Comput. Geosci. 1998; 2 :259–290), where the local position of flux continuity defines the quadrature point and each particular scheme. The subcell tensor approximation ensures that a symmetric positive‐definite (SPD) discretization matrix is obtained for the base member (q=1) of the formulation. The physical‐space schemes are shown to be non‐symmetric for general quadrilateral cells. Conditions for discrete ellipticity of the non‐symmetric schemes are derived with respect to the local symmetric part of the tensor. The relationship with the mixed finite element method is given for both the physical‐space and subcell‐space q‐families of schemes. M‐matrix monotonicity conditions for these schemes are summarized. A numerical convergence study of the schemes shows that while the physical‐space schemes are the most accurate, the subcell tensor approximation reduces solution errors when compared with earlier cell‐wise constant tensor schemes and that subcell tensor approximation using the control‐volume face geometry yields the best SPD scheme results. A particular quadrature point is found to improve numerical convergence of the subcell schemes for the cases tested. Copyright © 2007 John Wiley & Sons, Ltd.     

Equivalent mechanical properties of textile monolayers from discrete asymptotic homogenization

The determination of the effective mechanical moduli of textiles from mechanical measurements is usually difficult due to their discrete architecture, which makes micromechanical analyses a relevant alternative to access those properties. Micropolar continuum models describing the effective mechanical behavior of woven fabric monolayers are constructed from the homogenization of an identified repetitive pattern of the textile within a representative unit cell. The interwoven yarns within the textile are represented as a network of trusses connected by nodes at their crossover points. These trusses have extensional and bending rigidities to allow for yarn stretching and flexion, and a transverse shear deformation is additionally considered. Interactions between yarns at the crossover points are captured by beam segments connecting the nodes. The woven fabric is modeled after homogenization as an anisotropic planar continuum with two preferred material directions in the mean plane of the textile. Based on the developed methodology, the effective mechanical properties of plain weave and twill are evaluated, including their bending moduli and characteristic flexural lengths. A satisfactory agreement is obtained between the effective moduli obtained by homogenization and numerical values obtained by finite element simulations performed over periodic unit cells.     

Linear and higher order displacement theories for adhesively bonded lap joints

This work presents an adhesive model for stress analysis of bonded lap joints, which can be applied to model thin and thick adhesive layers. In this theory, linear variations of displacement components along the adhesive thickness are firstly assumed, and the longitudinal strain and the Poisson's effect of the adhesive are modeled. A differential form of the equilibrium equations for the adherends is analytically solved by means of compatible relations of the adhesive deformation. The derived shear and peel stresses are compared with the classical adhesive model of continuous springs with constant shear and peel stresses, and validated with two-dimensional finite element results of the geometrically nonlinear analysis using a commercial package. The numerical results show that the present linear displacement theory can be applied to both thin and moderately thick adhesive layers. The present formulation of the linear displacement theory is then extended to the higher order displacement theory for stress analysis of a thick adhesive, whose numerical results are also compared with those of the finite element computation.     

Effective thermoelastic properties of graded doublyperiodic particulate matrix composites in varying externalstress fields

We consider a linear elastic composite medium, which consists of a homogeneousmatrix containing aligned ellipsoidal uncoated or coated inclusions arranged in a doubly periodicarray and subjected to inhomogeneous boundary conditions. The hypothesis of effective fieldhomogeneity near the inclusions is used. The general integral equation obtained reduces theanalysis of infinite number of inclusion problems to the analysis of a finite number of inclusions insome representative volume element (RVE) . The integral equation is solved by a modifiedversion of the Neumann series; the fast convergence of this method is demonstrated for concreteexamples. The nonlocal macroscopic constitutive equation relating the cell averages of stress andstrain is derived in explicit iterative form of an integral equation. A doubly periodic inclusion fieldin a finite ply subjected to a stress gradient along the functionally graded direction is considered.The stresses averaged over the cell are explicitly represented as functions of the boundaryconditions. Finally, the employed of proposed explicit relations for numerical simulations oftensors describing the local and nonlocal effective elastic properties of finite inclusion pliescontaining a simple cubic lattice of rigid inclusions and voids are considered. The local andnonlocal parts of average strains are estimated for inclusion plies of different thickness. Theboundary layers and scale effects for effective local and nonlocal effective properties as well as foraverage stresses will be revealed.     

Linear and Nonlinear Algorithms for Stress Separation in Photoelasticity

An experimental-numerical hybrid method for the stress separation in photoelasticity is proposed in this study. In the proposed method, boundary conditions for a local finite element model, that is, tractions along boundaries are inversely determined from photoelastic fringes. Two algorithms are proposed for determining the boundary condition. One is a linear algorithm in which the tractions are obtained by the method of linear least-squares from both principal stress difference and principal direction. Another is the nonlinear algorithm in which the tractions are determined only from the principal stress difference. After determining the boundary conditions for the local finite element model, the stresses can be obtained by finite element direct analysis. The effectiveness is demonstrated by applying the proposed method to a perforated plate under tension and contact problems. Results show that the boundary conditions of the local finite element model can be determined from the photoelastic fringes and then the individual stresses can be obtained by the proposed method. Furthermore, the stresses can be evaluated even if the boundary condition is complicated such as at the contact surface. It is expected that the proposed method can be powerful tool for stress analysis.     

三维机织复合材料的弹性性能预报模型

建立了基于等效响应比拟技术的三维机织复合材料弹性性能预报模型．首先将三维机织物的结构单元分解为4个子元(经纱、纬纱、填充纱和接结纱),用几何模型去估算这些子元的体积分数．然后依据不同的外载形式,将复合材料的应力-应变关系等效地表达为3组诸子元所组成的三维弹簧网络．根据刚度系数的物理意义,采用不同的弹簧网络连接形式,并按体积平均化方法获得材料总体刚度矩阵中相应的刚度系数,进而计算得到三维机织复合材料的9个弹性系数．该模型考虑了层内交织经纱、层间交织接结纱的弯曲以及材料内部纯树脂区对三维机织复合材料弹性性能的影响．试验结果与模型的理论预测值进行比较,表明这个模型是有效的。     

Stress recovery algorithm for reduced order models of mechanical systems in nonlinear dynamic operative conditions

Nonlinear forced response analyses of mechanical systems in the presence of contact interfaces are usually performed in built-in numerical codes on reduced order models (ROM). Most of the cases these derive from complex finite element (FE) models, resulting from the high accuracy the designers require in modeling and meshing the components in commercial FE software. In the technical literature several numerical methods are proposed for the identification of the nonlinear forced response in terms of a kinematic quantity (i.e. displacement, velocity and acceleration) associated either to the master degrees-of-freedom retained in the ROM, or to the slave ones after having expanded the reduced response through the reduction matrix. In fact, the displacement is the quantity usually adopted to monitor the nonlinear response, and to evaluate the effectiveness of a partially loose friction interface in damping vibrations, with respect to a linear case where no friction interfaces exist and no energy dissipation can take place. However, when a ROM is used the engineering quantities directly involved in the mechanical design, i.e. the strains and stresses, cannot be retrieved without a further data processing. Moreover, in the case of a strong nonlinear behavior of the mechanical joints, the distributions of the nonlinear strains and stresses over the structure is likely different than the one obtained as a superposition of linear mode shapes whose definition require a-priori assumptions on the boundary conditions at the contact interface. This means that the mentioned approximation cannot be used to predict the safety margins of a structure working in real (nonlinear) operative conditions. This paper addresses this topic and presents a novel stress recovery algorithm for the identification of the strains and stresses resulting from a nonlinear forced response analysis on a ROM. The algorithm is applied to a bladed disk with friction contacts at the shroud joint, which make the behavior of the blades nonlinear and non-predictable by means of standard linear analyses in commercial FE software.

     

Deformation Characteristics of Laminar Composites under Nonlinear Strains

It is shown that the type of governing relations in a composite can change, namely, a laminar composite formed by layers of physically linear materials under nonlinear strains should be described by nonlinear Hooke's law. Local stresses can be not proportional to elastic constants of the layers under nonlinear strains.     

薄壁细胞受微吸管与探针共同作用接触模型及数值模拟

建立了以典型的薄壁球型植物细胞为原型的细胞、微吸管及探针接触模型.模型的细胞壁采用封闭球形薄膜,其本构关系为体积不可压超弹性,膜球内充满有压流体以模拟细胞质.应用轴对称几何非线性方法得出了基本微分方程组,并应用龙格-库塔法进行了求解;同时,应用流固耦合有限元法进行了数值模拟以资比较.两种方法得出了较为一致的变形和应力分...