Investigation of FRP and SFRC Technologies for Efficient Tunnel Reinforcement Using the Cohesive Zone Model

Amid urbanization and the continuous expansion of transportation networks,the necessity for tunnel construction and maintenance has become paramount.Addressing this need requires the investigation of efficient,economical,and robust tunnel reinforcement techniques.This paper explores fiber reinforced polymer(FRP)and steel fiber reinforced concrete(SFRC)technologies,which have emerged as viable solutions for enhancing tunnel structures.FRP is celebrated for its lightweight and high-strength attributes,effectively augmenting load-bearing capacity and seismic resistance,while SFRC’s notable crack resistance and longevity potentially enhance the performance of tunnel segments.Nonetheless,current research predominantly focuses on experimental analysis,lacking comprehensive theoretical models.To bridge this gap,the cohesive zone model(CZM),which utilizes cohesive elements to characterize the potential fracture surfaces of concrete/SFRC,the rebar-concrete interface,and the FRP-concrete interface,was employed.A modeling approach was subsequently proposed to construct a tunnel segment model reinforced with either SFRC or FRP.Moreover,the corresponding mixed-mode constitutive models,considering interfacial friction,were integrated into the proposed model.Experimental validation and numerical simulations corroborated the accuracy of the proposed model.Additionally,this study examined the reinforcement design of tunnel segments.Through a numerical evaluation,the effectiveness of innovative reinforcement schemes,such as substituting concrete with SFRC and externally bonding FRP sheets,was assessed utilizing a case study from the Fuzhou Metro Shield Tunnel Construction Project.

Extended finite element-based cohesive zone method for modeling simultaneous hydraulic fracture height growth in layered reservoirs

In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.

Numerical simulation of the mechanical behavior of superconducting tape in conductor on round core cable using the cohesive zone model

Cables composed of rare-earth barium copper oxide(REBCO)tapes have been extensively used in various superconducting devices.In recent years,conductor on round core(CORC)cable has drawn the attention of researchers with its outstanding current-carrying capacity and mechanical properties.The REBCO tapes are wound spirally on the surface of CORC cable.Under extreme loadings,the REBCO tapes with layered composite structures are vulnerable,which can lead to degradation of critical current and even quenching of superconducting devices.In this paper,we simulate the deformation of CORC cable under external loads,and analyze the damage inside the tape with the cohesive zone model(CZM).Firstly,the fabrication and cabling of CORC are simulated,and the stresses and strains generated in the tape are extracted as the initial condition of the next step.Then,the tension and bending loads are applied to CORC cable,and the damage distribution inside the tape is presented.In addition,the effects of some parameters on the damage are discussed during the bending simulations.

Numerical study of fatigue damage of asphalt concrete using cohesive zone model

In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model （CZM） is developed. Integrated with the CZM, a fatigue damage evolution model is established to indicate the gradual degradation of cohesive properties of asphalt concrete under cyclic loading. Then the model is implemented in the finite element software ABAQUS through a user-defined subroutine. Based on the proposed model, an indirect tensile fatigue test is finally simulated. The fatigue lives obtained through numerical analysis show good agreement with laboratory results. Fatigue damage accumulates in a nonlinear manner during the cyclic loading process and damage initiation phase is the major part of fatigue failure. As the stress ratio increases, the time of the steady damage growth stage decreases significantly. It is found that the proposed fatigue damage evolution model can serve as an accurate and efficient tool for the prediction of fatigue damage of asphalt concrete.

CRACK PROPAGATION IN POLYCRYSTALLINE ELASTIC-VISCOPLASTIC MATERIALS USING COHESIVE ZONE MODELS

