Objective: To compare the stress distribution in the periodontal ligament under different orthodontic forces during canine distalization using long-arm brackets, and to determine the optimal force value for this devic...Objective: To compare the stress distribution in the periodontal ligament under different orthodontic forces during canine distalization using long-arm brackets, and to determine the optimal force value for this device in orthodontic treatment. Methods: A finite element model was constructed after extracting the mandibular first premolar, and a long-arm bracket with a traction height of 6 mm was placed on the labial side of the mandibular canine. Three working conditions of 50 g, 100 g, and 150 g were simulated, and the magnitude and distribution of von Mises stress in the periodontal ligament were compared for each condition. Results: The maximum von Mises stress in the periodontal ligament was 0.013281 MPa in the 50 g condition, 0.02536 MPa in the 100 g condition, and 0.035549 MPa in the 150 g condition. As the orthodontic force increased, the stress distribution area in the periodontal ligament also expanded. Conclusion: A 100 g orthodontic force is the most suitable when using long-arm brackets, providing a relatively uniform stress distribution in the periodontal ligament and keeping the stress within a reasonable range.展开更多
Fractional-order time-delay differential equations can describe many complex physical phenomena with memory or delay effects, which are widely used in the fields of cell biology, control systems, signal processing, et...Fractional-order time-delay differential equations can describe many complex physical phenomena with memory or delay effects, which are widely used in the fields of cell biology, control systems, signal processing, etc. Therefore, it is of great significance to study fractional-order time-delay differential equations. In this paper, we discuss a finite volume element method for a class of fractional-order neutral time-delay differential equations. By introducing an intermediate variable, the fourth-order problem is transformed into a system of equations consisting of two second-order partial differential equations. The L1 formula is used to approximate the time fractional order derivative terms, and the finite volume element method is used in space. A fully discrete format of the equations is established, and we prove the existence, uniqueness, convergence and stability of the solution. Finally, the validity of the format is verified by numerical examples.展开更多
In this article, a finite volume element algorithm is presented and discussed for the numerical solutions of a time-fractional nonlinear fourth-order diffusion equation with time delay. By choosing the second-order sp...In this article, a finite volume element algorithm is presented and discussed for the numerical solutions of a time-fractional nonlinear fourth-order diffusion equation with time delay. By choosing the second-order spatial derivative of the original unknown as an additional variable, the fourth-order problem is transformed into a second-order system. Then the fully discrete finite volume element scheme is formulated by using L1approximation for temporal Caputo derivative and finite volume element method in spatial direction. The unique solvability and stable result of the proposed scheme are proved. A priori estimate of L2-norm with optimal order of convergence O(h2+τ2−α)where τand hare time step length and space mesh parameter, respectively, is obtained. The efficiency of the scheme is supported by some numerical experiments.展开更多
Total hip arthroplasty for adults with sequelae from childhood hip disorders poses significant challenges due to altered anatomy.The paper published by Oommen et al reviews the essential management strategies for thes...Total hip arthroplasty for adults with sequelae from childhood hip disorders poses significant challenges due to altered anatomy.The paper published by Oommen et al reviews the essential management strategies for these complex cases.This article explores the integration of finite element analysis(FEA)to enhance surgical precision and outcomes.FEA provides detailed biomechanical insights,aiding in preoperative planning,implant design,and surgical technique optimization.By simulating implant configurations and assessing bone quality,FEA helps in customizing implants and evaluating surgical techniques like subtrochanteric shortening osteotomy.Advanced imaging techniques,such as 3D printing,virtual reality,and augmented reality,further enhance total hip arthroplasty precision.Future research should focus on validating FEA models,developing patient-specific simulations,and promoting multidisciplinary collaboration.Integrating FEA and advanced technologies in total hip arthroplasty can improve functional outcomes,reduce complications,and enhance quality of life for patients with childhood hip disorder sequelae.展开更多
We present an assumed enhanced strain finite element framework for the simulation of tensile fracturing processes in transversely isotropic rocks.Fractures along the weak bedding planes and through the anisotropic roc...We present an assumed enhanced strain finite element framework for the simulation of tensile fracturing processes in transversely isotropic rocks.Fractures along the weak bedding planes and through the anisotropic rock matrix are treated with distinct enrichment,and a recently proposed dualmechanism tensile failure criterion for transversely isotropic rocks is adopted to determine crack initiation for the two failure modes.The cohesive crack model is adopted to characterize the response of embedded cracks.As for the numerical implementation of the proposed framework,both algorithms for the update of local history variables at Gauss points and of the global finite element system are derived.Four boundary-value problem simulations are carried out with the proposed framework,including uniaxial tension tests of Argillite,pre-notched square loaded in tension,three-point bending tests on Longmaxi shale,and simulations of tensile cracks induced by a strip load around a tunnel in transversely isotropic rocks.Simulation results reveal that the proposed framework can properly capture the tensile strength anisotropy and the anisotropic evolution of tensile cracks in transversely isotropic rocks.展开更多
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 goal of this paper is to represent a numerical study of magnetohydrodynamic mixed convection heat transfer in a lid-driven vertical wavy enclosure with a fin attached to the bottomwall.We use a finite elementmetho...The goal of this paper is to represent a numerical study of magnetohydrodynamic mixed convection heat transfer in a lid-driven vertical wavy enclosure with a fin attached to the bottomwall.We use a finite elementmethod based on Galerkin weighted residual(GWR)techniques to set up the appropriate governing equations for the present flow model.We have conducted a parametric investigation to examine the impact of Hartmann and Richardson numbers on the flow pattern and heat transmission features inside a wavy cavity.We graphically represent the numerical results,such as isotherms,streamlines,velocity profiles,local and mean Nusselt numbers,and average surface temperature.Comparisons between the results of this work and previously published work in a literature review have been produced to examine the reliability and consistency of the data.The different sizes of the fin surface significantly impact flow creation and temperature fields.Additionally,the long fin size is necessary to enhance the heat transfer rate on the right surface at large Richardson numbers and low Hartmann numbers.Fin surfaces can significantly increase the mixing of fluid inside the enclosure,which can mean reductions in reaction times and operating costs,along with increases in heat transfer and efficiency.展开更多
In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hy...In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.展开更多
The fatigue life and reliability of wrought carbon steel castings produced with an optimized mold design are predicted using a finite element method integrated with reliability calculations.The optimization of the mol...The fatigue life and reliability of wrought carbon steel castings produced with an optimized mold design are predicted using a finite element method integrated with reliability calculations.The optimization of the mold is carried out using MAGMASoft mainly based on porosity reduction as a response.After validating the initial mold design with experimental data,a spring flap,a common component of an automotive suspension system is designed and optimized followed by fatigue life prediction based on simulation using Fe-safe.By taking into consideration the variation in both stress and strength,the stress-strength model is used to predict the reliability of the component under fatigue loading.Under typical loading conditions of 70 kN,the analysis showed that 95%of the steel spring flaps achieve infinite life.However,under maximum loading conditions of 90 kN,reliability declined significantly,with only 65%of the spring flaps expected to withstand the stress without failure.The study also identified a safe load-induced stress of 95 MPa on the spring flap.The findings suggest that transitioning from forged to cast spring flaps is a promising option,particularly if further improvements in casting design reduce porosity to negligible levels,potentially achieving 100%reliability under typical loading conditions.This integrated approach of mold optimization coupled with reliability estimation under realistic service loading conditions offers significant potential for the casting industry to produce robust,cost-effective products.展开更多
The mechanical properties of an extruded Mg-10Gd sample, specifically designed for vascular stents, are crucial for predicting its behavior under service conditions. Achieving homogeneous stresses in the hoop directio...The mechanical properties of an extruded Mg-10Gd sample, specifically designed for vascular stents, are crucial for predicting its behavior under service conditions. Achieving homogeneous stresses in the hoop direction, essential for characterizing vascular stents, poses challenges in experimental testing based on standard specimens featuring a reduced cross section. This study utilizes an elasto-visco-plastic self-consistent polycrystal model(ΔEVPSC) with the predominant twinning reorientation(PTR) scheme as a numerical tool, offering an alternative to mechanical testing. For verification, various mechanical experiments, such as uniaxial tension, compression, notched-bar tension, three-point bending, and C-ring compression tests, were conducted. The resulting force vs. displacement curves and textures were then compared with those based on the ΔEVPSC model. The computational model's significance is highlighted by simulation results demonstrating that the differential hardening along with a weak strength differential effect observed in the Mg-10Gd sample is a result of the interplay between micromechanical deformation mechanisms and deformation-induced texture evolution. Furthermore, the study highlights that incorporating the axisymmetric texture from the as-received material incorporating the measured texture gradient significantly improves predictive accuracy on the strength in the hoop direction. Ultimately, the findings suggest that the ΔEVPSC model can effectively predict the mechanical behavior resulting from loading scenarios that are impossible to realize experimentally, emphasizing its valuable contribution as a digital twin.展开更多
A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines...A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines as trial functions to approximate the state and costate variables in two space dimensions.A Crank-Nicolson difference scheme is constructed for time discretization.The resulting numerical solutions belong to C2in space,and the order of the coefficient matrix is low.Moreover,the Bogner-Fox-Schmit element is considered for comparison.Two numerical experiments demonstrate the feasibility and effectiveness of the proposed method.展开更多
For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of ...For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.展开更多
Microstructures determine mechanical properties of steels,but in actual steel product process it is difficult to accurately control the microstructure to meet the requirements.General microstructure characterization m...Microstructures determine mechanical properties of steels,but in actual steel product process it is difficult to accurately control the microstructure to meet the requirements.General microstructure characterization methods are time consuming and results are not rep-resentative for overall quality level as only a fraction of steel sample was selected to be examined.In this paper,a macro and micro coupled 3D model was developed for nondestructively characterization of steel microstructures.For electromagnetic signals analysis,the relative permeability value computed by the micro cellular model can be used in the macro electromagnetic sensor model.The effects of different microstructure components on the relative permeability of duplex stainless steel(grain size,phase fraction,and phase distribu-tion)were discussed.The output inductance of an electromagnetic sensor was determined by relative permeability values and can be val-idated experimentally.The findings indicate that the inductance value of an electromagnetic sensor at low frequency can distinguish dif-ferent microstructures.This method can be applied to real-time on-line characterize steel microstructures in process of steel rolling.展开更多
A combined method of wave superposition and finite element is proposed to solve the radiation noise of targets in shallow sea.Taking the sound propagation of spherical sound source in shallow sea as an example,the rad...A combined method of wave superposition and finite element is proposed to solve the radiation noise of targets in shallow sea.Taking the sound propagation of spherical sound source in shallow sea as an example,the radiation sound field of the spherical sound source is equivalent to the linear superposition of the radiation sound field of several internal point sound sources,and then the radiated noise induced by spherical sound source can be predicted quickly.The accuracy and efficiency of the method are verified by comparing with the numerical results of finite element method,and the rapid prediction of underwater radiated noise of cylindrical shell is carried out based on the method.The results show that compared with the finite element method,the relative error of the calculation results under different simulation conditions does not exceed 0.1%,and the calculation time is about 1/10 of the finite element method,so this method can be used to solve the radiated noise of shallow underwater targets.展开更多
A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization...A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.展开更多
High-performance finite element research has always been a major focus of finite element method studies.This article introduces isogeometric analysis into the finite element method and proposes a new isogeometric fini...High-performance finite element research has always been a major focus of finite element method studies.This article introduces isogeometric analysis into the finite element method and proposes a new isogeometric finite element method.Firstly,the physical field is approximated by uniform B-spline interpolation,while geometry is represented by non-uniform rational B-spline interpolation.By introducing a transformation matrix,elements of types C^(0)and C^(1)are constructed in the isogeometric finite element method.Subsequently,the corresponding calculation formats for one-dimensional bars,beams,and two-dimensional linear elasticity in the isogeometric finite element method are derived through variational principles and parameter mapping.The proposed method combines element construction techniques of the finite element method with geometric construction techniques of isogeometric analysis,eliminating the need for mesh generation and maintaining flexibility in element construc-tion.Two elements with interpolation characteristics are constructed in the method so that boundary conditions and connections between elements can be processed like the finite element method.Finally,the test results of several examples show that:(1)Under the same degree and element node numbers,the constructed elements are almost consistent with the results obtained by traditional finite element method;(2)For bar problems with large local field variations and beam problems with variable cross-sections,high-degree and multi-nodes elements constructed can achieve high computational accuracy with fewer degrees of freedom than finite element method;(3)The computational efficiency of isogeometric finite element method is higher than finite element method under similar degrees of freedom,while as degrees of freedom increase,the computational efficiency between the two is similar.展开更多
The variation of the principal stress of formations with the working and geo-mechanical conditions can trigger wellbore instabilities and adversely affect the well completion.A finite element model,based on the theory...The variation of the principal stress of formations with the working and geo-mechanical conditions can trigger wellbore instabilities and adversely affect the well completion.A finite element model,based on the theory of poro-elasticity and the Mohr-Coulomb rock damage criterion,is used here to analyze such a risk.The changes in wellbore stability before and after reservoir acidification are simulated for different pressure differences.The results indicate that the risk of wellbore instability grows with an increase in the production-pressure difference regardless of whether acidification is completed or not;the same is true for the instability area.After acidizing,the changes in the main geomechanical parameters(i.e.,elastic modulus,Poisson’s ratio,and rock strength)cause the maximum wellbore instability coefficient to increase.展开更多
This study aimed to conduct finite element(FE)analysis matched with an in vitro experiment to analyze traumatic retinal detachments(TrRD)resulting from blunt trauma and provide stress and strain thresholds to predict ...This study aimed to conduct finite element(FE)analysis matched with an in vitro experiment to analyze traumatic retinal detachments(TrRD)resulting from blunt trauma and provide stress and strain thresholds to predict the occurrence of TrRD.The in vitro experiment was performed on forty-eight porcine eyes using a pendulum device.We examined dynamic mechanical responses at four energy levels.A FE model,based on experimental results and published data,was used to simulate TrRD.Fifty-one additional eyes underwent immediate pathological examination following blunt impact.A dynamic variation of velocities was observed post-impact,displaying an approximate cosine oscillation-attenuation profile.Energy absorption increased as the initial energy and differed significantly at four energy levels(p<0.001).FE simulation showed a peak strain of 0.462 in the anterior vitreous body and a peak stress of 1.408 MPa at the cornea at the high-energy level.During the energy transfer,the stress was initially observed in retinal region along the impact direction at the separation.TrRD were observed in injured eyes,where a few detachments were detected in control eyes.Correlations were performed between the proportion of pathological outcomes and FE results.In conclusion,this study suggests that stress contributes to the development of retinal detachment,providing an indicator to distinguish the occurrence of TrRD.展开更多
At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The who...At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The whole BFGSP skew-plates is placed on a variable visco-elastic foundation(VEF)in the hygro-thermal environment and subjected to the blast load.The BFGSP skew-plate thickness is permitted to vary non-linearly over both the length and width of the skew-plate,thereby faithfully representing the real behavior of the structure itself.The analysis is based on a four-node planar quadrilateral element with eight degrees of freedom per node,which is approximated using Lagrange Q_(4)shape function and C^(1)level non-conforming Hermite shape function based on refined higher-order shear deformation plate theory.The forced vibration parameters of the non-uniform thickness BFGSP skew-plate are fully determined using Hamilton's principle and the Newmark-βdirect integration technique.Accuracy of the calculation program is validated by comparing its numerical results with those from reputable sources.Furthermore,a thorough assessment is conducted to determine the impact of various parameters on the free and forced vibration responses of the non-uniform thickness BFGSP skew-plate.The findings of the paper may be used in the development of civil and military structures in situations that are prone to exceptional forces,such as explosions and impacts load.展开更多
文摘Objective: To compare the stress distribution in the periodontal ligament under different orthodontic forces during canine distalization using long-arm brackets, and to determine the optimal force value for this device in orthodontic treatment. Methods: A finite element model was constructed after extracting the mandibular first premolar, and a long-arm bracket with a traction height of 6 mm was placed on the labial side of the mandibular canine. Three working conditions of 50 g, 100 g, and 150 g were simulated, and the magnitude and distribution of von Mises stress in the periodontal ligament were compared for each condition. Results: The maximum von Mises stress in the periodontal ligament was 0.013281 MPa in the 50 g condition, 0.02536 MPa in the 100 g condition, and 0.035549 MPa in the 150 g condition. As the orthodontic force increased, the stress distribution area in the periodontal ligament also expanded. Conclusion: A 100 g orthodontic force is the most suitable when using long-arm brackets, providing a relatively uniform stress distribution in the periodontal ligament and keeping the stress within a reasonable range.
文摘Fractional-order time-delay differential equations can describe many complex physical phenomena with memory or delay effects, which are widely used in the fields of cell biology, control systems, signal processing, etc. Therefore, it is of great significance to study fractional-order time-delay differential equations. In this paper, we discuss a finite volume element method for a class of fractional-order neutral time-delay differential equations. By introducing an intermediate variable, the fourth-order problem is transformed into a system of equations consisting of two second-order partial differential equations. The L1 formula is used to approximate the time fractional order derivative terms, and the finite volume element method is used in space. A fully discrete format of the equations is established, and we prove the existence, uniqueness, convergence and stability of the solution. Finally, the validity of the format is verified by numerical examples.
文摘In this article, a finite volume element algorithm is presented and discussed for the numerical solutions of a time-fractional nonlinear fourth-order diffusion equation with time delay. By choosing the second-order spatial derivative of the original unknown as an additional variable, the fourth-order problem is transformed into a second-order system. Then the fully discrete finite volume element scheme is formulated by using L1approximation for temporal Caputo derivative and finite volume element method in spatial direction. The unique solvability and stable result of the proposed scheme are proved. A priori estimate of L2-norm with optimal order of convergence O(h2+τ2−α)where τand hare time step length and space mesh parameter, respectively, is obtained. The efficiency of the scheme is supported by some numerical experiments.
文摘Total hip arthroplasty for adults with sequelae from childhood hip disorders poses significant challenges due to altered anatomy.The paper published by Oommen et al reviews the essential management strategies for these complex cases.This article explores the integration of finite element analysis(FEA)to enhance surgical precision and outcomes.FEA provides detailed biomechanical insights,aiding in preoperative planning,implant design,and surgical technique optimization.By simulating implant configurations and assessing bone quality,FEA helps in customizing implants and evaluating surgical techniques like subtrochanteric shortening osteotomy.Advanced imaging techniques,such as 3D printing,virtual reality,and augmented reality,further enhance total hip arthroplasty precision.Future research should focus on validating FEA models,developing patient-specific simulations,and promoting multidisciplinary collaboration.Integrating FEA and advanced technologies in total hip arthroplasty can improve functional outcomes,reduce complications,and enhance quality of life for patients with childhood hip disorder sequelae.
基金supported by the National Natural Science Foundation of China(Grant Nos.52038005 and 52201326)the fellowship of China Postdoctoral Science Foundation(Grant No.2022M721883)Tsinghua University Initiative Scientific Research Program.
文摘We present an assumed enhanced strain finite element framework for the simulation of tensile fracturing processes in transversely isotropic rocks.Fractures along the weak bedding planes and through the anisotropic rock matrix are treated with distinct enrichment,and a recently proposed dualmechanism tensile failure criterion for transversely isotropic rocks is adopted to determine crack initiation for the two failure modes.The cohesive crack model is adopted to characterize the response of embedded cracks.As for the numerical implementation of the proposed framework,both algorithms for the update of local history variables at Gauss points and of the global finite element system are derived.Four boundary-value problem simulations are carried out with the proposed framework,including uniaxial tension tests of Argillite,pre-notched square loaded in tension,three-point bending tests on Longmaxi shale,and simulations of tensile cracks induced by a strip load around a tunnel in transversely isotropic rocks.Simulation results reveal that the proposed framework can properly capture the tensile strength anisotropy and the anisotropic evolution of tensile cracks in transversely isotropic rocks.
基金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.
基金the Deanship of Scientific Research at Umm Al-Qura University for supporting this work through Grant Code:22UQU4240002DSR19.
文摘The goal of this paper is to represent a numerical study of magnetohydrodynamic mixed convection heat transfer in a lid-driven vertical wavy enclosure with a fin attached to the bottomwall.We use a finite elementmethod based on Galerkin weighted residual(GWR)techniques to set up the appropriate governing equations for the present flow model.We have conducted a parametric investigation to examine the impact of Hartmann and Richardson numbers on the flow pattern and heat transmission features inside a wavy cavity.We graphically represent the numerical results,such as isotherms,streamlines,velocity profiles,local and mean Nusselt numbers,and average surface temperature.Comparisons between the results of this work and previously published work in a literature review have been produced to examine the reliability and consistency of the data.The different sizes of the fin surface significantly impact flow creation and temperature fields.Additionally,the long fin size is necessary to enhance the heat transfer rate on the right surface at large Richardson numbers and low Hartmann numbers.Fin surfaces can significantly increase the mixing of fluid inside the enclosure,which can mean reductions in reaction times and operating costs,along with increases in heat transfer and efficiency.
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
文摘In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.
基金funded by Interdisciplinary Research Center for Intelligent Manufacturing and Robotics at King Fahd University of Petroleum and Minerals (KFUPM),Dhahran,through Project#INMR2107.
文摘The fatigue life and reliability of wrought carbon steel castings produced with an optimized mold design are predicted using a finite element method integrated with reliability calculations.The optimization of the mold is carried out using MAGMASoft mainly based on porosity reduction as a response.After validating the initial mold design with experimental data,a spring flap,a common component of an automotive suspension system is designed and optimized followed by fatigue life prediction based on simulation using Fe-safe.By taking into consideration the variation in both stress and strength,the stress-strength model is used to predict the reliability of the component under fatigue loading.Under typical loading conditions of 70 kN,the analysis showed that 95%of the steel spring flaps achieve infinite life.However,under maximum loading conditions of 90 kN,reliability declined significantly,with only 65%of the spring flaps expected to withstand the stress without failure.The study also identified a safe load-induced stress of 95 MPa on the spring flap.The findings suggest that transitioning from forged to cast spring flaps is a promising option,particularly if further improvements in casting design reduce porosity to negligible levels,potentially achieving 100%reliability under typical loading conditions.This integrated approach of mold optimization coupled with reliability estimation under realistic service loading conditions offers significant potential for the casting industry to produce robust,cost-effective products.
基金supports from the National Research Foundation of Korea funded by the Ministry of Education (No. 2018R1A6A1A03024509, NRF-2023R1A2C1005121)
文摘The mechanical properties of an extruded Mg-10Gd sample, specifically designed for vascular stents, are crucial for predicting its behavior under service conditions. Achieving homogeneous stresses in the hoop direction, essential for characterizing vascular stents, poses challenges in experimental testing based on standard specimens featuring a reduced cross section. This study utilizes an elasto-visco-plastic self-consistent polycrystal model(ΔEVPSC) with the predominant twinning reorientation(PTR) scheme as a numerical tool, offering an alternative to mechanical testing. For verification, various mechanical experiments, such as uniaxial tension, compression, notched-bar tension, three-point bending, and C-ring compression tests, were conducted. The resulting force vs. displacement curves and textures were then compared with those based on the ΔEVPSC model. The computational model's significance is highlighted by simulation results demonstrating that the differential hardening along with a weak strength differential effect observed in the Mg-10Gd sample is a result of the interplay between micromechanical deformation mechanisms and deformation-induced texture evolution. Furthermore, the study highlights that incorporating the axisymmetric texture from the as-received material incorporating the measured texture gradient significantly improves predictive accuracy on the strength in the hoop direction. Ultimately, the findings suggest that the ΔEVPSC model can effectively predict the mechanical behavior resulting from loading scenarios that are impossible to realize experimentally, emphasizing its valuable contribution as a digital twin.
基金supported by the National Natural Science Foundation of China(11871312,12131014)the Natural Science Foundation of Shandong Province,China(ZR2023MA086)。
文摘A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines as trial functions to approximate the state and costate variables in two space dimensions.A Crank-Nicolson difference scheme is constructed for time discretization.The resulting numerical solutions belong to C2in space,and the order of the coefficient matrix is low.Moreover,the Bogner-Fox-Schmit element is considered for comparison.Two numerical experiments demonstrate the feasibility and effectiveness of the proposed method.
基金supported by National Natural Science Foundation of China(11771257)the Shandong Provincial Natural Science Foundation of China(ZR2023YQ002,ZR2023MA007,ZR2021MA004)。
文摘For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.
基金supported by the National Natural Science Foundation of China(No.52204340)the Natural Science Foundation of Guangxi,China(No.2022GXNSFBA035621)The authors wish to thank the Advanced Manufacturing and Materials Centre from Warwick Manufacturing Group(WMG),University of Warwick for the provision of facilities and equipment.
文摘Microstructures determine mechanical properties of steels,but in actual steel product process it is difficult to accurately control the microstructure to meet the requirements.General microstructure characterization methods are time consuming and results are not rep-resentative for overall quality level as only a fraction of steel sample was selected to be examined.In this paper,a macro and micro coupled 3D model was developed for nondestructively characterization of steel microstructures.For electromagnetic signals analysis,the relative permeability value computed by the micro cellular model can be used in the macro electromagnetic sensor model.The effects of different microstructure components on the relative permeability of duplex stainless steel(grain size,phase fraction,and phase distribu-tion)were discussed.The output inductance of an electromagnetic sensor was determined by relative permeability values and can be val-idated experimentally.The findings indicate that the inductance value of an electromagnetic sensor at low frequency can distinguish dif-ferent microstructures.This method can be applied to real-time on-line characterize steel microstructures in process of steel rolling.
基金Foundation item:This study was financially supported by the National Natural Science Foundation of China(Grant No.52101351)。
文摘A combined method of wave superposition and finite element is proposed to solve the radiation noise of targets in shallow sea.Taking the sound propagation of spherical sound source in shallow sea as an example,the radiation sound field of the spherical sound source is equivalent to the linear superposition of the radiation sound field of several internal point sound sources,and then the radiated noise induced by spherical sound source can be predicted quickly.The accuracy and efficiency of the method are verified by comparing with the numerical results of finite element method,and the rapid prediction of underwater radiated noise of cylindrical shell is carried out based on the method.The results show that compared with the finite element method,the relative error of the calculation results under different simulation conditions does not exceed 0.1%,and the calculation time is about 1/10 of the finite element method,so this method can be used to solve the radiated noise of shallow underwater targets.
基金supported by a Major Research Project in Higher Education Institutions in Henan Province,with Project Number 23A560015.
文摘A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.
基金funded by the Zhejiang Province Science and Technology Plan Project under grant number 2023C01069the Hebei Provincial Program on Key Basic Research Project under grant number 23311808Dthe Wenzhou Major Science and Technology Innovation Project of China under grant number ZG2022004。
文摘High-performance finite element research has always been a major focus of finite element method studies.This article introduces isogeometric analysis into the finite element method and proposes a new isogeometric finite element method.Firstly,the physical field is approximated by uniform B-spline interpolation,while geometry is represented by non-uniform rational B-spline interpolation.By introducing a transformation matrix,elements of types C^(0)and C^(1)are constructed in the isogeometric finite element method.Subsequently,the corresponding calculation formats for one-dimensional bars,beams,and two-dimensional linear elasticity in the isogeometric finite element method are derived through variational principles and parameter mapping.The proposed method combines element construction techniques of the finite element method with geometric construction techniques of isogeometric analysis,eliminating the need for mesh generation and maintaining flexibility in element construc-tion.Two elements with interpolation characteristics are constructed in the method so that boundary conditions and connections between elements can be processed like the finite element method.Finally,the test results of several examples show that:(1)Under the same degree and element node numbers,the constructed elements are almost consistent with the results obtained by traditional finite element method;(2)For bar problems with large local field variations and beam problems with variable cross-sections,high-degree and multi-nodes elements constructed can achieve high computational accuracy with fewer degrees of freedom than finite element method;(3)The computational efficiency of isogeometric finite element method is higher than finite element method under similar degrees of freedom,while as degrees of freedom increase,the computational efficiency between the two is similar.
基金This work is financially sponsored by Tarim Oilfield“Study on Adaptability Evaluation and Parameter Optimization of Completion Technology in Bozi Block,Tarim Oilfield”(Item Number:201021113436).
文摘The variation of the principal stress of formations with the working and geo-mechanical conditions can trigger wellbore instabilities and adversely affect the well completion.A finite element model,based on the theory of poro-elasticity and the Mohr-Coulomb rock damage criterion,is used here to analyze such a risk.The changes in wellbore stability before and after reservoir acidification are simulated for different pressure differences.The results indicate that the risk of wellbore instability grows with an increase in the production-pressure difference regardless of whether acidification is completed or not;the same is true for the instability area.After acidizing,the changes in the main geomechanical parameters(i.e.,elastic modulus,Poisson’s ratio,and rock strength)cause the maximum wellbore instability coefficient to increase.
基金supported by the National Nature Science Foundation of China(Grant Nos.11972066,U20A20390,and 11827803)the support of Open Fund of State Key Laboratory of Virtual Reality Technology and Systems。
文摘This study aimed to conduct finite element(FE)analysis matched with an in vitro experiment to analyze traumatic retinal detachments(TrRD)resulting from blunt trauma and provide stress and strain thresholds to predict the occurrence of TrRD.The in vitro experiment was performed on forty-eight porcine eyes using a pendulum device.We examined dynamic mechanical responses at four energy levels.A FE model,based on experimental results and published data,was used to simulate TrRD.Fifty-one additional eyes underwent immediate pathological examination following blunt impact.A dynamic variation of velocities was observed post-impact,displaying an approximate cosine oscillation-attenuation profile.Energy absorption increased as the initial energy and differed significantly at four energy levels(p<0.001).FE simulation showed a peak strain of 0.462 in the anterior vitreous body and a peak stress of 1.408 MPa at the cornea at the high-energy level.During the energy transfer,the stress was initially observed in retinal region along the impact direction at the separation.TrRD were observed in injured eyes,where a few detachments were detected in control eyes.Correlations were performed between the proportion of pathological outcomes and FE results.In conclusion,this study suggests that stress contributes to the development of retinal detachment,providing an indicator to distinguish the occurrence of TrRD.
文摘At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The whole BFGSP skew-plates is placed on a variable visco-elastic foundation(VEF)in the hygro-thermal environment and subjected to the blast load.The BFGSP skew-plate thickness is permitted to vary non-linearly over both the length and width of the skew-plate,thereby faithfully representing the real behavior of the structure itself.The analysis is based on a four-node planar quadrilateral element with eight degrees of freedom per node,which is approximated using Lagrange Q_(4)shape function and C^(1)level non-conforming Hermite shape function based on refined higher-order shear deformation plate theory.The forced vibration parameters of the non-uniform thickness BFGSP skew-plate are fully determined using Hamilton's principle and the Newmark-βdirect integration technique.Accuracy of the calculation program is validated by comparing its numerical results with those from reputable sources.Furthermore,a thorough assessment is conducted to determine the impact of various parameters on the free and forced vibration responses of the non-uniform thickness BFGSP skew-plate.The findings of the paper may be used in the development of civil and military structures in situations that are prone to exceptional forces,such as explosions and impacts load.
