A numerical model for reinforced concrete structures

This paper describes a two-dimensional approach to model fracture of reinforced concrete structures under (increasing) static loading conditions. The first part is dedicated to the concrete material. The concrete is described in compression by a non-local isotropic damage constitutive law. In tension, a fictitious crack/crack band model is proposed. The influence of biaxial stress states is incorporated in the constitutive relations. In the second part a bond model is described. It accounts for different failure mechanisms, a pullout failure and a splitting failure. This approach is applied to prestressed concrete beams with different failure mechanisms. The numerical results are compared to experimental data and show good agreement.     

准脆性材料的弹塑性各向异性损伤耦合本构模型研究

针对准脆性材料的非线性特征：强度软化和刚度退化、单边效应、侧限强化和拉压软化、不可恢复变形、剪胀及非弹性体胀,在热动力学框架内,建立了准脆性材料的弹塑性与各向异性损伤耦合的本构关系。对准脆性材料的变形机理和损伤诱发的各向异性进行了诠释,并给出了损伤构形和有效构形中各物理量之间的关系。在有效应力空间内,建立了塑性屈服准则、拉压不同的塑性随动强化法则和各向同性强化法则。在损伤构形中,采用应变能释放率,建立了拉压损伤准则、拉压不同的损伤随动强化法则和各向同性强化法则。基于塑性屈服准则和损伤准则,构建了塑性势泛函和损伤势泛函,并由正交性法则,给出了塑性和损伤强化效应内变量的演化规律,同时,联立塑性屈服面和损伤加载面,给出了塑性流动和损伤演化内变量的演化法则。将损伤力学和塑性力学结合起来,建立了应变驱动的应力-应变增量本构关系,给出了本构数值积分的要点。以单轴加载-卸载往复试验识别和校准了本构材料常数,并对单轴单调试验、单轴加载-卸载往复试验、二轴受压、二轴拉压试验和三轴受压试验进行了预测,并与试验结果作了比较,结果表明,所建本构模型对准脆性材料的非线性材料性能有良好的预测能力。     

Dissymétrie de comportement élastique anisotrope couplé ou non à l'endommagement

A thermodynamic framework is proposed to model anisotropic elasticity different in tension and in compression. Based on Kelvin decomposition of the compliance tensor, it applies to any 3D loading. Coupling with damage is made considering fourth- and second-order damage tensors. The proposed formulation automatically satisfies the continuity of the stress tensor and of the energy release rate tensor. The particular case of initially isotropic materials is exposed.     

Modelling the slight compressibility of anisotropic soft tissue

In order to avoid the numerical difficulties in locally enforcing the incompressibility constraint using the displacement formulation of the Finite Element Method, slight compressibility is typically assumed when simulating transversely isotropic, soft tissue. The current standard method of accounting for slight compressibility of hyperelastic soft tissue assumes an additive decomposition of the strain-energy function into a volumetric and a deviatoric part. This has been shown, however, to be inconsistent with the linear theory for anisotropic materials. It is further shown here that, under hydrostatic tension or compression, a transversely isotropic cube modelled using this additive split is simply deformed into another cube regardless of the size of the deformation, in contravention of the physics of the problem. A remedy for these defects is proposed here: the trace of the Cauchy stress is assumed linear in both volume change and fibre stretch. The general form of the strain-energy function consistent with this model is obtained and is shown to be a generalisation of the current standard model. A specific example is used to clearly demonstrate the differences in behaviour between the two models in hydrostatic tension and compression.     

基于临界状态的砂土本构模型研究

砂土孔隙比及所受压力是其力学特性的重要影响因素. 本文基于砂土临界状态线特性分析,采用以e-(p/p_a)ξ平面内的线性关系描述其等向压缩线. 通过对比分析两种不同压缩线函数 与临界状态线函数之间的关系提出更适合描述砂土在等向压缩下的参考压缩线,并给出了基于参考压缩线的等向硬化规律. 建议了适用于 描述砂土剪切特性的屈服面函数,并给出利用等向压缩和等p路径确定屈服面形状参数μ的方法. 将不同应力比对应的压缩线作为砂土状态参量参考线,以获取潜在强度M_f与特征状态应力比M_c,进而描述砂土压缩与剪切特性;基于等向压缩与等p路径建立了当前应力比与状态参量参考线之间的相关关系,从而实现了砂土状 态参量参考线由参考压缩线向临界状态线平稳过渡. 建立的砂土本构模型共11个参数,均能够通过常规土工试验或经验获取. 基于模型预测与Toyoura砂的等向压缩、三轴不排水剪切试验及排水剪切试验的对比结果,本文建立的砂土本构模型很好地描述了Toyoura 砂在不同孔隙比和不同压力下的压缩与剪切特性.      

