This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is e...This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.展开更多
A numerical method is put forward in this paper, using the boundary element method （BEM） to model 3D terrain effects on magnetotelluric （MT） surveys, Using vector integral theory and electromagnetic field boundary...A numerical method is put forward in this paper, using the boundary element method （BEM） to model 3D terrain effects on magnetotelluric （MT） surveys, Using vector integral theory and electromagnetic field boundary conditions, the boundary problem of two electromagnetic fields in the upper half space （air） and lower half space （earth medium） was transformed into two vector integral equations just related to the topography ： one magnetic equation for computing the magnetic field and the other electrical equation for computing the electrical field. The topography integral is decomposed into a series of integrals in a triangle element. For the integral in a triangle element, we suppose that the electromagnetic field in it is the stack of the electromagnetic field in the homogeneous earth and the topography response which is a constant; so the computation becomes simple, convenient and highly accurate. By decomposition and computation, each vector integral equation can be calculated by solving three linear equations that are related to the three Cartesian directions. The matrix of these linear equations is diagonally dominant and can be solved using the Symmetric Successive Over-Relaxation （SSOR） method. The apparent resistivity curve of MT on two 3D terrains calculated by BEM is shown in this paper.展开更多
The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical ...The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical profiles within the basin are used to analyze the generation and propagation of surface waves and their relation to the basin structures in both the time and frequency domain.The amplification effects are analyzed by the distribution of peak ground velocity(PGV)and cumulative kinetic energy(Ek) in the basin.The results show that in the 3D basin,the excitation of the fundamental and higher surface wave modes are similar to that of the 2D model.Small bowls in the basin have great influence on the amplification and distribution of strong ground motion,due to their lateral resonances when the wavelengths of the lateral surface waves are comparable to the size of the bowls.Obvious basin edge effects can be seen at the basin edge closer to the source for constructive interference between direct body waves and the basin-induced surface waves.The Ek distribution maps show very large values in small bowls and some corners in the basin due to the interference of waves propagating in different directions.A high impedance contrast model can excite more surface wave modes,resulting in longer shaking durations as well as more complex seismograms and PGV and Ek distributions.展开更多
For higher accuracy in simulating the transformation of three dimensional waves, in consideration of the advantages of constant panels and linear elements, a combined boundary elements is applied in this research. The...For higher accuracy in simulating the transformation of three dimensional waves, in consideration of the advantages of constant panels and linear elements, a combined boundary elements is applied in this research. The method can be used to remove the transverse vibration due to the accumulation of computational errors. A combined boundary condition of sponge layer and Sommerfeld radiation condition is used to remove the reflected waves from the computing domain. By following the water particle on the water surface, the third order Stokes wave transform is simulated by the numerical wave flume technique. The computed results are in good agreement with theoretical ones.展开更多
In order to study the failure patterns and strength of 3D braided composites from the microscopic view, the damage propagation under tensile loading steps in three kinds of unit cells is simulated. The homogenization ...In order to study the failure patterns and strength of 3D braided composites from the microscopic view, the damage propagation under tensile loading steps in three kinds of unit cells is simulated. The homogenization formula of micro-stress and the solving approach of finite element method are given firstly. A criterion is presented to determine the damage and its pattern of each element, and then the stiffness degradation method based on Murakami's geometric damage theory is used to simulate the status of damage under tensile loading steps for three kinds of unit cells. It can be seen that the damage percentage and damage pattern of damaged unit cell are totally different for different kind of unit cells. More damaged elements are observed for face cell and corner cell than that for body cell. It is also observed that the damage firstly occurs at the area of face cell, which agrees well with experimental results. It is verified that considering the effects of face and corner cells are important for the damage and strength analysis of 3D braided composites.展开更多
Applying stiffness migration method,a 3D finite element mechanical model is established to simulate the excavation and advance processes.By using 3D nonlinear finite element method,the tunnel boring machine(TBM) excav...Applying stiffness migration method,a 3D finite element mechanical model is established to simulate the excavation and advance processes.By using 3D nonlinear finite element method,the tunnel boring machine(TBM) excavation process is dynamically simulated to analyze the stress and strain field status of surrounding rock and segment.The maximum tensile stress of segment ring caused by tunnel construction mainly lies in arch bottom and presents zonal distribution.The stress increases slightly and limitedly in the course of excavation.The maximum and minimum displacements of segment,manifesting as zonal distribution,distribute in arch bottom and vault respectively.The displacements slightly increase with the advance of TBM and gradually tend to stability.展开更多
Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass ...Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.展开更多
In numerical simulation of the mechanical responses and acoustic emission(AE)characteristics of rocks under cyclic loading,the impacts of compositional heterogeneities of mineral grains have barely been considered.Thi...In numerical simulation of the mechanical responses and acoustic emission(AE)characteristics of rocks under cyclic loading,the impacts of compositional heterogeneities of mineral grains have barely been considered.This will lead to a poor reproduction of rock’s behaviors in terms of stress-strain relationship and micro-seismic characteristics in numerical simulation.This work aims to analyze and reveal the impact of parameter heterogeneity on the rock’s fatigue and micro-seismic properties based on PFC3D.Two distribution patterns(uniform and Weibull distributions),are implemented to assign four critical parameters(i.e.tensile strength,cohesion,parallel bond stiffness and linear stiffness)for 32 sets of numerical schemes.The results show that the models with high heterogeneity of tensile strength and cohesion can better reproduce the stress-strain relationship as well as the patterns of cumulative AE counts and energy magnitude.The evolution of the proportion of three-level AE events in the laboratory test is consistent with the numerical results when the highly heterogeneous tensile strength and cohesion are distributed.The numerical results can provide practical guidance to the PFC-based modeling of rock heterogeneity when exposed to multi-level cyclic loading and AE monitoring.展开更多
In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform ...In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.展开更多
The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good ac...The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.展开更多
On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate an...On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.展开更多
The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for th...The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for the solution of 3 D frictional contact problems, the key part is the determination of the tangential slip states at the contact points, and a great amount of computing work is needed for a high accuracy result. A new method based on a combination of programming and iteration methods, which are respectively known as two main kinds of methods for contact analysis, was put forward to deal with 3 D elastic-plastic contact problems. Numerical results demonstrate the efficiency of the algorithm illustrated here.展开更多
The current national criteria cannot accurately reflect the multi-crack interference effect of defective gas pipelines,and thus result in conservative assessment results.In order to improve the safety assessment accur...The current national criteria cannot accurately reflect the multi-crack interference effect of defective gas pipelines,and thus result in conservative assessment results.In order to improve the safety assessment accuracy of defective gas pipelines,we compared the singular element method with the three-dimensional virtual crack closure technology(3D-VCCT),then,selected the 3D-VCCT to build a multi-crack interference model of gas pipelines,and finally,studied the interference effect of subsidiary cracks on the front points of the main crack by analyzing the variation of the interference factors between cracks.The following research results were obtained:first,as the subsidiary crack size changes,the strongest interference effect of a parallel collinear crack occurs near the surface of the crack,exhibiting an enhances effect;second,under certain conditions,the interference effect of a parallel coaxial crack is more significant at the deepest point and the surface point of the crack,and exhibits an enhanced effect at the deepest point and a weakened effect at the surface point;third,when the horizontal spacing between cracks is 6 times greater than the major semi-axis of the main crack,the interference effect of parallel collinear cracks and parallel bias cracks is very little,and multi-cracks can be simplified as a single crack for analysis.The research results not only can provide technical references for fracture analysis of multi-crack gas pipelines,but also can provide theoretical basis for the safety assessment of multi-cracks of the gas pipelines.展开更多
Current developments in magnetohydrodynamic(MHD)convection and nanofluid engineering technology have have greatly enhanced heat transfer performance in process systems,particularly through the use of carbon nanotube(C...Current developments in magnetohydrodynamic(MHD)convection and nanofluid engineering technology have have greatly enhanced heat transfer performance in process systems,particularly through the use of carbon nanotube(CNT)–based fluids that offer exceptional thermal conductivity.Despite extensive research on MHD natural convection in enclosures,the combined effects of complex obstacle geometries,magnetic fields,and CNT nanofluids in three-dimensional configurations remain insufficiently explored.This research investigates MHD natural convection of carbon nanotube(CNT)-water nanofluid within a three-dimensional cavity.The study considers an inclined cross-shaped hot obstacle,a configuration not extensively explored in previous works.The work aims to elucidate the combined effects of CNT nanofluid concentration,magnetic field strength,and obstacle inclination on fluid flow patterns and heat transfer characteristics.Numerical simulations are performed using the finite element method(FEM)based on the Galerkin Weighted Residual approach.The analysis systematically considers variations in Rayleigh number(Ra),Hartmann number(Ha),nanoparticle volume fraction(Φ),and obstacle inclination angle(θ).Results show that increasing Ra from 103 to 106 enhances convective heat transfer by up to 228%,while raising the CNT volume fraction to 4.5%improves heat transfer by about 64%.In contrast,strengthening the magnetic field from Ha=0 to Ha=100 suppresses fluid motion and reduces heat transfer by nearly 67%,whereas varying the obstacle inclination from 0○to 45○leads to a 4.6%decrease in efficiency.The addition of nanoparticles slightly increases viscosity,reducing flow intensity by 8.3%when Ha=0.Furthermore,a novel multiparametric correlation is proposed,accurately predicting the average Nusselt number as a function of Ra,Ha,ϕ,andθ,with an R2 of 0.98.These findings provide new insights into the role of geometry,magnetic effects,and nanofluids in heat transfer enhancement,offering practical guidance for the design and optimization of advanced thermal systems.展开更多
The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy ...The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.展开更多
To study the influence of roll casting process parameters on temperature and thermal-stress fields for the AZ31 magnesium alloy sheets,three-dimensional geometric and 3D finite element models for roll casting were est...To study the influence of roll casting process parameters on temperature and thermal-stress fields for the AZ31 magnesium alloy sheets,three-dimensional geometric and 3D finite element models for roll casting were established based on the symmetry of roll casting by ANSYS software.Meshing method and smart-sizing algorithm were used to divide finite element mesh in ANSYS software.A series of researches on the temperature and stress distributions during solidification process with different process parameters were done by 3D finite element method.The temperatures of both the liquid-solid two-phase zone and liquid phase zone were elevated with increasing pouring temperature.With the heat transfer coefficient increasing,the two-phase region for liquid-solid becomes smaller.With the pouring temperature increasing and the increase of casting speed,the length of two-phase zone rises.The optimized of process parameters（casting speed 2 m/min,pouring temperature 640 ℃ and heat transfer coefficient 15 kW/（m2·℃） with the water pouring at roller exit was used to produce magnesium alloy AZ31 sheet,and equiaxed grains with the average grain size of 50 μm were achieved after roll casting.The simulation results give better understanding of the temperature variation in phase transformation zone and the formation mechanism of hot cracks in plates during roll casting and help to design the optimized process parameters of roll casting for Mg alloy.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.51890912,51979025 and 52011530189).
文摘This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.
基金This paper is supported by the National Natural Science Foundation ofChina (No .40344002) .
文摘A numerical method is put forward in this paper, using the boundary element method （BEM） to model 3D terrain effects on magnetotelluric （MT） surveys, Using vector integral theory and electromagnetic field boundary conditions, the boundary problem of two electromagnetic fields in the upper half space （air） and lower half space （earth medium） was transformed into two vector integral equations just related to the topography ： one magnetic equation for computing the magnetic field and the other electrical equation for computing the electrical field. The topography integral is decomposed into a series of integrals in a triangle element. For the integral in a triangle element, we suppose that the electromagnetic field in it is the stack of the electromagnetic field in the homogeneous earth and the topography response which is a constant; so the computation becomes simple, convenient and highly accurate. By decomposition and computation, each vector integral equation can be calculated by solving three linear equations that are related to the three Cartesian directions. The matrix of these linear equations is diagonally dominant and can be solved using the Symmetric Successive Over-Relaxation （SSOR） method. The apparent resistivity curve of MT on two 3D terrains calculated by BEM is shown in this paper.
基金National Natural Science Foundation of China under Grant No.51078337,No.51108431 and No.91315301
文摘The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical profiles within the basin are used to analyze the generation and propagation of surface waves and their relation to the basin structures in both the time and frequency domain.The amplification effects are analyzed by the distribution of peak ground velocity(PGV)and cumulative kinetic energy(Ek) in the basin.The results show that in the 3D basin,the excitation of the fundamental and higher surface wave modes are similar to that of the 2D model.Small bowls in the basin have great influence on the amplification and distribution of strong ground motion,due to their lateral resonances when the wavelengths of the lateral surface waves are comparable to the size of the bowls.Obvious basin edge effects can be seen at the basin edge closer to the source for constructive interference between direct body waves and the basin-induced surface waves.The Ek distribution maps show very large values in small bowls and some corners in the basin due to the interference of waves propagating in different directions.A high impedance contrast model can excite more surface wave modes,resulting in longer shaking durations as well as more complex seismograms and PGV and Ek distributions.
基金National Natural Science Foundation of China(No.49876026)
文摘For higher accuracy in simulating the transformation of three dimensional waves, in consideration of the advantages of constant panels and linear elements, a combined boundary elements is applied in this research. The method can be used to remove the transverse vibration due to the accumulation of computational errors. A combined boundary condition of sponge layer and Sommerfeld radiation condition is used to remove the reflected waves from the computing domain. By following the water particle on the water surface, the third order Stokes wave transform is simulated by the numerical wave flume technique. The computed results are in good agreement with theoretical ones.
基金National Natural Science Foundation of China (10772115)
文摘In order to study the failure patterns and strength of 3D braided composites from the microscopic view, the damage propagation under tensile loading steps in three kinds of unit cells is simulated. The homogenization formula of micro-stress and the solving approach of finite element method are given firstly. A criterion is presented to determine the damage and its pattern of each element, and then the stiffness degradation method based on Murakami's geometric damage theory is used to simulate the status of damage under tensile loading steps for three kinds of unit cells. It can be seen that the damage percentage and damage pattern of damaged unit cell are totally different for different kind of unit cells. More damaged elements are observed for face cell and corner cell than that for body cell. It is also observed that the damage firstly occurs at the area of face cell, which agrees well with experimental results. It is verified that considering the effects of face and corner cells are important for the damage and strength analysis of 3D braided composites.
基金Supported by National Natural Science Foundation of China(No.90815019)National Key Basic Research Program of China("973" Program,No.2007CB714101)Key Project in the National Science and Technology Pillar Program during the Eleventh Five-Year Plan Period(No.2006BAB04A13)
文摘Applying stiffness migration method,a 3D finite element mechanical model is established to simulate the excavation and advance processes.By using 3D nonlinear finite element method,the tunnel boring machine(TBM) excavation process is dynamically simulated to analyze the stress and strain field status of surrounding rock and segment.The maximum tensile stress of segment ring caused by tunnel construction mainly lies in arch bottom and presents zonal distribution.The stress increases slightly and limitedly in the course of excavation.The maximum and minimum displacements of segment,manifesting as zonal distribution,distribute in arch bottom and vault respectively.The displacements slightly increase with the advance of TBM and gradually tend to stability.
基金supported by the National Natural Science Foundation of China(Grant No.11304160)the Natural Science Foundation of Jiangsu Provincial Higher Education Institutions,China(Grant No.13KJB140008)the Foundation of Nanjing University of Posts and Telecommunications,China(Grant No.NY213018)
文摘Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.
基金funded by the Funds from Joint National-Local Engineering Research Center for Safe and Precise Coal Mining(Grant No.EC2021004).
文摘In numerical simulation of the mechanical responses and acoustic emission(AE)characteristics of rocks under cyclic loading,the impacts of compositional heterogeneities of mineral grains have barely been considered.This will lead to a poor reproduction of rock’s behaviors in terms of stress-strain relationship and micro-seismic characteristics in numerical simulation.This work aims to analyze and reveal the impact of parameter heterogeneity on the rock’s fatigue and micro-seismic properties based on PFC3D.Two distribution patterns(uniform and Weibull distributions),are implemented to assign four critical parameters(i.e.tensile strength,cohesion,parallel bond stiffness and linear stiffness)for 32 sets of numerical schemes.The results show that the models with high heterogeneity of tensile strength and cohesion can better reproduce the stress-strain relationship as well as the patterns of cumulative AE counts and energy magnitude.The evolution of the proportion of three-level AE events in the laboratory test is consistent with the numerical results when the highly heterogeneous tensile strength and cohesion are distributed.The numerical results can provide practical guidance to the PFC-based modeling of rock heterogeneity when exposed to multi-level cyclic loading and AE monitoring.
基金China Postdoctoral Science Foundation Under Grant No.20100480321National Basic Research Program of China Under Grant No. 2007CB714200
文摘In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.
文摘The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.
基金supported by the Key Projects of the National Science Foundation of China(Nos.52178190,52078250,11832013)
文摘On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.
基金theNationalKeyBasicResearchSpecialFoundation (G1 9990 3 2 80 5 ) the FoundationforUniversityKeyTeacherbytheMinistryofEducationo
文摘The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for the solution of 3 D frictional contact problems, the key part is the determination of the tangential slip states at the contact points, and a great amount of computing work is needed for a high accuracy result. A new method based on a combination of programming and iteration methods, which are respectively known as two main kinds of methods for contact analysis, was put forward to deal with 3 D elastic-plastic contact problems. Numerical results demonstrate the efficiency of the algorithm illustrated here.
基金Project supported by National Science and Technology Infrastructure Program“Risk-based key technologies for preventing pipeline accidents”(No.:2011BAK06B01-11).
文摘The current national criteria cannot accurately reflect the multi-crack interference effect of defective gas pipelines,and thus result in conservative assessment results.In order to improve the safety assessment accuracy of defective gas pipelines,we compared the singular element method with the three-dimensional virtual crack closure technology(3D-VCCT),then,selected the 3D-VCCT to build a multi-crack interference model of gas pipelines,and finally,studied the interference effect of subsidiary cracks on the front points of the main crack by analyzing the variation of the interference factors between cracks.The following research results were obtained:first,as the subsidiary crack size changes,the strongest interference effect of a parallel collinear crack occurs near the surface of the crack,exhibiting an enhances effect;second,under certain conditions,the interference effect of a parallel coaxial crack is more significant at the deepest point and the surface point of the crack,and exhibits an enhanced effect at the deepest point and a weakened effect at the surface point;third,when the horizontal spacing between cracks is 6 times greater than the major semi-axis of the main crack,the interference effect of parallel collinear cracks and parallel bias cracks is very little,and multi-cracks can be simplified as a single crack for analysis.The research results not only can provide technical references for fracture analysis of multi-crack gas pipelines,but also can provide theoretical basis for the safety assessment of multi-cracks of the gas pipelines.
基金Deputyship for Research&Innovation,Ministry of Education in Saudi Arabia for funding this research work through the project number RI-44-0451.
文摘Current developments in magnetohydrodynamic(MHD)convection and nanofluid engineering technology have have greatly enhanced heat transfer performance in process systems,particularly through the use of carbon nanotube(CNT)–based fluids that offer exceptional thermal conductivity.Despite extensive research on MHD natural convection in enclosures,the combined effects of complex obstacle geometries,magnetic fields,and CNT nanofluids in three-dimensional configurations remain insufficiently explored.This research investigates MHD natural convection of carbon nanotube(CNT)-water nanofluid within a three-dimensional cavity.The study considers an inclined cross-shaped hot obstacle,a configuration not extensively explored in previous works.The work aims to elucidate the combined effects of CNT nanofluid concentration,magnetic field strength,and obstacle inclination on fluid flow patterns and heat transfer characteristics.Numerical simulations are performed using the finite element method(FEM)based on the Galerkin Weighted Residual approach.The analysis systematically considers variations in Rayleigh number(Ra),Hartmann number(Ha),nanoparticle volume fraction(Φ),and obstacle inclination angle(θ).Results show that increasing Ra from 103 to 106 enhances convective heat transfer by up to 228%,while raising the CNT volume fraction to 4.5%improves heat transfer by about 64%.In contrast,strengthening the magnetic field from Ha=0 to Ha=100 suppresses fluid motion and reduces heat transfer by nearly 67%,whereas varying the obstacle inclination from 0○to 45○leads to a 4.6%decrease in efficiency.The addition of nanoparticles slightly increases viscosity,reducing flow intensity by 8.3%when Ha=0.Furthermore,a novel multiparametric correlation is proposed,accurately predicting the average Nusselt number as a function of Ra,Ha,ϕ,andθ,with an R2 of 0.98.These findings provide new insights into the role of geometry,magnetic effects,and nanofluids in heat transfer enhancement,offering practical guidance for the design and optimization of advanced thermal systems.
基金Projects(41074085,41374118)supported by the National Natural Science Foundation of ChinaProject(20120162110015)supported by Doctoral Fund of Ministry of Education of ChinaProject(NCET-12-0551)supported by Program for New Century Excellent Talents in University,China
文摘The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.
基金Project(CSTC 2010BB4301) supported by Natural Science Foundation Project of Chongqing,ChinaProject supported by the Open Fund for Key Laboratory of Manufacture and Test Techniques for Automobile Parts of Ministry of Education Chongqing University of Technology,2003,China
文摘To study the influence of roll casting process parameters on temperature and thermal-stress fields for the AZ31 magnesium alloy sheets,three-dimensional geometric and 3D finite element models for roll casting were established based on the symmetry of roll casting by ANSYS software.Meshing method and smart-sizing algorithm were used to divide finite element mesh in ANSYS software.A series of researches on the temperature and stress distributions during solidification process with different process parameters were done by 3D finite element method.The temperatures of both the liquid-solid two-phase zone and liquid phase zone were elevated with increasing pouring temperature.With the heat transfer coefficient increasing,the two-phase region for liquid-solid becomes smaller.With the pouring temperature increasing and the increase of casting speed,the length of two-phase zone rises.The optimized of process parameters（casting speed 2 m/min,pouring temperature 640 ℃ and heat transfer coefficient 15 kW/（m2·℃） with the water pouring at roller exit was used to produce magnesium alloy AZ31 sheet,and equiaxed grains with the average grain size of 50 μm were achieved after roll casting.The simulation results give better understanding of the temperature variation in phase transformation zone and the formation mechanism of hot cracks in plates during roll casting and help to design the optimized process parameters of roll casting for Mg alloy.