Cohesive zone model was used to simulate two-dimensional plane strain crack propagation at the grain level model including grain boundary zones. Simulated results show that the original crack-tip may not be separated firstly in an elastic-viscoplastic polycrystals. The grain interior＇s material properties （e.g. strain rate sensitivity） characterize the competitions between plastic and cohesive energy dissipation mechanisms. The higher the strain rate sensitivity is, the larger amount of the external work is transformed into plastic dissipation energy than into cohesive energy, which delays the cohesive zone rupturing. With the strain rate sensitivity decreased, the material property tends to approach the elastic-plastic responses. In this case, the plastic dissipation energy decreases and the cohesive dissipation energy increases which accelerates the cohesive zones debonding. Increasing the cohesive strength or the critical separation displacement will reduce the stress triaxiality at grain interiors and grain boundaries. Enhancing the cohesive zones ductility can improve the matrix materials resistance to void damage.

Finite Element Simulations of the Localized Failure and Fracture Propagation in Cohesive Materials with Friction

Strain localization frequently occurs in cohesive materials with friction(e.g.,composites,soils,rocks)and is widely recognized as a fundamental cause of progressive structural failure.Nonetheless,achieving high-fidelity simulation for this issue,particularly concerning strong discontinuities and tension-compression-shear behaviors within localized zones,remains significantly constrained.In response,this study introduces an integrated algorithmwithin the finite element framework,merging a coupled cohesive zone model(CZM)with the nonlinear augmented finite elementmethod(N-AFEM).The coupledCZMcomprehensively describes tension-compression and compressionshear failure behaviors in cohesive,frictional materials,while the N-AFEM allows nonlinear coupled intraelement discontinuities without necessitating extra nodes or nodal DoFs.Following CZM validation using existing experimental data,this integrated algorithm was utilized to analyze soil slope failure mechanisms involving a specific tensile strength and to assess the impact of mechanical parameters(e.g.,tensile strength,weighting factor,modulus)in soils.

基于Voronoi建模和Cohesive单元的PBX细观力学模型

针对PBX的细观力学模型存在模型结构与真实PBX细观形貌不符、破坏行为描述不全面(仅考虑界面脱粘)、研究对象单一(大多针对PBX-9501)等问题,建立了基于Voronoi建模和Cohesive单元的PBX细观力学模型。首先,基于光学显微镜下的细观结构图,采用Voronoi方法构建了更贴近实际形貌的PBX炸药细观模型;然后,采用内聚力有限元方法在炸药晶体内部、黏结剂内部及炸药颗粒/黏结剂界面处都引入了Cohesive单元,以在二维尺度上实现任意路径破坏行为的模拟;最后,针对某HMX基-F2311黏结剂炸药和某HMX基-F2313黏结剂炸药标定了模型参数,模拟了这两种炸药的准静态拉伸行为。结果表明,两种炸药数值模拟的拉伸变形曲线、破坏应力、破坏应变都与试验数据吻合良好,表明该模型不仅适用于不同PBX炸药材料(线弹性PBX、非线性PBX)的变形行为描述,还可以分析PBX炸药在拉伸破坏过程中的细观机理,刻画不同PBX炸药的裂纹萌生、扩展和贯穿的过程。

Multi-scenario Simulation and Spatial-temporal Analysis of LUCC in China's Coastal Zone Based on Coupled SD-FLUS Model

Increased human activities in China's coastal zone have resulted in the depletion of ecological land resources.Thus,conducting current and future multi-scenario simulation research on land use and land cover change(LUCC)is crucial for guiding the healthy and sustainable development of coastal zones.System dynamic(SD)-future land use simulation(FLUS)model,a coupled simulation model,was developed to analyze land use dynamics in China's coastal zone.This model encompasses five scenarios,namely,SSP1-RCP2.6(A),SSP2-RCP4.5(B),SSP3-RCP4.5(C),SSP4-RCP4.5(D),and SSP5-RCP8.5(E).The SD model simulates land use demand on an annual basis up to the year 2100.Subsequently,the FLUS model determines the spatial distribution of land use for the near term(2035),medium term(2050),and long term(2100).Results reveal a slowing trend in land use changes in China's coastal zone from 2000–2020.Among these changes,the expansion rate of construction land was the highest and exhibited an annual decrease.By 2100,land use predictions exhibit high accuracy,and notable differences are observed in trends across scenarios.In summary,the expansion of production,living,and ecological spaces toward the sea remains prominent.Scenario A emphasizes reduced land resource dependence,benefiting ecological land protection.Scenario B witnesses an intensified expansion of artificial wetlands.Scenario C sees substantial land needs for living and production,while Scenario D shows coastal forest and grassland shrinkage.Lastly,in Scenario E,the conflict between humans and land intensifies.This study presents pertinent recommendations for the future development,utilization,and management of coastal areas in China.The research contributes valuable scientific support for informed,long-term strategic decision making within coastal regions.

