Numerical analysis of projectile penetration and perforation of plain and fiber reinforced concrete slabs

The research presented in this paper deals with the numerical analysis of projectile impact on regular strength concrete (RSC), high‐strength concrete (HSC), and engineered cementitious composites (ECC) using the Lattice Discrete Particle Model (LDPM). The LDPM is chosen in this study as it naturally captures the failure mechanisms at the length scale of coarse aggregate of concrete, and its capabilities include the accurate depiction of both intrinsic and apparent rate effects in concrete, as well as fiber reinforcement effects. The model is used to predict the experimental impact response performed by four independent testing laboratories, and for each data set the model parameters are calibrated and validated using a combination of uniaxial compression, triaxial compression, uniaxial strain compression, and dogbone tests. In the first study, perforation experiments on RSC and HSC for varied impact velocities are carried out, and the exit velocity is compared with the available experimental data. The second study focuses on ECC, where multiple impact of steel and plastic fiber reinforced concrete panels are explored. A third investigation is performed on four RSC panels with varied thicknesses and subjected to the same impact velocity. In this instance, the model is used to predict the penetration depths for the different cases. Finally, in the last study, the response of large‐thickness infinite panels of sizes ranging from 300 mm to 700 mm under projectile impact is considered. Copyright © 2016 John Wiley & Sons, Ltd.     

Multiscale modeling of a sensitive marine clay

This paper examines the mechanical behavior of a sensitive marine clay. Various laboratory tests on intact and reconstituted samples of Guinea Gulf marine clay were performed under isotropic compression and drained triaxial compression at constant confining stresses. Microstructure analysis on intact and reconstituted samples was also carried out under different loading conditions. The effect of inter‐aggregates bonding on mechanical properties is discussed. Based on experimental analysis, a new modeling method is proposed. In this approach, the clay is regarded as an assembly of aggregates of clay particles. An inter‐aggregate contact law is introduced relating contact forces to aggregates relative displacements. The deformation of the assembly can be obtained by integrating the movement of the inter‐aggregate contacts in all orientations. Thus, the effect of inter‐aggregates bonds and debonding is considered in a direct way. The model is evaluated through comparisons between the predicted and measured results on Guinea Gulf marine clay. The evolutions of local stresses, strains, and bonds in inter‐aggregates planes are discussed to explain the anisotropy induced by the applied loading. Copyright © 2010 John Wiley & Sons, Ltd.     

Microplane damage model for jointed rock masses

The paper presents a new microplane constitutive model for the inelastic behavior of jointed rock masses that takes into account the mechanical behavior and geometric characteristics of cracks and joints. The basic idea is that the microplane modeling of rock masses under general triaxial loading, including compression, requires the isotropic rock matrix and the joints to be considered as two distinct phases coupled in parallel. A joint continuity factor is defined as a microplane damage variable to represent the stress‐carrying area fraction of the joint phase. Based on the assumption of parallel coupling between the rock joint and the rock matrix, the overall mechanical behavior of the rock is characterized by microplane constitutive laws for the rock matrix and for the rock joints, along with an evolution law for the microplane joint continuity factor. The inelastic response of the rock matrix and the rock joints is controlled on the microplane level by the stress–strain boundaries. Based on the arguments enunciated in developing the new microplane model M7 for concrete, the previously used volumetric–deviatoric splits of the elastic strains and of the tensile boundary are avoided. The boundaries are tensile normal, compressive normal, and shear. The numerical simulations demonstrate satisfactory fits of published triaxial test data on sandstone and on jointed plaster mortar, including quintessential features such as the strain softening and dilatancy under low confining pressure, as well as the brittle–ductile transition under higher confining pressure, and the decrease of jointed rock strength and Young's modulus with an increasing dip angle of the joint. Copyright © 2014 John Wiley & Sons, Ltd.     

Meso‐scale particle modeling of concrete deterioration caused by alkali‐aggregate reaction

