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.

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.

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.

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.

Building fire zone model with symbolic mathematics

To apply the fire modelling for the fire engineer with symbolic mathematics,the key equations of a zone model were demonstrated. There were thirteen variables with nine constraints,so only four ordinary differential equations (ODEs) were required to solve. A typical fire modelling with two-room structure was studied. Accordingly,the source terms included in the ODEs were simplified and modelled,and the fourth Runge-Kutta method was used to solve the ordinary differential equations (ODEs) with symbolic mathematics. Then a zone model could be used with symbolic mathematics. It is proposed that symbolic mathematics is possible for use by fire engineer.

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

一个改进接口连贯的地区模型为接口接触的模拟被开发，在混合模式 loading.A 下面新 debonding 开始标准和 debonding 法律的繁殖，在接触上考虑压力应力影响砍强度， proposed.The 模型用骆驼毛的织物 QUS Explicit.An 边槽口的子程序 VUINTER 在一个有限元素的程序被实现四点的弯曲过程和表穿孔机形成测试的把压成薄片的颤动抑制钢与improve

Progressive fragmentation of granular assemblies within rockslides: Insights from discrete-continuous numerical modeling

Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive disintegration and kinematics of multi-deformable rock blocks during rockslides.The present study proposes a discrete-continuous numerical model,based on a cohesive zone model,to explicitly incorporate the progressive fragmentation and intricate interparticle interactions inherent in rockslides.Breakable rock granular assemblies are released along an inclined plane and flow onto a horizontal plane.The numerical scenarios are established to incorporate variations in slope angle,initial height,friction coefficient,and particle number.The evolutions of fragmentation,kinematic,runout and depositional characteristics are quantitatively analyzed and compared with experimental and field data.A positive linear relationship between the equivalent friction coefficient and the apparent friction coefficient is identified.In general,the granular mass predominantly exhibits characteristics of a dense granular flow,with the Savage number exhibiting a decreasing trend as the volume of mass increases.The process of particle breakage gradually occurs in a bottom-up manner,leading to a significant increase in the angular velocities of the rock blocks with increasing depth.The simulation results reproduce the field observations of inverse grading and source stratigraphy preservation in the deposit.We propose a disintegration index that incorporates factors such as drop height,rock mass volume,and rock strength.Our findings demonstrate a consistent linear relationship between this index and the fragmentation degree in all tested scenarios.

Isogeometric cohesive zone model for thin shell delamination analysis based on Kirchhoff-Love shell model

We present a cohesive zone model for delamination in thin shells and composite structures.The isogeometric(IGA)thin shell model is based on Kirchhoff-Love theory.Non-Uniform Rational B-Splines(NURBS)are used to discretize the exact mid-surface of the shell geometry exploiting their C 1-continuity property which avoids rotational degrees of freedom.The fracture process zone is modeled by interface elements with a cohesive law.Two numerical examples are presented to test and validate the proposed formulation in predicting the delamination behavior of composite structures.

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 effcient 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.

Genesis of continental seismogenic zone and a new fault zone model

Experiments were conducted repeatedly on Mannari granite under different temperature and confining pressure conditions. Systematic micro- and submicro-structural and mechanical analyses of granite samples deformed under 1.5 GPa (confining pressure), at 25 ℃-650 ℃ temperatures and at 2×10-6s-1 strain rate show the brittle-ductile deformation microstructures and microstructural associations similar to those observed in naturally deformed crustal rocks and minerals. Brittle fracturing and crystalline plasticity co-exist and react with each other in the brittle-ductile transition domain of the continental lithosphere. The interaction between the different mechanisms in the transitional domain results in the variation of anomalous strength values, which may best explain the genesis of the continental seismogenic zone. A new fault zone model is proposed on the basis of detailed micromechanical and microstructural analyses.

EVOLUTIONARY MODEL OF FREE ECONOMIC ZONES——Different Generations and Structural Features

Free economic zone (FEZ) has a long history and plays a more and more important role in the world economy. Most studies, however, focused on the theoretical analysis of benefit and cost as well as the economic role of FEZ in the less developed countries and little attention has been paid to the evolution of FEZ. This paper will improve the above-mentioned studies and put forward the structural and spatial evolutionary model of FEZ by analyzing the development of objectives, preferential policy, governance structure, industrial sectors and location of FEZs based on the international economic and political development. FEZs develop towards: 1) more comprehensive and macro objectives, 2) more industry-oriented and multi-preferential policies, 3) more cross-national and combination zones with administrative areas, 4) more technology-intensive and multi-industries, 5) more flexible location and larger spatial dimensions, 6) more rapid evolution and typologies, and 7) more economic integration to the host economy.

FLUID-SOLID COUPLING MATHEMATICAL MODEL OF CONTAMINANT TRANSPORT IN UNSATURATED ZONE AND ITS ASYMPTOTICAL SOLUTION

