基于对混凝土破损过程的模拟研究

从细观力学层次出发,把混凝土看作是由骨料、砂浆基质及两者之间的界面过渡区组成的三相复合材料,通过引入细观损伤变量及损伤变量的演化规律模拟了混凝土的破坏过程。首先生成骨料位置随机分布且颗粒级配符合Fuller曲线的数值试件,然后对每一相材料分别赋以服从Weibull分布的随机力学参数(弹模、强度和泊松比等)并规定相应的破坏准则,将通过上述途径获得的混凝土数值试件导入有限元软件Abaqus进行数值模拟。与实验结果的对比显示,提出的随机损伤数值模型能够较好地模拟混凝土板式试件在单拉、单压下的应力-应变关系曲线以及混凝土中裂纹的萌生、开裂过程。     

混凝土材料与结构破坏过程模拟分析

采用离散单元法对混凝土材料和混凝土结构破坏机理进行分析。在细观尺度上将混凝土材料视为由粗骨料、水泥砂浆及界面过渡区三相组成,建立了混凝土材料的离散元模型;在宏观尺度上将混凝土视为均质材料建立了混凝土结构离散单元模型。计算分析结果表明:细观尺度上的二维离散单元模型可以用来很好地模拟混凝土材料的单轴受力破坏过程,但不能很好地模拟复合受力状态下的混凝土材料的破坏;宏观尺度上的离散单元模型可以很好地模拟钢筋混凝土构件的破坏过程,但模拟结果对单元的形状有较大的依赖性;宏观尺度上的离散单元模型可以很好地模拟结构的倒塌过程,但计算效率有待提高。     

考虑过渡区界面影响的混凝土宏观力学性质研究

混凝土材料的宏观力学特性及破坏机理由其细观组分来决定,界面过渡区是影响混凝土断裂破坏路径及宏观力学特性的重要因素。认为界面过渡区是区别于远处砂浆基质的一层含较高孔隙率的近场砂浆材料,采用“两步等效法”得到了混凝土细观单元的等效本构关系模型。最后基于细观单元等效化方法分析了在单轴拉伸、单轴压缩及弯拉载荷条件下混凝土试件的破坏过程及宏观力学性质,探讨了界面过渡区对混凝土力学特性的影响,并与随机骨料模型分析结果进行了对比。结果表明:界面相的存在对混凝土的弹性模量、强度及残余强度等力学性质有很大影响,在对混凝土宏观力学特性及细观断裂破坏过程进行研究时不可忽略其影响。     

力学损伤混凝土热传导行为细观研究

为研究力学损伤对温度传导行为的影响,考虑混凝土细观结构的非均质性,将其视为由骨料、砂浆和界面过渡区组成的三相复合材料,提出并建立了力学-热学单向耦合作用研究的细观数值方法与模型。方法的思想是:首先对非均质混凝土在荷载作用下开裂损伤力学行为模拟分析,获得其内部的损伤分布情况;然后,将力学分析结果作为初始输入条件,结合复合材料均匀化理论,获得损伤混凝土各细观单元的导热性能,从而实现了力学-热学单向耦合作用的数值模拟。以单轴压缩加载为例,通过与试验结果的对比,验证细观方法的合理性与有效性。基于该套细观数值模拟方法,对荷载作用后混凝土试件的宏观有效导热系数与温度场分布进行了细观计算分析,对比了不同复合材料力学模型的影响,探讨分析了压缩荷载水平的影响规律。     

钢筋混凝土柱偏心受压力学性能的细观数值研究

钢筋混凝土构件的宏观力学性能由其组分-钢筋和混凝土两部分的力学性能决定。结合混凝土细观结构形式,认为混凝土是由骨料颗粒、砂浆基质及界面过渡区组成的复合材料,假定钢筋与混凝土之间完好粘结,基于钢筋混凝土柱偏心受压试验,建立了钢筋混凝土柱偏心受压加载下力学特性及破坏行为研究的细观尺度力学分析模型。通过对混凝土方形和矩形试件进行受压力学特性模拟,采用反演法确定了界面的力学参数,进而模拟了钢筋混凝土柱偏心受压加载下的宏观力学性能。结果表明,相比于宏观尺度模型,细观数值分析模型能够充分体现材料的非均质性,能够较好的模拟试件的宏观力学性能,并且能够细致的描述裂缝发展及试件破坏过程,与试验结果吻合良好。该文建立的细观尺度分析模型与方法,为钢筋混凝土构件层次宏观力学非线性及其尺寸效应研究提供了理论支持。     