A meso‐scale particle model is presented to simulate the expansion of concrete subjected to alkali‐aggregate reaction (AAR) and to analyze the AAR‐induced degradation of the mechanical properties. It is the first attempt to evaluate the deterioration mechanism due to AAR using the discrete‐element method. A three‐phase meso‐scale model for concrete composed of aggregates, mortar and the interface is established with the combination of a pre‐processing approach and the particle flow code, PFC2D. A homogeneous aggregate expansion approach is applied to model the AAR expansion. Uniaxial compression tests are conducted for the AAR‐affected concrete to examine the effects on the mechanical properties. Two specimens with different aggregate sizes are analyzed to consider the effects of aggregate size on AAR. The results show that the meso‐scale particle model is valid to predict the expansion and the internal micro‐cracking patterns caused by AAR. The two different specimens exhibit similar behavior. The Young's modulus and compressive strength are significantly reduced with the increase of AAR expansion. The shape of the stress–strain curves obtained from the compression tests clearly reflects the influence of internal micro‐cracks: an increased nonlinearity before the peak loading and a more gradual softening for more severely affected specimens. Similar macroscopic failure patterns of the specimens under compression are observed in terms of diagonal macroscopic cracks splitting the specimen into several triangular pieces, whereas localized micro‐cracks forming in slightly affected specimens are different from branching and diffusing cracks in severely affected ones, demonstrating different failure mechanisms. Copyright © 2012 John Wiley & Sons, Ltd.     

混凝土二维骨料随机分布对物质扩散的影响

将混凝土视为骨料和砂浆组成的二相复合材料。首先建立二维骨料的随机分布模型,并改进模型骨料重叠的判断方法,在此基础上以混凝土的中性化为例,通过计算机模拟骨料随机分布对混凝土中性化时CO2扩散的延缓作用。结果表明,单个骨料的形状、大小和摆放位置都会对CO2扩散路径产生影响,而骨料集合体通过扩散路径的长度、迂曲度的变化影响扩散速度,表现为混凝土的扩散系数Da与骨料圆度?、平均粒径d和骨料含量ra具有良好的对应关系。最后,通过数值模拟确定了二维扩散中Da-?、Da-d和Da-ra的关系式,可供工程应用参考。     

Numerical modeling of the elastoplastic damage behavior of dry and saturated concrete targets subjected to rigid projectile penetration

The objective of the present paper is to present a numerical study on the penetration performance of concrete targets with 2 different water contents. Numerical analysis has been performed by using the finite element code Abaqus/Explicit, in which a coupled elastoplastic damage model has been developed for saturated/unsaturated concrete under a wide range of confining pressures. The performance of proposed model has been firstly verified by simulating the triaxial compression tests and penetration tests realized with saturated/dry concretes. Comparisons of available experimental results and numerical simulations show that the proposed model is able to reproduce satisfactorily the mechanical behavior of saturated and dry concretes. A higher failure stress and a more important pores closing are generally obtained in dry concrete samples with respect to saturated ones. Furthermore, the main observed patterns of penetration test realized with saturated concrete targets are also satisfactorily simulated by the numerical results. Therefore, the proposed model is used to numerically predict the penetration performance of dry concrete target, and the penetration performance of dry/saturated concrete target is discussed. We observe that in dry concrete target, the penetration of projectile is strongly declined, and a smaller damage zone is created. The numerical predictions and discussions can help engineers to enhance their understandings on the influence of hydraulic conditions on structural vulnerability of concrete structures subjected to near‐field detonations or impacts.     

A macroscopic damage‐plastic constitutive law for modeling quasi‐brittle fracture and ductile behavior of concrete

A new phenomenological macroscopic constitutive model for the numerical simulation of quasi‐brittle fracture and ductile concrete behavior, under general triaxial stress conditions, is presented. The model is particularly addressed to simulate a wide range of confinement stress states, as also, to capture the strong influence of the mean stress value in the concrete failure mechanisms. The model is based on a two‐surface damage‐plastic formulation. The mechanical behavior in different domains of the stress space is separately described by means of a quasi‐brittle or ductile material response:

相似文献   

Experimental study and constitutive modelling of elasto‐plastic damage in heat‐treated mortar

