The finite element approach is used for the first time to simulate and examine the free oscillation and transient response of a visco-elastic multi-directional functionally graded porous(MFGP)skew-nanoplate,taking int...The finite element approach is used for the first time to simulate and examine the free oscillation and transient response of a visco-elastic multi-directional functionally graded porous(MFGP)skew-nanoplate,taking into account surface effects using nonlocal strain gradient hypothesis.The mechanical characteristics of the materials vary in all three directions of length,width,and thickness of the plate according to the exponential law.Additionally,it has viscoelastic behavior according to the Kelvin-Voigt model.The novelty of this paper lies in the incorporation of the spatial variability of nonlocal and lengthscale factors as additional mechanical characteristics of the material.The overall equation of motion for the plate is derived by including the classical plate hypothesis and Hamilton’s principle.A quadrilateral plate element with four nodes and six degrees of freedom is created using a non-conforming C2-level Hermitian function.This function offers precise results and rapid convergence for various forms and boundary conditions(BCs)that low-order elements cannot accomplish.The Newmark-beta direct integration technique is used to calculate the transient responses of the visco-elastic MFGP skew-nanoplate under various BCs.Furthermore,a thorough assessment of the impacts of several factors such as residual surface stress,grading indices,elastic foundation stiffness,skew angle,other geometrical parameters,and BCs on the transient responses of the viscoelastic MFGP skew-nanoplate has been uncovered.展开更多
This study develops a surrogate super-resolution(SR)framework that accelerates finite element method(FEM)-based computational fluid dynamics(CFD)using deep learning.High-resolution(HR)FEM-based CFDremains computationa...This study develops a surrogate super-resolution(SR)framework that accelerates finite element method(FEM)-based computational fluid dynamics(CFD)using deep learning.High-resolution(HR)FEM-based CFDremains computationally prohibitive for time-sensitive applications,including patient-specific aneurysm hemodynamics where rapid turnaround is valuable.The proposed pipeline learns to reconstruct HR velocity-magnitude fields fromlow-resolution(LR)FEM solutions generated under the same governing equations and boundary conditions.It consistsof three modules:(i)offline pre-training of a residual network on representative vascular geometries;(ii)lightweightfine-tuning to adapt the pretrained model to geometric variability,including patient-specific aneurysm morphologies;and(iii)an unstructured-to-structured sampling strategy with region-of-interest upsampling that concentrates resolution in flow-critical zones(e.g.,the aneurysm sac)rather than the full domain.This targeted reconstruction substantiallyreduces inference and post-processing cost while preserving key HR flow features.Experiments on cerebral aneurysmmodels show that HR velocity-magnitude fields can be recovered with accuracy comparable to direct HR simulationsat less than 1%of the direct HR simulation cost per analysis(LR simulation and SR inference),while adaptation to newgeometries requires only lightweight fine-tuning with limited target-specific HR data.While clinical endpoints andadditional variables(e.g.,pressure or wall-based metrics)are left for future work,the results indicate that the proposedsurrogate SR approach can streamline FEM-based CFD workflows toward near real-time hemodynamic analysis acrossmorphologically similar vascular models.展开更多
This paper proposes a state-of-the-art three-dimensional Voronoi cell finite element method(3D VCFEM)aimed at investigating the mechanical properties of particle-reinforced composites(PRCs)in space under different mic...This paper proposes a state-of-the-art three-dimensional Voronoi cell finite element method(3D VCFEM)aimed at investigating the mechanical properties of particle-reinforced composites(PRCs)in space under different microstructural properties.Firstly,the modified residual energy generalized function of 3D VCFEM was proposed by applying the hybrid stress element method,and the element format of the 3D Voronoi element was constructed.On this basis,the interaction between the matrix and the inclusions was considered,and the higher-order stress function including the interaction stress term was constructed.Secondly,to solve the difficulty of integrating easily due to the complexity and irregularity of the integration region in space,Delaunay tetrahedra were introduced within the 3D Voronoi element for mesh refinement.It simplified the integration process.Finally,to verify the accuracy and efficiency of the 3D VCFEM model,comparative models of 3D VCFENM and FEM were established for analysis and discussion.The stress field and strain field were compared and analyzed for the first time.An example was also given for the presence of a large number of randomly distributed inclusion particles.The results showed that under the same accuracy,3D VCFEM had the advantages of convenient mesh delineation and high computational efficiency compared with FEM,which provided a new way of thinking to analyze the actual PCRs.展开更多
In this study,a collaborative optimization framework combining artificial fish swarm algorithm(AFSA)and Voronoi cell finite element method(VCFEM)is proposed to solve the problem of the influence of microstructure dist...In this study,a collaborative optimization framework combining artificial fish swarm algorithm(AFSA)and Voronoi cell finite element method(VCFEM)is proposed to solve the problem of the influence of microstructure distribution on the macroscopic mechanical properties of particle-reinforced composites.Compared with the traditional displacement-based finite element method,VCFEM significantly improves the computational efficiency of multi-inclusion problems by means of generalized stress function and element local adaptation technology,and accurately captures the stress field discontinuity and interface stress concentration phenomenon.At the same time,AFSA abandons the dependence of traditional optimization methods on simplified models,and directly searches for the global optimal solution through discrete topological variables(particle positions).Its natural selection mechanism and swarm intelligence characteristics effectively overcome the local optimal trap.The numerical simulation results show that the proposed method exhibits high accuracy in both single-and multi-inclusion models(the error with the commercial software MARC is less than 5%),and in the complex model with 100 inclusions,the maximum Mises stress is reduced by 32.6%after optimization.The synergistic effect of the two breaks through the trade-off between efficiency and accuracy in traditional optimization,and provides the ability of both computational efficiency and global optimization for multi-scale composite material design.展开更多
Slope stability issues in geotechnical engineering significantly impact both structural safety and economic efficiency, making stability analysis a critical component of design and construction processes. The finite e...Slope stability issues in geotechnical engineering significantly impact both structural safety and economic efficiency, making stability analysis a critical component of design and construction processes. The finite element method is employed to develop numerical models of geotechnical slopes, incorporating appropriate boundary conditions and material parameters to conduct in-depth investigations into slope behavior under various loading conditions. During model development, comprehensive consideration is given to the heterogeneity of each rock layer and detailed internal structural characteristics, with meticulous examination of stress distribution patterns and deformation accumulation mechanisms under different stress environments. Research findings demonstrate that the finite element method effectively captures internal stress distributions and identifies precise locations of potential failure surfaces, providing robust technical support for early detection of slope instability risks. The study also elucidates the underlying mechanisms governing slope stability under multiple influencing factors, clarifying the operational mechanisms and interrelationships of key control parameters. Based on these findings, practical engineering mitigation measures are proposed, offering solid theoretical foundations and operational guidance for safe slope design and risk management in geotechnical engineering. These advancements enhance overall project safety, improve comprehensive disaster prevention and response systems, and contribute significantly to industry excellence.展开更多
To address challenges in architectural extensibility and cross-module collaboration of CAE software,this study proposes OPFEM(open-source Python-based finite element modeling)—an open-source framework featuring a uni...To address challenges in architectural extensibility and cross-module collaboration of CAE software,this study proposes OPFEM(open-source Python-based finite element modeling)—an open-source framework featuring a unified four-layer architecture.The geometric modeling framework achieves plug-in support for geometric kernels through an interface abstraction layer and adapter patterns,decoupling kernel-specific implementations while enabling state machine-driven interaction design and parametric sketching.The pre-processing modules establish multi-level associations among materials,sections,and geometric entities using Composite and Factory patterns,while implementing Observer pattern to ensure geometric-mesh consistency and employing finite-state machines to optimize boundary workflows.The computational modules implement a modular finite element library that decouples topology from element attributes,along with a surface boundary element technique for load conversion and task-scheduling management,validated through a cantilever beam large-deformation case.The post-processing module facilitates standardized data storage and dynamic field visualization through architecture-level standardized interface definitions and hierarchical component design.Collectively,OPFEM achieves full-process integration from parametric modeling to nonlinear solving and visualization,enhancing configuration efficiency and providing an extensible,pattern-driven solution for complex CAE challenges.展开更多
A fully nonlinear potential flow(FNPF)solver has been developed using the Finite Element Method(FEM)to simulate time-domain interactions between free-surface waves and marine structures.The ALE framework is implemente...A fully nonlinear potential flow(FNPF)solver has been developed using the Finite Element Method(FEM)to simulate time-domain interactions between free-surface waves and marine structures.The ALE framework is implemented alongside a segment spring analogy-based moving mesh strategy to accurately track evolving free surfaces and moving boundaries of floating bodies.The solver employs a preconditioned conjugate gradient method to efficiently resolve the resulting sparse,symmetric linear system at each time step.Temporal evolution is managed through a standard fourth-order Runge-Kutta scheme,while Chebyshev 5-point smoothing suppresses non-physical saw-tooth instabilities.The solver’s performance and reliability are verified through comprehensive benchmark tests,including freesurface sloshing,nonlinear wave propagation,and wave-structure interactions with submerged or floating bodies.Furthermore,the study explores a modified potential flow model incorporating a quadratic damping term to address viscous effects in gap/moonpool resonance problems.展开更多
This study presents and verifies a hybrid methodology for reliable determination of parameters in structural rheological models(Zener,Burgers,and Maxwell)describing the viscoelastic behavior of polyurethane specimens ...This study presents and verifies a hybrid methodology for reliable determination of parameters in structural rheological models(Zener,Burgers,and Maxwell)describing the viscoelastic behavior of polyurethane specimens manufactured using extrusion-based 3D printing.Through comprehensive testing,including cyclic compression at strain rates ranging from 0.12 to 120 mm/min(0%-15%strain)and creepelaxation experiments(10%-30%strain),the lumped parameters were independently determined using both analytical and numerical solutions of the models’differential equations,followed by cross-verification in additional experiments.Numerical solutions for creep and relaxation problems were obtained using finite element analysis,with the three-parameter Mooney-Rivlin model and Prony series employed to simulate elastic and viscous stress components,respectively.Energy dissipation per cycle was quantified during cyclic compression tests.The results demonstrate that all three models adequately describe material behavior within the 0%-15%strain range across various strain rates.Comparative analysis revealed the Burgers model’s superior performance in characterizing creep and stress relaxation at low strain levels.While Zener and Burgers model parameters from uniaxial compression showed limited applicability for energy dissipation calculations,the generalized Maxwell model effectively captured viscoelastic properties across different strain rates.Notably,parameters derived from creep tests provided a more universal assessment of dissipative properties due to optimization based on characteristic curve regions.Both parameter sets described polyurethane’s elastic-hysteretic behavior with approximately 20%error,proving significantly more accurate than the linear strain-time dependence hypothesis.Finite element analysis(FEA)complemented numerical modeling by demonstrating that while the generalized Maxwell model effectively describes initial rapid stress-strain changes,FEA provides superior characterization of steady-state processes.This computational approach yields more physically representative results compared to simplified analytical solutions,despite certain limitations in transient analysis.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solv...This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solving a specific numerical problem under the scope of the linear finite element method(LFEM),so the method is termed computational method for analytical solutions with finite elements(CMAS-FE).The primary objective of the CMAS-FE is to construct analytical expressions for displacements and reaction forces at nodes,as well as for strains and stresses at elemental quadrature points,all of which are formulated as infinite series solutions of various orders of Poisson’s ratios.Like the conventional LFEM,the CMAS-FE forms global sparse linear equations,but the Young’s modulus and Poisson’s ratio remain variables(or symbols).By employing a direct inverse method to solve these symbolic linear systems,an analytical expression of the displacement field can be constructed.The CMAS-FE is validated via patch and bending tests,which demonstrate convergence with mesh and term refine-ment.Furthermore,the CMAS-FE is applied to obtain the bending stiffness of a beam structure and to estimate an approximate stress intensity factor for a straight crack within a square-shaped plate.展开更多
In modern ZnO varistors,traditional aging mechanisms based on increased power consumption are no longer relevant due to reduced power consumption during DC aging.Prolonged exposure to both AC and DC voltages results i...In modern ZnO varistors,traditional aging mechanisms based on increased power consumption are no longer relevant due to reduced power consumption during DC aging.Prolonged exposure to both AC and DC voltages results in increased leakage current,decreased breakdown voltage,and lower nonlinearity,ultimately compromising their protective performance.To investigate the evolution in electrical properties during DC aging,this work developed a finite element model based on Voronoi networks and conducted accelerated aging tests on commercial varistors.Throughout the aging process,current-voltage characteristics and Schottky barrier parameters were measured and analyzed.The results indicate that when subjected to constant voltage,current flows through regions with larger grain sizes,forming discharge channels.As aging progresses,the current focus increases on these channels,leading to a decline in the varistor’s overall performance.Furthermore,analysis of the Schottky barrier parameters shows that the changes in electrical performance during aging are non-monotonic.These findings offer theoretical support for understanding the aging mechanisms and condition assessment of modern stable ZnO varistors.展开更多
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.展开更多
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.展开更多
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.展开更多
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(Pf)by considering seismic stability with spatial heterogeneity.The Pfwas 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 finite element approach is used for the first time to simulate and examine the free oscillation and transient response of a visco-elastic multi-directional functionally graded porous(MFGP)skew-nanoplate,taking into account surface effects using nonlocal strain gradient hypothesis.The mechanical characteristics of the materials vary in all three directions of length,width,and thickness of the plate according to the exponential law.Additionally,it has viscoelastic behavior according to the Kelvin-Voigt model.The novelty of this paper lies in the incorporation of the spatial variability of nonlocal and lengthscale factors as additional mechanical characteristics of the material.The overall equation of motion for the plate is derived by including the classical plate hypothesis and Hamilton’s principle.A quadrilateral plate element with four nodes and six degrees of freedom is created using a non-conforming C2-level Hermitian function.This function offers precise results and rapid convergence for various forms and boundary conditions(BCs)that low-order elements cannot accomplish.The Newmark-beta direct integration technique is used to calculate the transient responses of the visco-elastic MFGP skew-nanoplate under various BCs.Furthermore,a thorough assessment of the impacts of several factors such as residual surface stress,grading indices,elastic foundation stiffness,skew angle,other geometrical parameters,and BCs on the transient responses of the viscoelastic MFGP skew-nanoplate has been uncovered.
摘要This study develops a surrogate super-resolution(SR)framework that accelerates finite element method(FEM)-based computational fluid dynamics(CFD)using deep learning.High-resolution(HR)FEM-based CFDremains computationally prohibitive for time-sensitive applications,including patient-specific aneurysm hemodynamics where rapid turnaround is valuable.The proposed pipeline learns to reconstruct HR velocity-magnitude fields fromlow-resolution(LR)FEM solutions generated under the same governing equations and boundary conditions.It consistsof three modules:(i)offline pre-training of a residual network on representative vascular geometries;(ii)lightweightfine-tuning to adapt the pretrained model to geometric variability,including patient-specific aneurysm morphologies;and(iii)an unstructured-to-structured sampling strategy with region-of-interest upsampling that concentrates resolution in flow-critical zones(e.g.,the aneurysm sac)rather than the full domain.This targeted reconstruction substantiallyreduces inference and post-processing cost while preserving key HR flow features.Experiments on cerebral aneurysmmodels show that HR velocity-magnitude fields can be recovered with accuracy comparable to direct HR simulationsat less than 1%of the direct HR simulation cost per analysis(LR simulation and SR inference),while adaptation to newgeometries requires only lightweight fine-tuning with limited target-specific HR data.While clinical endpoints andadditional variables(e.g.,pressure or wall-based metrics)are left for future work,the results indicate that the proposedsurrogate SR approach can streamline FEM-based CFD workflows toward near real-time hemodynamic analysis acrossmorphologically similar vascular models.
基金funded by the National Natural Science Foundation of China(Grant No.12227801).
摘要This paper proposes a state-of-the-art three-dimensional Voronoi cell finite element method(3D VCFEM)aimed at investigating the mechanical properties of particle-reinforced composites(PRCs)in space under different microstructural properties.Firstly,the modified residual energy generalized function of 3D VCFEM was proposed by applying the hybrid stress element method,and the element format of the 3D Voronoi element was constructed.On this basis,the interaction between the matrix and the inclusions was considered,and the higher-order stress function including the interaction stress term was constructed.Secondly,to solve the difficulty of integrating easily due to the complexity and irregularity of the integration region in space,Delaunay tetrahedra were introduced within the 3D Voronoi element for mesh refinement.It simplified the integration process.Finally,to verify the accuracy and efficiency of the 3D VCFEM model,comparative models of 3D VCFENM and FEM were established for analysis and discussion.The stress field and strain field were compared and analyzed for the first time.An example was also given for the presence of a large number of randomly distributed inclusion particles.The results showed that under the same accuracy,3D VCFEM had the advantages of convenient mesh delineation and high computational efficiency compared with FEM,which provided a new way of thinking to analyze the actual PCRs.
基金funded by the National Natural Science Foundation of China(Grant No.12227801).
摘要In this study,a collaborative optimization framework combining artificial fish swarm algorithm(AFSA)and Voronoi cell finite element method(VCFEM)is proposed to solve the problem of the influence of microstructure distribution on the macroscopic mechanical properties of particle-reinforced composites.Compared with the traditional displacement-based finite element method,VCFEM significantly improves the computational efficiency of multi-inclusion problems by means of generalized stress function and element local adaptation technology,and accurately captures the stress field discontinuity and interface stress concentration phenomenon.At the same time,AFSA abandons the dependence of traditional optimization methods on simplified models,and directly searches for the global optimal solution through discrete topological variables(particle positions).Its natural selection mechanism and swarm intelligence characteristics effectively overcome the local optimal trap.The numerical simulation results show that the proposed method exhibits high accuracy in both single-and multi-inclusion models(the error with the commercial software MARC is less than 5%),and in the complex model with 100 inclusions,the maximum Mises stress is reduced by 32.6%after optimization.The synergistic effect of the two breaks through the trade-off between efficiency and accuracy in traditional optimization,and provides the ability of both computational efficiency and global optimization for multi-scale composite material design.
摘要Slope stability issues in geotechnical engineering significantly impact both structural safety and economic efficiency, making stability analysis a critical component of design and construction processes. The finite element method is employed to develop numerical models of geotechnical slopes, incorporating appropriate boundary conditions and material parameters to conduct in-depth investigations into slope behavior under various loading conditions. During model development, comprehensive consideration is given to the heterogeneity of each rock layer and detailed internal structural characteristics, with meticulous examination of stress distribution patterns and deformation accumulation mechanisms under different stress environments. Research findings demonstrate that the finite element method effectively captures internal stress distributions and identifies precise locations of potential failure surfaces, providing robust technical support for early detection of slope instability risks. The study also elucidates the underlying mechanisms governing slope stability under multiple influencing factors, clarifying the operational mechanisms and interrelationships of key control parameters. Based on these findings, practical engineering mitigation measures are proposed, offering solid theoretical foundations and operational guidance for safe slope design and risk management in geotechnical engineering. These advancements enhance overall project safety, improve comprehensive disaster prevention and response systems, and contribute significantly to industry excellence.
基金supported by the National Natural Science Foundation of China(Grant Nos.12202321,12432009,12172262,and 12202322)the National Key R&D Program of China(Grant No.2022YFE0113100)。
摘要To address challenges in architectural extensibility and cross-module collaboration of CAE software,this study proposes OPFEM(open-source Python-based finite element modeling)—an open-source framework featuring a unified four-layer architecture.The geometric modeling framework achieves plug-in support for geometric kernels through an interface abstraction layer and adapter patterns,decoupling kernel-specific implementations while enabling state machine-driven interaction design and parametric sketching.The pre-processing modules establish multi-level associations among materials,sections,and geometric entities using Composite and Factory patterns,while implementing Observer pattern to ensure geometric-mesh consistency and employing finite-state machines to optimize boundary workflows.The computational modules implement a modular finite element library that decouples topology from element attributes,along with a surface boundary element technique for load conversion and task-scheduling management,validated through a cantilever beam large-deformation case.The post-processing module facilitates standardized data storage and dynamic field visualization through architecture-level standardized interface definitions and hierarchical component design.Collectively,OPFEM achieves full-process integration from parametric modeling to nonlinear solving and visualization,enhancing configuration efficiency and providing an extensible,pattern-driven solution for complex CAE challenges.
基金financially supported by the Innovation Fund for Young Scholars of State Key Laboratory of Coastal and Offshore Engineering(Dalian University of Technology)(Grant No.LY2304)the National Natural Science Foundation of China(Grant No.52371277)+2 种基金the Guangdong Basic and Applied Basic Research Foundation(Grant No.2023A1515010890)the State Key Laboratory of Hydraulic Engineering Intelligent Construction and Operation(Tianjin University)(Grant No.HESS-2323)the Collaborative Research Program of the Research Institute for Applied Mechanics,Kyushu University(Grant No.25RE-1).
摘要A fully nonlinear potential flow(FNPF)solver has been developed using the Finite Element Method(FEM)to simulate time-domain interactions between free-surface waves and marine structures.The ALE framework is implemented alongside a segment spring analogy-based moving mesh strategy to accurately track evolving free surfaces and moving boundaries of floating bodies.The solver employs a preconditioned conjugate gradient method to efficiently resolve the resulting sparse,symmetric linear system at each time step.Temporal evolution is managed through a standard fourth-order Runge-Kutta scheme,while Chebyshev 5-point smoothing suppresses non-physical saw-tooth instabilities.The solver’s performance and reliability are verified through comprehensive benchmark tests,including freesurface sloshing,nonlinear wave propagation,and wave-structure interactions with submerged or floating bodies.Furthermore,the study explores a modified potential flow model incorporating a quadratic damping term to address viscous effects in gap/moonpool resonance problems.
摘要This study presents and verifies a hybrid methodology for reliable determination of parameters in structural rheological models(Zener,Burgers,and Maxwell)describing the viscoelastic behavior of polyurethane specimens manufactured using extrusion-based 3D printing.Through comprehensive testing,including cyclic compression at strain rates ranging from 0.12 to 120 mm/min(0%-15%strain)and creepelaxation experiments(10%-30%strain),the lumped parameters were independently determined using both analytical and numerical solutions of the models’differential equations,followed by cross-verification in additional experiments.Numerical solutions for creep and relaxation problems were obtained using finite element analysis,with the three-parameter Mooney-Rivlin model and Prony series employed to simulate elastic and viscous stress components,respectively.Energy dissipation per cycle was quantified during cyclic compression tests.The results demonstrate that all three models adequately describe material behavior within the 0%-15%strain range across various strain rates.Comparative analysis revealed the Burgers model’s superior performance in characterizing creep and stress relaxation at low strain levels.While Zener and Burgers model parameters from uniaxial compression showed limited applicability for energy dissipation calculations,the generalized Maxwell model effectively captured viscoelastic properties across different strain rates.Notably,parameters derived from creep tests provided a more universal assessment of dissipative properties due to optimization based on characteristic curve regions.Both parameter sets described polyurethane’s elastic-hysteretic behavior with approximately 20%error,proving significantly more accurate than the linear strain-time dependence hypothesis.Finite element analysis(FEA)complemented numerical modeling by demonstrating that while the generalized Maxwell model effectively describes initial rapid stress-strain changes,FEA provides superior characterization of steady-state processes.This computational approach yields more physically representative results compared to simplified analytical solutions,despite certain limitations in transient analysis.
基金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 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.
基金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.
基金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.
基金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 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.
基金supported by the National Natural Science Foundation of China Excellence Research Group Program for“Multiscale Problems in Nonlinear Mechanics”(Grant No.12588201)the National Key R&D Program of China(Grant No.2023YFA1008901)+1 种基金the National Nat-ural Science Foundation of China(Grant No.12172009)supported by“The Fundamental Research Funds for the Central Universities,Peking University”.
摘要This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solving a specific numerical problem under the scope of the linear finite element method(LFEM),so the method is termed computational method for analytical solutions with finite elements(CMAS-FE).The primary objective of the CMAS-FE is to construct analytical expressions for displacements and reaction forces at nodes,as well as for strains and stresses at elemental quadrature points,all of which are formulated as infinite series solutions of various orders of Poisson’s ratios.Like the conventional LFEM,the CMAS-FE forms global sparse linear equations,but the Young’s modulus and Poisson’s ratio remain variables(or symbols).By employing a direct inverse method to solve these symbolic linear systems,an analytical expression of the displacement field can be constructed.The CMAS-FE is validated via patch and bending tests,which demonstrate convergence with mesh and term refine-ment.Furthermore,the CMAS-FE is applied to obtain the bending stiffness of a beam structure and to estimate an approximate stress intensity factor for a straight crack within a square-shaped plate.
摘要In modern ZnO varistors,traditional aging mechanisms based on increased power consumption are no longer relevant due to reduced power consumption during DC aging.Prolonged exposure to both AC and DC voltages results in increased leakage current,decreased breakdown voltage,and lower nonlinearity,ultimately compromising their protective performance.To investigate the evolution in electrical properties during DC aging,this work developed a finite element model based on Voronoi networks and conducted accelerated aging tests on commercial varistors.Throughout the aging process,current-voltage characteristics and Schottky barrier parameters were measured and analyzed.The results indicate that when subjected to constant voltage,current flows through regions with larger grain sizes,forming discharge channels.As aging progresses,the current focus increases on these channels,leading to a decline in the varistor’s overall performance.Furthermore,analysis of the Schottky barrier parameters shows that the changes in electrical performance during aging are non-monotonic.These findings offer theoretical support for understanding the aging mechanisms and condition assessment of modern stable ZnO varistors.
基金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.
基金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.
基金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.
基金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(Pf)by considering seismic stability with spatial heterogeneity.The Pfwas 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.