基于细观单元等效化方法的混凝土动态破坏行为分析

混凝土材料具有明显的应变率效应,对其力学性质增强机理的认识还不统一。在细观随机骨料模型基础上,采用特征单元尺度划分试件网格,推导了考虑材料拉/压强度应变率效应的细观单元等效本构关系,建立了非均质混凝土材料的细观单元等效化数值模型。基于二维模型对Dilger等混凝土动态压缩试验进行了数值模拟,获得的数值结果与试验数据及随机骨料模型结果吻合良好,证明了细观单元等效化方法的准确性;进而对三维混凝土试件动态单轴拉伸和压缩破坏模式及宏观力学性质的加载速率效应进行了研究。数值结果表明:随着加载速率的增加,混凝土裂纹(损伤)数量增大,混凝土破坏将耗散更多的能量,是混凝土动态强度提高的主要原因。     

骨料级配对二维模型混凝土单轴抗拉强度影响的理论研究

混凝土尺寸效应及其宏观力学非线性根源于其材料细观组成的非均质性。结合混凝土细观结构形式,将混凝土看作由骨料颗粒、砂浆基质及界面过渡区组成的复合材料。采用双线性弹性损伤模型来描述砂浆基质及界面过渡区的力学行为,假定骨料颗粒为弹性体而不发生破坏,进而推导并获得了单轴拉伸条件下不同骨料颗粒级配混凝土断裂裂缝扩展路径长度及其抗拉强度的理论解。最后,对比了建立的理论公式结果与细观尺度数值模拟结果,验证了构建的关于裂缝长度及抗拉强度理论解的准确性和合理性。     

大坝混凝土试件三维细观力学并行计算研究

基于混凝土材料基本特性及其有关试验研究结果,给出了混凝土材料的损伤演化关系、及其抗拉强度和弹性模量的应变率强化定律,建立了既考虑细观损伤弱化又计及应变率强化效应的混凝土非线性有限元静动力学虚功方程。根据FEPG有限元语言的规则,编写了混凝土试件细观静、动力学方程的描述文件及位移、应力场的算法文件,并基于PFEPG软件平台开发了混凝土试件的三维细观数值并行计算程序。利用三维随机骨料随机参数模型对混凝土试件弯折损伤破坏过程进行了细观数值模拟。其计算结果与试验结果相吻合,显示了该算法及所编并行程序的有效性。     

非均质混凝土材料破坏的三维细观数值模拟

在随机骨料模型的基础上,采用特征单元尺度对网格进行剖分,建立了非均质混凝土材料损伤破坏及宏观力学特性研究的三维细观单元等效化模型。对单轴拉伸、单轴压缩条件下湿筛混凝土试件的破坏过程及宏观力学性能进行了模拟分析;研究了混凝土梁的三分点弯拉力学特性,并与平面模型结果作了对比。研究表明:1) 与平面计算模型相比,三维模型更真实地模拟混凝土材料在外荷载作用下的损伤破坏过程,更准确的描述非均质混凝土材料的宏观力学性能,且与实验结果吻合;2) 骨料空间分布形式基本不影响混凝土材料的宏观弹性模量及强度,但影响其破坏过程和破损路径;3) 与随机骨料模型等细观力学方法相比,该方法具有高效性。     

钢筋混凝土梁抗冲击性能研究的细观数值模拟方法

为更真实地描述钢筋混凝土梁在冲击荷载作用下复杂的破坏过程,结合混凝土细观结构非均质性,假定混凝土是由骨料颗粒、砂浆基质以及界面过渡区(ITZ)组成的三相复合材料,考虑混凝土细观组分的率效应,将钢筋网嵌入至素混凝土梁,建立了钢筋混凝土梁抗冲击力学行为研究的三维细观尺度数值模型。基于该数值分析模型与方法,研究了落锤速度对钢筋混凝土梁抗冲击性能的影响规律。将细观数值模拟结果与宏观模型结果以及已有试验结果进行了对比分析,包括破坏模式、冲击力、跨中位移以及支座反力等,发现模拟结果与试验结果吻合良好,验证了细观数值分析模型在钢筋混凝土梁抗冲击性能模拟中运用的可行性。     

Micromechanical investigation on size effect of tensile strength for recycled aggregate concrete using BFEM