The process of contaminant transport is a problem of multicomponent and multiphase flow in unsaturated zone. Under the presupposition that gas existence affects water transport, a coupled mathematical model of contaminant transport in unsaturated zone has been established based on fluid_solid interaction mechanics theory. The asymptotical solutions to the nonlinear coupling mathematical model were accomplished by the perturbation and integral transformation method. The distribution law of pore pressure, pore water velocity and contaminant concentration in unsaturated zone has been presented under the conditions of with coupling and without coupling gas phase. An example problem was used to provide a quantitative verification and validation of the model. The asymptotical solution was compared with Faust model solution. The comparison results show reasonable agreement between asymptotical solution and Faust solution, and the gas effect and media deformation has a large impact on the contaminant transport. The theoretical basis is provided for forecasting contaminant transport and the determination of the relationship among pressure_saturation_permeability in laboratory.

Evaluation of mobility impact on urban work zones using statistical models

This work correlated the detailed work zone location and time data from the Wis LCS system with the five-min inductive loop detector data. One-sample percentile value test and two-sample Kolmogorov-Smirnov(K-S) test were applied to compare the speed and flow characteristics between work zone and non-work zone conditions. Furthermore, we analyzed the mobility characteristics of freeway work zones within the urban area of Milwaukee, WI, USA. More than 50% of investigated work zones have experienced speed reduction and 15%-30% is necessary reduced volumes. Speed reduction was more significant within and at the downstream of work zones than at the upstream.

Modelling of ultrasonic motor with dead-zone based on Hammerstein model structure

The ultrasonic motor (USM) possesses heavy nonlinearities which vary with driving conditions and load-dependent characteristics such as the dead-zone. In this paper, an identification method for the rotary travelling-wave type ultrasonic motor (RTWUSM) with dead-zone is proposed based on a modified Hammerstein model structure. The driving voltage contributing effect on the nonlinearities of the RTWUSM was transformed to the change of dynamic parameters against the driving voltage. The dead-zone of the RTWUSM is identified based upon the above transformation. Experiment results showed good agreement be- tween the output of the proposed model and actual measured output.

Finite element modelling of the geodynamic processes of the Central Andes subduction zone:A Reference Model

This paper presents preliminary results of three-dimensional thermo mechanical finite-element models of a parameter study to compute the current temperature and stress distribution in the subduction zone of the central Andes(16°S-26°S) up to a depth of 400 km, the bottom of the asthenosphere. For this purpose a simulation running over c. 50,000 years will be realized based on the geometry of a generic subduction zone and an elasto-viscoplastic Drucker-Prager rheology. The kinematic and thermal boundary conditions as well as the rheological parameters represent the current state of the study area.In future works the model will be refined using a systematic study of physical parameters in order to estimate the influence of the main parameters(e.g. viscosity, fault friction, velocity, shear heating) on the results of the reference model presented here. The reference model is kept as simple as possible to be able to estimate the influence of the parameters in future studies in the best possible way, whilst minimizing computational time.

Simulation of a compartment fire using a zone model

This paper presents the zone modeling analysis of a single compartment flashover fire. Two criteria are ap- plied in the model to judge the onset of ignition for different combustibles. By calculating the total received energy through radiation or the surface temperature of the combus- tible, the fire growth can be quantitatively determined. The improved zone fire model shows the influence of different combustibles upon the fire growth. This model is better than the traditional zone model because the common criteria of flashover, i.e. an upper layer temperature of 600℃ and the heat radiation intensity received by the floor of 20 kW/m2, have not been applied in it.

Prediction of mushy zone permeability of Al-4.5wt%Cu alloy during solidification by phase field model and CFD simulation

Liquid permeability of the mushy zone is important for porosity formation during the solidification process. In order to investigate the permeability of the mushy zone, an integrated model was developed by combining the phase field model and computational fluid dynamics (CFD) model. The three-dimensional multigrain dendrite morphology was obtained by using the phase field model. Subsequently, the computer-aided design (CAD) geometry and mesh were generated based on calculated dendrite morphologies. Finally, the permeability of the dendritic mushy zone was obtained by solving the Navier-Stokes and continuity equations in ANSYS Fluent software. As an example, the dendritic mushy zone permeability of Al-4.5wt%Cu alloy and its relationship with the solid fractions were studied in detail. The predicted permeability data can be input to the solidification model on a greater length scale for macro segregation and porosity simulations.

Parametric modeling on spatial effect of excavation-damaged zone of underground cavern

Regarding excavation-damaged zone (EDZ) around underground opening as non-homogeneous rockmass with spatial deterioration effect on stuffiness and strength, a parametric model of EDZ using radius-displacement-dependent deformation modulus (RDDM) was proposed. Considering the nonlinearity characteristic of deformation and locality otherness of surrounding rock, deterioration parameter field of deformation modulus of rockmass around opening was quantitatively calculated through a given function. Applicability for multi-cavern condition and parameter sensibility of the model was analyzed by numerical experiments using synthetic data. Furthermore, the model was applied to identify EDZ of underground caverns of Pubugou hydropower station by calculating deterioration parameter field. Based on the parametric analysis of spatial effect and geological investigation, it is recognized that large radial deformation of deep fractured rock at the spandrel position and insufficient supporting bolts mainly result in great deformation pressure to act on the shotcrete and cause partial crack and spalling. It is shown that deterioration parameter field along the longitudinal axis of main powerhouse is evidently non-homogeneous in space and distributes exponentially along the radius from the opening. The model provides a simple and convenient way to identify the EDZ in the working state for rapid construction feedback analysis and support optimization of underground cavern from quantitative point of view and also aids in interpreting monitoring displacements and estimating support requirements.