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.展开更多
The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling ...The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.展开更多
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.展开更多
The present paper first investigates the collapse behavior of a conventional pipe-framed greenhouse under snow loading based on a 3-D finite element analysis,in which both geometrical and material non-linearities are ...The present paper first investigates the collapse behavior of a conventional pipe-framed greenhouse under snow loading based on a 3-D finite element analysis,in which both geometrical and material non-linearities are considered.Three snow load distribution patterns related to the wind-driven snow particle movement are used in the analysis.It is found that snow load distribution affects the deformation and collapse behavior of the pipe-framed greenhouse significantly.The results obtained in this study are consistent with the actual damage observed.Next,discussion is made of the effects of reinforcements by adding members to the basic frame on the strength of the whole structure,in which seven kinds of reinforcement methods are examined.A buckling analysis is also carried out.The results indicate that the most effective reinforcement method depends on the snow load distribution pattern.展开更多
The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent ...The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.展开更多
The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the tw...The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the two identical and coaxial half stators. The calculation of the field with or without current in the windings (respectively with or without permanent magnet) is done using a mixed formulation with strong coupling. In addition, the local high saturation of the ferromagnetic material and the radial and axial components of the magnetic flux are taken into account. The results obtained make it possible to clearly observe, as a function of the intensity of the bus current or the remanent induction, the saturation zones, the lines, the orientations and the magnetic flux densities. 3D finite element modelling provide more accurate numerical data on the magnetic field through multiphysics analysis. This analysis considers the actual operating conditions and leads to the design of an optimized machine structure, with or without current in the windings and/or permanent magnet.展开更多
Percutaneous electrical nerve stimulation of an injured nerve can promote and accelerate peripheral nerve regeneration and improve function.When performing acupuncture and moxibustion,locating the injured nerve using ...Percutaneous electrical nerve stimulation of an injured nerve can promote and accelerate peripheral nerve regeneration and improve function.When performing acupuncture and moxibustion,locating the injured nerve using ultrasound before percutaneous nerve stimulation can help prevent further injury to an already injured nerve.However,stimulation parameters have not been standardized.In this study,we constructed a multi-layer human forearm model using finite element modeling.Taking current density and activated function as optimization indicators,the optimal percutaneous nerve stimulation parameters were established.The optimal parameters were parallel placement located 3 cm apart with the injury site at the midpoint between the needles.To validate the efficacy of this regimen,we performed a randomized controlled trial in 23 patients with median nerve transection who underwent neurorrhaphy.Patients who received conventional rehabilitation combined with percutaneous electrical nerve stimulation experienced greater improvement in sensory function,motor function,and grip strength than those who received conventional rehabilitation combined with transcutaneous electrical nerve stimulation.These findings suggest that the percutaneous electrical nerve stimulation regimen established in this study can improve global median nerve function in patients with median nerve transection.展开更多
Isothermal hot compression experiments were conducted on homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy to investigate hot deformation behavior at the temperature range of 673-773 K and the strain rate range of 0.001-1 s...Isothermal hot compression experiments were conducted on homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy to investigate hot deformation behavior at the temperature range of 673-773 K and the strain rate range of 0.001-1 s^(-1)by using a Gleeble-1500D thermo mechanical simulator.Metallographic characterization on samples deformed to true strain of 0.70 illustrates the occurrence of flow localization and/or microcrack at deformation conditions of 673 K/0.01 s^(-1),673 K/1 s^(-1)and 698 K/1 s^(-1),indicating that these three deformation conditions should be excluded during hot working of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy.Based on the measured true stress-strain data,the strain-compensated Arrhenius constitutive model was constructed and then incorporated into UHARD subroutine of ABAQUS software to study hot deformation process of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy.By comparison with measured force-displacement curves,the predicted results can describe well the rheological behavior of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy,verifying the validity of finite element simulation of hot compression process with this complicated constitutive model.Numerical results demonstrate that the distribution of values of material parameters(α,n,Q and ln A)within deformed sample is inhomogeneous.This issue is directly correlated to the uneven distribution of equivalent plastic strain due to the friction effect.Moreover,at a given temperature the increase of strain rate would result in the decrease of equivalent plastic strain within the central region of deformed sample,which hinders the occurrence of dynamic recrystallization(DRX).展开更多
The comprehensive tire building and shaping processes are investigated through the finite element method(FEM)in this article.The mechanical properties of the uncured rubber from different tire components are investiga...The comprehensive tire building and shaping processes are investigated through the finite element method(FEM)in this article.The mechanical properties of the uncured rubber from different tire components are investigated through cyclic loading-unloading experiments under different strain rates.Based on the experiments,an elastoviscoplastic constitutive model is adopted to describe themechanical behaviors of the uncured rubber.The distinct mechanical properties,including the stress level,hysteresis and residual strain,of the uncured rubber can all be well characterized.The whole tire building process(including component winding,rubber bladder inflation,component stitching and carcass band folding-back)and the shaping process are simulated using this constitutive model.The simulated green tire profile is in good agreement with the actual profile obtained through 3D scanning.The deformation and stress of the rubber components and the cord reinforcements during production can be obtained fromthe FE simulation,which is helpful for judging the rationality of the tire construction design.Finally,the influence of the parameter“drum width”is investigated,and the simulated result is found to be consistent with the experimental observations,which verifies the effectiveness of the simulation.The established simulation strategy provides some guiding significance for the improvement of tire design parameters and the elimination of tire production defects.展开更多
In this study,an iterative-based three-dimensional finite element lower bound in association with the second-order cone programming method is adopted to evaluate the limit load of a single pile embedded in cross-aniso...In this study,an iterative-based three-dimensional finite element lower bound in association with the second-order cone programming method is adopted to evaluate the limit load of a single pile embedded in cross-anisotropic soils under general loading condition.The lower bound solutions of the pile embedded in an anisotropic soil deposit can be found by formulating the element equilibrium,equilibrium of shear and normal stresses along discontinuities,boundary conditions,yield function,and optimizing the objective function through the second-order cone programming method in conjunction with an iterative-based update procedure.A general loading condition is considered to profile the expansion of the safe load in the vertical-horizontal-moment(V-H-M)space.The results of this study are compared and validated against three different cases including an isotropic lateral loading,anisotropic end bearing capacity,and a pile embedded in an isotropic soil deposit under general loading condition.A parametric study is conducted to evaluate the impact of different influencing factors.It was found that the effect of anisotropy on the variation of lateral limit load of a single pile is more pronounced than the corresponding vertical and bending moment limit loads,whereas the interface properties have more significant effects on the vertical and bending moment limit loads in comparison to the lateral limit load.展开更多
This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their...This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their expressions and asymptotical stability criteria.Second,for the semi-discrete and one-parameter fully-discrete finite element methods solving the above equations,we work out the sufficient conditions for assuring that the finite element solutions are asymptotically stable.Finally,with a typical example with numerical experiments,we illustrate the applicability of the obtained theoretical results.展开更多
In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avo...In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avoid mesh distortion.Additional mesh smoothing and boundary node smoothing techniques are incorporated to improve the mesh quality and solution accuracy.The field variables are mapped from the old to the new mesh using the closest point projection method to minimize the mapping error.The procedures of the proposed Abaqus-based dynamic PFEM(Abaqus-DPFEM)analysis and its implementation in Abaqus are detailed.The accuracy and robustness of the proposed approach are examined via four illustrative numerical examples.The numerical results show a satisfactory agreement with published results and further confirm the applicability of the Abaqus-DPFEM to solving dynamic large-deformation problems in geotechnical engineering.展开更多
The large storage requirement is a critical issue in cross-correlation imaging-condition based reverse time migration(RTM),because it requires the operation of the source and receiver wavefields at the same time.The b...The large storage requirement is a critical issue in cross-correlation imaging-condition based reverse time migration(RTM),because it requires the operation of the source and receiver wavefields at the same time.The boundary value method(BVM),based on the finite difference method(FDM),can be used to reconstruct the source wavefield in the reverse time propagation in the same way as the receiver wavefield,which can reduce the storage burden of the RTM data.Considering that the FDM cannot well handle models with discontinuous material properties and rough interfaces,we develop a source wavefield reconstruction strategy based on the finite element method(FEM),using proper orthogonal decomposition(POD)to enhance computational efficiency.In this method,we divide the whole time period into several segments,and construct the POD basis functions to get a reduced order model(ROM)for the source wavefield reconstruction in each segment.We show the corresponding quantitative analysis of the storage requirement of the POD-FEM.Numerical tests on the homogeneous model show the effectiveness of the proposed method,while the layered model and part of the Marmousi model tests indicate that the POD-FEM can keep an excellent balance between computational efficiency and memory usage compared with the full-stored method(FSM)and the BVM,and can be effectively applied in imaging.展开更多
As key components connecting offshore floating production platforms and subsea imports, offshore flexible pipes play significant roles in oil, natural gas, and water injection. It is found that torsional failure is on...As key components connecting offshore floating production platforms and subsea imports, offshore flexible pipes play significant roles in oil, natural gas, and water injection. It is found that torsional failure is one of the failure modes of flexible pipes during transportation and laying. In this paper, a theoretical model(TM) of a flexible pipe section mechanics is established, in which the carcass layer and the pressure armor layer are equivalent to the orthogonal anisotropic layers due to its complex cross-section structure. The calculation results of the TM are consistent with those of a finite element model(FEM), which can accurately describe the torsional response of the flexible pipe.Subsequently, the TM and FEM are used to discuss the influence of boundary conditions on the torsional response.The structure of the flexible pipe is stable when twisted counterclockwise. However, limiting the top axial displacement can improve the axial and radial instability of the tensile armor layer when twisted clockwise. Finally, it is recommended that the flexible pipe can be kept under top fixation during service or installation to avoid torsional failure.展开更多
We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is ...We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.展开更多
The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can ov...The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.展开更多
This paper deals with the numerical implementation of the exponential Drucker-Parger plasticitymodel in the commercial finite element software,ABAQUS,via user subroutine UMAT for adhesive joint simulations.The influen...This paper deals with the numerical implementation of the exponential Drucker-Parger plasticitymodel in the commercial finite element software,ABAQUS,via user subroutine UMAT for adhesive joint simulations.The influence of hydrostatic pressure on adhesive strength was investigated by a modified Arcan fixture designed particularly to induce a different state of hydrostatic pressure within an adhesive layer.The developed user subroutine UMAT,which utilizes an associated plastic flow during a plastic deformation,can provide a good agreement between the simulations and the experimental data.Better numerical stability at highly positive hydrostatic pressure loads for a very high order of exponential function can also be achieved compared to when a non-associated flow is used.展开更多
Rail vehicles generate huge longitudinal impact loads in collisions.If unreasonable matching exists between the compressive strength of the intermediate coupler and the structural strength of the car body,the risk of ...Rail vehicles generate huge longitudinal impact loads in collisions.If unreasonable matching exists between the compressive strength of the intermediate coupler and the structural strength of the car body,the risk of car body structure damage and train derailment will increase.Herein,a four-stage rigid-flexible coupling finite element model of the coupler is established considering the coupler buckling load.The influence of the coupler buckling load on the train longitudinal-vertical-hori-zontal buckling behavior was studied,and the mechanism of the train horizontal buckling instability in train collisions was revealed.Analysis results show that an intermediate coupler should be designed to ensure that the actual buckling load is less than the compressive load when the car body structure begins to deform plastically.The actual buckling load of the coupler and the asymmetry of the structural strength of the car body in the lateral direction are two important influencing factors for the lateral buckling of a train collision.If the strength of the two sides of the car body structure in the lateral direction is asymmetrical,the deformation on the weaker side will be larger,and the end of the car body will begin to deflect under the action of the coupler force,which in turn causes the train to undergo sawtooth buckling.展开更多
The family of Falk-Neilan P_(k)finite elements,combined with the Argyris P_(k+1)finite elements,solves the Reissner-Mindlin plate equation quasi-optimally and locking-free,on triangular meshes.The method is truly conf...The family of Falk-Neilan P_(k)finite elements,combined with the Argyris P_(k+1)finite elements,solves the Reissner-Mindlin plate equation quasi-optimally and locking-free,on triangular meshes.The method is truly conforming or consistent in the sense that no projection/reduction is introduced.Theoretical proof and numerical confirmation are presented.展开更多
To prevent dense chromium oxide covers from cracking during firing,such as internal cracking,surface cracking and fracture,the firing process of dense chromium oxide products was simulated,and the performance paramete...To prevent dense chromium oxide covers from cracking during firing,such as internal cracking,surface cracking and fracture,the firing process of dense chromium oxide products was simulated,and the performance parameters at several specific temperatures before and after densification of chromium oxide green bodies were tested.The temperature difference and the stress difference of chromium oxide green bodies were calculated by Finite Element Software at heating and cooling rates of 10,20,30 and 40℃·h^(-1),respectively.The results show that large temperature difference and thermal stress difference are the main causes of internal cracking,surface cracking and fracture of the brick.Suitable heating rates reduce or avoid internal cracking of the brick.When the heating rate is 30℃·h^(-1) below 1450℃and 10℃·h^(-1) above 1450℃,the temperature difference and the thermal stress difference in the brick are below 4℃and 4 MPa,respectively,and there is no internal crack in the brick after firing.The initial cooling stage is the key stage that causes cracking or fracture of bricks.When the cooling rate is 20℃·h^(-1) above 1450℃,the thermal stress difference in the brick is less than 4 MPa,which can significantly reduce the surface cracks and fracture of the brick.展开更多
基金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.
基金The Construction S&T Project of the Department of Transportation of Sichuan Province(Grant No.2023A02)the National Natural Science Foundation of China(No.52109135).
文摘The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.
基金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.
基金financially supported by the Steel Structure Research and Education Promotion Project of the Japan Iron and Steel Federation in FY2016.
文摘The present paper first investigates the collapse behavior of a conventional pipe-framed greenhouse under snow loading based on a 3-D finite element analysis,in which both geometrical and material non-linearities are considered.Three snow load distribution patterns related to the wind-driven snow particle movement are used in the analysis.It is found that snow load distribution affects the deformation and collapse behavior of the pipe-framed greenhouse significantly.The results obtained in this study are consistent with the actual damage observed.Next,discussion is made of the effects of reinforcements by adding members to the basic frame on the strength of the whole structure,in which seven kinds of reinforcement methods are examined.A buckling analysis is also carried out.The results indicate that the most effective reinforcement method depends on the snow load distribution pattern.
文摘The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.
文摘The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the two identical and coaxial half stators. The calculation of the field with or without current in the windings (respectively with or without permanent magnet) is done using a mixed formulation with strong coupling. In addition, the local high saturation of the ferromagnetic material and the radial and axial components of the magnetic flux are taken into account. The results obtained make it possible to clearly observe, as a function of the intensity of the bus current or the remanent induction, the saturation zones, the lines, the orientations and the magnetic flux densities. 3D finite element modelling provide more accurate numerical data on the magnetic field through multiphysics analysis. This analysis considers the actual operating conditions and leads to the design of an optimized machine structure, with or without current in the windings and/or permanent magnet.
基金supported by the National Natural Science Foundation of China,No.81801787(to XZS)China Postdoctoral Science Foundation,No.2018M640238(to XZS)the Natural Science Foundation of Tianjin,No.20JCQNJC01690(to XLC)。
文摘Percutaneous electrical nerve stimulation of an injured nerve can promote and accelerate peripheral nerve regeneration and improve function.When performing acupuncture and moxibustion,locating the injured nerve using ultrasound before percutaneous nerve stimulation can help prevent further injury to an already injured nerve.However,stimulation parameters have not been standardized.In this study,we constructed a multi-layer human forearm model using finite element modeling.Taking current density and activated function as optimization indicators,the optimal percutaneous nerve stimulation parameters were established.The optimal parameters were parallel placement located 3 cm apart with the injury site at the midpoint between the needles.To validate the efficacy of this regimen,we performed a randomized controlled trial in 23 patients with median nerve transection who underwent neurorrhaphy.Patients who received conventional rehabilitation combined with percutaneous electrical nerve stimulation experienced greater improvement in sensory function,motor function,and grip strength than those who received conventional rehabilitation combined with transcutaneous electrical nerve stimulation.These findings suggest that the percutaneous electrical nerve stimulation regimen established in this study can improve global median nerve function in patients with median nerve transection.
基金supported by the National Natural Science Foundation of China(Grant Nos.51805064,51701034)the Scientific and Technological Research Program of Chongqing Municipal Education Commission(Grant Nos.KJQN201801137,KJ1600922)+1 种基金the Basic and Advanced Research Project of Chongqing Science and Technology Commission(Grant Nos.cstc2017jcyj AX0062,cstc2018jcyj AX0035)the Chongqing University Key Laboratory of Micro/Nano Materials Engineering and Technology(Grant Nos.KFJJ2003)
文摘Isothermal hot compression experiments were conducted on homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy to investigate hot deformation behavior at the temperature range of 673-773 K and the strain rate range of 0.001-1 s^(-1)by using a Gleeble-1500D thermo mechanical simulator.Metallographic characterization on samples deformed to true strain of 0.70 illustrates the occurrence of flow localization and/or microcrack at deformation conditions of 673 K/0.01 s^(-1),673 K/1 s^(-1)and 698 K/1 s^(-1),indicating that these three deformation conditions should be excluded during hot working of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy.Based on the measured true stress-strain data,the strain-compensated Arrhenius constitutive model was constructed and then incorporated into UHARD subroutine of ABAQUS software to study hot deformation process of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy.By comparison with measured force-displacement curves,the predicted results can describe well the rheological behavior of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy,verifying the validity of finite element simulation of hot compression process with this complicated constitutive model.Numerical results demonstrate that the distribution of values of material parameters(α,n,Q and ln A)within deformed sample is inhomogeneous.This issue is directly correlated to the uneven distribution of equivalent plastic strain due to the friction effect.Moreover,at a given temperature the increase of strain rate would result in the decrease of equivalent plastic strain within the central region of deformed sample,which hinders the occurrence of dynamic recrystallization(DRX).
基金funded by the NationalNatural Science Foundation of China (Nos.11902229,11502181)the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant Nos.XDB22040502,XDC06030200).
文摘The comprehensive tire building and shaping processes are investigated through the finite element method(FEM)in this article.The mechanical properties of the uncured rubber from different tire components are investigated through cyclic loading-unloading experiments under different strain rates.Based on the experiments,an elastoviscoplastic constitutive model is adopted to describe themechanical behaviors of the uncured rubber.The distinct mechanical properties,including the stress level,hysteresis and residual strain,of the uncured rubber can all be well characterized.The whole tire building process(including component winding,rubber bladder inflation,component stitching and carcass band folding-back)and the shaping process are simulated using this constitutive model.The simulated green tire profile is in good agreement with the actual profile obtained through 3D scanning.The deformation and stress of the rubber components and the cord reinforcements during production can be obtained fromthe FE simulation,which is helpful for judging the rationality of the tire construction design.Finally,the influence of the parameter“drum width”is investigated,and the simulated result is found to be consistent with the experimental observations,which verifies the effectiveness of the simulation.The established simulation strategy provides some guiding significance for the improvement of tire design parameters and the elimination of tire production defects.
文摘In this study,an iterative-based three-dimensional finite element lower bound in association with the second-order cone programming method is adopted to evaluate the limit load of a single pile embedded in cross-anisotropic soils under general loading condition.The lower bound solutions of the pile embedded in an anisotropic soil deposit can be found by formulating the element equilibrium,equilibrium of shear and normal stresses along discontinuities,boundary conditions,yield function,and optimizing the objective function through the second-order cone programming method in conjunction with an iterative-based update procedure.A general loading condition is considered to profile the expansion of the safe load in the vertical-horizontal-moment(V-H-M)space.The results of this study are compared and validated against three different cases including an isotropic lateral loading,anisotropic end bearing capacity,and a pile embedded in an isotropic soil deposit under general loading condition.A parametric study is conducted to evaluate the impact of different influencing factors.It was found that the effect of anisotropy on the variation of lateral limit load of a single pile is more pronounced than the corresponding vertical and bending moment limit loads,whereas the interface properties have more significant effects on the vertical and bending moment limit loads in comparison to the lateral limit load.
文摘This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their expressions and asymptotical stability criteria.Second,for the semi-discrete and one-parameter fully-discrete finite element methods solving the above equations,we work out the sufficient conditions for assuring that the finite element solutions are asymptotically stable.Finally,with a typical example with numerical experiments,we illustrate the applicability of the obtained theoretical results.
基金the National Natural Science Foundation of China(Grant No.41807223)the Fundamental Research Funds for the Central Universities(Grant No.B210202096)the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDA 23090202).
文摘In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avoid mesh distortion.Additional mesh smoothing and boundary node smoothing techniques are incorporated to improve the mesh quality and solution accuracy.The field variables are mapped from the old to the new mesh using the closest point projection method to minimize the mapping error.The procedures of the proposed Abaqus-based dynamic PFEM(Abaqus-DPFEM)analysis and its implementation in Abaqus are detailed.The accuracy and robustness of the proposed approach are examined via four illustrative numerical examples.The numerical results show a satisfactory agreement with published results and further confirm the applicability of the Abaqus-DPFEM to solving dynamic large-deformation problems in geotechnical engineering.
基金This work was supported by Natural Science Basic Research Program of Shaanxi(Program No.2023-JC-YB-269)the National Natural Science Foundation of China(Grant No.41974122).
文摘The large storage requirement is a critical issue in cross-correlation imaging-condition based reverse time migration(RTM),because it requires the operation of the source and receiver wavefields at the same time.The boundary value method(BVM),based on the finite difference method(FDM),can be used to reconstruct the source wavefield in the reverse time propagation in the same way as the receiver wavefield,which can reduce the storage burden of the RTM data.Considering that the FDM cannot well handle models with discontinuous material properties and rough interfaces,we develop a source wavefield reconstruction strategy based on the finite element method(FEM),using proper orthogonal decomposition(POD)to enhance computational efficiency.In this method,we divide the whole time period into several segments,and construct the POD basis functions to get a reduced order model(ROM)for the source wavefield reconstruction in each segment.We show the corresponding quantitative analysis of the storage requirement of the POD-FEM.Numerical tests on the homogeneous model show the effectiveness of the proposed method,while the layered model and part of the Marmousi model tests indicate that the POD-FEM can keep an excellent balance between computational efficiency and memory usage compared with the full-stored method(FSM)and the BVM,and can be effectively applied in imaging.
基金financially supported by the Natural Science Starting Project of SWPU (Grant No. 2022QHZ002)Sichuan Natural Science Foundation Youth Fund Project (Grant No. 2023NSFC0918)。
文摘As key components connecting offshore floating production platforms and subsea imports, offshore flexible pipes play significant roles in oil, natural gas, and water injection. It is found that torsional failure is one of the failure modes of flexible pipes during transportation and laying. In this paper, a theoretical model(TM) of a flexible pipe section mechanics is established, in which the carcass layer and the pressure armor layer are equivalent to the orthogonal anisotropic layers due to its complex cross-section structure. The calculation results of the TM are consistent with those of a finite element model(FEM), which can accurately describe the torsional response of the flexible pipe.Subsequently, the TM and FEM are used to discuss the influence of boundary conditions on the torsional response.The structure of the flexible pipe is stable when twisted counterclockwise. However, limiting the top axial displacement can improve the axial and radial instability of the tensile armor layer when twisted clockwise. Finally, it is recommended that the flexible pipe can be kept under top fixation during service or installation to avoid torsional failure.
基金supported by the National Natural Science Foundation of China(Nos.12132001 and 52192632)。
文摘We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.
文摘The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.
基金funded by King Mongkut’s University of Technology North Bangkok.Contract No.KMUTNB-PHD-62-07.
文摘This paper deals with the numerical implementation of the exponential Drucker-Parger plasticitymodel in the commercial finite element software,ABAQUS,via user subroutine UMAT for adhesive joint simulations.The influence of hydrostatic pressure on adhesive strength was investigated by a modified Arcan fixture designed particularly to induce a different state of hydrostatic pressure within an adhesive layer.The developed user subroutine UMAT,which utilizes an associated plastic flow during a plastic deformation,can provide a good agreement between the simulations and the experimental data.Better numerical stability at highly positive hydrostatic pressure loads for a very high order of exponential function can also be achieved compared to when a non-associated flow is used.
基金This work was supported by the National Natural Science Foundation of China(No.52172409)Sichuan Outstanding Youth Fund(No.2022JDJQ0025).
文摘Rail vehicles generate huge longitudinal impact loads in collisions.If unreasonable matching exists between the compressive strength of the intermediate coupler and the structural strength of the car body,the risk of car body structure damage and train derailment will increase.Herein,a four-stage rigid-flexible coupling finite element model of the coupler is established considering the coupler buckling load.The influence of the coupler buckling load on the train longitudinal-vertical-hori-zontal buckling behavior was studied,and the mechanism of the train horizontal buckling instability in train collisions was revealed.Analysis results show that an intermediate coupler should be designed to ensure that the actual buckling load is less than the compressive load when the car body structure begins to deform plastically.The actual buckling load of the coupler and the asymmetry of the structural strength of the car body in the lateral direction are two important influencing factors for the lateral buckling of a train collision.If the strength of the two sides of the car body structure in the lateral direction is asymmetrical,the deformation on the weaker side will be larger,and the end of the car body will begin to deflect under the action of the coupler force,which in turn causes the train to undergo sawtooth buckling.
文摘The family of Falk-Neilan P_(k)finite elements,combined with the Argyris P_(k+1)finite elements,solves the Reissner-Mindlin plate equation quasi-optimally and locking-free,on triangular meshes.The method is truly conforming or consistent in the sense that no projection/reduction is introduced.Theoretical proof and numerical confirmation are presented.
文摘To prevent dense chromium oxide covers from cracking during firing,such as internal cracking,surface cracking and fracture,the firing process of dense chromium oxide products was simulated,and the performance parameters at several specific temperatures before and after densification of chromium oxide green bodies were tested.The temperature difference and the stress difference of chromium oxide green bodies were calculated by Finite Element Software at heating and cooling rates of 10,20,30 and 40℃·h^(-1),respectively.The results show that large temperature difference and thermal stress difference are the main causes of internal cracking,surface cracking and fracture of the brick.Suitable heating rates reduce or avoid internal cracking of the brick.When the heating rate is 30℃·h^(-1) below 1450℃and 10℃·h^(-1) above 1450℃,the temperature difference and the thermal stress difference in the brick are below 4℃and 4 MPa,respectively,and there is no internal crack in the brick after firing.The initial cooling stage is the key stage that causes cracking or fracture of bricks.When the cooling rate is 20℃·h^(-1) above 1450℃,the thermal stress difference in the brick is less than 4 MPa,which can significantly reduce the surface cracks and fracture of the brick.
