复合材料蠕变本构关系及实验测定Ⅰ.本构关系

蠕变是复合材料最重要的力学性能之一,实验表明:复合材料在蠕变条件下的变形可以分为弹性变形、粘弹性变形和粘塑性变形.应用不可逆过程的热力学和广义变量的概念可以分析材料的蠕变变形.本文首先回顾了热力学的基本方程;基于Schapery本构关系的假设和思路推导了蠕变本构关系的一般形式,其中包括弹性变形、粘弹性变形和粘塑性变形;考虑到广义力选取的不唯一性,本文提出了广义力选取的原则以使得到的本构关系尽可能地简单;由此本文给出了复合材料的一维蠕变,各向同性复合材料的二维蠕变和纤维增强复合材料平面内的蠕变的本构关系.     

Einfache Stoffgleichungen für mehrachsiges Kriechen nichtlinearviskoelastischer Werkstoffe

Simple Constitutive Equations for Multiaxial Creep of Non-linear Viscoelastic Materials Isochronous creep data in uniaxial tension, compression and shear were used to formulate and test constitutive equations for nonlinear viscoelastic materials with different creep response in simple tension and compression.     

A nonlinear constitutive relationship for asphalt binders

A nonlinear constitutive relationship was developed for asphalt binders. Two binders, one polymer modified and one unmodified, were tested in shear using creep and recovery loading. Five different stress levels and four loading times were considered, to capture the response of the binders in the linear and nonlinear viscoelastic range. The creep response of the binders was successfully modeled with a nonlinear power law function. The modified superposition principle was unable to predict the recovery phase of the testing data. A nonlinear constitutive relationship composed of a nonlinear viscous part plus a linear viscoelastic part was developed. The constitutive relationships successfully predicted the binders’ response in creep and recovery. The predictions of the constitutive relationships matched accurately the response of the binders subjected to the Multiple Stress Creep and Recovery loading pattern.     

Micro-Macro Analysis of Viscoelastic Unidirectional Laminated Composite Plates Using DR Method

The Dynamic Relaxation (DR) technique together with finite difference discritization is used to study the bending behavior of Mindlin composite plate including geometric nonlinearity. The overall behavior of the unidirectional composite is obtained from a three-dimensional (3D) micromechanical model, in any combination of normal and shear loading conditions, based on the assumptions of Simplified Unit Cell Method (SUCM). The composite system consists of nonlinear viscoelastic matrix reinforced by transversely isotropic elastic fibers. A recursive formulation for the hereditary integral of the Schapery viscoelastic constitutive equation in multiaxial stress state is used to model the nonlinear viscoelastic matrix material in the material level. The creep tests data is used for verification of the predicted response of the current approach. Under uniform lateral pressure, the laminated plate deformation with clamped and hinged edged constraints is predicted for various time steps.     

An elastic-viscoplastic model for large deformation

It is shown that the theory of an elastic-viscoplastic work hardening material proposed by Bodner and Partom for small deformations may be generalized for large deformations by reformulating the equations using Lagrangian quantities. Restrictions on the general constitutive equations were obtained using the thermodynamic procedures proposed by Green and Naghdi. In this formulation the stress is determined directly from deformation quantities and in particular is not calculated using a hypoelastic type equation for a stress rate. Also, since Lagrangian quantities are used there is no need to introduce special rates like the Jaumann rate in the evolution equations. Specific constitutive equations were proposed for a material exhibiting isotropic-elastic response in its reference configuration, strain-rate and temperature dependent plastic flow with isotropic and directional hardening, and thermal recovery of hardening. These specific equations use only the material constants obtained from the corresponding small deformation theory. Examples of simple tension and simple shear show that these equations predict physically plausible material response for large deformations.     

Viscoplastic response and creep failure time prediction of polymers based on the transient network model