Cameroon Climate Predictions Using the SARIMA-LSTM Machine Learning Model: Adjustment of a Climate Model for the Sudano-Sahelian Zone of Cameroon

It is acknowledged today within the scientific community that two types of actions must be considered to limit global warming: mitigation actions by reducing GHG emissions, to contain the rate of global warming, and adaptation actions to adapt societies to Climate Change, to limit losses and damages [1] [2]. As far as adaptation actions are concerned, numerical simulation, due to its results, its costs which require less investment than tests carried out on complex mechanical structures, and its implementation facilities, appears to be a major step in the design and prediction of complex mechanical systems. However, despite the quality of the results obtained, biases and inaccuracies related to the structure of the models do exist. Therefore, there is a need to validate the results of this SARIMA-LSTM-digital learning model adjusted by a matching approach, “calculating-test”, in order to assess the quality of the results and the performance of the model. The methodology consists of exploiting two climatic databases (temperature and precipitation), one of which is in-situ and the other spatial, all derived from grid points. Data from the dot grids are processed and stored in specific formats and, through machine learning approaches, complex mathematical equations are worked out and interconnections within the climate system established. Through this mathematical approach, it is possible to predict the future climate of the Sudano-Sahelian zone of Cameroon and to propose adaptation strategies.

COHESIVE ZONE FINITE ELEMENT-BASED MODELING OF HYDRAULIC FRACTURES

Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.

Failure simulation in resistance spot-welded lap-joints using cohesive zone modeling

This paper concentrates on simulating fracture in thin walled single-lap joints connected by resistance spot-welding(RSW)process which were subjected to tensile loading.For this purpose,three sets of lap-joints with different spot configurations were tested to achieve the joints’tensile behavior.To simulate the joints tensile behavior,firstly a 2D axisymmetric finite element(FE)model was used to calculate residual stresses induced during the welding process.Then the results were transferred to 3D models as pre-stress.In this step,cohesive zone model(CZM)technique was used to simulate fracture in the models under tensile load.Cohesive zone parameters were extracted using coach-peel and shear lap specimens.The results were employed to simulate deformation and failure in single lap spot weld samples.It has been shown that considering the residual stresses in simulating deformation and fracture load enables quite accurate predictions.

Delamination analysis of woven fabrication laminates using cohesive zone model

A new test method was proposed to evaluate the cohesive strength of composite laminates. Cohesive strength and the critical strain energy for Mode-II interlamiar fracture of E-glass/epoxy woven fabrication were determined from the single lap joint(SLJ) and end notch flexure(ENF) test, respectively. In order to verify their adequacy, a cohesive zone model simulation based on interface finite elements was performed. A closed form solution for determination of the penalty stiffness parameter was proposed. Modified form of Park-Paulino-Roesler traction-separation law was provided and conducted altogether with trapezoidal and bilinear mixed-mode damage models to simulate damage using Abaqus cohesive elements. It was observed that accurate damage prediction and numerical convergence were obtained using the proposed penalty stiffness. Comparison between three damage models reveals that good simulation of fracture process zone and delamination prediction were obtained using the modified PPR model as damage model. Cohesive zone length as a material property was determined. To ensure the sufficient dissipation of energy, it was recommended that at least 4 elements should span cohesive zone length.

A Modified Cohesive Zone Model for Simulation of Delamination Behavior in Laminated Composites

