This paper proposes a case study in the control of a heavy oil pyrolysis/cracking furnace with a newly extended U-model based pole placement controller(U-PPC). The major work of the paper includes: 1) establishing a c...This paper proposes a case study in the control of a heavy oil pyrolysis/cracking furnace with a newly extended U-model based pole placement controller(U-PPC). The major work of the paper includes: 1) establishing a control oriented nonlinear dynamic model with Naphtha cracking and thermal dynamics; 2) analysing a U-model(i.e., control oriented prototype) representation of various popular process model sets; 3)designing the new U-PPC to enhance the control performance in pole placement and stabilisation; 4) taking computational bench tests to demonstrate the control system design and performance with a user-friendly step by step procedure.展开更多
The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to appl...The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to apply the calibrated model to simulating the formation of excavation damage zone(EDZ)around underground excavations.Firstly,a new cohesive crack model is implemented into the universal distinct element code(UDEC)to control the fracturing behaviour of materials under various loading modes.Next,a methodology for calibration of the components of the UDEC-Voronoi model is discussed.The role of connectivity of induced microcracks on increasing the permeability of laboratory-scale samples is investigated.The calibrated samples are used to investigate the influence of pore fluid pressure on weakening the drained strength of the laboratory-scale rock.The validity of the Terzaghi’s effective stress law for the drained peak strength of low-porosity rock is tested by performing a series of biaxial compression test simulations.Finally,the evolution of damage and pore pressure around two unsupported circular tunnels in crystalline granitic rock is studied.展开更多
A grain-based distinct element model featuring three-dimensional (3D) Voronoi tessellations (randompoly-crystals) is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in pol...A grain-based distinct element model featuring three-dimensional (3D) Voronoi tessellations (randompoly-crystals) is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional (2D) codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength (UCS) tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model's macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped (elongated) grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic (containing a fabric) rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
It is well known that quantitative estimation of slip distributions on fault plane is one of the most important issues for earthquake source inversion related to the fault rupture process. The characteristics of slip ...It is well known that quantitative estimation of slip distributions on fault plane is one of the most important issues for earthquake source inversion related to the fault rupture process. The characteristics of slip distribution on the main fault play a fundamental role to control strong ground motion pattern. A large amount of works have also suggested that variable slip models inverted from longer period ground motion recordings are relevant for the prediction of higher frequency ground motions. Zhang et al. (Chin J Geophys 56:1412-1417, 2013) and Wang et al. (Chin J Geophys 56:1408-1411,2013) published their source inversions for the fault rupturing process soon after the April 20, 2013 Lushan earthquake in Sichuan, China. In this study, first, we synthesize two forward source slip models: the value of maximum slip, fault dimension, size, and dimension of major asperities, and comer wave number obtained from Wang's model is adopted to constrain the gen- eration of k-2 model and crack model. Next, both inverted and synthetic slip models are employed to simulate the ground motions for the Lushan earthquake based on the stochastic finite-fault method. In addition, for a comparison purpose, a stochastic slip model and another k-2 model (k 2 model II) with 2 times value of comer wave number of the original k-2 model (k 2 model I) are also employed for simulation for Lushan event. The simulated results characterized by Modified Mer- calli Intensity (MMI) show that the source slip models based on the inverted and synthetic slip distributions could capture many basic features associated with the ground motion patterns. Moreover, the simulated MMI distributions reflect the rupture directivity effect and the influence of the shallow velocity structure well. On the other hand, the simulated MMI bystochastic slip model and k 2 model II is apparently higher than observed intensity. By contrast, our simulation results show that the higher frequency ground motion is sensitive to the degree of slip roughness; therefore, we suggest that, for realistic ground- motion simulations due to future earthquake, it is imperative to properly estimate the slip roughness distribution.展开更多
With the theory of subcritical crack growth, we can deduce the fundamental equation of regional seismicity acceleration model. Applying this model to intraplate earthquake regions, we select three earthquake subplates...With the theory of subcritical crack growth, we can deduce the fundamental equation of regional seismicity acceleration model. Applying this model to intraplate earthquake regions, we select three earthquake subplates: North China Subplate, Chuan Dian Block and Xinjiang Subplate, and divide the three subplates into seven researched regions by the difference of seismicity and tectonic conditions. With the modified equation given by Sornette and Sammis (1995), we analysis the seismicity of each region. To those strong earthquakes already occurred in these region, the model can give close fitting of magnitude and occurrence time, and the result in this article indicates that the seismicity acceleration model can also be used for describing the seismicity of intraplate. In the article, we give the magnitude and occurrence time of possible strong earthquakes in Shanxi, Ordos, Bole Tuokexun, Ayinke Wuqia earthquake regions. In the same subplate or block, the earthquake periods for each earthquake region are similar in time interval. The constant α in model can be used to describe the intensity of regional seismicity, and for the Chinese Mainland, α is 0.4 generally. To the seismicity in Taiwan and other regions with complex tectonic conditions, the model does not fit well at present.展开更多
Based on the Hamiltonian governing equations of plane elasticity for sectorial domain, the variable separation and eigenfunction expansion techniques were employed to develop a novel analytical finite element for the ...Based on the Hamiltonian governing equations of plane elasticity for sectorial domain, the variable separation and eigenfunction expansion techniques were employed to develop a novel analytical finite element for the fictitious crack model in fracture mechanics of concrete. The new analytical element can be implemented into FEM program systems to solve fictitious crack propagation problems for concrete cracked plates with arbitrary shapes and loads. Numerical results indicate that the method is more efficient and accurate than ordinary finite element method.展开更多
The propagation for the model I crack in aluminum single crystal has been directly studied by in-situ TEM observation.The equation of energy barrier of the dislocation building-up and emission at the model I crack tip...The propagation for the model I crack in aluminum single crystal has been directly studied by in-situ TEM observation.The equation of energy barrier of the dislocation building-up and emission at the model I crack tip has been established by means of Peierls-Nabarro dislocation model and starting from angle of energy.By means of calculation,the critical value of spontaneous emission of the dislocations from tip of the model I crack was obtained.展开更多
In the last decades,a new type of distress has been observed more and more frequently on asphalt pavements.This distress,ascribable to fatigue failure,has been named top-down cracking(TDC) because it consists in longi...In the last decades,a new type of distress has been observed more and more frequently on asphalt pavements.This distress,ascribable to fatigue failure,has been named top-down cracking(TDC) because it consists in longitudinal cracks that initiate on the pavement surface and then propagate downwards.A series of surveys recently carried out on Italian motorways highlighted that TDC can affect up to 20%-30% of the slow traffic lane.Therefore,in order to achieve a better understanding of such distress,this paper reviews causes,models and experimental tools and highlights future challenges for TDC.The literature review indicates that TDC can evolve on the pavement surface in three stages(i.e.,single crack,sister cracks,alligator cracking) and,below a certain depth,the cracks can form angles of 20°-40° with respect to the vertical plane.Even though multiple factors contribute to TDC development,thick pavements are more likely to fail due to TDC induced by tire-pavement contact stresses,especially in the presence of open-graded friction courses(OGFCs).Moreover,in literature there are several TDC models based on mechanics(e.g.,fracture mechanics or continuum damage mechanics) which allow a rigorous study of crack initiation and propagation.Future challenges include the identification of a reliable and feasible test method,among those proposed in literature,to study the TDC performance of asphalt mixtures and the implementation of TDC in pavement management systems(PMSs) through the definition of criteria for TDC recognition in the field as well as for the rehabilitation depth evaluation.Finally,more research is needed for open-graded asphalt mixtures,which present critical drawbacks in terms of TDC.展开更多
The initiation and propagation of failure in intact rock are a matter of fundamental importance in rock engineering. At low confining pressures, tensile fracturing initiates in samples at 40%-60% of the uniaxial compr...The initiation and propagation of failure in intact rock are a matter of fundamental importance in rock engineering. At low confining pressures, tensile fracturing initiates in samples at 40%-60% of the uniaxial compressive strength and as loading continues, and these tensile fractures increase in density, ultimately coalescing and leading to strain localization and macro-scale shear failure of the samples. The Griffith theory of brittle failure provides a simplified model and a useful basis for discussion of this process. The Hoek-Brown failure criterion provides an acceptable estimate of the peak strength for shear failure but a cutoff has been added for tensile conditions. However, neither of these criteria adequately explains the progressive coalition of tensile cracks and the final shearing of the specimens at higher confining stresses. Grain-based numerical models, in which the grain size distributions as well as the physical properties of the component grains of the rock are incorporated, have proved to be very useful in studying these more complex fracture processes.展开更多
In situ tensile tests in a transmission electronmicroscope(TEM)show that dislocations emitted from a modelⅡcrack tipwill form a inverse piled-up group after equilibrium or a doublepiled-up group when they meet a obst...In situ tensile tests in a transmission electronmicroscope(TEM)show that dislocations emitted from a modelⅡcrack tipwill form a inverse piled-up group after equilibrium or a doublepiled-up group when they meet a obstruction, e.g., grain boundary orsecond phase. A microcrack can initiates in front of the piled-upgroup of dislocations. Micromechanics analysis shows thatdislocations emitted from a modelⅡblunt crack tip can form a inversepiled-up or double piled-up group, depending upon the applied stressintensity factor K_Ⅱa, lattice friction stressτ_f and the distanceof the obstruction from the crack tip L.展开更多
A number of aftershocks of the May 10th 1997, Zirkuh (Ghaen-Birjand) destructive earthquake have been used to investigate the anisotropy in the upper crust by observing shear wave splitting. Particle motion diagram an...A number of aftershocks of the May 10th 1997, Zirkuh (Ghaen-Birjand) destructive earthquake have been used to investigate the anisotropy in the upper crust by observing shear wave splitting. Particle motion diagram and aspect ratio methods were used as two different approaches to obtain splitting parameters. Clear shear wave splitting was observed on the records of the selected aftershocks, indicating that the media in the region was highly anisotropic. By using particle motion method, the direction of fast shear wave was found 22°N±19°E, while the delay time between the fast and slow shear waves was obtained to be (65±16) ms. By aspect ratio method, the direction of fast shear wave was determined to be 35°N±18°E and the delay time between fast and slow shear waves was found to be (49±10) ms. For a simple horizontal layer with a thickness about 5 km and uniformly distributed anisotropy, a stress aligned cracks model was used and the result was interpreted in terms of vertical aligned cracks in the direction of N22°E, having a density about 0.01. It is assumed that cracks are fluid-filled since they are located in the upper crust. Finally, by using Hudson cracks model for three crack densities 0.005, 0.01, 0.03, the velocity curves of shear wave were plotted as a function of angle between the symmetrical axis of cracks and the azimuth of source to receiver. It was concluded that when shear wave was polarized parallel to the crack surface, the velocity was uniform, but the velocity curve varied clearly if shear wave was polarized perpendicular to the crack surface.展开更多
To determine the solutions of the well-known problem of a finite width strip with single edge crack,some results on elasto-plastic fracture analysis for metallic foams are reported.Meanwhile,in order to discuss and pu...To determine the solutions of the well-known problem of a finite width strip with single edge crack,some results on elasto-plastic fracture analysis for metallic foams are reported.Meanwhile,in order to discuss and put an insight into the nonlinear fracture analysis,the Dugdale model for plastic deformation of this configuration for metallic foams is recommended and solved.Combining the asymptotic solution with the Dugdale model and elastic solution,the stress field in the plastic zone and the size of the plastic zone are expressed as analytical forms.Based on Williams expansion method,the estimate of the scale factor is also completed and analyzed.In view of these analytical solutions,the results show the scale factor is a useful parameter for the fracture theory of metallic foams.展开更多
The numerical simulation of internal features,such as inclusions and voids,is important to analyze their impact on the performance of composite materials.However,the complex geometries of internal features and the ind...The numerical simulation of internal features,such as inclusions and voids,is important to analyze their impact on the performance of composite materials.However,the complex geometries of internal features and the induced continuous-discontinuous(C-D)deformation fields are challenges to their numerical simulation.In this study,a 3D approach using a simple mesh to simulate irregular internal geometries is developed for the first time.With the help of a developed voxel crack model,image models that are efficient when recording complex geometries are directly imported into the simulation.Surface reconstructions,which are usually labor-intensive,are excluded from this approach.Moreover,using image models as the geometric input,image processing techniques are applied to detect material interfaces and develop contact pairs.Then,the C-D deformations of the complex internal features are directly calculated based on the numerical manifold method.The accuracy and convergence of the developed3D approach are examined based on multiple benchmarks.Successful 3D C-D simulation of sandstones with naturally formed complex microfeatures demonstrates the capability of the developed approach.展开更多
This study proposes an algorithm of embedding cohesive elements in Abaqus and develops the computer code to model 3D complex cragk propagation in quasi-brittle materials in a relatively easy and efficient manner. The ...This study proposes an algorithm of embedding cohesive elements in Abaqus and develops the computer code to model 3D complex cragk propagation in quasi-brittle materials in a relatively easy and efficient manner. The cohesive elements with softening traction-separation relations and damage initiation and evolution laws are embedded between solid elements in regions of interest in the initial mesh to model potential cracks. The initial mesh can consist of tetrahedrons, wedges, bricks or a mixture of these elements. Neither remeshing nor objective crack propagation criteria are needed. Four examples of concrete specimens, including a wedgesplitting test, a notched beam under torsion, a pull-out test of an anchored cylinder and a notched beam under impact, were modelled and analysed. The simulated crack propagation processes and load-displacement curves agreed well with test results or other numerical simulations for all the examples using initial meshes with reasonable densities. Making use of Abaqus's rich pre/post- processing functionalities and powerful standard/explicit solvers, the developed method offers a practical tool for engineering analysts to model complex 3D fracture problems.展开更多
This paper investigates the influence of crack orientation on damage quantification using Lamb wave in plate structures. Finite element simulation is performed to acquire Lamb wave signal responses for different confi...This paper investigates the influence of crack orientation on damage quantification using Lamb wave in plate structures. Finite element simulation is performed to acquire Lamb wave signal responses for different configurations of crack orientations and crack lengths. Two Lamb wave features, namely the normalized amplitude and the phase change, are used as damage sensitive features to develop a crack size quantification model. A hypothesis based on the geometrical influence on signal features is proposed, and the crack size quantification model incorporating the orientation angle is established using the hypothesis. An index of Probability of Reliable Quantification(PRQ) is proposed to evaluate the performance of the model. The index can be used to determine the sizing risk in terms of probabilities. A realistic aluminum plate is used to obtain the experimental data using piezoelectric(PZT) wafer-type sensors around a center through crack. The experimental data are used to validate the overall method. Results indicate that the proposed model can yield reliable results for size quantification of inclined cracks.展开更多
An extended dielectric crack model is proposed to capture the effects of the physical properties of crack interior on crack-tip thermoelectroelastic fields.The typical crack-face boundary conditions can be retrieved b...An extended dielectric crack model is proposed to capture the effects of the physical properties of crack interior on crack-tip thermoelectroelastic fields.The typical crack-face boundary conditions can be retrieved by considering the limiting cases of electrical permeability and thermal conductivity inside a crack.Making use of the Fourier transform technique,the problem of a thermopiezoelectric strip weakened by a Griffith crack is investigated and transformed to solve the system of the second kind Fredholm integral equations with Cauchy kernel.The Lobatto-Chebyshev collocation method is used to form a nonlinear system of algebraic equations,which is solved by elaborating on an algorithm.The crack-tip thermoelectroelastic fields are determined by using the approximate solutions.Numerical simulations are carried out to show the variations of the fracture parameters of concern under applied thermoelectromechanical loads,the physical properties of the dielectric medium inside the crack and the geometry of the cracked thermopiezoelectric strip.Some comparisons with the experimental results are reported to reveal the effectiveness of the extended dielectric crack model.展开更多
Three-dimensional crack closure correction methods are investigated in this paper.The fatigue crack growth tests of surface cracks in 14MnNbq steel for bridge plate subjected to tensile and bending loadings are system...Three-dimensional crack closure correction methods are investigated in this paper.The fatigue crack growth tests of surface cracks in 14MnNbq steel for bridge plate subjected to tensile and bending loadings are systematically conducted.The experimentally measured fatigue crack growth rates of surface cracks are compared with those of through-thickness cracks in detail.It is found that the crack growth rates of surface cracks are lower than those of through-thickness cracks.In order to correct their differences in fatigue crack growth rates, a dimensionless crack closure correction model is proposed.Although this correction model is determined only by the experimental data of surface cracks under tensile loading with a constant ratio R=0.05, it can correlate the surface crack growth rates with reasonable accuracy under tensile and bending loadings with various stress ratios ranging from 0 to 0.5.Furthermore, predictions of fatigue life and crack aspect ratio for surface cracks are discussed, and the predicted results are also compared with those obtained from other prediction approaches.Comparison results show that the proposed crack closure correction model gives better prediction of fatigue life than other models.展开更多
The formation and propagation of the popular off-corner subsurface cracks in bloom continuous casting were investigated through thermo-mechanical analysis using three coupled thermo-mechanical models. A two-dimensiona...The formation and propagation of the popular off-corner subsurface cracks in bloom continuous casting were investigated through thermo-mechanical analysis using three coupled thermo-mechanical models. A two-dimensional thermo-elasto-visco-plastic finite element model was developed to predict the mould gap evolution, temperature profiles and deformation behavior of the solidified shell in the mould region. Then, a three-dimensional model was adopted to calculate the shell growth, tempera- ture history and the development of stresses and strains of the shell in the following secondary cooling zones. Finally, another three-dimensional model was used to analyze the stress distributions in the straightening region, The results showed that the off-corner cracks in the shell originated from the mould owing to the tensile strain developed in the crack sensitive regions of the solidification front, and they could be driven deeper by the possible severe surface temperature rebound and the extensive tensile stress in the secondary cooling zone, especially upon the straightening operation of the bloom casting. It is revealed that more homogenous shell temperature and thickness can be obtained through optimization of mould corner radius, casting speed and secondary cooling scheme, which help to decrease stress and strain concentration and therefore prevent the initiation of the cracks.展开更多
The well-known fatigue crack growth (FCG) curves are two-parameter dependents. The range of the stress intensity factor AK and the load ratio R are the parameters normally used for describing these curves. For engin...The well-known fatigue crack growth (FCG) curves are two-parameter dependents. The range of the stress intensity factor AK and the load ratio R are the parameters normally used for describing these curves. For engineering purposes, the mathematical representation of these curves should be integrated between the initial and final crack sizes in order to obtain the safety factors for stresses and life. First of all, it is necessary to reduce the dependence of the FCG curves to only one parameter. AK is almost always selected and, in these conditions, considered as the crack driving force. Using experimental data from literature, the present paper shows how to perform multiple regression analyses using the traditional Walker approach and the more recent unified approach. The correlations so obtained are graphically analyzed in three dimensions. Numerical examples of crack growth analysis for cracks growing under nominal stresses of constant amplitude in smooth and notched geometries are performed, assuming an identical material component as that of the available experimental data. The resulting curves of crack size versus number of cycles (a vs. N) are then compared. The two models give approximately the same (a vs. N) curves in both geometries. Differences between the behaviors of the (a vs. N) curves in smooth and notched geometries are highlighted, and the reasons for these particular behaviors are discussed.展开更多
基金partially supported by the National Natural Science Foundation of China(61273188,61473312)Taishan Scholar Construction Engineering Special Funding of Shandong
文摘This paper proposes a case study in the control of a heavy oil pyrolysis/cracking furnace with a newly extended U-model based pole placement controller(U-PPC). The major work of the paper includes: 1) establishing a control oriented nonlinear dynamic model with Naphtha cracking and thermal dynamics; 2) analysing a U-model(i.e., control oriented prototype) representation of various popular process model sets; 3)designing the new U-PPC to enhance the control performance in pole placement and stabilisation; 4) taking computational bench tests to demonstrate the control system design and performance with a user-friendly step by step procedure.
文摘The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to apply the calibrated model to simulating the formation of excavation damage zone(EDZ)around underground excavations.Firstly,a new cohesive crack model is implemented into the universal distinct element code(UDEC)to control the fracturing behaviour of materials under various loading modes.Next,a methodology for calibration of the components of the UDEC-Voronoi model is discussed.The role of connectivity of induced microcracks on increasing the permeability of laboratory-scale samples is investigated.The calibrated samples are used to investigate the influence of pore fluid pressure on weakening the drained strength of the laboratory-scale rock.The validity of the Terzaghi’s effective stress law for the drained peak strength of low-porosity rock is tested by performing a series of biaxial compression test simulations.Finally,the evolution of damage and pore pressure around two unsupported circular tunnels in crystalline granitic rock is studied.
文摘A grain-based distinct element model featuring three-dimensional (3D) Voronoi tessellations (randompoly-crystals) is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional (2D) codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength (UCS) tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model's macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped (elongated) grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic (containing a fabric) rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
基金supported by the CAS/SAFEA International Partnership Program for Creative Research Teams(KZZD-EW-TZ-19)
文摘It is well known that quantitative estimation of slip distributions on fault plane is one of the most important issues for earthquake source inversion related to the fault rupture process. The characteristics of slip distribution on the main fault play a fundamental role to control strong ground motion pattern. A large amount of works have also suggested that variable slip models inverted from longer period ground motion recordings are relevant for the prediction of higher frequency ground motions. Zhang et al. (Chin J Geophys 56:1412-1417, 2013) and Wang et al. (Chin J Geophys 56:1408-1411,2013) published their source inversions for the fault rupturing process soon after the April 20, 2013 Lushan earthquake in Sichuan, China. In this study, first, we synthesize two forward source slip models: the value of maximum slip, fault dimension, size, and dimension of major asperities, and comer wave number obtained from Wang's model is adopted to constrain the gen- eration of k-2 model and crack model. Next, both inverted and synthetic slip models are employed to simulate the ground motions for the Lushan earthquake based on the stochastic finite-fault method. In addition, for a comparison purpose, a stochastic slip model and another k-2 model (k 2 model II) with 2 times value of comer wave number of the original k-2 model (k 2 model I) are also employed for simulation for Lushan event. The simulated results characterized by Modified Mer- calli Intensity (MMI) show that the source slip models based on the inverted and synthetic slip distributions could capture many basic features associated with the ground motion patterns. Moreover, the simulated MMI distributions reflect the rupture directivity effect and the influence of the shallow velocity structure well. On the other hand, the simulated MMI bystochastic slip model and k 2 model II is apparently higher than observed intensity. By contrast, our simulation results show that the higher frequency ground motion is sensitive to the degree of slip roughness; therefore, we suggest that, for realistic ground- motion simulations due to future earthquake, it is imperative to properly estimate the slip roughness distribution.
文摘With the theory of subcritical crack growth, we can deduce the fundamental equation of regional seismicity acceleration model. Applying this model to intraplate earthquake regions, we select three earthquake subplates: North China Subplate, Chuan Dian Block and Xinjiang Subplate, and divide the three subplates into seven researched regions by the difference of seismicity and tectonic conditions. With the modified equation given by Sornette and Sammis (1995), we analysis the seismicity of each region. To those strong earthquakes already occurred in these region, the model can give close fitting of magnitude and occurrence time, and the result in this article indicates that the seismicity acceleration model can also be used for describing the seismicity of intraplate. In the article, we give the magnitude and occurrence time of possible strong earthquakes in Shanxi, Ordos, Bole Tuokexun, Ayinke Wuqia earthquake regions. In the same subplate or block, the earthquake periods for each earthquake region are similar in time interval. The constant α in model can be used to describe the intensity of regional seismicity, and for the Chinese Mainland, α is 0.4 generally. To the seismicity in Taiwan and other regions with complex tectonic conditions, the model does not fit well at present.
文摘Based on the Hamiltonian governing equations of plane elasticity for sectorial domain, the variable separation and eigenfunction expansion techniques were employed to develop a novel analytical finite element for the fictitious crack model in fracture mechanics of concrete. The new analytical element can be implemented into FEM program systems to solve fictitious crack propagation problems for concrete cracked plates with arbitrary shapes and loads. Numerical results indicate that the method is more efficient and accurate than ordinary finite element method.
文摘The propagation for the model I crack in aluminum single crystal has been directly studied by in-situ TEM observation.The equation of energy barrier of the dislocation building-up and emission at the model I crack tip has been established by means of Peierls-Nabarro dislocation model and starting from angle of energy.By means of calculation,the critical value of spontaneous emission of the dislocations from tip of the model I crack was obtained.
基金sponsored by Autostrade per l’Italia S.p.A. (Italy)which gave both financial and technical support within the framework of the Highway Pavement Evolutive Research (HiPER) project。
文摘In the last decades,a new type of distress has been observed more and more frequently on asphalt pavements.This distress,ascribable to fatigue failure,has been named top-down cracking(TDC) because it consists in longitudinal cracks that initiate on the pavement surface and then propagate downwards.A series of surveys recently carried out on Italian motorways highlighted that TDC can affect up to 20%-30% of the slow traffic lane.Therefore,in order to achieve a better understanding of such distress,this paper reviews causes,models and experimental tools and highlights future challenges for TDC.The literature review indicates that TDC can evolve on the pavement surface in three stages(i.e.,single crack,sister cracks,alligator cracking) and,below a certain depth,the cracks can form angles of 20°-40° with respect to the vertical plane.Even though multiple factors contribute to TDC development,thick pavements are more likely to fail due to TDC induced by tire-pavement contact stresses,especially in the presence of open-graded friction courses(OGFCs).Moreover,in literature there are several TDC models based on mechanics(e.g.,fracture mechanics or continuum damage mechanics) which allow a rigorous study of crack initiation and propagation.Future challenges include the identification of a reliable and feasible test method,among those proposed in literature,to study the TDC performance of asphalt mixtures and the implementation of TDC in pavement management systems(PMSs) through the definition of criteria for TDC recognition in the field as well as for the rehabilitation depth evaluation.Finally,more research is needed for open-graded asphalt mixtures,which present critical drawbacks in terms of TDC.
文摘The initiation and propagation of failure in intact rock are a matter of fundamental importance in rock engineering. At low confining pressures, tensile fracturing initiates in samples at 40%-60% of the uniaxial compressive strength and as loading continues, and these tensile fractures increase in density, ultimately coalescing and leading to strain localization and macro-scale shear failure of the samples. The Griffith theory of brittle failure provides a simplified model and a useful basis for discussion of this process. The Hoek-Brown failure criterion provides an acceptable estimate of the peak strength for shear failure but a cutoff has been added for tensile conditions. However, neither of these criteria adequately explains the progressive coalition of tensile cracks and the final shearing of the specimens at higher confining stresses. Grain-based numerical models, in which the grain size distributions as well as the physical properties of the component grains of the rock are incorporated, have proved to be very useful in studying these more complex fracture processes.
基金[This project was financially supported by the special funds for the Major State Basic Research Projects (G19990650), Visiting Scholar Foundation of Key Laboratory of the Ministry of Education PRC and by the NNSF of China.]
文摘In situ tensile tests in a transmission electronmicroscope(TEM)show that dislocations emitted from a modelⅡcrack tipwill form a inverse piled-up group after equilibrium or a doublepiled-up group when they meet a obstruction, e.g., grain boundary orsecond phase. A microcrack can initiates in front of the piled-upgroup of dislocations. Micromechanics analysis shows thatdislocations emitted from a modelⅡblunt crack tip can form a inversepiled-up or double piled-up group, depending upon the applied stressintensity factor K_Ⅱa, lattice friction stressτ_f and the distanceof the obstruction from the crack tip L.
文摘A number of aftershocks of the May 10th 1997, Zirkuh (Ghaen-Birjand) destructive earthquake have been used to investigate the anisotropy in the upper crust by observing shear wave splitting. Particle motion diagram and aspect ratio methods were used as two different approaches to obtain splitting parameters. Clear shear wave splitting was observed on the records of the selected aftershocks, indicating that the media in the region was highly anisotropic. By using particle motion method, the direction of fast shear wave was found 22°N±19°E, while the delay time between the fast and slow shear waves was obtained to be (65±16) ms. By aspect ratio method, the direction of fast shear wave was determined to be 35°N±18°E and the delay time between fast and slow shear waves was found to be (49±10) ms. For a simple horizontal layer with a thickness about 5 km and uniformly distributed anisotropy, a stress aligned cracks model was used and the result was interpreted in terms of vertical aligned cracks in the direction of N22°E, having a density about 0.01. It is assumed that cracks are fluid-filled since they are located in the upper crust. Finally, by using Hudson cracks model for three crack densities 0.005, 0.01, 0.03, the velocity curves of shear wave were plotted as a function of angle between the symmetrical axis of cracks and the azimuth of source to receiver. It was concluded that when shear wave was polarized parallel to the crack surface, the velocity was uniform, but the velocity curve varied clearly if shear wave was polarized perpendicular to the crack surface.
基金Supported by the National Natural Science Foundation of China(10972035)
文摘To determine the solutions of the well-known problem of a finite width strip with single edge crack,some results on elasto-plastic fracture analysis for metallic foams are reported.Meanwhile,in order to discuss and put an insight into the nonlinear fracture analysis,the Dugdale model for plastic deformation of this configuration for metallic foams is recommended and solved.Combining the asymptotic solution with the Dugdale model and elastic solution,the stress field in the plastic zone and the size of the plastic zone are expressed as analytical forms.Based on Williams expansion method,the estimate of the scale factor is also completed and analyzed.In view of these analytical solutions,the results show the scale factor is a useful parameter for the fracture theory of metallic foams.
基金supported by the National Natural Science Foundation of China(Grant Nos.41807277,42172306,and U1965204)the Natural Science Foundation of Hebei Province(Grant No.D2019202440)。
文摘The numerical simulation of internal features,such as inclusions and voids,is important to analyze their impact on the performance of composite materials.However,the complex geometries of internal features and the induced continuous-discontinuous(C-D)deformation fields are challenges to their numerical simulation.In this study,a 3D approach using a simple mesh to simulate irregular internal geometries is developed for the first time.With the help of a developed voxel crack model,image models that are efficient when recording complex geometries are directly imported into the simulation.Surface reconstructions,which are usually labor-intensive,are excluded from this approach.Moreover,using image models as the geometric input,image processing techniques are applied to detect material interfaces and develop contact pairs.Then,the C-D deformations of the complex internal features are directly calculated based on the numerical manifold method.The accuracy and convergence of the developed3D approach are examined based on multiple benchmarks.Successful 3D C-D simulation of sandstones with naturally formed complex microfeatures demonstrates the capability of the developed approach.
基金supported by EPSRC UK(No.EP/F00656X/1)Xiangting Su's one-year visit to the University of Liverpoosupported by the China Scholarship Council and the National Natural Science Foundation of China(No.50579081).
文摘This study proposes an algorithm of embedding cohesive elements in Abaqus and develops the computer code to model 3D complex cragk propagation in quasi-brittle materials in a relatively easy and efficient manner. The cohesive elements with softening traction-separation relations and damage initiation and evolution laws are embedded between solid elements in regions of interest in the initial mesh to model potential cracks. The initial mesh can consist of tetrahedrons, wedges, bricks or a mixture of these elements. Neither remeshing nor objective crack propagation criteria are needed. Four examples of concrete specimens, including a wedgesplitting test, a notched beam under torsion, a pull-out test of an anchored cylinder and a notched beam under impact, were modelled and analysed. The simulated crack propagation processes and load-displacement curves agreed well with test results or other numerical simulations for all the examples using initial meshes with reasonable densities. Making use of Abaqus's rich pre/post- processing functionalities and powerful standard/explicit solvers, the developed method offers a practical tool for engineering analysts to model complex 3D fracture problems.
基金supported by Science Challenge Project of China (No. TZ2018007)National Natural Science Foundation of China (Nos. 11872088 and51975546)。
文摘This paper investigates the influence of crack orientation on damage quantification using Lamb wave in plate structures. Finite element simulation is performed to acquire Lamb wave signal responses for different configurations of crack orientations and crack lengths. Two Lamb wave features, namely the normalized amplitude and the phase change, are used as damage sensitive features to develop a crack size quantification model. A hypothesis based on the geometrical influence on signal features is proposed, and the crack size quantification model incorporating the orientation angle is established using the hypothesis. An index of Probability of Reliable Quantification(PRQ) is proposed to evaluate the performance of the model. The index can be used to determine the sizing risk in terms of probabilities. A realistic aluminum plate is used to obtain the experimental data using piezoelectric(PZT) wafer-type sensors around a center through crack. The experimental data are used to validate the overall method. Results indicate that the proposed model can yield reliable results for size quantification of inclined cracks.
基金The work was supported by the National Natural Science Foundation of China(Nos.11872155 and 11362002)the Guangxi Natural Science Foundation(No.2016GXNSFAA380261)the innovation project of Guangxi Graduate Education(YCSW2019045).
文摘An extended dielectric crack model is proposed to capture the effects of the physical properties of crack interior on crack-tip thermoelectroelastic fields.The typical crack-face boundary conditions can be retrieved by considering the limiting cases of electrical permeability and thermal conductivity inside a crack.Making use of the Fourier transform technique,the problem of a thermopiezoelectric strip weakened by a Griffith crack is investigated and transformed to solve the system of the second kind Fredholm integral equations with Cauchy kernel.The Lobatto-Chebyshev collocation method is used to form a nonlinear system of algebraic equations,which is solved by elaborating on an algorithm.The crack-tip thermoelectroelastic fields are determined by using the approximate solutions.Numerical simulations are carried out to show the variations of the fracture parameters of concern under applied thermoelectromechanical loads,the physical properties of the dielectric medium inside the crack and the geometry of the cracked thermopiezoelectric strip.Some comparisons with the experimental results are reported to reveal the effectiveness of the extended dielectric crack model.
文摘Three-dimensional crack closure correction methods are investigated in this paper.The fatigue crack growth tests of surface cracks in 14MnNbq steel for bridge plate subjected to tensile and bending loadings are systematically conducted.The experimentally measured fatigue crack growth rates of surface cracks are compared with those of through-thickness cracks in detail.It is found that the crack growth rates of surface cracks are lower than those of through-thickness cracks.In order to correct their differences in fatigue crack growth rates, a dimensionless crack closure correction model is proposed.Although this correction model is determined only by the experimental data of surface cracks under tensile loading with a constant ratio R=0.05, it can correlate the surface crack growth rates with reasonable accuracy under tensile and bending loadings with various stress ratios ranging from 0 to 0.5.Furthermore, predictions of fatigue life and crack aspect ratio for surface cracks are discussed, and the predicted results are also compared with those obtained from other prediction approaches.Comparison results show that the proposed crack closure correction model gives better prediction of fatigue life than other models.
文摘The formation and propagation of the popular off-corner subsurface cracks in bloom continuous casting were investigated through thermo-mechanical analysis using three coupled thermo-mechanical models. A two-dimensional thermo-elasto-visco-plastic finite element model was developed to predict the mould gap evolution, temperature profiles and deformation behavior of the solidified shell in the mould region. Then, a three-dimensional model was adopted to calculate the shell growth, tempera- ture history and the development of stresses and strains of the shell in the following secondary cooling zones. Finally, another three-dimensional model was used to analyze the stress distributions in the straightening region, The results showed that the off-corner cracks in the shell originated from the mould owing to the tensile strain developed in the crack sensitive regions of the solidification front, and they could be driven deeper by the possible severe surface temperature rebound and the extensive tensile stress in the secondary cooling zone, especially upon the straightening operation of the bloom casting. It is revealed that more homogenous shell temperature and thickness can be obtained through optimization of mould corner radius, casting speed and secondary cooling scheme, which help to decrease stress and strain concentration and therefore prevent the initiation of the cracks.
文摘The well-known fatigue crack growth (FCG) curves are two-parameter dependents. The range of the stress intensity factor AK and the load ratio R are the parameters normally used for describing these curves. For engineering purposes, the mathematical representation of these curves should be integrated between the initial and final crack sizes in order to obtain the safety factors for stresses and life. First of all, it is necessary to reduce the dependence of the FCG curves to only one parameter. AK is almost always selected and, in these conditions, considered as the crack driving force. Using experimental data from literature, the present paper shows how to perform multiple regression analyses using the traditional Walker approach and the more recent unified approach. The correlations so obtained are graphically analyzed in three dimensions. Numerical examples of crack growth analysis for cracks growing under nominal stresses of constant amplitude in smooth and notched geometries are performed, assuming an identical material component as that of the available experimental data. The resulting curves of crack size versus number of cycles (a vs. N) are then compared. The two models give approximately the same (a vs. N) curves in both geometries. Differences between the behaviors of the (a vs. N) curves in smooth and notched geometries are highlighted, and the reasons for these particular behaviors are discussed.