In this paper, the validity and performance of base force element method (BFEM) based on potential energy principle was studied by some numerical examples. And the BFEM on damage mechanics is used to analyze the size effect on tensile strength for recycled aggregate concrete (RAC) at meso-level. The recycled aggregate concrete is taken as five-phase composites consisting of natural coarse aggregate, new mortar, new interfacial transition zone (ITZ), old mortar and old ITZ on meso-level. The random aggregate model is used to simulate the meso-structure of recycled aggregate concrete. The size effects of mechanical properties of RAC under uniaxial tensile loading are simulated using the BFEM on damage mechanics. The simulation results agree with the test results. This analysis method is the new way for investigating fracture mechanism and numerical simulation of mechanical properties for RAC.     

Modelling cohesive crack growth using a two-step finite element-scaled boundary finite element coupled method

A two-step method, coupling the finite element method (FEM) and the scaled boundary finite element method (SBFEM), is developed in this paper for modelling cohesive crack growth in quasi-brittle normal-sized structures such as concrete beams. In the first step, the crack trajectory is fully automatically predicted by a recently-developed simple remeshing procedure using the SBFEM based on the linear elastic fracture mechanics theory. In the second step, interfacial finite elements with tension-softening constitutive laws are inserted into the crack path to model gradual energy dissipation in the fracture process zone, while the elastic bulk material is modelled by the SBFEM. The resultant nonlinear equation system is solved by a local arc-length controlled solver. Two concrete beams subjected to mode-I and mixed-mode fracture respectively are modelled to validate the proposed method. The numerical results demonstrate that this two-step SBFEM-FEM coupled method can predict both satisfactory crack trajectories and accurate load-displacement relations with a small number of degrees of freedom, even for crack growth problems with strong snap-back phenomenon. The effects of the tensile strength, the mode-I and mode-II fracture energies on the predicted load-displacement relations are also discussed.     

Nonlinear fracture mechanics and plasticity of the split cylinder test

The split cylinder test is subjected to an analysis combining nonlinear fracture mechanics and plasticity. The fictitious crack model is applied for the analysis of splitting tensile fracture, and the Mohr-Coulomb yield criterion is adopted for modelling the compressive crushing/sliding failure. Two models are presented, a simple semi-analytical model based on analytical solutions for the crack propagation in a rectangular prismatic body, and a finite element model including plasticity in bulk material as well as crack propagation in interface elements. A numerical study applying these models demonstrates the influence of varying geometry or constitutive properties. For a split cylinder test in load control it is shown how the ultimate load is either plasticity dominated or fracture mechanics dominated. The transition between the two modes is related to changes in geometry or constitutive properties. This implies that the linear elastic interpretation of the ultimate splitting force in term of the uniaxial tensile strength of the material is only valid for special situations, e.g. for very large cylinders. Furthermore, the numerical analysis suggests that the split cylinder test is not well suited for determining the tensile strength of early age or fibre reinforced concrete.     

CDM based modelling of damage and fracture mechanisms in concrete under tension and compression

Anisotropic damage evolution and crack propagation in the elastic–brittle materials is analysed by the concepts of continuum damage mechanics (CDM) and finite element method (FEM). The modified Murakami–Kamiya (MMK) model of elastic-damage material is used to describe damage anisotropy in concrete. The Helmholtz free energy representation is discussed. The unilateral crack opening/closure effect is incorporated in such a way that the continuity requirement during unloading holds. The incremental form of the stress–strain equations is developed. The general failure criterion is proposed by checking the positive definiteness of the Hessian matrix of the free energy function. The local approach to fracture (LAF) by FEM is applied to the pre-critical damage evolution that precedes the crack initiation, and the post-critical damage/fracture interaction. Crack is modelled as the assembly of failed finite elements in the mesh, the stiffness of which is reduced to zero when the critical points at stress–strain curves are reached. A concrete specimen with the pre-load, inclined crack is analysed in order to simulate different fracture mechanisms in tension or compression. The constitutive model is capable of predicting the kinked-type crack under tension and the wing-type crack under compression.     

FRP-混凝土三点受弯梁损伤粘结模型有限元分析

该文采用双线形损伤粘结模型研究带切口FRP-混凝土三点受弯梁(3PBB)I型加载下的界面断裂性能。通过有限元参数分析,详细讨论了界面粘结强度、界面粘结能、混凝土抗拉强度、混凝土断裂能对3PBB受力性能的影响。数值模拟表明,FRP-混凝土界面有两种破坏形式,包括FRP-混凝土界面的损伤脱粘和界面混凝土的损伤脱粘破坏,与实验所观察到的现象一致。两种破坏形式尽管在宏观上均表现为界面脱粘,但破坏机制却不同。FRP-混凝土界面的损伤粘结模型与混凝土的拉伸塑性损伤模型相结合,不但再现了3PBB的宏观力学性能,数值分析得到的荷载-位移曲线接近实验结果,而且还能详细展示FRP-混凝土界面的损伤、断裂破坏过程以及损伤在FRP-混凝土界面和界面混凝土之间的转移,能够预测构件的承载力,有助于界面优化设计,这是单纯以能量判据预测裂纹发展的经典断裂力学方法所无法做到的。     