Considering the promotion effect of interlaminar normal tensile stress and the inhibition effect of interlaminar normal compressive stress,two kinds of elimination initial criteria were proposed in this paper.Based on these two delamination initial criteria,a modified cohesive zone model(CZM)was established to simulate the delamination behavior in laminated composites.Numerical simulations of double cantilever beam(DCB),mixed-mode bending(MMB)and end notched flexure(ENF)tests were conducted.The results show that the proposed model can do a better job than common ones when it is used to predict laminates’delamination under interlaminar compression stress.Moreover,a factor r,named cohesive strength coefficient,was defined in this paper on account of the difference between cohesive strength and interlaminar fracture strength.With changing factor r,it shows that a moderate variation of cohesive strength will not cause significant influences on global load-displacement responses.Besides,in order to obtain a good balance between prediction accuracy and computational efficiency,there shall be two or three numerical elements within the cohesive zone.

Determination of key parameters of Al–Li alloy adhesively bonded joints using cohesive zone model

The key parameters of the adhesive layer of a reinforcing patch are of great significance and affect the ability to suppress crack propagation in an Al–Li alloy patch-reinforced structure.This paper proposes a method to determine the key parameters of the adhesive layer of adhesively bonded joints in the Al–Li alloy patch-reinforced structure.A zero-thickness cohesive zone model(CZM)was selected to simulate the adhesive layer’s fracture process,and an orthogonal simulation was designed to compare against the test results.A three-dimensional progressive damage model of an Al–Li alloy patch-reinforced structure with single-lap adhesively bonded joints was developed.The simulation’s results closely agree with the test results,demonstrating that this method of determining the key parameters is likely accurate.The results also verify the correctness of the cohesive strength and fracture energy,the two key parameters of the cohesive zone model.The model can accurately predict the strength and fracture process of adhesively bonded joints,and can be used in research to suppress crack propagation in Al–Li alloy patch-reinforced structures.

Determination of Key Cohesive Zone Model’s Parameters for Orthotropic Paper and Its Static Fracture Simulation

Investigation of paper cutting process is vital for the design of cutting tools,but the fracture mechanism of paper cutting is still unclear.Here,we focus on the cutting process of paper,including the key parameters of cohesive zone model(CZM)for the orthotropic paper,to simulate the shear fracture process.Firstly,the material constants of the orthotropic paper are determined by longitudinal and transverse tensile test.Secondly,based on the tensile stressstrain curves,combined with damage theory and numerical simulations,the key parameters of the CZM for the orthotropic paper are obtained.Finally,a model III fracture is simulated to verify the accuracy of the model.Results show that the load-displacement curves obtained by the simulation is consistent with the test results.

Finite element simulation of the micromachining of nanosized-silicon-carbide-particle reinforced composite materials based on the cohesive zone model

A finite element method based on the cohesive zone model was used to study the micromachining process of nanosized silicon-carbide-particle(SiCp) reinforced aluminum matrix composites. As a hierarchical multiscale simulation method, the parameters for the cohesive zone model were obtained from the stress-displacement curves of the molecular dynamics simulation. The model considers the random properties of the siliconcarbide-particle distribution and the interface of bonding between the silicon carbide particles and the matrix.The machining mechanics was analyzed according to the chip morphology, stress distribution, cutting temperature, and cutting force. The simulation results revealed that the random distribution of nanosized SiCp causes non-uniform interaction between the tool and the reinforcement particles. This deformation mechanics leads to inhomogeneous stress distribution and irregular cutting force variation.

IMPROVED COHESIVE ZONE MODEL AND ITS APPLICATION IN INTERFACE CONTACT ANALYSIS

An improved interface cohesive zone model is developed for the simulation of interface contact, under mixed-mode loading. A new debonding initiation criterion and propagation of debonding law, taking into account the pressure stress influence on contact shear strength, is proposed. The model is implemented in a finite-element program using subroutine VUINTER of ABAQUS Explicit. An edge-notch four-point bending process and laminated vibration damping steel sheet punch forming test are simulated with the improved model in ABAQUS Explicit. The numerical predictions agree satisfactorily with the corresponding experimental results.

Predicting the Thickness of an Excavation Damaged Zone around the Roadway Using the DA-RF Hybrid Model