混凝土的弹塑性损伤双面本构模型

针对混凝土材料拉压应变空间下损伤机制的不同,结合连续损伤力学和塑性理论建立了一个全新的本构模型。该模型中损伤和塑性变形的演变由应变空间的同一个非弹性曲面来控制,但对拉压应变空间中非弹性曲面的演变分别采用了随动强化法则和各向同性演化规律。计算结果表明,该模型能较好地描述混凝土材料在单轴及多轴单调加载和低周反复荷载下的典型非线性特征。     

Formulation and numerical implementation of a constitutive law for concrete with strain-based damage and plasticity

A triaxial constitutive law for concrete within the framework of isotropic damage combined with plasticity is proposed in this paper. It covers typical characteristics of concrete like non-linear uniaxial compression and tension, bi- and triaxial failure criteria and dilatancy with a unified strain-based approach. Thus, this model is quite simple and especially suitable for strain-driven methods like common finite elements. It is complemented with a regularization method based on the crack band approach. A further issue is discussed with procedures for the model parameter determination for a wide range of concrete grades. The application of the model is demonstrated with typical benchmark tests for plain concrete.     

Multiaxial yield surface of transversely isotropic foams: Part I—Modeling

A new yield criterion is proposed for transversely isotropic solid foams. Its derivation is based on the hypothesis that the yielding in foams is driven by the total strain energy density, rather than a completely phenomenological approach. This allows defining the yield surface with minimal number of parameters and does not require complex experiments. The general framework used leads to the introduction of new scalar measures of stress and strain (characteristic stress and strain) for transversely isotropic foams. Furthermore, the central hypothesis that the total strain energy density drives yielding in foams ascribes to the characteristic stress an analogous role of von Mises stress in metal plasticity. Unlike the overwhelming majority of yield models in literature the proposed model recognizes the tension–compression difference in yield behavior of foams through a linear mean stress term. Predictions of the proposed yield model are in excellent agreement with the results of uniaxial, biaxial and triaxial FE analyses implemented on both isotropic and transversely isotropic Kelvin foam models.     

混凝土单轴受拉的非局部本构模型

混凝土受拉本构行为存在很强的局部软化现象,使得单轴受拉试验无法给出应力-应变关系,而只能给出应力-位移关系。本文根据内变量理论和等效应变假设建立了基于真实应变的混凝土单轴受力本构方程,并根据Weibull分布可以描述混凝土等脆性材料断裂过程的试验现象,建立了关于弹性应变的损伤演化规律。然后,通过假设平均应变与真实弹性应变的函数关系,在应力-平均应变的本构关系中采用平均弹性应变以描述其非局部行为,而在材料的损伤演化规律中采用真实弹性应变以描述其局部行为,由此建立了单轴受拉荷载条件下的非局部本构模型。最后,对一个单调受拉试验和一个反复受拉试验的仿真结果表明所提出的非局部本构模型可以准确地模拟试验结果。     

Analysis of creep deformation and creep damage in thin-walled branched shells from materials with different behavior in tension and compression

A constitutive model for describing the creep and creep damage in initially isotropic materials with different properties in tension and compression has been applied to the modeling of creep deformation and creep damage growth in thin-walled shells of revolution with the branched meridian. The approach of establishing the basic equations for axisymmetrically loaded branched shells under creep deformation and creep damage conditions has been introduced. To solve the initial/boundary-value problem, the fourth-order Runge–Kutta–Merson’s method of time integration with the combination of the numerically stable Godunov’s method of discrete orthogonalization is used. The solution of the boundary value problem for the branched shell at each time instant is reduced to integration of the series of systems of ordinary differential equations describing the deformation of each branch and the shell with basic meridian. Some numerical examples are considered, and the processes of creep deformation and creep damage growth in a shell with non-branched meridian as well as in a branched shell are analyzed. The influence of the tension–compression asymmetry on the stress–strain state and damage evolution in a shell with non-branched meridian as well as in a branched shell with time are discussed.     

Modeling stress state dependent damage evolution in a cast Al–Si–Mg aluminum alloy

