To compare finite element analysis(FEA)predictions and stereovision digital image correlation(StereoDIC)strain measurements at the same spatial positions throughout a region of interest,a field comparison procedure is...To compare finite element analysis(FEA)predictions and stereovision digital image correlation(StereoDIC)strain measurements at the same spatial positions throughout a region of interest,a field comparison procedure is developed.The procedure includes(a)conversion of the finite element data into a triangular mesh,(b)selection of a common coordinate system,(c)determination of the rigid body transformation to place both measurements and FEA data in the same system and(d)interpolation of the FEA nodal information to the same spatial locations as the StereoDIC measurements using barycentric coordinates.For an aluminum Al-6061 double edge notched tensile specimen,FEA results are obtained using both the von Mises isotropic yield criterion and Hill’s quadratic anisotropic yield criterion,with the unknown Hill model parameters determined using full-field specimen strain measurements for the nominally plane stress specimen.Using Hill’s quadratic anisotropic yield criterion,the point-by-point comparison of experimentally based full-field strains and stresses to finite element predictions are shown to be in excellent agreement,confirming the effectiveness of the field comparison process.展开更多
This paper describes the structure of the base,on which two working arms are installed simultaneously.To ensure structura safety,the fatigue failure analysis and statics analysis are finished using the finite element ...This paper describes the structure of the base,on which two working arms are installed simultaneously.To ensure structura safety,the fatigue failure analysis and statics analysis are finished using the finite element method.The calculation can make sure that the structure of the base meets the design standard,and the material can be reduced one grade.展开更多
Materials with a negative Poisson's ratio(PR)are called auxetics;they are characterized by expansion/contraction when tensioned/compressed.Given this counterintuitive behavior,they present very particular character...Materials with a negative Poisson's ratio(PR)are called auxetics;they are characterized by expansion/contraction when tensioned/compressed.Given this counterintuitive behavior,they present very particular characteristics and mechanical behavior.Geometrical models have been developed to justify and artificiall reproduce such materials' auxetic behavior.The focus of this study is the exploration of a reentrant model by analyzing the variation in the PR of reentrant structures as a function of geometrical and base material parameters.It is shown that,even in the presence of protruding ribs,there may not be auxetic behavior,and this depends on the geometry of each reentrant structure.Values determined for these parameters can be helpful as approximate reference data in the design and fabrication of auxetic lattices using reentrant geometries.展开更多
The alternating method based on the fundamental solutions of the infinite domain containing a crack,namely Muskhelishvili’s solutions,divides the complex structure with a crack into a simple model without crack which...The alternating method based on the fundamental solutions of the infinite domain containing a crack,namely Muskhelishvili’s solutions,divides the complex structure with a crack into a simple model without crack which can be solved by traditional numerical methods and an infinite domain with a crack which can be solved by Muskhelishvili’s solutions.However,this alternating method cannot be directly applied to the edge crack problems since partial crack surface of Muskhelishvili’s solutions is located outside the computational domain.In this paper,an improved alternating method,the spline fictitious boundary element alternating method(SFBEAM),based on infinite domain with the combination of spline fictitious boundary element method(SFBEM)and Muskhelishvili’s solutions is proposed to solve the edge crack problems.Since the SFBEM and Muskhelishvili’s solutions are obtained in the framework of infinite domain,no special treatment is needed for solving the problem of edge cracks.Different mixed boundary conditions edge crack problems with varies of computational parameters are given to certify the high precision,efficiency and applicability of the proposed method compared with other alternating methods and extend finite element method.展开更多
In recent years,finite element analysis is increasingly being proposed in slope stability problems as a competitive method to traditional limit equilibrium methods(LEMs)which are known for their inherent deficiencies....In recent years,finite element analysis is increasingly being proposed in slope stability problems as a competitive method to traditional limit equilibrium methods(LEMs)which are known for their inherent deficiencies.However,the application of finite element method(FEM)to slope stability as a strength reduction method(SRM)or as finite element limit analysis(FELA)is not always a success for the drawbacks that characterize both methods.To increase the performance of finite element analysis in this problem,a new approach is proposed in this paper.It consists in gradually expanding the mobilized stress Mohr’s circles until the soil failure occurs according to a prescribed non-convergence criterion.The present approach called stress deviator increasing method(SDIM)is considered rigorous for three main reasons.Firstly,it preserves the definition of the factor of safety(FOS)as the ratio of soil shear strength to the mobilized shear stress.Secondly,it maintains the progressive development of shear stress resulting from the increase in the principal stress deviator on the same plane,on which the shear strength takes place.Thirdly,by introducing the concept of equivalent stress loading,the resulting trial stresses are checked against the violation of the actual yield criterion formed with the real strength parameters rather than those reduced by a trial factor.The new numerical procedure was encoded in a Fortran computer code called S^(4)DINA and verified by several examples.Comparisons with other numerical methods such as the SRM,gravity increasing method(GIM)or even FELA by assessing both the FOS and contours of equivalent plastic strains showed promising results.展开更多
Each surface of roughness has different shape of asperity which is modeled with various shapes of analytical models. In this paper, the differences among various models of shape of asperity investigate using the Finit...Each surface of roughness has different shape of asperity which is modeled with various shapes of analytical models. In this paper, the differences among various models of shape of asperity investigate using the Finite Element Method (FEM) and various analytical models. The contact stresses in rough surfaces are calculated analytically using various asperity shape models. Finite element analysis is also carried out assuming three types of material properties namely, the linear, the elastic-perfect plastic and the elastic-nonlinear hardening. The analytical results are compared with the results obtained by the finite element method. The results illustrate for using a deterministic approach which the numerical models are suitable. In hertz model, the result of force is very big in interface of causing deformation plastic, while Model Zhao has almost same result with FEM nonlinear property model. It is observed that the results obtained from Zhao’s model are generally in a better agreement with the results obtained from various finite element models especially in elastic-plastic and plastic zones, hence it may be concluded that Zhao’s model can be used for analyzing the rough surfaces in contact mechanics.展开更多
The FEM analysis of stress and deformation of steel base of the special diamond saw blades are carried out.The small carbide round segments are welded on the side of blades in order to increase the side wear re sistan...The FEM analysis of stress and deformation of steel base of the special diamond saw blades are carried out.The small carbide round segments are welded on the side of blades in order to increase the side wear re sistance of blade.Comparing with conventional saw blade,the maximum stress val ues of reasonable special saw blade are reduced respectively about 17%20% 33%,and the maximum deformation values are reduced respectively about 26%22 %44.7%in thetangential(X),radial(Y)and axial(Z)direction.The stress conce ntration zone is decreased for the special structure diamond saw blade.The stru ctures of diamond saw blade with different number of hard material pellet are an alyzed and optimized.展开更多
Chemo-mechanical coupling behavior of materials is a transformation process between mechanical and chemical energy.In this paper,based on the coupled chemo-mechanical constitutive equations and governing equations dur...Chemo-mechanical coupling behavior of materials is a transformation process between mechanical and chemical energy.In this paper,based on the coupled chemo-mechanical constitutive equations and governing equations during isothermal process,the equivalent integral forms of chemo-mechanical coupling governing equations and corresponding finite element procedure are obtained by using Hamilton’s principle.An isoparametric plane element for chemo-mechanical coupling is associated into ABAQUS finite element package through user element subroutine UEL.The numerical examples exhibit that the ionic concentration variation can cause mechanical deformation and mechanical action can produce redistribution of ionic concentration for hydrogels.It is proved that the present developed chemo-mechanical coupling finite element procedure can be utilized to model the coupling behavior of hydrogels effectively.展开更多
In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split i...In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton’s method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.展开更多
In this paper, we present a theoretical analysis for linear finite element superconvergent gradient recovery on Par6 mesh, the dual of which is centroidal Voronoi tessellations with the lowest energy per unit volume a...In this paper, we present a theoretical analysis for linear finite element superconvergent gradient recovery on Par6 mesh, the dual of which is centroidal Voronoi tessellations with the lowest energy per unit volume and is the congruent cell predicted by the three-dimensional Gersho's conjecture. We show that the linear finite element solution uh and the linear interpolation uI have superclose gradient on Par6 meshes. Consequently, the gradient recovered from the finite element solution by using the superconvergence patch recovery method is superconvergent to Vu. A numerical example is presented to verify the theoretical result.展开更多
A stabilized and convergent finite element formulation for the generalized Stokes problem is proposed and a posteriori analysis is performed to produce an error indicator. On this basis adaptive numerical method for s...A stabilized and convergent finite element formulation for the generalized Stokes problem is proposed and a posteriori analysis is performed to produce an error indicator. On this basis adaptive numerical method for solying the problem is developed . Numerical calculations are performed to confirm the reliability and effectiveness of the method.展开更多
Lateral closing wedge osteotomy is a treatment for Kienbock’s disease. It is one of the most frequently used treatment options, which has been reported with relatively good long-term results. However, the results abo...Lateral closing wedge osteotomy is a treatment for Kienbock’s disease. It is one of the most frequently used treatment options, which has been reported with relatively good long-term results. However, the results about the treatment are still controversial in some literatures and some key mechanisms are still not clear. The objective of the cur-rent study was to study the biomedical mechanism of the treatment. A finite element model was developed based on the geometry of carpal bones. Various situations in-cluding inclination angle changes by cutting the radial with 0?, 5?, 10? and 15? osteot-omy angles were studied. The effectiveness of the treatment was also studied for the carpal structure with abnormal positions of the lunate bone. The results show that the effectiveness of the stress reduction with the angle depends on many situations such as the initial morphology of the carpal structure.展开更多
In this article, we study the explicit expressions of the constants in the error estimate of the nonconforming finite element method. We explicitly obtain the approximation error estimate and the consistency error est...In this article, we study the explicit expressions of the constants in the error estimate of the nonconforming finite element method. We explicitly obtain the approximation error estimate and the consistency error estimate for the Wilson's element without the regular assumption, respectively, which implies the final finite element error estimate. Such explicit a priori error estimates can be used as computable error bounds.展开更多
基金Financial support provided by Correlated Solutions Incorporated to perform StereoDIC experimentsthe Department of Mechanical Engineering at the University of South Carolina for simulation studies is deeply appreciated.
文摘To compare finite element analysis(FEA)predictions and stereovision digital image correlation(StereoDIC)strain measurements at the same spatial positions throughout a region of interest,a field comparison procedure is developed.The procedure includes(a)conversion of the finite element data into a triangular mesh,(b)selection of a common coordinate system,(c)determination of the rigid body transformation to place both measurements and FEA data in the same system and(d)interpolation of the FEA nodal information to the same spatial locations as the StereoDIC measurements using barycentric coordinates.For an aluminum Al-6061 double edge notched tensile specimen,FEA results are obtained using both the von Mises isotropic yield criterion and Hill’s quadratic anisotropic yield criterion,with the unknown Hill model parameters determined using full-field specimen strain measurements for the nominally plane stress specimen.Using Hill’s quadratic anisotropic yield criterion,the point-by-point comparison of experimentally based full-field strains and stresses to finite element predictions are shown to be in excellent agreement,confirming the effectiveness of the field comparison process.
文摘This paper describes the structure of the base,on which two working arms are installed simultaneously.To ensure structura safety,the fatigue failure analysis and statics analysis are finished using the finite element method.The calculation can make sure that the structure of the base meets the design standard,and the material can be reduced one grade.
文摘Materials with a negative Poisson's ratio(PR)are called auxetics;they are characterized by expansion/contraction when tensioned/compressed.Given this counterintuitive behavior,they present very particular characteristics and mechanical behavior.Geometrical models have been developed to justify and artificiall reproduce such materials' auxetic behavior.The focus of this study is the exploration of a reentrant model by analyzing the variation in the PR of reentrant structures as a function of geometrical and base material parameters.It is shown that,even in the presence of protruding ribs,there may not be auxetic behavior,and this depends on the geometry of each reentrant structure.Values determined for these parameters can be helpful as approximate reference data in the design and fabrication of auxetic lattices using reentrant geometries.
基金supported by the National Natural Science Foundation of China(51078150)the National Natural Science Foundation of China(11602087)+1 种基金the State Key Laboratory of Subtropical Building Science,South China University of Technology(2017ZB32)National Undergraduate Innovative and Entrepreneurial Training Program(201810561180).
文摘The alternating method based on the fundamental solutions of the infinite domain containing a crack,namely Muskhelishvili’s solutions,divides the complex structure with a crack into a simple model without crack which can be solved by traditional numerical methods and an infinite domain with a crack which can be solved by Muskhelishvili’s solutions.However,this alternating method cannot be directly applied to the edge crack problems since partial crack surface of Muskhelishvili’s solutions is located outside the computational domain.In this paper,an improved alternating method,the spline fictitious boundary element alternating method(SFBEAM),based on infinite domain with the combination of spline fictitious boundary element method(SFBEM)and Muskhelishvili’s solutions is proposed to solve the edge crack problems.Since the SFBEM and Muskhelishvili’s solutions are obtained in the framework of infinite domain,no special treatment is needed for solving the problem of edge cracks.Different mixed boundary conditions edge crack problems with varies of computational parameters are given to certify the high precision,efficiency and applicability of the proposed method compared with other alternating methods and extend finite element method.
文摘In recent years,finite element analysis is increasingly being proposed in slope stability problems as a competitive method to traditional limit equilibrium methods(LEMs)which are known for their inherent deficiencies.However,the application of finite element method(FEM)to slope stability as a strength reduction method(SRM)or as finite element limit analysis(FELA)is not always a success for the drawbacks that characterize both methods.To increase the performance of finite element analysis in this problem,a new approach is proposed in this paper.It consists in gradually expanding the mobilized stress Mohr’s circles until the soil failure occurs according to a prescribed non-convergence criterion.The present approach called stress deviator increasing method(SDIM)is considered rigorous for three main reasons.Firstly,it preserves the definition of the factor of safety(FOS)as the ratio of soil shear strength to the mobilized shear stress.Secondly,it maintains the progressive development of shear stress resulting from the increase in the principal stress deviator on the same plane,on which the shear strength takes place.Thirdly,by introducing the concept of equivalent stress loading,the resulting trial stresses are checked against the violation of the actual yield criterion formed with the real strength parameters rather than those reduced by a trial factor.The new numerical procedure was encoded in a Fortran computer code called S^(4)DINA and verified by several examples.Comparisons with other numerical methods such as the SRM,gravity increasing method(GIM)or even FELA by assessing both the FOS and contours of equivalent plastic strains showed promising results.
文摘Each surface of roughness has different shape of asperity which is modeled with various shapes of analytical models. In this paper, the differences among various models of shape of asperity investigate using the Finite Element Method (FEM) and various analytical models. The contact stresses in rough surfaces are calculated analytically using various asperity shape models. Finite element analysis is also carried out assuming three types of material properties namely, the linear, the elastic-perfect plastic and the elastic-nonlinear hardening. The analytical results are compared with the results obtained by the finite element method. The results illustrate for using a deterministic approach which the numerical models are suitable. In hertz model, the result of force is very big in interface of causing deformation plastic, while Model Zhao has almost same result with FEM nonlinear property model. It is observed that the results obtained from Zhao’s model are generally in a better agreement with the results obtained from various finite element models especially in elastic-plastic and plastic zones, hence it may be concluded that Zhao’s model can be used for analyzing the rough surfaces in contact mechanics.
文摘The FEM analysis of stress and deformation of steel base of the special diamond saw blades are carried out.The small carbide round segments are welded on the side of blades in order to increase the side wear re sistance of blade.Comparing with conventional saw blade,the maximum stress val ues of reasonable special saw blade are reduced respectively about 17%20% 33%,and the maximum deformation values are reduced respectively about 26%22 %44.7%in thetangential(X),radial(Y)and axial(Z)direction.The stress conce ntration zone is decreased for the special structure diamond saw blade.The stru ctures of diamond saw blade with different number of hard material pellet are an alyzed and optimized.
基金The financial support from the National Natural Science Foundation of China under grants#11172012,#11472020 is gratefully acknowledged.
文摘Chemo-mechanical coupling behavior of materials is a transformation process between mechanical and chemical energy.In this paper,based on the coupled chemo-mechanical constitutive equations and governing equations during isothermal process,the equivalent integral forms of chemo-mechanical coupling governing equations and corresponding finite element procedure are obtained by using Hamilton’s principle.An isoparametric plane element for chemo-mechanical coupling is associated into ABAQUS finite element package through user element subroutine UEL.The numerical examples exhibit that the ionic concentration variation can cause mechanical deformation and mechanical action can produce redistribution of ionic concentration for hydrogels.It is proved that the present developed chemo-mechanical coupling finite element procedure can be utilized to model the coupling behavior of hydrogels effectively.
基金financially supported by the National Natural Science Foundation of China(Grant No.51349011)the Foundation of Si’chuan Educational Committee(Grant No.17ZB0452)+1 种基金the Innovation Team Project of Si’chuan Educational Committee(Grant No.18TD0019)the Longshan Academic Talent Research Support Program of the Southwest of Science and Technology(Grant Nos.18LZX715 and 18LZX410)
文摘In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton’s method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.
基金supported by Singapore AcRF RG59/08 (M52110092)Singapore NRF 2007 IDM-IDM002-010.
文摘In this paper, we present a theoretical analysis for linear finite element superconvergent gradient recovery on Par6 mesh, the dual of which is centroidal Voronoi tessellations with the lowest energy per unit volume and is the congruent cell predicted by the three-dimensional Gersho's conjecture. We show that the linear finite element solution uh and the linear interpolation uI have superclose gradient on Par6 meshes. Consequently, the gradient recovered from the finite element solution by using the superconvergence patch recovery method is superconvergent to Vu. A numerical example is presented to verify the theoretical result.
文摘A stabilized and convergent finite element formulation for the generalized Stokes problem is proposed and a posteriori analysis is performed to produce an error indicator. On this basis adaptive numerical method for solying the problem is developed . Numerical calculations are performed to confirm the reliability and effectiveness of the method.
文摘Lateral closing wedge osteotomy is a treatment for Kienbock’s disease. It is one of the most frequently used treatment options, which has been reported with relatively good long-term results. However, the results about the treatment are still controversial in some literatures and some key mechanisms are still not clear. The objective of the cur-rent study was to study the biomedical mechanism of the treatment. A finite element model was developed based on the geometry of carpal bones. Various situations in-cluding inclination angle changes by cutting the radial with 0?, 5?, 10? and 15? osteot-omy angles were studied. The effectiveness of the treatment was also studied for the carpal structure with abnormal positions of the lunate bone. The results show that the effectiveness of the stress reduction with the angle depends on many situations such as the initial morphology of the carpal structure.
基金supported by National Natural Science Foundation of China (11071226 11201122)
文摘In this article, we study the explicit expressions of the constants in the error estimate of the nonconforming finite element method. We explicitly obtain the approximation error estimate and the consistency error estimate for the Wilson's element without the regular assumption, respectively, which implies the final finite element error estimate. Such explicit a priori error estimates can be used as computable error bounds.