This study investigates the effect of a heat‐treatment upon the thermo‐mechanical behaviour of a model cement‐based material, i.e. a normalized mortar, with a (w/c) ratio of 0.5. First, a whole set of varied experimental results is provided, in order to either identify or validate a thermo‐mechanical constitutive model, presented in the second paper part. Experimental responses of both hydraulic and mechanical behaviour are given after different heating/cooling cycling levels (105, 200, 300, 400^°C). The reference state, used for comparison purposes, is taken after mass stabilization at 60^°C. Typical uniaxial compression tests are provided, and original triaxial deviatoric compressive test responses are also given. Hydraulic behaviour is identified simultaneously to triaxial deviatoric compressive loading through gas permeability K_gas assessment. K_gas is well correlated with volumetric strain evolution: gas permeability increases hugely when ε_v testifies of a dilatant material behaviour, instead of contractile from the test start. Finally, the thermo‐mechanical model, based on a thermodynamics approach, is identified using the experimental results on uniaxial and triaxial deviatoric compression. It is also positively validated at residual state for triaxial deviatoric compression, but also by using a different stress path in lateral extension, which is at the origin of noticeable plasticity. Copyright © 2009 John Wiley & Sons, Ltd.     

基于细观力学的混凝土劈拉破坏率相关特性研究

提出了一种高效的混凝土细观力学预处理方法,其中包括基于背景网格的骨料投放与细观三组分的识别;建立了颗粒离散元的混凝土细观力学模型。利用该模型对不同应变率下混凝土的劈拉试验进行了数值仿真,仿真结果显示出的混凝土动力强度提高的规律以及不同应变率下混凝土的破坏模式与试验结果吻合较好。研究表明,随着应变率的提高,混凝土破坏时产生更多弥散状分布的裂纹,同时有些裂纹穿过强度较高的骨料,导致混凝土在高应变率破坏过程中消耗更多的能量,是混凝土的动力强度提高的重要原因。此外,在高应变率下,力链产生更多的分叉,空间分布的不均匀性和随机性增加,因此,在高应变率下混凝土动力强度表现出更大的离散性。     

Numerical evaluation of soft inter‐slab joint in concrete‐faced rockfill dam with dual mortar finite element method

Recently constructed concrete‐faced rockfill dams (CFRDs) often use soft inter‐slab joints to prevent axial compression‐induced extrusion damage in the concrete face. Due to the complexity of the multibody contact and the lack of information on the actual behavior of soft joints, it is highly challenging to numerically assess the effect of soft joints in CFRDs. In this paper, we present a numerical approach for the three‐dimensional modeling of CFRDs with hard and soft joints. A dual mortar finite element method with Lagrange multiplier is developed to treat the multibody contact in hard joints with impenetrability condition. The soft joint slab‐filler‐slab contact system is modeled using an equivalent contact interface approach, where the soft contact constraints are imposed using a perturbed Lagrange formulation. Through a series of laboratory tests, the mechanical behavior of soft joint is investigated. An extrusion model for the soft joint is presented and implemented in the dual mortar finite element method. The proposed numerical method is applied to the three‐dimensional analysis of Tianshengqiao‐1 CFRD. Despite the complex multibody contact and strong material and geometry nonlinearities in the CFRD, the proposed method is stable and capable of capturing salient characteristics of the CFRD. Numerical results show that in Tianshengqiao‐1, the employment of soft joints can effectively reduce the axial compression stress, thus greatly alleviating the risk of extrusion damage in the concrete face.     

Three‐dimensional elasto‐plastic model for unsaturated compacted soils with different initial densities

This paper presents an elasto‐plastic model for unsaturated compacted soils and experimental results obtained from a series of suction‐controlled triaxial tests on unsaturated compacted clay with different initial densities. The initial density dependency of the compacted soil behaviour is modelled by establishing experimental relationships between the initial density and the corresponding yield stress and thereby between the initial density and the location and slope of normal compression line. The model is generalized to three‐dimensional stress states by assuming that the shapes of the failure surface and the yield surface in the deviatoric plane are given by the extended SMP criterion. A considerable number of the isotropic compression, triaxial compression and extension tests on unsaturated compacted clay with different initial densities were performed using a suction‐controllable triaxial apparatus, to measure the stress–strain–volume change in different stress paths and wetting paths. The model has well‐predicting capabilities to reproduce the mechanical behaviour of specimens compacted under different conditions not only in isotropic compression but also in triaxial compression and triaxial extension. Copyright © 2003 John Wiley & Sons, Ltd.     

Optimization framework for calibration of constitutive models enhanced by neural networks