The nonlinear viscoelastic/viscoplastic response of polymeric materials is described by a new model based on previous works in terms of monotonic loading, stress–relaxation, and creep. In the proposed analysis, following a constitutive equation of viscoelasticity, based on the transient network theory, essential modifications are introduced, which account for the nonlinearity and viscoplasticity at small elastic and finite plastic strain regime. In addition, viscoplastic response is successfully analyzed by a proper kinematic formulation, which is combined with a functional form of the rate of plastic deformation. A three-dimensional constitutive equation is then derived for an isotropic incompressible medium. This analysis is capable of capturing the main aspects of inelastic response and the instability stage taking place at the tertiary creep, related to the creep failure. Model simulations described successfully the experimental data of polypropylene, which were performed elsewhere.     

Branching of strain histories for nonlinear viscoelastic solids with a strain clock

Summary An important class of constitutive equations for nonlinear viscoelastic response utilizes the concept of a strain clock. The clock takes the form of a material time variable which is defined in terms of the strain history and which increases faster than physical time. Important consequences of the strain clock are that stress relaxation and creep occur faster as strain increases, and the stress may not increase monotonically with time. In this work, we discuss whether this non-monotonic response implies that strain histories may branch into multiple histories.     

基于不可逆内变量热力学的岩石材料粘弹-粘塑性本构方程

岩石材料的粘弹性和粘塑性变形是与时间相关的能量耗散行为。在Rice不可逆内变量热力学框架下,引入两组内变量分别用来描述在粘弹性和粘塑性变形过程中材料的内部结构调整。通过给定比余能的具体形式和内变量的演化方程,推导出内变量粘弹-粘塑性本构方程。粘弹性本构方程具有普遍性,能涵盖Kelvin-Voigt和Poynting-Thomson在内的经典粘弹性模型的本构方程。并指出热力学力与应力呈线性关系是组合元件模型为线性模型的根本原因。粘塑性本构方程能较好地刻画岩石材料在粘塑性变形过程中的硬化现象。对模拟岩石的模型相似材料进行单轴加卸载蠕变试验,将蠕变过程中的粘弹性和粘塑性变形分离并根据试验数据对本构方程的材料参数进行辨识。试验数据和理论曲线对比结果表明该文提出的本构方程能很好地模拟材料的蠕变行为。该类型的本构方程能为岩石工程的长期稳定性的预测、评价以及加固分析提供基础。     

Creep buckling and post-buckling analyses for a hybrid laminated viscoelastic FGM cylindrical shell under in-plane loading

In this paper, creep buckling and post-buckling of a hybrid laminated viscoelastic functionally graded material (FGM) cylindrical shell under in-plane loading are investigated. Considering the high-order transverse shear deformation and geometric nonlinear theory, the von Karman geometric relation of the hybrid laminated viscoelastic FGM cylindrical shell with initial deflection is established. Based on the Donnell theory, elastic piezoelectric theory and Boltzmann superposition principle, nonlinear creep governing equations of the hybrid laminated viscoelastic FGM cylindrical shell under in-plane loading are derived. By means of the finite difference method and the Newton–Newmark method, the problem for creep buckling and post-buckling of the laminated shell’s structure is solved. Numerical results are presented to show effects of geometric parameters, power law index and loading on creep buckling and post-buckling of the hybrid laminated viscoelastic FGM cylindrical shell.     

Zur exakten und näherungsweisen Berechnung unidirektionalverstärkter Kunststoffe

Exact and Approximate Calculations of UD-Layers It is shown that the quasi-elastic solution of layers nearly almost used in practical viscoelastic calculations is an excellent approach. The exact solution achieved with simple integral equations are compared with the quasi-elastic solution and the occuring deviations are described. The creep functions of layers are given the form of the uniaxial creep function of the isotropic viscoelastic material. This makes the creep behaviour of layers very clear and symplifies practical calculations. The given solutions are under no restriction with regard to the chosen form of the uni-axial creep function.     

Nonlinear viscoelastic analysis of orthotropic beams using a general third-order theory