After the excavation of the roadway,the original stress balance is destroyed,resulting in the redistribution of stress and the formation of an excavation damaged zone(EDZ)around the roadway.The thickness of EDZ is the key basis for roadway stability discrimination and support structure design,and it is of great engineering significance to accurately predict the thickness of EDZ.Considering the advantages of machine learning(ML)in dealing with high-dimensional,nonlinear problems,a hybrid prediction model based on the random forest(RF)algorithm is developed in this paper.The model used the dragonfly algorithm(DA)to optimize two hyperparameters in RF,namely mtry and ntree,and used mean absolute error(MAE),rootmean square error(RMSE),determination coefficient(R^(2)),and variance accounted for(VAF)to evaluatemodel prediction performance.A database containing 217 sets of data was collected,with embedding depth(ED),drift span(DS),surrounding rock mass strength(RMS),joint index(JI)as input variables,and the excavation damaged zone thickness(EDZT)as output variable.In addition,four classic models,back propagation neural network(BPNN),extreme learning machine(ELM),radial basis function network(RBF),and RF were compared with the DA-RF model.The results showed that the DARF mold had the best prediction performance(training set:MAE=0.1036,RMSE=0.1514,R^(2)=0.9577,VAF=94.2645;test set:MAE=0.1115,RMSE=0.1417,R^(2)=0.9423,VAF=94.0836).The results of the sensitivity analysis showed that the relative importance of each input variable was DS,ED,RMS,and JI from low to high.

A novel method for modeling the phase change of iron ore particles in the cohesive zone of a blast furnace

Cohesive zone plays a vital role in the stable operation of a blast furnace(BF),yet the complex phase change process of iron ore particles in this zone is still not well understood.In this study,a novel one-dimensional(1D)unsteady phase change model was developed to elucidate the heat transfer and melting mechanisms of iron ore particles.After model validation,the effects of several key operating parameters(e.g.,particle diameter,gas velocity,initial temperature)on the phase change behavior of iron ore particles were analyzed,and the joint effect of multiple parameters was discussed.The results show that larger-sized iron ore particles possess lower specific surface areas,which in turn reduces their convective heat absorption capacity.Consequently,the distance from the solid-liquid phase interface to the particle surface increases,thereby slowing down the movement of the phase interface and pro-longing the melting duration of the particles.Increasing the gas velocity and the initial temperature does not have a significant impact on reducing the duration of the complete melting process.Under the specified conditions,it is observed that increasing the gas velocity by 3-fold and 9-fold results in a reduction of the melting duration by 2.4%and 8.3%,respectively.Elevating the initial temperature of iron ore particles results in a decrease in the core-to-surface temperature difference,a slower heating rate,and a shorter duration to achieve melting.Among the factors affecting the melting process,the particle diameter is found to be the most significant in terms of the liquid phase precipitation,mushy zone thickness,and core-to-surface temperature difference of iron ore particles.

Numerical model for homogeneous cohesive dam breaching due to overtopping failure

Based on the large-scale model tests, a simplified dam breach model for homogeneous cohesive dam due to overtopping failure is put forward. The model considers headcut erosion as one of the key homogeneous cohesive dam breaching mechanisms and we calculate the time-averaged headcut migration rate using an energy-based empirical formula. A numerical method is adopted to determine the initial scour position at the downstream slope in terms of the water head and dam height, and the broad-crested weir equation is utilized to simulate the breach flow. The limit equilibrium method is used to analyze the stability of breach slope during the breach process. An iterative method is developed to simulate the coupling process of soil and water at each time step. The calculated results of three dam breach cases testify the reasonability of the model, and the sensitivity studies of soil erodibility show that sensitivity is dependent on each test case's soil conditions. In addition, three typical dam breach models, NWS BREACH, WinDAM B, and HR BREACH, are also chosen to compare with the proposed model. It is found that NWS BREACH may have large errors for cohesive dams, since it uses a noncohesive sediment transport model and does notconsider headcut erosion, WinDAM B and HR BREACH consider headcut erosion as the breaching mechanism and handle well homogeneous cohesive dam overtopping failure, but overall, the proposed model has the best performance.