In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
The design of finite element analysis program using object-oriented programming (OOP) techniques is presented. The objects, classes and the subclasses used in the programming are explained. The system of classes libra...The design of finite element analysis program using object-oriented programming (OOP) techniques is presented. The objects, classes and the subclasses used in the programming are explained. The system of classes library of finite element analysis program and Windows-type Graphical User Interfaces by VC + + and its MFC are developed. The reliability, reusability and extensibility of program are enhanced. It is a reference to develop the large-scale, versatile and powerful systems of object-oriented finite element software.展开更多
In order to investigate the mechanical behavior of systems with complex architecture and a large number of contacting bodies,a finite element software,named LiToTac,has been developed by using the object-oriented prog...In order to investigate the mechanical behavior of systems with complex architecture and a large number of contacting bodies,a finite element software,named LiToTac,has been developed by using the object-oriented programming technique.This software,with an interactive graphical user interface,is able to handle highly non-linear problems including multiple contacts and large deformation.More importantly,the contact detection based on a hybrid three-stages methodology can be performed automatically,which is more efficient than the common strategies of pre-defining contact zones in commercial FEM software like ANSYS,ABAQUS,etc.In addition,the contact solver in LiToTac is portable between dynamic and quasi-static codes and can accurately solve contact coupled with friction in a reduced system.Several numerical examples are carried out to illustrate the functionality and capacity of the software package.展开更多
The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess v...The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.展开更多
This paper describes the object-oriented implementational method of finite element structural analysis, gives the basic concepts of the object-oriented method and objectoriented programming, develops a complete class ...This paper describes the object-oriented implementational method of finite element structural analysis, gives the basic concepts of the object-oriented method and objectoriented programming, develops a complete class hierarchy structure of object-oriented finite element structural analysis, and gives a part C+ + code description.展开更多
The Voronoi cell finite element method (VCFEM) is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric va...The Voronoi cell finite element method (VCFEM) is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric variational principle and quadratic programming method are developed for elastic-plastic Voronoi finite element analysis of two-dimensional problems. Finite element formulations are derived and a standard quadratic programming model is deduced from the elastic-plastic equations. Influence of microscopic heterogeneities on the overall mechanical response of heterogeneous materials is studied in detail. The overall properties of heterogeneous materials depend mostly on the size, shape and distribution of the material phases of the microstructure. Numerical examples are presented to demonstrate the validity and effectiveness of the method developed.展开更多
Recently an object-oriented approach has been applied in the fields of finite element analysis with a view to treating the various complexities within these. It has been demonstrated that finite element software desig...Recently an object-oriented approach has been applied in the fields of finite element analysis with a view to treating the various complexities within these. It has been demonstrated that finite element software designed using an object-oriented approach can be significantly more robust than traditional codes. This paper describes a special kind of implementation of object-oriented programming which is rather hybrid in nature, in the development of a finite element code for engineering analysis of metal working problems using C++, and discusses the advantages of this approach.展开更多
In the present study a Genetic Programing model (GP) proposed for the prediction of relative crest settlement of concrete faced rock fill dams. To this end information of 30 large dams constructed in seven countries a...In the present study a Genetic Programing model (GP) proposed for the prediction of relative crest settlement of concrete faced rock fill dams. To this end information of 30 large dams constructed in seven countries across the world is gathered with their reported settlements. The results showed that the GP model is able to estimate the dam settlement properly based on four properties, void ratio of dam’s body (e), height (H), vertical deformation modulus (Ev) and shape factor (Sc) of the dam. For verification of the model applicability, obtained results compared with other research methods such as Clements’s formula and the finite element model. The comparison showed that in all cases the GP model led to be more accurate than those of performed in literature. Also a proper compatibility between the GP model and the finite element model was perceived.展开更多
Bilinear theological lubrication mechanics provides an important basis for the designs of re- cently developed electrorheological(ER)'smart'journal bearings and those lubricated by mixed fluid-solid lubri- can...Bilinear theological lubrication mechanics provides an important basis for the designs of re- cently developed electrorheological(ER)'smart'journal bearings and those lubricated by mixed fluid-solid lubri- cants.But there is not yet a reliable and efficient numerical method for such a problem of non-Newtonian flu- id mechanics.In the present paper,a finite element method(FEM)together with mat hematical programming solution is successfully used to solve such a problem.A reliable and generalized numerical method for the designs of electrorheological 'smart' journal bearings and the bearings lubricated by mixed fluid- solid lubri- cant is presented.展开更多
In this paper, natural frequencies and vibration modes of spindl D1203 are calculated by us-ing structural analysis program SAP84. It is found that the second natural frequency of D1203spindle is too close to the rang...In this paper, natural frequencies and vibration modes of spindl D1203 are calculated by us-ing structural analysis program SAP84. It is found that the second natural frequency of D1203spindle is too close to the range of working speed. In order to raise the second natural frequencyand to reduce the vibration noise, simple optimal design of several structural parameters are alsodiscussed by means of the orthogonal design method. It is shown that the better design is toshorten the length of the spindle blade appropriately and to increase the elasticity of its uppersupport.展开更多
The rigid-plastic analysis of mental forming simulation is formulated as a discrete nonlinear mathematical programming problem with equality and inequality constraints by means of the finite element technique. An iter...The rigid-plastic analysis of mental forming simulation is formulated as a discrete nonlinear mathematical programming problem with equality and inequality constraints by means of the finite element technique. An iteration algorithm is used to solve this formulation, which distinguishes the integration points of the rigid zones and the plastic zones and solves a series of the quadratic programming to overcome the difficulties caused by the nonsmoothness and the nonlinearity of the objective function. This method has been used to carry out the rigid-plastic FEM analysis. An example is given to demonstrate the effectiveness of this method.展开更多
This paper deals with the limit analyses of perfect rigid-plastic continua.Based on the kinematic theorem of the limit analysis theory,a mathematical programming finite element formula for determining the upper bound ...This paper deals with the limit analyses of perfect rigid-plastic continua.Based on the kinematic theorem of the limit analysis theory,a mathematical programming finite element formula for determining the upper bound load multiplier has been established,and an iteration algorithm proposed accordingly.In this algorithm the plastic and rigid zones are distinguished for every iteration step,and the goal function is modified gradually.The difficulties caused by the nonsmoothness of the goal function are over- come.Some examples solved by this algorithm are presented.展开更多
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
文摘The design of finite element analysis program using object-oriented programming (OOP) techniques is presented. The objects, classes and the subclasses used in the programming are explained. The system of classes library of finite element analysis program and Windows-type Graphical User Interfaces by VC + + and its MFC are developed. The reliability, reusability and extensibility of program are enhanced. It is a reference to develop the large-scale, versatile and powerful systems of object-oriented finite element software.
基金We gratefully acknowledge the financial support of the National Key R&D Program of China(Grant No.2017YFB0703200)the National Natural Science Foundation of China(Grant No.11772274).
文摘In order to investigate the mechanical behavior of systems with complex architecture and a large number of contacting bodies,a finite element software,named LiToTac,has been developed by using the object-oriented programming technique.This software,with an interactive graphical user interface,is able to handle highly non-linear problems including multiple contacts and large deformation.More importantly,the contact detection based on a hybrid three-stages methodology can be performed automatically,which is more efficient than the common strategies of pre-defining contact zones in commercial FEM software like ANSYS,ABAQUS,etc.In addition,the contact solver in LiToTac is portable between dynamic and quasi-static codes and can accurately solve contact coupled with friction in a reduced system.Several numerical examples are carried out to illustrate the functionality and capacity of the software package.
基金part of the TPS projecta Vied-Newton PhD scholarship+1 种基金a Dixon scholarship from Imperial College London,UKthe Dean’s Fund from Imperial College London for financial support(2017-2020)。
文摘The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.
文摘This paper describes the object-oriented implementational method of finite element structural analysis, gives the basic concepts of the object-oriented method and objectoriented programming, develops a complete class hierarchy structure of object-oriented finite element structural analysis, and gives a part C+ + code description.
基金Project supported by the National Natural Science Foundation of China(Nos.10225212, 10421002 and 10332010)the NCET Program provided by the Ministry of Education and the National Key Basic Research Special Foundation of China (No.2005CB321704)
文摘The Voronoi cell finite element method (VCFEM) is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric variational principle and quadratic programming method are developed for elastic-plastic Voronoi finite element analysis of two-dimensional problems. Finite element formulations are derived and a standard quadratic programming model is deduced from the elastic-plastic equations. Influence of microscopic heterogeneities on the overall mechanical response of heterogeneous materials is studied in detail. The overall properties of heterogeneous materials depend mostly on the size, shape and distribution of the material phases of the microstructure. Numerical examples are presented to demonstrate the validity and effectiveness of the method developed.
文摘Recently an object-oriented approach has been applied in the fields of finite element analysis with a view to treating the various complexities within these. It has been demonstrated that finite element software designed using an object-oriented approach can be significantly more robust than traditional codes. This paper describes a special kind of implementation of object-oriented programming which is rather hybrid in nature, in the development of a finite element code for engineering analysis of metal working problems using C++, and discusses the advantages of this approach.
文摘In the present study a Genetic Programing model (GP) proposed for the prediction of relative crest settlement of concrete faced rock fill dams. To this end information of 30 large dams constructed in seven countries across the world is gathered with their reported settlements. The results showed that the GP model is able to estimate the dam settlement properly based on four properties, void ratio of dam’s body (e), height (H), vertical deformation modulus (Ev) and shape factor (Sc) of the dam. For verification of the model applicability, obtained results compared with other research methods such as Clements’s formula and the finite element model. The comparison showed that in all cases the GP model led to be more accurate than those of performed in literature. Also a proper compatibility between the GP model and the finite element model was perceived.
文摘Bilinear theological lubrication mechanics provides an important basis for the designs of re- cently developed electrorheological(ER)'smart'journal bearings and those lubricated by mixed fluid-solid lubri- cants.But there is not yet a reliable and efficient numerical method for such a problem of non-Newtonian flu- id mechanics.In the present paper,a finite element method(FEM)together with mat hematical programming solution is successfully used to solve such a problem.A reliable and generalized numerical method for the designs of electrorheological 'smart' journal bearings and the bearings lubricated by mixed fluid- solid lubri- cant is presented.
文摘In this paper, natural frequencies and vibration modes of spindl D1203 are calculated by us-ing structural analysis program SAP84. It is found that the second natural frequency of D1203spindle is too close to the range of working speed. In order to raise the second natural frequencyand to reduce the vibration noise, simple optimal design of several structural parameters are alsodiscussed by means of the orthogonal design method. It is shown that the better design is toshorten the length of the spindle blade appropriately and to increase the elasticity of its uppersupport.
文摘The rigid-plastic analysis of mental forming simulation is formulated as a discrete nonlinear mathematical programming problem with equality and inequality constraints by means of the finite element technique. An iteration algorithm is used to solve this formulation, which distinguishes the integration points of the rigid zones and the plastic zones and solves a series of the quadratic programming to overcome the difficulties caused by the nonsmoothness and the nonlinearity of the objective function. This method has been used to carry out the rigid-plastic FEM analysis. An example is given to demonstrate the effectiveness of this method.
基金The project supported by National Natural Science Foundation of China.
文摘This paper deals with the limit analyses of perfect rigid-plastic continua.Based on the kinematic theorem of the limit analysis theory,a mathematical programming finite element formula for determining the upper bound load multiplier has been established,and an iteration algorithm proposed accordingly.In this algorithm the plastic and rigid zones are distinguished for every iteration step,and the goal function is modified gradually.The difficulties caused by the nonsmoothness of the goal function are over- come.Some examples solved by this algorithm are presented.