Computational applications of a coupled plasticity-damage constitutive model for simulating plain concrete fracture

A coupled plasticity-damage model for plain concrete is presented in this paper. Based on continuum damage mechanics (CDM), an isotropic and anisotropic damage model coupled with a plasticity model is proposed in order to effectively predict and simulate plain concrete fracture. Two different damage evolution laws for both tension and compression are formulated for a more accurate prediction of the plain concrete behavior. In order to derive the constitutive equations and for the easiness in the numerical implementation, in the CDM framework the strain equivalence hypothesis is adopted such that the strain in the effective (undamaged) configuration is equivalent to the strain in the nominal (damaged) configuration. The proposed constitutive model has been shown to satisfy the thermodynamics requirements. Detailed numerical algorithms are developed for the finite element implementation of the proposed coupled plasticity-damage model. The numerical algorithm is coded using the user subroutine UMAT and then implemented in the commercial finite element analysis program Abaqus. Special emphasis is placed on identifying the plasticity and damage model material parameters from loading-unloading uniaxial test results. The overall performance of the proposed model is verified by comparing the model predictions to various experimental data, such as monotonic uniaxial tension and compression tests, monotonic biaxial compression test, loading-unloading uniaxial tensile and compressive tests, and mixed-mode fracture tests.     

纤维增强复合材料界面脱层和基体裂纹的模拟分析

基于Ghosh提出的Voronoi单元有限元方法,构造能同时反映纤维增强复合材料界面脱层和基体裂纹扩展的单元(X-VCFEM单元);应用界面力学理论和断裂力学理论,建立界面脱层、界面裂纹扩展方向和基体裂纹扩展的判断准则;结合网格重划分技术,模拟分析了只含有一个夹杂时界面脱层和基体裂纹扩展的过程,并通过与传统有限元计算结果的比较,验证X-VCFEM单元的可靠性和有效性;同时,模拟分析含任意随机分布夹杂的纤维增强复合材料界面脱层和基体裂纹的产生和扩展过程。结果表明：应用该方法模拟复杂多相复合材料裂纹问题具有计算速度快和精度高的优越性。     

Combined continuum damage‐embedded discontinuity model for explicit dynamic fracture analyses of quasi‐brittle materials

In this paper, a novel constitutive model combining continuum damage with embedded discontinuity is developed for explicit dynamic analyses of quasi‐brittle failure phenomena. The model is capable of describing the rate‐dependent behavior in dynamics and the three phases in failure of quasi‐brittle materials. The first phase is always linear elastic, followed by the second phase corresponding to fracture‐process zone creation, represented with rate‐dependent continuum damage with isotropic hardening formulated by utilizing consistency approach. The third and final phase, involving nonlinear softening, is formulated by using an embedded displacement discontinuity model with constant displacement jumps both in normal and tangential directions. The proposed model is capable of describing the rate‐dependent ductile to brittle transition typical of cohesive materials (e.g., rocks and ice). The model is implemented in the finite element setting by using the CST elements. The displacement jump vector is solved for implicitly at the local (finite element) level along with a viscoplastic return mapping algorithm, whereas the global equations of motion are solved with explicit time‐stepping scheme. The model performance is illustrated by several numerical simulations, including both material point and structural tests. The final validation example concerns the dynamic Brazilian disc test on rock material under plane stress assumption. Copyright © 2014 John Wiley & Sons, Ltd.     

Irreversible constitutive law for modeling the delamination process using interfacial surface discontinuities

An irreversible cohesive–decohesive constitutive law is postulated for interfacial surface discontinuities to predict initiation and progression of delamination. An exponential function is used for the constitutive law that naturally satisfies a multi-axial stress criterion for the onset of delamination and a mixed mode fracture criterion for the progression of delamination. The constitutive equations are made thermomechanically consistent by including a damage parameter to prevent the restoration of the previous cohesive state between the interfacial surfaces. To demonstrate the capability to predict delamination and the irreversibility capability of the constitutive law, steady-state delamination growth is simulated for quasi-static loading–unloading cycle of various fracture test specimens. The finite element results are in good agreement with either experimental data available in the literature or with linear elastic fracture mechanics analytical solutions which are valid for sufficiently long cracks.