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.展开更多
Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump mate...Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.展开更多
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.展开更多
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.展开更多
Due to their superior properties, the interest in nanostructures is increasing today in engineering. This study presents a new two-noded curved finite element for analyzing the in-plane static behaviors of curved nano...Due to their superior properties, the interest in nanostructures is increasing today in engineering. This study presents a new two-noded curved finite element for analyzing the in-plane static behaviors of curved nanobeams. Opposite to traditional curved finite elements developed by using approximate interpolation functions, the proposed curved finite element is developed by using exact analytical solutions. Although this approach was first introduced for analyzing the mechanical behaviors of macro-scale curved beams by adopting the local theory of elasticity, the exact analytical expressions used in this study were obtained from the solutions of governing equations that were expressed via the differential form of the nonlocal theory of elasticity. Therefore, the effects of shear strain and axial extension included in the analytical formulation are also inherited by the curved finite element developed here. The rigidity matrix and the consistent force vector are developed for a circular finite element. To demonstrate the applicability of the method, static analyses of various curved nanobeams subjected to different boundary conditions and loading scenarios are performed, and the obtained results are compared with the exact analytical ones. The presented study provides an accurate and low computational cost method for researchers to investigate the in-plane static behavior of curved nanobeams.展开更多
Purpose–The brake pipe system was an essential braking component of the railway freight trains,but the existing E-type sealing rings had problems such as insufficient low-temperature resistance,poor heat stability an...Purpose–The brake pipe system was an essential braking component of the railway freight trains,but the existing E-type sealing rings had problems such as insufficient low-temperature resistance,poor heat stability and short service life.To address these issues,low-phenyl silicone rubber was prepared and tested,and the finite element analysis and experimental studies on the sealing performance of its sealing rings were carried out.Design/methodology/approach–The low-temperature resistance and thermal stability of the prepared lowphenyl silicone rubber were studied using low-temperature tensile testing,differential scanning calorimetry,dynamic thermomechanical analysis and thermogravimetric analysis.The sealing performance of the lowphenyl silicone rubber sealing ring was studied by using finite element analysis software abaqus and experiments.Findings–The prepared low-phenyl silicone rubber sealing ring possessed excellent low-temperature resistance and thermal stability.According to the finite element analysis results,the finish of the flange sealing surface and groove outer edge should be ensured,and extrusion damage should be avoided.The sealing rings were more susceptible to damage in high compression ratio and/or low-temperature environments.When the sealing effect was ensured,a small compression ratio should be selected,and rubbers with hardness and elasticity less affected by temperature should be selected.The prepared low-phenyl silicone rubber sealing ring had zero leakage at both room temperature(RT)and�508C.Originality/value–The innovation of this study is that it provides valuable data and experience for the future development of the sealing rings used in the brake pipe flange joints of the railway freight cars in China.展开更多
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 work,the microstructure evolution and mechanical behavior of extruded SiC/ZA63 Mg matrix composites are investigated via combined experimental study and three-dimensionalfinite element modelling(3D FEM)based on...In this work,the microstructure evolution and mechanical behavior of extruded SiC/ZA63 Mg matrix composites are investigated via combined experimental study and three-dimensionalfinite element modelling(3D FEM)based on the actual 3D microstructure achieved by synchrotron tomography.The results show that the average grain size of composite increases from 0.57μm of 8μm-SiC/ZA63 to 8.73μm of 50μm-SiC/ZA63.The type of texture transforms from the typicalfiber texture in 8μm-SiC/ZA63 to intense basal texture in 50μm-SiC/ZA63 composite and the intensity of texture increases sharply with increase of SiC particle size.The dynamic recrystallization(DRX)mechanism is also changed with increasing SiC particle size.Experimental and simulation results verify that the strength and elongation both decrease with increase of SiC particle size.The 8μm-SiC/ZA63 composite possesses the optimal mechanical property with yield strength(YS)of 383 MPa,ultimate tensile strength(UTS)of 424 MPa and elongation of 6.3%.The outstanding mechanical property is attributed to the ultrafine grain size,high-density precipitates and dislocation,good loading transfer effect and the interface bonding between SiC and matrix,as well as the weakened basal texture.The simulation results reveal that the micro-cracks tend to initiate at the interface between SiC and matrix,and then propagate along the interface between particle and Mg matrix or at the high strain and stress regions,and further connect with other micro-cracks.The main fracture mechanism in 8μm-SiC/ZA63 composite is ductile damage of matrix and interfacial debonding.With the increase of particle size,interface strength and particle strength decrease,and interface debonding and particle rupture become the main fracture mechanism in the 30μm-and 50μm-SiC/ZA63 composites.展开更多
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.展开更多
Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis r...Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.展开更多
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 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.展开更多
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.展开更多
In generator design field,waveform total harmonic distortion(THD)and telephone harmonic factor(THF)are parameters commonly used to measure the impact of generator no-load voltage harmonics on the power communication q...In generator design field,waveform total harmonic distortion(THD)and telephone harmonic factor(THF)are parameters commonly used to measure the impact of generator no-load voltage harmonics on the power communication quality.Tubular hydrogenerators are considered the optimal generator for exploiting low-head,high-flow hydro resources,and they have seen increasingly widespread application in China's power systems recent years.However,owing to the compact and constrained internal space of such generators,their internal magnetic-field harmonics are pronounced.Therefore,accurate calculation of their THD and THF is crucial during the analysis and design stages to ensure the quality of power communication.Especially in the electromagnetic field finite element modeling analysis of such generators,the type and order of the finite element meshes may have a significant impact on the THD and THF calculation results,which warrants in-depth research.To address this,this study takes a real 34 MW large tubular hydrogenerator as an example,and establishes its electromagnetic field finite element model under no-load conditions.Two types of meshes,five mesh densities,and two mesh orders are analyzed to reveal the effect of electromagnetic field finite element mesh types and orders on the calculation results of THD and THF for such generators.展开更多
Coronal shear fractures of the femoral neck (CSFF) are the most challenging to treat among proximal femur fractures, directly affecting the life expectancy of patients with osteoporosis. However, an adequate osteosynt...Coronal shear fractures of the femoral neck (CSFF) are the most challenging to treat among proximal femur fractures, directly affecting the life expectancy of patients with osteoporosis. However, an adequate osteosynthesis method has not been elucidated yet. This study investigated the displacement direction of the femoral head fragment and its effect on the bone using finite element method. A finite element model for CSFF was developed from CT image data of a patient with osteoporosis using Mechanical Finder (ver. 11). Subsequently, finite element analyses were performed on six osteosynthesis models under maximum load applied during walking. The compressive stresses, tensile stresses, and compressive strains of each model were examined. The results suggested that the compressive and tensile stress distributions were concentrated on the anterior side of the femoral neck. Compressive strain distribution in the femoral head and neck was concentrated in four areas: at the tip of the blade or lag screw, the anteroinferior side of the blade or lag screw near the fracture site, and the upper right and lower left near the junction of the blade or lag screw and nail. Thus, the distribution of both these stresses revealed that the femoral head fragment was prone to anterior and inferior displacement. Distribution of compressive strains revealed the direction of the stress exerted by the osteosynthetic implant on the bone. The same results were observed in all osteosynthetic implants;thus, the findings could lay the foundation for developing methods for placing osteosynthetic implants less prone to displacement and the osteosynthetic implants themselves. In particular, the study provides insight into the optimal treatment of CSFF.展开更多
In this paper, the optimal control problem of parabolic integro-differential equations is solved by gradient recovery based two-grid finite element method. Piecewise linear functions are used to approximate state and ...In this paper, the optimal control problem of parabolic integro-differential equations is solved by gradient recovery based two-grid finite element method. Piecewise linear functions are used to approximate state and co-state variables, and piecewise constant function is used to approximate control variables. Generally, the optimal conditions for the problem are solved iteratively until the control variable reaches error tolerance. In order to calculate all the variables individually and parallelly, we introduce a gradient recovery based two-grid method. First, we solve the small scaled optimal control problem on coarse grids. Next, we use the gradient recovery technique to recover the gradients of state and co-state variables. Finally, using the recovered variables, we solve the large scaled optimal control problem for all variables independently. Moreover, we estimate priori error for the proposed scheme, and use an example to validate the theoretical results.展开更多
The effect of structure,elastic modulus and thickness of lower modulus layer in porous titanium implants on the stress distribution at the implant-bone interface was investigated.Three-dimensional finite element model...The effect of structure,elastic modulus and thickness of lower modulus layer in porous titanium implants on the stress distribution at the implant-bone interface was investigated.Three-dimensional finite element models of different titanium implants were constructed.The structures of the implants included the whole lower modulus style （No.1）,bio-mimetic style （No.2）,the whole lower modulus style in cancellous bone （No.3） and the whole dense style No.4.The stress distributions at bone-implant interface under static loading were analyzed using Ansys Workbench 10.0 software.The results indicated that the distribution of interface stress is strongly depended on the structure of the implants.The maximum stresses in cancellous bone and root region of implant No.2 are lower than those in the other three implants.A decrease in the modulus of the low modulus layer facilitates the interface stress transferring.Increasing the thickness of the low modulus layer can reduce the stress and induce a more uniform stress distribution at the interface.Among the four implants,biomimetic style implant No.2 is superior in transferring implant-bone interface stress to surrounding bones.展开更多
The grain statistics effect was investigated through asymmetric rolling of pure copper foil by a realistic polycrystalline aggregates model and crystal plasticity element finite model.A polycrystalline aggregate model...The grain statistics effect was investigated through asymmetric rolling of pure copper foil by a realistic polycrystalline aggregates model and crystal plasticity element finite model.A polycrystalline aggregate model was generated and a crystal plasticity-based finite element model was developed for each grain and the specimen as a whole.The crystal plasticity model itself is rate dependent and accounts for local dissipative hardening effects and the original orientation of each grain was generated based on the orientation distribution function（ODF）.The deformation behaviors,including inhomogeneous material flow,decrease of contact press and roll force with the increase of grain size for the constant size of specimens,were studied.It is revealed that when the specimens are composed of only a few grains across thickness,the grains with different sizes,shapes and orientations are unevenly distributed in the specimen and each grain plays a significant role in micro-scale plastic deformation and leads to inhomogeneous deformation and the scatter of experimental and simulation results.The slip system activity was examined and the predicted results are consistent with the surface layer model.The slip band is strictly influenced by the misorientation of neighbor grain with consideration of slip system activity.Furthermore,it is found that the decrease of roll force and the most active of slip system in surface grains are caused by the increase of free surface grain effect when the grain size is increased.The results of the physical experiment and simulation provide a basic understanding of micro-scaled plastic deformation behavior in asymmetric foil rolling.展开更多
文摘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.
基金the financial support provided by MHRD,Govt.of IndiaCoal India Limited for providing financial assistance for the research(Project No.CIL/R&D/01/73/2021)the partial financial support provided by the Ministry of Education,Government of India,under SPARC project(Project No.P1207)。
文摘Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.
文摘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.
文摘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.
基金supported by Scientific Research Projects Department of Istanbul Technical University.Project Number:MGA-2018-41546.Grant receiver:E.T.
文摘Due to their superior properties, the interest in nanostructures is increasing today in engineering. This study presents a new two-noded curved finite element for analyzing the in-plane static behaviors of curved nanobeams. Opposite to traditional curved finite elements developed by using approximate interpolation functions, the proposed curved finite element is developed by using exact analytical solutions. Although this approach was first introduced for analyzing the mechanical behaviors of macro-scale curved beams by adopting the local theory of elasticity, the exact analytical expressions used in this study were obtained from the solutions of governing equations that were expressed via the differential form of the nonlocal theory of elasticity. Therefore, the effects of shear strain and axial extension included in the analytical formulation are also inherited by the curved finite element developed here. The rigidity matrix and the consistent force vector are developed for a circular finite element. To demonstrate the applicability of the method, static analyses of various curved nanobeams subjected to different boundary conditions and loading scenarios are performed, and the obtained results are compared with the exact analytical ones. The presented study provides an accurate and low computational cost method for researchers to investigate the in-plane static behavior of curved nanobeams.
基金supported by the Science and Technology Research and Development Plan of the China State Railway Group Company Limited(No.Q2023J012).
文摘Purpose–The brake pipe system was an essential braking component of the railway freight trains,but the existing E-type sealing rings had problems such as insufficient low-temperature resistance,poor heat stability and short service life.To address these issues,low-phenyl silicone rubber was prepared and tested,and the finite element analysis and experimental studies on the sealing performance of its sealing rings were carried out.Design/methodology/approach–The low-temperature resistance and thermal stability of the prepared lowphenyl silicone rubber were studied using low-temperature tensile testing,differential scanning calorimetry,dynamic thermomechanical analysis and thermogravimetric analysis.The sealing performance of the lowphenyl silicone rubber sealing ring was studied by using finite element analysis software abaqus and experiments.Findings–The prepared low-phenyl silicone rubber sealing ring possessed excellent low-temperature resistance and thermal stability.According to the finite element analysis results,the finish of the flange sealing surface and groove outer edge should be ensured,and extrusion damage should be avoided.The sealing rings were more susceptible to damage in high compression ratio and/or low-temperature environments.When the sealing effect was ensured,a small compression ratio should be selected,and rubbers with hardness and elasticity less affected by temperature should be selected.The prepared low-phenyl silicone rubber sealing ring had zero leakage at both room temperature(RT)and�508C.Originality/value–The innovation of this study is that it provides valuable data and experience for the future development of the sealing rings used in the brake pipe flange joints of the railway freight cars in China.
基金supported by the National Natural Science Foundation of China(Grant Nos.42102346,42172301).
文摘Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant importance.The traditional finite element method(FEM)remains one of the primary approaches for addressing such issues.However,the application of FEM typically necessitates the use of a fine finite element mesh to accurately capture the heterogeneous properties of the materials and meet the required computational precision,which inevitably leads to a reduction in computational efficiency.To enhance the computational accuracy and efficiency of the FEM for heterogeneous multi-field coupling problems,this study presents the coupling magneto-electro-elastic multiscale finite element method(CM-MsFEM)for heterogeneous MEE structures.Unlike the conventional multiscale FEM(MsFEM),the proposed algorithm simultaneously constructs displacement,electric,and magnetic potential multiscale basis functions to address the heterogeneity of the corresponding parameters.The macroscale formulation of CM-MsFEM was derived,and the macroscale/microscale responses of the problems were obtained through up/downscaling calculations.Evaluation using numerical examples analyzing the transient behavior of heterogeneous MEE structures demonstrated that the proposed method outperforms traditional FEM in terms of both accuracy and computational efficiency,making it an appropriate choice for numerically modeling the dynamics of heterogeneous MEE structures.
基金supported by the National Natural Science Foundation of China[51974058,52371005,52022017,51927801]the Fundamental Research Funds for the Central Universities(DUT23YG104).
文摘In this work,the microstructure evolution and mechanical behavior of extruded SiC/ZA63 Mg matrix composites are investigated via combined experimental study and three-dimensionalfinite element modelling(3D FEM)based on the actual 3D microstructure achieved by synchrotron tomography.The results show that the average grain size of composite increases from 0.57μm of 8μm-SiC/ZA63 to 8.73μm of 50μm-SiC/ZA63.The type of texture transforms from the typicalfiber texture in 8μm-SiC/ZA63 to intense basal texture in 50μm-SiC/ZA63 composite and the intensity of texture increases sharply with increase of SiC particle size.The dynamic recrystallization(DRX)mechanism is also changed with increasing SiC particle size.Experimental and simulation results verify that the strength and elongation both decrease with increase of SiC particle size.The 8μm-SiC/ZA63 composite possesses the optimal mechanical property with yield strength(YS)of 383 MPa,ultimate tensile strength(UTS)of 424 MPa and elongation of 6.3%.The outstanding mechanical property is attributed to the ultrafine grain size,high-density precipitates and dislocation,good loading transfer effect and the interface bonding between SiC and matrix,as well as the weakened basal texture.The simulation results reveal that the micro-cracks tend to initiate at the interface between SiC and matrix,and then propagate along the interface between particle and Mg matrix or at the high strain and stress regions,and further connect with other micro-cracks.The main fracture mechanism in 8μm-SiC/ZA63 composite is ductile damage of matrix and interfacial debonding.With the increase of particle size,interface strength and particle strength decrease,and interface debonding and particle rupture become the main fracture mechanism in the 30μm-and 50μm-SiC/ZA63 composites.
基金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.
文摘Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.
基金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(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.
基金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.
基金sponsored by the National Natural Science Foundation,Youth Foundation of China,Grant/Award Number:51607146Sichuan Natural Sciences Fund,Grant/Award Number:2023NSFSC0295。
文摘In generator design field,waveform total harmonic distortion(THD)and telephone harmonic factor(THF)are parameters commonly used to measure the impact of generator no-load voltage harmonics on the power communication quality.Tubular hydrogenerators are considered the optimal generator for exploiting low-head,high-flow hydro resources,and they have seen increasingly widespread application in China's power systems recent years.However,owing to the compact and constrained internal space of such generators,their internal magnetic-field harmonics are pronounced.Therefore,accurate calculation of their THD and THF is crucial during the analysis and design stages to ensure the quality of power communication.Especially in the electromagnetic field finite element modeling analysis of such generators,the type and order of the finite element meshes may have a significant impact on the THD and THF calculation results,which warrants in-depth research.To address this,this study takes a real 34 MW large tubular hydrogenerator as an example,and establishes its electromagnetic field finite element model under no-load conditions.Two types of meshes,five mesh densities,and two mesh orders are analyzed to reveal the effect of electromagnetic field finite element mesh types and orders on the calculation results of THD and THF for such generators.
文摘Coronal shear fractures of the femoral neck (CSFF) are the most challenging to treat among proximal femur fractures, directly affecting the life expectancy of patients with osteoporosis. However, an adequate osteosynthesis method has not been elucidated yet. This study investigated the displacement direction of the femoral head fragment and its effect on the bone using finite element method. A finite element model for CSFF was developed from CT image data of a patient with osteoporosis using Mechanical Finder (ver. 11). Subsequently, finite element analyses were performed on six osteosynthesis models under maximum load applied during walking. The compressive stresses, tensile stresses, and compressive strains of each model were examined. The results suggested that the compressive and tensile stress distributions were concentrated on the anterior side of the femoral neck. Compressive strain distribution in the femoral head and neck was concentrated in four areas: at the tip of the blade or lag screw, the anteroinferior side of the blade or lag screw near the fracture site, and the upper right and lower left near the junction of the blade or lag screw and nail. Thus, the distribution of both these stresses revealed that the femoral head fragment was prone to anterior and inferior displacement. Distribution of compressive strains revealed the direction of the stress exerted by the osteosynthetic implant on the bone. The same results were observed in all osteosynthetic implants;thus, the findings could lay the foundation for developing methods for placing osteosynthetic implants less prone to displacement and the osteosynthetic implants themselves. In particular, the study provides insight into the optimal treatment of CSFF.
文摘In this paper, the optimal control problem of parabolic integro-differential equations is solved by gradient recovery based two-grid finite element method. Piecewise linear functions are used to approximate state and co-state variables, and piecewise constant function is used to approximate control variables. Generally, the optimal conditions for the problem are solved iteratively until the control variable reaches error tolerance. In order to calculate all the variables individually and parallelly, we introduce a gradient recovery based two-grid method. First, we solve the small scaled optimal control problem on coarse grids. Next, we use the gradient recovery technique to recover the gradients of state and co-state variables. Finally, using the recovered variables, we solve the large scaled optimal control problem for all variables independently. Moreover, we estimate priori error for the proposed scheme, and use an example to validate the theoretical results.
基金Project(30770576) supported by the National Natural Science Foundation of ChinaProject(2007AA03Z114) supported by Hi-tech Research and Development Program of ChinaProject supported by State Key Laboratory of Powder Metallurgy,China
文摘The effect of structure,elastic modulus and thickness of lower modulus layer in porous titanium implants on the stress distribution at the implant-bone interface was investigated.Three-dimensional finite element models of different titanium implants were constructed.The structures of the implants included the whole lower modulus style （No.1）,bio-mimetic style （No.2）,the whole lower modulus style in cancellous bone （No.3） and the whole dense style No.4.The stress distributions at bone-implant interface under static loading were analyzed using Ansys Workbench 10.0 software.The results indicated that the distribution of interface stress is strongly depended on the structure of the implants.The maximum stresses in cancellous bone and root region of implant No.2 are lower than those in the other three implants.A decrease in the modulus of the low modulus layer facilitates the interface stress transferring.Increasing the thickness of the low modulus layer can reduce the stress and induce a more uniform stress distribution at the interface.Among the four implants,biomimetic style implant No.2 is superior in transferring implant-bone interface stress to surrounding bones.
基金Project(51374069)supported by the National Natural Science Foundation of ChinaProject(U1460107)supported by the Joint Fund of the National Natural Science Foundation of China
文摘The grain statistics effect was investigated through asymmetric rolling of pure copper foil by a realistic polycrystalline aggregates model and crystal plasticity element finite model.A polycrystalline aggregate model was generated and a crystal plasticity-based finite element model was developed for each grain and the specimen as a whole.The crystal plasticity model itself is rate dependent and accounts for local dissipative hardening effects and the original orientation of each grain was generated based on the orientation distribution function（ODF）.The deformation behaviors,including inhomogeneous material flow,decrease of contact press and roll force with the increase of grain size for the constant size of specimens,were studied.It is revealed that when the specimens are composed of only a few grains across thickness,the grains with different sizes,shapes and orientations are unevenly distributed in the specimen and each grain plays a significant role in micro-scale plastic deformation and leads to inhomogeneous deformation and the scatter of experimental and simulation results.The slip system activity was examined and the predicted results are consistent with the surface layer model.The slip band is strictly influenced by the misorientation of neighbor grain with consideration of slip system activity.Furthermore,it is found that the decrease of roll force and the most active of slip system in surface grains are caused by the increase of free surface grain effect when the grain size is increased.The results of the physical experiment and simulation provide a basic understanding of micro-scaled plastic deformation behavior in asymmetric foil rolling.
