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.展开更多
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.展开更多
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.展开更多
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.展开更多
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 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.展开更多
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.展开更多
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 Finite Element Method of Lines （FEMOL） is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to...The Finite Element Method of Lines （FEMOL） is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to a set of examples such as edge-crack plate, the central-crack plate, the plate with cracks emanating from a hole under tensile or under combination loads of tensile and bending. Their dimensionless stress distribution, the stress intensify factor （SIF） and crack opening displacement （COD） are obtained, and comparison with known solutions by other methods are reported. It is found that a good accuracy is achieved by FEMOL. The method is successfully modified to remarkably increase the accuracy and reduce convergence difficulties. So it is a very useful and new tool in studying fracture mechanics problems.展开更多
A 3D femoral model was built to obtain the three-dimensional temperature distribution of femur and its surrounding tissues and provide references for clinical applications. According to the relationship between gray-v...A 3D femoral model was built to obtain the three-dimensional temperature distribution of femur and its surrounding tissues and provide references for clinical applications. According to the relationship between gray-value and material properties,the model was assigned with various materials to make sure that it is more similar to the real femur in geometry and physical properties. 3D temperature distribution is obtained by using finite element analysis software ANSYS 11. 0 on the basis of heat conduction theory,Laplace equation,Pennes bio-heat transfer equation,thermo physical parameters of bone tissues,the boundary condition,and initial conditions. Taken the asymmetry of the 3D distribution of temperature into account,it is necessary to adopt the heating method with multiple heat sources. This method can ensure that the temperature fields match well with the tumor tissues and kill the tumor cells efficiently under the condition of protecting the normal tissues from damage. The analysis results supply important guidance for determining the needle position and the needle number and controlling the intensity of heating.展开更多
Precisely quantifying the strength of the proximal femur and accurately assessing hip fracture risk would enable those at high risk to be identified so that preventive interventions could be taken.Development of bette...Precisely quantifying the strength of the proximal femur and accurately assessing hip fracture risk would enable those at high risk to be identified so that preventive interventions could be taken.Development of better measures of femoral strength using the clinically展开更多
In order to analyze the load carrying capacity of prestressed concrete box girders, failure behaviors of in-situ deteriorated continuous prestressed concrete box girders under loading are experimentally observed and a...In order to analyze the load carrying capacity of prestressed concrete box girders, failure behaviors of in-situ deteriorated continuous prestressed concrete box girders under loading are experimentally observed and a finite failure analysis method for predicting behaviors of box girders is developed. A degenerated solid shell element is used to simulate box girders and material nonlinearity is considered. Since pre-stressed concrete box girders usually have a large number of curve prestressed tendons, a type of combined element is presented to simulate the prestressed tendons of box girders, and then the number of elements can be significantly reduced. The analytical results are compared with full-scale failure test results. The comparison shows that the presented method can be effectively applied to the failure analysis of in-situ continuous prestressed concrete box girders, and it also shows that the studied old bridge still has enough load carrying capacity.展开更多
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.展开更多
Objective： The biomechanical characters of the bone fracture of the man femoral hip joint under impact loads are explored. Methods ：A biosystem model of the man femoral hip joint by using the GE （ General Electric...Objective： The biomechanical characters of the bone fracture of the man femoral hip joint under impact loads are explored. Methods ：A biosystem model of the man femoral hip joint by using the GE （ General Electric） lightspeed multi-lay spiral CT is conducted. A 3D finite element model is established by employing the finite element software ANSYS. The FE analysis mainly concentrates on the effects of the impact directions arising from intense movements and the parenchyma on the femoral hip joint on the stress distributions of the proximal femur. Results：The parenchyma on the hip joint has relatively large relaxation effect on the impact loads. Conclusion：Effects of the angle δ of the impact load to the anterior direction and the angle γ of the impact load to the femur shaft on the bone fracture are given;δ has larger effect on the stress and strain distributions than the angle γ,which mainly represents the fracture of the upper femur including the femoral neck fracture when the posterolateral femur is impacted, consistent with the clinical resuits.展开更多
To analyze the stress state of steel orthotropic deck pavement and provide reference for the design of the overlay, the inner stress state and strain distribution of surfacing under the load of the deformation of the ...To analyze the stress state of steel orthotropic deck pavement and provide reference for the design of the overlay, the inner stress state and strain distribution of surfacing under the load of the deformation of the whole bridge structure and tyre load are analyzed by the finite element method of submodeling. Influence of surfacing modulus on the strain state of the overlay is analyzed for the purpose of the optimal design of the overlay structure. Analysis results show that the deformation of the whole bridge structure has no evident influence on the stress state of the overlay. The key factor of the overlay design is the transverse tensile strain in the overlay above the upper edge of web plate of rib. The stress state of the overlay is influenced evidently by the modulus of rigidity transform overlay. And the stress state of the overlay can be optimized and lowered by increasing the modulus and thickness of rigidity transform overlay, The fatigue test has been done to evaluate the fatigue performance and modulus of different deck pavement materials such as epoxy asphalt, SBS modified asphalt, rosphalt asphalt which can provide reference for deck pavement structure design.展开更多
The Taizhou Water Station Site is an ancient masonry structure ruin built in the Southern Song Dynasty. The main structure was severely damaged. In order to understand its current structural properties and provide a s...The Taizhou Water Station Site is an ancient masonry structure ruin built in the Southern Song Dynasty. The main structure was severely damaged. In order to understand its current structural properties and provide a scientific basis for protection design nonlinear finite element analysis and parameter analysis are carried out.The crack patterns deformations and stresses of the main structure under four load cases are analyzed by nonlinear finite element analysis and the effect of the backfill bulk density and modulus on the maximum principal tensile stress and maximum compressive stress are studied by parameter analysis.The results show that the most unfavorable condition for the foundation is the combination of weight＋backfill soil pressure＋additional load the most unfavorable load case to the main structure is weight ＋backfill soil pressure ＋water pressure＋additional heap load the maximum principal tensile stress of the main structure is very sensitive to the changes in the bulk density of the backfill soil.展开更多
In order to analyze the electrostatic field concerned with electrostatic proximity fuze problem using the available finite analysis software package, the technology to model the problem with a scale reduction object a...In order to analyze the electrostatic field concerned with electrostatic proximity fuze problem using the available finite analysis software package, the technology to model the problem with a scale reduction object and boundary was presented. The boundary is determined by the maximum distance the sensor can detect. The object model is obtained by multiplying the terms in Poisson's equation with a scale reduction factor and the real value can be reconstructed with the same reverse process after software calculation. Using the finite element analysis program, the simulation value is close to the theoretical value with a little error. The boundary determination and scale reduction method is suitable to modeling the irregular electrostatic field around air targets, such as airplane, missile and so on, which is based on commonly used personal computer (PC). The technology reduces the calculation and storage cost greatly.展开更多
Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh...Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh independently of f luid motion, and the container fixed noninertial coordinate system is employed to establish the governing equations so that the mesh is needed to be updated in this coordinate system only. This leads to a very simple mesh moving algorithm which makes it easy to trace the motion of the moving boundaries and the free su rface without producing undesirable distortion of the computational mesh. The fi nite element method and finite difference method are used spacewise and timewise , respectively. A numerical example involving either forced horizontal oscillati on or forced pitching oscillation of the fluid filled container is presented to illustrate the effectiveness and the robustness of the method. In additi on, this work can be extended for the fluid structure interaction problems.展开更多
In conjunction with ANSYS,we use the finite element method to analyze the bonding stresses of Si/ GaAs. We also apply a numerical model to investigate a contour map and the distribution of normal stress,shearing stres...In conjunction with ANSYS,we use the finite element method to analyze the bonding stresses of Si/ GaAs. We also apply a numerical model to investigate a contour map and the distribution of normal stress,shearing stress, and peeling stress, taking into full consideration the thermal expansion coefficient as a function of temperature. Novel bonding structures are proposed for reducing the effect of thermal stress as compared with conventional structures. Calculations show the validity of this new structure.展开更多
The difficulty in establishing the finite element model of head and cervical spine is interpreted in the study. A head-neck 3D model is constructed accurately and quickly by the technology of CT scan,the automatically...The difficulty in establishing the finite element model of head and cervical spine is interpreted in the study. A head-neck 3D model is constructed accurately and quickly by the technology of CT scan,the automatically modeling of Mimics software and the RE technology of Geomagic software. Then the finite element model of the head-neck which is close to the real one is set up by the preprocessor of the FEM soft ware ANSYS. After the transient finite element analysis is performed on the model,the historical response of the displacement of the head is obtained. The result is validated by the result of the existed experiment. The stress,as well as the deformation,of nodes in the head and the cervical spine at any time benefits a lot to the clinic study on the injure to the head and neck caused by the impacts. And all the analysis is done by limited computer available.展开更多
文摘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.
基金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 National Natural Science Foundation of China(Grant Nos.51890912,51979025 and 52011530189).
文摘This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.
基金supported by the 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.
文摘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 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 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.
基金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 Finite Element Method of Lines （FEMOL） is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to a set of examples such as edge-crack plate, the central-crack plate, the plate with cracks emanating from a hole under tensile or under combination loads of tensile and bending. Their dimensionless stress distribution, the stress intensify factor （SIF） and crack opening displacement （COD） are obtained, and comparison with known solutions by other methods are reported. It is found that a good accuracy is achieved by FEMOL. The method is successfully modified to remarkably increase the accuracy and reduce convergence difficulties. So it is a very useful and new tool in studying fracture mechanics problems.
基金Sponsored by the National Natural Science Foundation of China(Grant No.61272387)the Program for New Century Excellent Talents in University(Grant No.NCET-13-0756)the Distinguished Young Scientists Funds of Heilongjiang Province(Grant No.JC201302)
文摘A 3D femoral model was built to obtain the three-dimensional temperature distribution of femur and its surrounding tissues and provide references for clinical applications. According to the relationship between gray-value and material properties,the model was assigned with various materials to make sure that it is more similar to the real femur in geometry and physical properties. 3D temperature distribution is obtained by using finite element analysis software ANSYS 11. 0 on the basis of heat conduction theory,Laplace equation,Pennes bio-heat transfer equation,thermo physical parameters of bone tissues,the boundary condition,and initial conditions. Taken the asymmetry of the 3D distribution of temperature into account,it is necessary to adopt the heating method with multiple heat sources. This method can ensure that the temperature fields match well with the tumor tissues and kill the tumor cells efficiently under the condition of protecting the normal tissues from damage. The analysis results supply important guidance for determining the needle position and the needle number and controlling the intensity of heating.
基金supported by The HongKong Polytechnic University Research Grants(No.1-BB81)grants from National Natural Science Foundation of China,Nos.10872078 and 10832012
文摘Precisely quantifying the strength of the proximal femur and accurately assessing hip fracture risk would enable those at high risk to be identified so that preventive interventions could be taken.Development of better measures of femoral strength using the clinically
基金Post-Doctoral Innovative Projects of Shandong Province(No.200703072)the National Natural Science Foundation of China(No.50574053)
文摘In order to analyze the load carrying capacity of prestressed concrete box girders, failure behaviors of in-situ deteriorated continuous prestressed concrete box girders under loading are experimentally observed and a finite failure analysis method for predicting behaviors of box girders is developed. A degenerated solid shell element is used to simulate box girders and material nonlinearity is considered. Since pre-stressed concrete box girders usually have a large number of curve prestressed tendons, a type of combined element is presented to simulate the prestressed tendons of box girders, and then the number of elements can be significantly reduced. The analytical results are compared with full-scale failure test results. The comparison shows that the presented method can be effectively applied to the failure analysis of in-situ continuous prestressed concrete box girders, and it also shows that the studied old bridge still has enough load carrying capacity.
基金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.
文摘Objective： The biomechanical characters of the bone fracture of the man femoral hip joint under impact loads are explored. Methods ：A biosystem model of the man femoral hip joint by using the GE （ General Electric） lightspeed multi-lay spiral CT is conducted. A 3D finite element model is established by employing the finite element software ANSYS. The FE analysis mainly concentrates on the effects of the impact directions arising from intense movements and the parenchyma on the femoral hip joint on the stress distributions of the proximal femur. Results：The parenchyma on the hip joint has relatively large relaxation effect on the impact loads. Conclusion：Effects of the angle δ of the impact load to the anterior direction and the angle γ of the impact load to the femur shaft on the bone fracture are given;δ has larger effect on the stress and strain distributions than the angle γ,which mainly represents the fracture of the upper femur including the femoral neck fracture when the posterolateral femur is impacted, consistent with the clinical resuits.
文摘To analyze the stress state of steel orthotropic deck pavement and provide reference for the design of the overlay, the inner stress state and strain distribution of surfacing under the load of the deformation of the whole bridge structure and tyre load are analyzed by the finite element method of submodeling. Influence of surfacing modulus on the strain state of the overlay is analyzed for the purpose of the optimal design of the overlay structure. Analysis results show that the deformation of the whole bridge structure has no evident influence on the stress state of the overlay. The key factor of the overlay design is the transverse tensile strain in the overlay above the upper edge of web plate of rib. The stress state of the overlay is influenced evidently by the modulus of rigidity transform overlay. And the stress state of the overlay can be optimized and lowered by increasing the modulus and thickness of rigidity transform overlay, The fatigue test has been done to evaluate the fatigue performance and modulus of different deck pavement materials such as epoxy asphalt, SBS modified asphalt, rosphalt asphalt which can provide reference for deck pavement structure design.
基金The National Natural Science Foundation of China(No.51108238)
文摘The Taizhou Water Station Site is an ancient masonry structure ruin built in the Southern Song Dynasty. The main structure was severely damaged. In order to understand its current structural properties and provide a scientific basis for protection design nonlinear finite element analysis and parameter analysis are carried out.The crack patterns deformations and stresses of the main structure under four load cases are analyzed by nonlinear finite element analysis and the effect of the backfill bulk density and modulus on the maximum principal tensile stress and maximum compressive stress are studied by parameter analysis.The results show that the most unfavorable condition for the foundation is the combination of weight＋backfill soil pressure＋additional load the most unfavorable load case to the main structure is weight ＋backfill soil pressure ＋water pressure＋additional heap load the maximum principal tensile stress of the main structure is very sensitive to the changes in the bulk density of the backfill soil.
文摘In order to analyze the electrostatic field concerned with electrostatic proximity fuze problem using the available finite analysis software package, the technology to model the problem with a scale reduction object and boundary was presented. The boundary is determined by the maximum distance the sensor can detect. The object model is obtained by multiplying the terms in Poisson's equation with a scale reduction factor and the real value can be reconstructed with the same reverse process after software calculation. Using the finite element analysis program, the simulation value is close to the theoretical value with a little error. The boundary determination and scale reduction method is suitable to modeling the irregular electrostatic field around air targets, such as airplane, missile and so on, which is based on commonly used personal computer (PC). The technology reduces the calculation and storage cost greatly.
文摘Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh independently of f luid motion, and the container fixed noninertial coordinate system is employed to establish the governing equations so that the mesh is needed to be updated in this coordinate system only. This leads to a very simple mesh moving algorithm which makes it easy to trace the motion of the moving boundaries and the free su rface without producing undesirable distortion of the computational mesh. The fi nite element method and finite difference method are used spacewise and timewise , respectively. A numerical example involving either forced horizontal oscillati on or forced pitching oscillation of the fluid filled container is presented to illustrate the effectiveness and the robustness of the method. In additi on, this work can be extended for the fluid structure interaction problems.
文摘In conjunction with ANSYS,we use the finite element method to analyze the bonding stresses of Si/ GaAs. We also apply a numerical model to investigate a contour map and the distribution of normal stress,shearing stress, and peeling stress, taking into full consideration the thermal expansion coefficient as a function of temperature. Novel bonding structures are proposed for reducing the effect of thermal stress as compared with conventional structures. Calculations show the validity of this new structure.
文摘The difficulty in establishing the finite element model of head and cervical spine is interpreted in the study. A head-neck 3D model is constructed accurately and quickly by the technology of CT scan,the automatically modeling of Mimics software and the RE technology of Geomagic software. Then the finite element model of the head-neck which is close to the real one is set up by the preprocessor of the FEM soft ware ANSYS. After the transient finite element analysis is performed on the model,the historical response of the displacement of the head is obtained. The result is validated by the result of the existed experiment. The stress,as well as the deformation,of nodes in the head and the cervical spine at any time benefits a lot to the clinic study on the injure to the head and neck caused by the impacts. And all the analysis is done by limited computer available.