Internal state variable rate equations are cast in a continuum framework to model void nucleation, growth, and coalescence in a cast Al–Si–Mg aluminum alloy. The kinematics and constitutive relations for damage resulting from void nucleation, growth, and coalescence are discussed. Because damage evolution is intimately coupled with the stress state, internal state variable hardening rate equations are developed to distinguish between compression, tension, and torsion straining conditions. The scalar isotropic hardening equation and second rank tensorial kinematic hardening equation from the Bammann–Chiesa–Johnson (BCJ) Plasticity model are modified to account for hardening rate differences under tension, compression, and torsion. A method for determining the material constants for the plasticity and damage equations is presented. Parameter determination for the proposed phenomenological nucleation rate equation, motivated from fracture mechanics and microscale physical observations, involves counting nucleation sites as a function of strain from optical micrographs. Although different void growth models can be included, the McClintock void growth model is used in this study. A coalescence model is also introduced. The damage framework is then evaluated with respect to experimental tensile data of notched Al–Si–Mg cast aluminum alloy specimens. Finite element results employing the damage framework are shown to illustrate its usefulness.     

复合材料层合板面内渐进损伤分析的CDM模型

基于连续介质损伤力学,提出了一个预测复合材料层合板面内渐进损伤分析的模型,它包括损伤表征、损伤判定和损伤演化3 部分. 模型能够区分纤维拉伸断裂、纤维压缩断裂、纤维间拉伸损伤和纤维间压缩损伤4 种损伤模式,定义了与4 个损伤模式对应的损伤状态变量,导出了材料主轴系下损伤前后材料本构之间的关系. 损伤起始采用Puck 准则判定,损伤演化由特征长度内应变能释放密度控制. 假定材料服从线性应变软化行为,建立了损伤状态变量关于断裂面上等效应变的渐进损伤演化法则. 模型涵盖了复合材料面内损伤起始、演化直至最终失效的全过程. 完成了含孔[45/0/-45/90]_2S 层合板在拉伸和压缩载荷下失效分析,结果表明该模型能合理进行层合板的强度预测和损伤失效分析.      

复合材料层合板面内渐进损伤分析的CDM模型简

基于连续介质损伤力学,提出了一个预测复合材料层合板面内渐进损伤分析的模型,它包括损伤表征、损伤判定和损伤演化3 部分. 模型能够区分纤维拉伸断裂、纤维压缩断裂、纤维间拉伸损伤和纤维间压缩损伤4 种损伤模式,定义了与4 个损伤模式对应的损伤状态变量,导出了材料主轴系下损伤前后材料本构之间的关系. 损伤起始采用Puck 准则判定,损伤演化由特征长度内应变能释放密度控制. 假定材料服从线性应变软化行为,建立了损伤状态变量关于断裂面上等效应变的渐进损伤演化法则. 模型涵盖了复合材料面内损伤起始、演化直至最终失效的全过程. 完成了含孔[45/0/-45/90]_2S 层合板在拉伸和压缩载荷下失效分析,结果表明该模型能合理进行层合板的强度预测和损伤失效分析.     

Effect of a characteristic length on crack spacing in a reinforced concrete bar under tension

The paper presents quasi-static FE-simulations of the crack formation in a reinforced concrete bar without stirrups subject to tension. The material was modeled with a continuum smeared crack model using an elasto-plastic constitutive law. A linear Rankine criterion with isotropic softening and associated flow rule was adopted in a tensile regime. To ensure the mesh-independency, the softening parameter was enhanced by a characteristic length of micro-structure by means of a non-local theory. Attention was laid to the effect of a different characteristic length of micro-structure and initial bond-slip stiffness on the spacing of localized zones.     

Application of some anisotropic damage model to the prediction of the failure of some complex industrial concrete structure

The aim of this paper is to present a new model for the anisotropic damage of concrete and apply it to the numerical simulation of the failure of some complex structure. This model considers damage as anisotropic in tension but isotropic in compression, and incorporates asymmetry between tension and compression. In spite of its complexity, it is expressed in a format fit for numerical calculations by the finite element method, and involves only six material parameters (in addition to the usual elastic constants), which makes it suitable for applications to large industrial structures. The material parameters can be identified from simple uniaxial and biaxial experiments. The application envisaged relates to the failure of some cylindrical nuclear containment vessel subjected to some excessive internal pressure. This structure is made of reinforced concrete containing both passive (initially stress-free) steel armatures and pre-stressed cables. The numerical results illustrate the importance of incorporation of both asymmetry between tension and compression and anisotropy of damage in the simulation to accurately predict the resistance of the structure to fracture.     

