A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FE...A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FEM is described and compared with ES-FEM.A practical modification of IES-FEM is then introduced that used the technique employed by ES-FEM for the nodal strain calculation.The differences in the strain computation among ES-FEM,IES-FEM,and FEM are then discussed.The modified IES-FEM exhibited superior performance in displacement and a slight advantage in stress compared to FEM using the same mesh according to the results obtained from both the regular and irregular elements.The robustness of the IES-FEM to severely deformed meshes was also verified.展开更多
The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES...The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.展开更多
We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to...We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to achieve exceptional computational efficiency and accuracy.The workflow is demonstrated through the modeling of wireline electromagnetic propagation resistivity logs,where the measured responses exhibit a highly nonlinear relationship with formation properties.The motivation for this research is the need for advanced modeling al-gorithms that are fast enough for use in modern quantitative interpretation tools,where thousands of simulations may be required in iterative inversion processes.The proposed algorithm achieves a remarkable enhancement in performance,being up to 3000 times faster than the finite element method alone when utilizing a GPU.While still ensuring high accuracy,this makes it well-suited for practical applications when reliable payzone assessment is needed in complex environmental scenarios.Furthermore,the algorithm’s efficiency positions it as a promising tool for stochastic Bayesian inversion,facilitating reliable uncertainty quantification in subsurface property estimation.展开更多
Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant i...Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant importance.The traditional finite element method(FEM)remains one of the primary approaches for addressing such issues.However,the application of FEM typically necessitates the use of a fine finite element mesh to accurately capture the heterogeneous properties of the materials and meet the required computational precision,which inevitably leads to a reduction in computational efficiency.To enhance the computational accuracy and efficiency of the FEM for heterogeneous multi-field coupling problems,this study presents the coupling magneto-electro-elastic multiscale finite element method(CM-MsFEM)for heterogeneous MEE structures.Unlike the conventional multiscale FEM(MsFEM),the proposed algorithm simultaneously constructs displacement,electric,and magnetic potential multiscale basis functions to address the heterogeneity of the corresponding parameters.The macroscale formulation of CM-MsFEM was derived,and the macroscale/microscale responses of the problems were obtained through up/downscaling calculations.Evaluation using numerical examples analyzing the transient behavior of heterogeneous MEE structures demonstrated that the proposed method outperforms traditional FEM in terms of both accuracy and computational efficiency,making it an appropriate choice for numerically modeling the dynamics of heterogeneous MEE structures.展开更多
The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX a...The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX are determined by the grain level deformation and sequentially coupled.In order to fully capture the microstructure and mechanical property evolution,a crystal plasticity finite element based modelling method for DRX and SRX is proposed in the current work.The grain level deformation is calculated with crystal plasticity which is coupled with the recrystallization model straightforwardly,and both the grain deformation and microstructure evolution are updated simultaneously.The proposed method is validated with discontinuous DRX experiments and the effects of initial deformation conditions are well-captured.Two controversial mechanisms for recrystallization microstructure evolution,i.e.oriented nucleation and growth selection,are discussed in the current framework with the advantages of accurate grain level deformation and interaction predictions.Furthermore,the sequentially coupled DRX and SRX are modelled seamlessly in the current work which provides a critical method for fully integrated thermo-mechanical processes analysis.展开更多
In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a ge...In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a generalized nonlinear Stokes problem of displacement vector field related to pseudo pressure and a diffusion problem of other pseudo pressure fields.Secondly,a fully discrete multiphysics finite element method is performed to solve the reformulated system numerically.Thirdly,existence and uniqueness of the weak solution of the reformulated model and stability analysis and optimal convergence order for the multiphysics finite element method are proven theoretically.Lastly,numerical tests are given to verify the theoretical results.展开更多
Skin panels on supersonic vehicles are subjected to aero-thermo-acoustic loads,resulting in a well-known multi-physics dynamic problem.The high-frequency dynamic response of these panels significantly impacts the stru...Skin panels on supersonic vehicles are subjected to aero-thermo-acoustic loads,resulting in a well-known multi-physics dynamic problem.The high-frequency dynamic response of these panels significantly impacts the structural safety of supersonic vehicles,but it has been rarely investigated.Given that existing methods are inefficient for high-frequency dynamic analysis in multi-physics fields,the present work addresses this challenge by proposing a Stochastic Energy Finite Element Method(SEFEM).SEFEM uses energy density instead of displacement to describe the dynamic response,thereby significantly enhancing its efficiency.In SEFEM,the effects of aerodynamic and thermal loads on the energy propagation characteristics are studied analytically and incorporated into the energy density governing equation.These effects are also considered when calculating the input power generated by the acoustic load,and two effective approaches named Frequency Response Function Method(FRFM)and Mechanical Impedance Method(MIM)are developed accordingly and integrated into SEFEM.The good accuracy,applicability,and high efficiency of the proposed SEFEM are demonstrated through numerical simulations performed on a two-dimensional panel under aero-thermoacoustic loads.Additionally,the effects and underlying mechanisms of aero-thermo-acoustic loads on the high-frequency response are explored.This work not only presents an efficient approach for predicting high-frequency dynamic response of panels subjected to aero-thermo-acoustic loads,but also provides insights into the high-frequency dynamic characteristics in multi-physics fields.展开更多
Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump mate...Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.展开更多
Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres wi...Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres with artificial inherent pores are an important type of lightweight proppant,enabling their transport to distant fracture extremities and enhancing fracture conductivity.However,the focus frequently gravitates towards the low-density advantage,often overlooking the pore geometry impacts on compressive strength by traditional strength evaluation.This paper numerically bypasses such limitations by using a combined finite and discrete element method(FDEM)considering experimental results.The mesh size of the model undergoes validation,followed by the calibration of cohesive element parameters via the single particle compression test.The stimulation elucidates that proppants with a smaller pore size(40μm)manifest crack propagation evolution at a more rapid pace in comparison to their larger-pore counterparts,though the influence of pore diameter on overall strength is subtle.The inception of pores not only alters the trajectory of crack progression but also,with an increase in porosity,leads to a discernible decline in proppant compressive strength.Intriguingly,upon crossing a porosity threshold of 10%,the decrement in strength becomes more gradual.A denser congregation of pores accelerates crack propagation,undermining proppant robustness,suggesting that under analogous conditions,hollow proppants might not match the strength of their porous counterparts.This exploration elucidates the underlying mechanisms of proppant failure from a microstructural perspective,furnishing pivotal insights that may guide future refinements in the architectural design of porous proppant.展开更多
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.展开更多
In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow e...In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.展开更多
Controlled nuclear fusion represents a significant solution for future clean energy,with ion cyclotron range of frequency(ICRF)heating emerging as one of the most promising technologies for heating the fusion plasma.T...Controlled nuclear fusion represents a significant solution for future clean energy,with ion cyclotron range of frequency(ICRF)heating emerging as one of the most promising technologies for heating the fusion plasma.This study primarily presents a self-developed 2D ion cyclotron resonance antenna electromagnetic field solver(ICRAEMS)code implemented on the MATLAB platform,which solves the electric field wave equation by using the finite element method,establishing perfectly matched layer(PML)boundary conditions,and post-processing the electromagnetic field data.This code can be utilized to facilitate the design and optimization processes of antennas for ICRF heating technology.Furthermore,this study examines the electric field distribution and power spectrum associated with various antenna phases to investigate how different antenna configurations affect the electromagnetic field propagation and coupling characteristics.展开更多
Within the roadway advanced detection methods, DC resistivity method has an extensive application because of its simple principle and operation. Numerical simulation of the effect of focusing current on advanced detec...Within the roadway advanced detection methods, DC resistivity method has an extensive application because of its simple principle and operation. Numerical simulation of the effect of focusing current on advanced detection was carried out using a three-dimensional finite element method （FEM）, meanwhile the electric-field distribution of the point source and nine-point power source were calculated and analyzed with the same electric charges. The results show that the nine-point power source array has a very good ability to focus, and the DC focus method can be used to predict the aquifer abnormality body precisely. By comparing the FEM modelling results with physical simulation results from soil sink, it is shown that the accuracy of forward simulation meets the requirement and the artificial disturbance from roadway has no impact on the DC focus method.展开更多
The equilibrium equations and the functional for tubing buckling in arbitrary straight wells are derived. The entire buckling process of tubing in deviated wells is analyzed for the first time by utilizing the finite ...The equilibrium equations and the functional for tubing buckling in arbitrary straight wells are derived. The entire buckling process of tubing in deviated wells is analyzed for the first time by utilizing the finite element method. The effects of gravity and torques on the buckling are included in the analyses and the calculated results are well compared with existing solutions. It is shown that the buckling only occurs at the lower portion of the tubing where the axial load is the largest, and the contact force of the well, the bending moment of the tubing and the buckling displacement of this portion vary periodically. The buckling spreads upwards from the bit with the increase of axial load. There is no buckling at the upper portion of the tubing where the bending moment is zero. And the contact force of this section increases only slightly with the increase of the axial load. With the increase of the deviation angle, the length of buckling portion and buckling displacement amplitude decrease, the contact force increases with the increase of load at the upper portion and its amplitude decreases at the lower buckling section, the bending moment remains zero at the upper portion and its amplitude decreases at the lower buckling portion. The buckling displacement increases with the increase of the torque, but the increment is very small.展开更多
A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underex...A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.展开更多
A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic gr...A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.展开更多
A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations....A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.展开更多
A combined method for the fast 3-D finite element modeling of defect responses in nondestructive testing of electromagnetics is presented. The method consists of three numerical techniques: zoom-in technique, differen...A combined method for the fast 3-D finite element modeling of defect responses in nondestructive testing of electromagnetics is presented. The method consists of three numerical techniques: zoom-in technique, difference field technique and iterative solution technique. Utilizing the zoom-in technique, the computational zone focuses on a relatively small domain around the defect. Employing the difference field technique, the axisymmetrical field solution corresponding to the case with no defect can be used to simplify the mesh generation and obtain the modeling results quickly. Using the iterative solution technique, the matrix equation system in the 3-D finite element modeling of nondestructive probe signals can easily be solved. The sample calculation shows that the presented method is highly effective and can consequently save significant computer resources.展开更多
A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element...A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.展开更多
The Runge-Kutta discontinuous Galerkin finite element method （RK-DGFEM） is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and ...The Runge-Kutta discontinuous Galerkin finite element method （RK-DGFEM） is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and it is explicit in the time domain. Consequently it is a best mixture of FEM and finite volume method （FVM）. RK-DGFEM can obtain local high-order accuracy by using high-order polynomial basis. Numerical experiments of transverse magnetic （TM） wave propagation in a 2-D resonator are performed. A high-order Lagrange polynomial basis is adopted. Numerical results agree well with analytical solution. And different order Lagrange interpolation polynomial basis impacts on simulation result accuracy are discussed. Computational results indicate that the accuracy is evidently improved when the order of interpolation basis is increased. Finally, L^2 errors of different order polynomial basis in RK-DGFEM are presented. Computational results show that L^2 error declines exponentially as the order of basis increases.展开更多
基金the National Natural Science Foundation of China(No.11672238)the 111 Project(No.BP0719007)the Shaanxi Province Natural Science Foundation(No.2020JZ-06)for the financial support.
文摘A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FEM is described and compared with ES-FEM.A practical modification of IES-FEM is then introduced that used the technique employed by ES-FEM for the nodal strain calculation.The differences in the strain computation among ES-FEM,IES-FEM,and FEM are then discussed.The modified IES-FEM exhibited superior performance in displacement and a slight advantage in stress compared to FEM using the same mesh according to the results obtained from both the regular and irregular elements.The robustness of the IES-FEM to severely deformed meshes was also verified.
基金funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under Grant number 107.02-2019.330。
文摘The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.
基金financially supported by the Russian federal research project No.FWZZ-2022-0026“Innovative aspects of electro-dynamics in problems of exploration and oilfield geophysics”.
文摘We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to achieve exceptional computational efficiency and accuracy.The workflow is demonstrated through the modeling of wireline electromagnetic propagation resistivity logs,where the measured responses exhibit a highly nonlinear relationship with formation properties.The motivation for this research is the need for advanced modeling al-gorithms that are fast enough for use in modern quantitative interpretation tools,where thousands of simulations may be required in iterative inversion processes.The proposed algorithm achieves a remarkable enhancement in performance,being up to 3000 times faster than the finite element method alone when utilizing a GPU.While still ensuring high accuracy,this makes it well-suited for practical applications when reliable payzone assessment is needed in complex environmental scenarios.Furthermore,the algorithm’s efficiency positions it as a promising tool for stochastic Bayesian inversion,facilitating reliable uncertainty quantification in subsurface property estimation.
基金supported by the National Natural Science Foundation of China(Grant Nos.42102346,42172301).
文摘Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant importance.The traditional finite element method(FEM)remains one of the primary approaches for addressing such issues.However,the application of FEM typically necessitates the use of a fine finite element mesh to accurately capture the heterogeneous properties of the materials and meet the required computational precision,which inevitably leads to a reduction in computational efficiency.To enhance the computational accuracy and efficiency of the FEM for heterogeneous multi-field coupling problems,this study presents the coupling magneto-electro-elastic multiscale finite element method(CM-MsFEM)for heterogeneous MEE structures.Unlike the conventional multiscale FEM(MsFEM),the proposed algorithm simultaneously constructs displacement,electric,and magnetic potential multiscale basis functions to address the heterogeneity of the corresponding parameters.The macroscale formulation of CM-MsFEM was derived,and the macroscale/microscale responses of the problems were obtained through up/downscaling calculations.Evaluation using numerical examples analyzing the transient behavior of heterogeneous MEE structures demonstrated that the proposed method outperforms traditional FEM in terms of both accuracy and computational efficiency,making it an appropriate choice for numerically modeling the dynamics of heterogeneous MEE structures.
基金supported by the National Natural Science Foundation of China(Nos.52105384 and U2141215).
文摘The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX are determined by the grain level deformation and sequentially coupled.In order to fully capture the microstructure and mechanical property evolution,a crystal plasticity finite element based modelling method for DRX and SRX is proposed in the current work.The grain level deformation is calculated with crystal plasticity which is coupled with the recrystallization model straightforwardly,and both the grain deformation and microstructure evolution are updated simultaneously.The proposed method is validated with discontinuous DRX experiments and the effects of initial deformation conditions are well-captured.Two controversial mechanisms for recrystallization microstructure evolution,i.e.oriented nucleation and growth selection,are discussed in the current framework with the advantages of accurate grain level deformation and interaction predictions.Furthermore,the sequentially coupled DRX and SRX are modelled seamlessly in the current work which provides a critical method for fully integrated thermo-mechanical processes analysis.
基金Supported by the National Natural Science Foundation of China(Grant Nos.12371393,11971150 and 11801143)Natural Science Foundation of Henan Province(Grant No.242300421047).
文摘In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a generalized nonlinear Stokes problem of displacement vector field related to pseudo pressure and a diffusion problem of other pseudo pressure fields.Secondly,a fully discrete multiphysics finite element method is performed to solve the reformulated system numerically.Thirdly,existence and uniqueness of the weak solution of the reformulated model and stability analysis and optimal convergence order for the multiphysics finite element method are proven theoretically.Lastly,numerical tests are given to verify the theoretical results.
基金financially supported by the National Natural Science Foundation of China(Nos.12302228 and 12372170)。
文摘Skin panels on supersonic vehicles are subjected to aero-thermo-acoustic loads,resulting in a well-known multi-physics dynamic problem.The high-frequency dynamic response of these panels significantly impacts the structural safety of supersonic vehicles,but it has been rarely investigated.Given that existing methods are inefficient for high-frequency dynamic analysis in multi-physics fields,the present work addresses this challenge by proposing a Stochastic Energy Finite Element Method(SEFEM).SEFEM uses energy density instead of displacement to describe the dynamic response,thereby significantly enhancing its efficiency.In SEFEM,the effects of aerodynamic and thermal loads on the energy propagation characteristics are studied analytically and incorporated into the energy density governing equation.These effects are also considered when calculating the input power generated by the acoustic load,and two effective approaches named Frequency Response Function Method(FRFM)and Mechanical Impedance Method(MIM)are developed accordingly and integrated into SEFEM.The good accuracy,applicability,and high efficiency of the proposed SEFEM are demonstrated through numerical simulations performed on a two-dimensional panel under aero-thermoacoustic loads.Additionally,the effects and underlying mechanisms of aero-thermo-acoustic loads on the high-frequency response are explored.This work not only presents an efficient approach for predicting high-frequency dynamic response of panels subjected to aero-thermo-acoustic loads,but also provides insights into the high-frequency dynamic characteristics in multi-physics fields.
基金the financial support provided by MHRD,Govt.of IndiaCoal India Limited for providing financial assistance for the research(Project No.CIL/R&D/01/73/2021)the partial financial support provided by the Ministry of Education,Government of India,under SPARC project(Project No.P1207)。
文摘Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.
基金the financial support provided by Tianfu Yongxing Laboratory Organized Research Project Funding(No.2023CXXM01)the ARC linkage program(No.LP200100420).
文摘Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres with artificial inherent pores are an important type of lightweight proppant,enabling their transport to distant fracture extremities and enhancing fracture conductivity.However,the focus frequently gravitates towards the low-density advantage,often overlooking the pore geometry impacts on compressive strength by traditional strength evaluation.This paper numerically bypasses such limitations by using a combined finite and discrete element method(FDEM)considering experimental results.The mesh size of the model undergoes validation,followed by the calibration of cohesive element parameters via the single particle compression test.The stimulation elucidates that proppants with a smaller pore size(40μm)manifest crack propagation evolution at a more rapid pace in comparison to their larger-pore counterparts,though the influence of pore diameter on overall strength is subtle.The inception of pores not only alters the trajectory of crack progression but also,with an increase in porosity,leads to a discernible decline in proppant compressive strength.Intriguingly,upon crossing a porosity threshold of 10%,the decrement in strength becomes more gradual.A denser congregation of pores accelerates crack propagation,undermining proppant robustness,suggesting that under analogous conditions,hollow proppants might not match the strength of their porous counterparts.This exploration elucidates the underlying mechanisms of proppant failure from a microstructural perspective,furnishing pivotal insights that may guide future refinements in the architectural design of porous proppant.
基金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.
基金supported by the Natural Science Foundation of Shandong Province(ZR2021MA019)the National Natural Science Foundation of China(11871312)。
文摘In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.
基金Project supported by the National MCF Energy R&D Program(Grant No.2022YFE03190100)the National Natural Science Foundation of China(Grant Nos.12422513,12105035,and U21A20438)the Xiaomi Young Talents Program.
文摘Controlled nuclear fusion represents a significant solution for future clean energy,with ion cyclotron range of frequency(ICRF)heating emerging as one of the most promising technologies for heating the fusion plasma.This study primarily presents a self-developed 2D ion cyclotron resonance antenna electromagnetic field solver(ICRAEMS)code implemented on the MATLAB platform,which solves the electric field wave equation by using the finite element method,establishing perfectly matched layer(PML)boundary conditions,and post-processing the electromagnetic field data.This code can be utilized to facilitate the design and optimization processes of antennas for ICRF heating technology.Furthermore,this study examines the electric field distribution and power spectrum associated with various antenna phases to investigate how different antenna configurations affect the electromagnetic field propagation and coupling characteristics.
基金Project(41174103)supported by the National Natural Science Foundation of ChinaProject(20110162130008)supported by the PhD Program Foundation of Ministry of Education of ChinaProject(2011BAB04B08)supported by the National Key Technology R&D Program during the 12th Five-Year Plan of China
文摘Within the roadway advanced detection methods, DC resistivity method has an extensive application because of its simple principle and operation. Numerical simulation of the effect of focusing current on advanced detection was carried out using a three-dimensional finite element method （FEM）, meanwhile the electric-field distribution of the point source and nine-point power source were calculated and analyzed with the same electric charges. The results show that the nine-point power source array has a very good ability to focus, and the DC focus method can be used to predict the aquifer abnormality body precisely. By comparing the FEM modelling results with physical simulation results from soil sink, it is shown that the accuracy of forward simulation meets the requirement and the artificial disturbance from roadway has no impact on the DC focus method.
文摘The equilibrium equations and the functional for tubing buckling in arbitrary straight wells are derived. The entire buckling process of tubing in deviated wells is analyzed for the first time by utilizing the finite element method. The effects of gravity and torques on the buckling are included in the analyses and the calculated results are well compared with existing solutions. It is shown that the buckling only occurs at the lower portion of the tubing where the axial load is the largest, and the contact force of the well, the bending moment of the tubing and the buckling displacement of this portion vary periodically. The buckling spreads upwards from the bit with the increase of axial load. There is no buckling at the upper portion of the tubing where the bending moment is zero. And the contact force of this section increases only slightly with the increase of the axial load. With the increase of the deviation angle, the length of buckling portion and buckling displacement amplitude decrease, the contact force increases with the increase of load at the upper portion and its amplitude decreases at the lower buckling section, the bending moment remains zero at the upper portion and its amplitude decreases at the lower buckling portion. The buckling displacement increases with the increase of the torque, but the increment is very small.
文摘A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.
基金Projects(51161011,11364024)supported by the National Natural Science Foundation of ChinaProject(1204GKCA065)supported by the Key Technology R&D Program of Gansu Province,China+1 种基金Project(201210)supported by the Fundamental Research Funds for the Universities of Gansu Province,ChinaProject(J201304)supported by the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China
文摘A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.
基金Supported by the National Natural Science Foundation of China(50976072,51106099,10902070)the Leading Academic Discipline Project of Shanghai Municipal Education Commission(J50501)the Science Foundation for the Excellent Youth Scholar of Higher Education of Shanghai(slg09003)~~
文摘A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.
文摘A combined method for the fast 3-D finite element modeling of defect responses in nondestructive testing of electromagnetics is presented. The method consists of three numerical techniques: zoom-in technique, difference field technique and iterative solution technique. Utilizing the zoom-in technique, the computational zone focuses on a relatively small domain around the defect. Employing the difference field technique, the axisymmetrical field solution corresponding to the case with no defect can be used to simplify the mesh generation and obtain the modeling results quickly. Using the iterative solution technique, the matrix equation system in the 3-D finite element modeling of nondestructive probe signals can easily be solved. The sample calculation shows that the presented method is highly effective and can consequently save significant computer resources.
基金Project supported by the National Basic Research Program of China (973Project) (No.2002CB412709) and the National Natural Science Foundation of China (Nos.50278012,10272027,19832010)
文摘A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.
文摘The Runge-Kutta discontinuous Galerkin finite element method （RK-DGFEM） is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and it is explicit in the time domain. Consequently it is a best mixture of FEM and finite volume method （FVM）. RK-DGFEM can obtain local high-order accuracy by using high-order polynomial basis. Numerical experiments of transverse magnetic （TM） wave propagation in a 2-D resonator are performed. A high-order Lagrange polynomial basis is adopted. Numerical results agree well with analytical solution. And different order Lagrange interpolation polynomial basis impacts on simulation result accuracy are discussed. Computational results indicate that the accuracy is evidently improved when the order of interpolation basis is increased. Finally, L^2 errors of different order polynomial basis in RK-DGFEM are presented. Computational results show that L^2 error declines exponentially as the order of basis increases.