A two‐level procedure designed for the estimation of constitutive model parameters is presented in this paper. The neural network (NN) approach at the first level is applied to achieve the first approximation of parameters. This technique is used to avoid potential pitfalls related to the conventional gradient‐based optimization techniques, considered here as a corrector that improves predicted parameters. The feed‐forward NN (FFNN) and the modified Gauss–Newton algorithms are briefly presented. The proposed framework is verified for the elasto‐plastic modified Cam Clay model that can be calibrated based on standard triaxial laboratory tests, i.e. the isotropic consolidation test and the drained compression test. Two different formulations of the input data to the NN, enhanced by a dimensional reduction of experimental data using principal component analysis, are presented. The determination of model characteristics is demonstrated, first on numerical pseudo‐experiments and then on the experimental data. The efficiency of the proposed approach by means of accuracy and computational effort is also discussed. Copyright © 2008 John Wiley & Sons, Ltd.     

岩盐弹塑性损伤耦合模型研究

岩盐力学模型是进行能源岩盐储存工程稳定性分析的基础,而损伤和塑性机制并存且相互耦合是岩盐力学行为的基本特点。采用云应岩盐,进行了多组围压条件下的三轴压缩试验,分析了不同围压下岩盐的变形特征。在试验分析的基础上,提出了一种能够描述岩盐特性的弹塑性损伤耦合的模型,该模型描述了岩盐损伤的演化和塑性变形的耦合关系,并引入了一种非关联的塑性流动法则来描述岩盐从塑性体积压缩到膨胀的转化。采用该模型对在三轴压缩下的岩盐应力-应变关系进行了模拟分析,并与试验数据进行了对比,结果表明该模型能够较好地描述岩盐的主要力学和变形特性。     

Experimental and theoretical investigation of the high‐pressure,undrained response of a cohesionless sand

Key in predicting stability in sands during dynamic events is gaining a fundamental understanding of the physics of its deformation and failure at high pressures. In this paper are reported results of an experimental investigation into the high‐pressure (up to 700 MPa) mechanical response of Quikrete sand. During all triaxial compression tests, the material exhibited hardening up to failure while both compressibility and dilatancy regimes of the volumetric response were observed. Furthermore, the transition from compressibility to dilatancy was found to be highly dependent on confining pressure. By performing triaxial compression tests with several creep stages followed by unloading–reloading cycles, the time influence of the overall response was detected. Using the experimental data, a new model that captures both compressibility and dilatancy has been developed. Comparison between model predictions and data showed that the proposed model describes well the main characteristics of the high‐pressure response of sand. Copyright © 2012 John Wiley & Sons, Ltd.     

Effect of the water/cement ratio on concrete behavior under extreme loading

This study focuses on identifying concrete behavior under severe triaxial loadings (near field detonation or ballistic impacts). In order to reproduce high stress levels with well‐controlled loading paths, static tests have been carried out on concrete samples by mean of a very high‐capacity triaxial press (stress levels on the order of 1 GPa). It is a longstanding fact that the water/cement ratio (W/C), upon entering the concrete composition, is a major parameter affecting the porosity and strength of the cement matrix of hardened concrete. The objective of this article is to quantify the effect of this ratio on concrete behavior under conditions of high confinement. From the composition of a reference ‘ordinary’ concrete (i.e. W/C=0.6), two other concretes have been produced with W/C ratios equal to 0.4 and 0.8, respectively. This article presents experimental results and their analysis regarding the effect of water/cement ratio (W/C) on concrete behavior under high confinement. It shows that when placed under high confinement, concrete behaves like a granular stacking composed of concrete without any influence from the level of cement matrix strength. Copyright © 2009 John Wiley & Sons, Ltd.     

Influence of size on the constitutive equations of concrete or rock dowels

The numerical fracture analysis of non‐homogeneous rock or concrete dowels subjected to shear and compression is described in detail. The method of analysis allows the consideration of scale and rate effects due to material non‐homogeneity and fracture. The proposed approach is verified by comparing numerical predictions with experimental results reported in the literature for a series of small rock samples, since experimental evidence for large bodies is not yet available (2007). Results generated by Monte Carlo simulation using the so‐called discrete element method to model the dowels suggest that a simple three parameters law can be used to predict the relationship between tangential stress at the base and lateral distortion. It is observed that the larger the size of the cubes, the smaller both the peak tangential stress and the rupture distortion. Size effects are also evaluated in samples with vertical restraint. The influence of loading rate is likewise numerically assessed for two sample sizes. The effect is compatible with experimental evidence available for concrete using small samples. Copyright © 2008 John Wiley & Sons, Ltd.     