The onset and progression of damage in isotropic paper sheets

An experimental investigation is performed and analyzed in order to examine the onset and evolution of damage processes in thin isotropic paper sheets made of mechanical pulp. A microscopy technique has been used to estimate the relative fraction of bond and fibre breaks. It has been found that the active damage mechanism is bond failure, hence supporting the assumption of an isotropic scalar valued damage variable.All experiments have been performed by simultaneous with the mechanical loading monitoring the acoustic emission activity. Three different experimental setups have been designed offering the possibility to analyze the influence of stress gradients, as well as different levels of the ratios between the in-plane normal stresses, on the onset of damage. It is concluded that stress gradients in the paper specimens have a large influence on the onset of damage. When stress gradients are present a non-local theory has to be used in the analysis. In this way compliance with an isotropic damage criterion is achieved. The characteristic length, determining the gradient sensitivity, has been found to be of the same order of magnitude as some average fibre length.To study the evolution of the damage processes, wide and short specimens have been loaded in tension resulting in stable damage processes. With the assumptions made regarding the mechanical behavior of the paper material after onset of damage, the damage and the cumulative number of acoustic events curve correlates very well. The experimentally obtained data is used to determine material parameters in a proposed damage evolution law. It is found that the assumed damage evolution law can, for isotropic paper materials with bond rupture as the prevalent failure mechanism, be further simplified as only one specific material dependent damage evolution parameter has to be determined in experiments.     

A PVC-foam material model based on a thermodynamically elasto-plastic-damage framework exhibiting failure and crushing

A well known foam for naval sandwiches is PVC (polyvinyl chloride) foam. This foam exhibits elasto-damage behavior under tension and elasto-plastic behavior under compression. A proper material model is required for the prediction of the failure and post-failure behavior of these sandwiches during (indirect) underwater explosion loading and blast loading. As material models available did not have combined elasto-plastic and elasto-damage behaviors, a new model needed to be developed. To this end a general thermodynamically consistent framework was used and damage and plastic functional were derived describing the PVC foam behavior. The damage functional contains fracture mechanics properties. Special attention was paid to the compression–tension load reversal of the material model, that may occur during underwater shock loading.     

On non-local and non-homogeneous elastic continua

A thermodynamic framework endowed with the concept of non-locality residual is adopted to derive non-local models of integral-type for non-homogeneous linear elastic materials. Two expressions of the free energy are considered: the former yields a one-component non-local stress, the latter leads to a two-component local–non-local stress since the stress is expressed as the sum of the classical local stress and of a non-local component identically vanishing in the case of constant strains. The attenuation effects are accounted for by a symmetric space weight function which guarantees the constant strain requirement as well as the dual constant stress condition everywhere in the body. The non-local and non-homogeneous elastic structural boundary-value problem under quasi-static loads is addressed in a geometrically linear range. The complete set of variational formulations for the structural problem is then provided in a unitary framework. The solution uniqueness of the non-local structural model is proved and the non-local FEM is addressed starting from the non-local counterpart of the total potential energy. Numerical applications are provided with reference to a non-homogeneous bar in tension using the Fredholm integral equation and the non-local FEM. The solutions show no pathological features such as numerical instability and mesh sensitivity for degraded bar conditions.     

Theory of creep deformation with kinematic hardening for materials with different properties in tension and compression

A constitutive model for creep deformation that describes the loading-history-dependent behavior of initially isotropic materials with different properties in tension and compression under stress vector rotations limited by 50–60° is presented within a thermodynamic framework. In the proposed constitutive model a kinematic hardening rule is adopted. This model also introduces an effective equivalent stress in the creep potential that is based on the first and second invariants of the effective stress tensor, and on the joint invariant of the effective stress tensor and eigenvector associated with the maximum principal Cauchy stress. The formulation of the kinematic hardening rule is presented and discussed. All the material parameters in the model have been obtained from a series of proposed basic experiments with constant stresses. These model parameters are then used to predict the creep deformation of the aluminum alloy under multiaxial loading with constant stresses, and under non-proportional uniaxial and non-proportional multiaxial loadings for both isothermal and nonisothermal processes.