The displacement based finite element model of a general third-order beam theory is developed to study the quasi-static behavior of viscoelastic rectangular orthotropic beams. The mechanical properties are considered to be linear viscoelastic in nature with a scope to undergo von Kármán nonlinear geometric deformations. A differential constitutive law is developed for an orthotropic linear viscoelastic beam under the assumptions of plane-stress. The fully discretized finite element equations are obtained by approximating the convolution integrals using a trapezoidal rule. A two-point recurrence scheme is developed that necessitates storage of data from the previous time step only, and not from the entire deformation history. Full integration is used to evaluate all the stiffness terms using spectral/hp lagrange polynomials. The Newton iterative scheme is employed to enhance the rate of convergence of the nonlinear finite element equations. Numerical examples are presented to study the viscoelastic phenomena like creep, cyclic creep and recovery for thick and thin beams using classical mechanical analogues like generalized n-parameter Kelvin-Voigt solids and Maxwell solids.     

Viscoelastic response of thin membranes with application to red blood cells

The rate-type constitutive analysis of viscoelastic response of thin membranes, which includes an instantaneous elastic response and viscous behavior in both shear and dilatation, is developed with the aim to study the mechanical response of red blood cells. A convenient set of generalized stress and strain variables is introduced, which facilitates the derivation and integration of the governing differential equations. Gradual or sudden loading and stepwise unloading histories are considered. The performed parametric study of the mechanical response illustrates the effects of the introduced material parameters on the coefficient of viscoelastic lateral contraction and the overall membrane deformation. A closed form solution to the problem of radial stretching of a viscoelastic hollow circular membrane is derived without referral to the correspondence principle, which is of interest for the micropipette aspiration experiment of the red blood cell. The effects of the material parameters on the instantaneous elastic response and the subsequent rate of creep are discussed.     

Applications of tensor functions to the formulation of constitutive equations involving damage and initial anisotropy

In this paper some results of the tensor function theory are applied to the formulation of constitutive equations of isotropic and anisotropic materials in the secondary and tertiary creep stage. The creep process, in its tertiary phase, is characterized by a damage tensor. Because of its microscopic nature, damage has, in general, an anisotropic character even in cases where the material was originally isotropic, i.e. isotropic in its virgin state. Fissure orientation and length cause anisotropic macroscopic behaviour. In the first part of the paper some possible ways of representing constitutive equations involving (initial) anisotropy of the material (e.g. from rolling) and involving anisotropic creep-damage are dealt with. The formulations of such equations are based upon theorems concerning tensor-valued functions. Furthermore, some simplified constitutive equations for more practical use are discussed. The main problem of this part is: to find an irreducible set of tensor generators. Besides the problem of finding such tensor generators it is very important to determine the scalar coefficients in constitutive equations as functions of the invariants and experimental data. The second part of the paper is concerned with the determination of the scalar functions. This can be done by using tensorial interpolation methods as pointed out in detail.     

Constitutive equations for ligament and other soft tissue: evaluation by experiment

Ligaments, tendons and other soft tissues are nonlinearly viscoelastic. To discriminate among various constitutive equations which may be used to describe the tissue, appropriate experimental modalities are requisite. Ideally, testing should span physiologic ranges for load (or strain), load history (recovery and reloading), and load onset and duration, and a robust model will fit all data. Methods to expand the experimental window of time for relaxation and creep are presented and evaluated. The role of ramp, relaxation and recovery protocols is studied in the context of viscoelasticity describable by linear, quasi-linear, nonlinear superposition, Schapery, and multiple integral formulations. The advantages associated with testing protocols that expand the time windows for creep or relaxation are presented.     

A ten node tetrahedral Cosserat Point Element (CPE) for nonlinear isotropic elastic materials

It is known that the standard full integration ten node tetrahedral element is inaccurate for thin nearly incompressible structures. Also, the commercial code ABAQUS recommends replacing this standard element (C3D10) with an undocumented patented modified element (C3D10M) for contact problems. The objective of this work is to develop a ten node tetrahedral Cosserat Point Element (CPE) for nonlinear isotropic hyperelastic materials. Hyperelastic constitutive equations for the CPE are developed by treating the element as a structure with a strain energy function that is restricted to satisfy a nonlinear form of the patch test. A number of examples are considered which demonstrate that the resulting CPE is accurate and robust, it does not exhibit the numerical stiffness for nearly incompressible materials observed for (C3D10) nor the unphysical instabilities observed for (C3D10M). Moreover, the CPE can be used for thin structures and three-dimensional bodies with a smooth transition from compressible to nearly incompressible material behavior.     