Experimental and theoretical investigation of the high‐pressure behavior of concrete

The results of an experimental study aimed at characterizing the behavior of concrete for high confining pressures (up to 500 MPa) are reported. The main characteristics of the response under deviatoric conditions are quasi‐linearity in the elastic regime, stress‐path dependency, and gradual change from compressibility to dilatancy under increasing deviatoric stress. By performing the cyclic triaxial compression tests with several load–creep–unload and reload cycles, the time influence on the overall behavior was detected. Further, whether the main features of the observed behavior can be described within the framework of elastic/viscoplasticity theory was investigated. For this purpose, Cristescu's (Rock Rheology. Kluwer Academic Publishers: The Netherlands, 1989) approach was used. It was shown that the proposed elastic/viscoplastic model captures the main features of concrete behavior at high pressures. Copyright © 2008 John Wiley & Sons, Ltd.     

隧道衬砌结构三轴压缩蠕变特性试验研究

使用MTS-815液压伺服系统对隧道衬砌的混凝土试件进行了三轴压缩蠕变试验,当压缩蠕变强度达到混凝土试件抗压强度的85%时进入加速蠕变阶段。分析了隧道衬砌混凝土三轴压缩蠕变变形的特点,确定采用摩尔-库仑准则和Burgers模型串联组合而成的黏弹塑性模型可以合理地描述衬砌混凝土的蠕变特性。组合模型的黏弹性状态与Burgers模型一致,塑性状态与Mohr-Coulomb模型一致,它不仅能够反映衬砌结构蠕变的全过程,且能模拟材料的黏弹塑性变形特性,解决了典型Burgers蠕变模型只能描述黏弹性蠕变的缺陷。同时也研究了三轴试验条件下组合蠕变模型参数的求解方法,得到隧道衬砌混凝土的蠕变参数。     

Predictions of large stress reversals in true triaxial tests on cross‐anisotropic sand

A rotational kinematic hardening constitutive model with the capability of predicting the behavior of soil during three‐dimensional stress reversals has been developed. An existing elasto‐plastic constitutive model, the Single Hardening Model, utilizing isotropic hardening serves as the basic framework in these formulations. The framework of the kinematic hardening model was discussed in a companion paper. The previously proposed cross‐anisotropic Single Hardening Model is added to the present kinematic hardening mechanism to capture inherent anisotropy of sands in addition to the stress reversals. This model involves 13 parameters, which can be determined from simple laboratory experiments, such as isotropic compression, drained triaxial compression and triaxial extension tests. The results from a series of true triaxial tests with large three‐dimensional stress reversals performed on medium dense cross‐anisotropic Santa Monica Beach sand are employed for comparison with predictions. Copyright © 2008 John Wiley & Sons, Ltd.     

Microplane constitutive model for porous isotropic rocks

The paper deals with constitutive modelling of contiguous rock located between rock joints. A fully explicit kinematically constrained microplane‐type constitutive model for hardening and softening non‐linear triaxial behaviour of isotropic porous rock is developed. The microplane framework, in which the constitutive relation is expressed in terms of stress and strain vectors rather than tensors, makes it possible to model various microstructural physical mechanisms associated with oriented internal surfaces, such as cracking, slip, friction and splitting of a particular orientation. Formulation of the constitutive relation is facilitated by the fact that it is decoupled from the tensorial invariance restrictions, which are satisfied automatically. In its basic features, the present model is similar to the recently developed microplane model M4 for concrete, but there are significant improvements and modifications. They include a realistic simulation of (1) the effects of pore collapse on the volume changes during triaxial loading and on the reduction of frictional strength, (2) recovery of frictional strength during shearing, and (3) the shear‐enhanced compaction in triaxial tests, manifested by a deviation from the hydrostatic stress–strain curve. The model is calibrated by optimal fitting of extensive triaxial test data for Salem limestone, and good fits are demonstrated. Although these data do not cover the entire range of behaviour, credence in broad capabilities of the model is lend by its similarity to model M4 for concrete—an artificial rock. The model is intended for large explicit finite‐element programs. Copyright © 2002 John Wiley & Sons, Ltd.