Bending and creep buckling response of viscoelastic functionally graded beam-columns

The viscoelastic creep response of flexural beams and beam-columns made with functionally graded materials is numerically investigated. The paper highlights the challenges associated with the modeling and analysis of such structures, and presents a nonlinear theoretical model for their bending and creep buckling analysis. The model accounts for the viscoelasticity of the materials using differential-type constitutive relations that are based on the linear Boltzmann’s principle of superposition. The model is general in terms of its ability to deal with any material volume faction distribution through the depth of the beam, and with different linear viscoelastic laws, boundary conditions, and loading schemes. The governing equations are solved through time stepping numerical integration, which yields an exponential algorithm following the expansion of the relaxation function into a Dirichlet series. A numerical study that examines the capabilities of the model and quantifies the creep response of functionally graded beam-columns is presented, with special focus on the stresses and strains redistribution over time and on the creep buckling response. The results show that the creep response of such structures can be strongly nonlinear due to the variation of the viscoelastic properties through the depth, along with unique phenomena that are not observed in homogenous structures.     

Boundary element method for finite elasticity

This paper describes some integral formulations and implementations of a Boundary Element Method to solve two- and three-dimensional finite deformation problems of rubber-like materials. The integral equations are formulated in terms of unknown incremental displacement and total boundary traction fields, or alternatively in terms of the incremental displacement and incremental boundary traction fields. The elastic material is either compressible or incompressible with given constitutive equations. Both formulations are implemented and tested. The uniaxial elongation and simple shear deformations of a model material are successfully simulated by both formulations. Some non-trivial examples are performed using the first formulation.     

Natural frequency and critical speed determination of an axially moving viscoelastic beam

In this paper, the natural frequency and critical speed of an axially moving viscoelastic beam with clamped and simple supports are calculated analytically based on the Euler–Bernoulli and Timoshenko theories. The beam is incompressible in bulk and viscoelastic in shear, which obeys the linear standard solid model with material time derivative. The axial speed is characterized by a simple harmonic variation about a constant mean speed. By defining some dimensionless parameters, the governing equations are derived from the Newton’s second law. They contain two coupled partial differential equations with time depended coefficients. The straightforward method in perturbation theory is used to solve these equations. By considering the homogeneous equation, the natural frequencies are calculated. The critical speed is determined by a constant speed assumption. By a parametric study, the effects of mechanical and geometrical parameters on the natural frequency and critical speed are investigated.     

Higher-order fractional constitutive equations of viscoelastic materials involving three different parameters and their relaxation and creep functions

Two higher-order fractional viscoelastic material models consisting of the fractional Voigt model (FVM) and the fractional Maxwell model (FMM) are considered. Their higher-order fractional constitutive equations are derived due to the models’ constructions. We call them the higher-order fractional constitutive equations because they contain three different fractional parameters and the maximum order of equations is more than one. The relaxation and creep functions of the higher-order fractional constitutive equations are obtained by Laplace transform method. As particular cases, the analytical solutions of standard (integer-order) quadratic constitutive equations are contained. The generalized Mittag–Leffler function and H-Fox function play an important role in the solutions of the higher-order fractional constitutive equations. Finally, experimental data of human cranial bone are used to fit with the models given by this paper. The fitting plots show that the models given in the paper are efficient in describing the property of viscoelastic materials.     

Modelling creep behaviour of orthotropic composites by the coincident element method

In this paper, the coincident method proposed previously is applied to model the four-point-bending creep experiments conducted at the Cooperative Research Centre for Advanced Composite Structures (CRC-ACS) on carbon-epoxy composite laminates. A parameter identification methodology is first developed to determine the elastic and viscoelastic material models to be used for a coincident element. Simulations are then conducted to model the flexural creep response of the composite laminates under different loading and temperature conditions. The predicted results are in reasonably good agreement with those obtained by experiments. It is demonstrated that the coincident element method is a relatively simple and useful tool for modelling orthotropic and viscoelastic response of laminated composites by using a finite element package that only supports isotropic viscoelastic material models.