On the basis of concept of element node topological analysis, the normalized element node topological matrices for finite element meshes are presented in the paper, including 3-node and 6-node triangular element, 4-n...On the basis of concept of element node topological analysis, the normalized element node topological matrices for finite element meshes are presented in the paper, including 3-node and 6-node triangular element, 4-node and 8-node quadrilateral element, 8-node and 20-node hexahedral element. It is beneficial to further analyzing topological characteristics of finite element models and automatic generation of meshes展开更多
It is known in the computational electromagnetics （CEM） that the node element has a relative wellconditioned matrix, but suffers from the spurious solution problem; whereas the edge element has no spurious solutions...It is known in the computational electromagnetics （CEM） that the node element has a relative wellconditioned matrix, but suffers from the spurious solution problem; whereas the edge element has no spurious solutions, but usually produces an ill-conditioned matrix. Particularly, when the mesh is over dense, the iterative solution of the matrix equation from edge element converges very slowly. Based on the node element and edge element, a node-edge element is presented, which has no spurious solutions and better-conditioned matrix. Numerical experiments demonstrate that the proposed node-edge element is more efficient than now-widely used edge element.展开更多
An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary ...An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.展开更多
Formulation and numerical evaluation of a novel four-node quadrilateral element with continuous nodal stress（Q4-CNS）are presented.Q4-CNS can be regarded as an improved hybrid FE-meshless four-node quadrilateral elem...Formulation and numerical evaluation of a novel four-node quadrilateral element with continuous nodal stress（Q4-CNS）are presented.Q4-CNS can be regarded as an improved hybrid FE-meshless four-node quadrilateral element（FE-LSPIM QUAD4）, which is a hybrid FE-meshless method.Derivatives of Q4-CNS are continuous at nodes, so the continuous nodal stress can be obtained without any smoothing operation.It is found that,compared with the standard four-node quadrilateral element（QUAD4）,Q4- CNS can achieve significantly better accuracy and higher convergence rate.It is also found that Q4-CNS exhibits high tolerance to mesh distortion.Moreover,since derivatives of Q4-CNS shape functions are continuous at nodes,Q4-CNS is potentially useful for the problem of bending plate and shell models.展开更多
How to automatically generate three-dimensional finite element Delaunay mesh by a peifected node connection method is introduced, where nodes are generated based on existing elements, instead of independence of node c...How to automatically generate three-dimensional finite element Delaunay mesh by a peifected node connection method is introduced, where nodes are generated based on existing elements, instead of independence of node creation and elements generation in traditional node connection method. Therefore, Ihe the difficulty about how to automatically create nodes in the traditional method is overcome.展开更多
A novel polygonal finite element method （PFEM） based on partition of unity is proposed, termed the virtual node method （VNM）. To test the performance of the present method, numerical examples are given for solid m...A novel polygonal finite element method （PFEM） based on partition of unity is proposed, termed the virtual node method （VNM）. To test the performance of the present method, numerical examples are given for solid mechanics problems. With a polynomial form, the VNM achieves better results than those of traditional PFEMs, including the Wachspress method and the mean value method in standard patch tests. Compared with the standard triangular FEM, the VNM can achieve better accuracy. With the ability to construct shape functions on polygonal elements, the VNM provides greater flexibility in mesh generation. Therefore, several fracture problems are studied to demonstrate the potential implementation. With the advantage of the VNM, the convenient refinement and remeshing strategy are applied.展开更多
Isopaxametric quadrilateral elements are widely used in the finite element method. However, they have a disadvantage of accuracy loss when elements are distorted. Spline functions have properties of simpleness and con...Isopaxametric quadrilateral elements are widely used in the finite element method. However, they have a disadvantage of accuracy loss when elements are distorted. Spline functions have properties of simpleness and conformality. A 17onode quadrilateral element has been developed using the bivaxiate quaxtic spline interpolation basis and the triangular area coordinates, which can exactly model the quartic displacement fields. Some appropriate examples are employed to illustrate that the element possesses high precision and is insensitive to mesh distortions.展开更多
A new method called node dynamic relaxation is proposed to simulate multilayer welding. A two dimensional plane strain model for multilayer welding is simulated and the results show that mesh distortion can be decreas...A new method called node dynamic relaxation is proposed to simulate multilayer welding. A two dimensional plane strain model for multilayer welding is simulated and the results show that mesh distortion can be decreased, and it is also found that the node dynamic relaxation is a kind of method to calculate welding deformation accurately by comparing experiment results with simulation results.展开更多
A locking-free rectangular Mindlin plate element with a new multi-resolution analysis (MRA) is proposed and a new finite element method is hence presented. The MRA framework is formulated out of a mutually nesting dis...A locking-free rectangular Mindlin plate element with a new multi-resolution analysis (MRA) is proposed and a new finite element method is hence presented. The MRA framework is formulated out of a mutually nesting displacement subspace sequence whose basis functions are constructed of scaling and shifting on the element domain of basic full node shape function. The basic full node shape function is constructed by extending the split node shape function of a traditional Mindlin plate element to other three quadrants around the coordinate zero point. As a result, a new rational MRA concept together with the resolution level (RL) is constituted for the element. The traditional 4-node rectangular Mindlin plate element and method is a mono-resolution one and also a special case of the proposed element and method. The meshing for the monoresolution plate element model is based on the empiricism while the RL adjusting for the multiresolution is laid on the rigorous mathematical basis. The analysis clarity of a plate structure is actually determined by the RL, not by the mesh. Thus, the accuracy of a plate structural analysis is replaced by the clarity, the irrational MRA by the rational and the mesh model by the RL that is the discretized model by the integrated.展开更多
This paper proposes a new, simple and efficient method for nonlinear simulation of arch dam cracking from the construction period to the operation period, which takes into account the arch dam construction process and...This paper proposes a new, simple and efficient method for nonlinear simulation of arch dam cracking from the construction period to the operation period, which takes into account the arch dam construction process and temperature loads. In the calculation mesh, the contact surface of pair nodes is located at places on the arch dam where cracking is possible. A new effective iterative method, the mixed finite element method for friction-contact problems, is improved and used for nonlinear simulation of the cracking process. The forces acting on the structure are divided into two parts： external forces and contact forces. The displacement of the structure is chosen as the basic variable and the nodal contact force in the possible contact region of the local coordinate system is chosen as the iterative variable, so that the nonlinear iterative process is only limited within the possible contact surface and is much more economical. This method was used to simulate the cracking process of the Shuanghe Arch Dam in Southwest China. In order to prove the validity and accuracy of this method and to study the effect of thermal stress on arch dam cracking, three schemes were designed for calculation. Numerical results agree with actual measured data, proving that it is feasible to use this method to simulate the entire process of nonlinear arch dam cracking.展开更多
A multi-resolution rectangular shell element with membrane-bending based on the Kirchhoff-Love theory is proposed. The multi-resolution analysis （MRA） framework is formulated out of a mutually nesting displacement s...A multi-resolution rectangular shell element with membrane-bending based on the Kirchhoff-Love theory is proposed. The multi-resolution analysis （MRA） framework is formulated out of a mutually nesting displacement subspace sequence, whose basis functions are constructed of scaling and shifting on the element domain of basic node shape functions. The basic node shape functions are constructed from shifting to other three quadrants around a specific node of a basic element in one quadrant and joining the corresponding node shape functions of four elements at the specific node. The MRA endows the proposed element with the resolution level （RL） to adjust the element node number, thus modulating structural analysis accuracy accordingly. The node shape functions of Kronecker delta property make the treatment of element boundary condition quite convenient and enable the stiffness matrix and the loading column vectors of the proposed element to be automatically acquired through quadraturing around nodes in RL adjusting. As a result, the traditional 4-node rectangular shell element is a mono-resolution one and also a special case of the proposed element. The accuracy of a structural analysis is actually determined by the RL, not by the mesh. The simplicity and clarity of node shape function formulation with the Kronecker delta property, and the rational MRA enable the proposed element method to be implemented more rationally, easily and efficiently than the conventional mono-resolution rectangular shell element method or other corresponding MRA methods.展开更多
文摘On the basis of concept of element node topological analysis, the normalized element node topological matrices for finite element meshes are presented in the paper, including 3-node and 6-node triangular element, 4-node and 8-node quadrilateral element, 8-node and 20-node hexahedral element. It is beneficial to further analyzing topological characteristics of finite element models and automatic generation of meshes
文摘It is known in the computational electromagnetics （CEM） that the node element has a relative wellconditioned matrix, but suffers from the spurious solution problem; whereas the edge element has no spurious solutions, but usually produces an ill-conditioned matrix. Particularly, when the mesh is over dense, the iterative solution of the matrix equation from edge element converges very slowly. Based on the node element and edge element, a node-edge element is presented, which has no spurious solutions and better-conditioned matrix. Numerical experiments demonstrate that the proposed node-edge element is more efficient than now-widely used edge element.
文摘An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.
文摘Formulation and numerical evaluation of a novel four-node quadrilateral element with continuous nodal stress（Q4-CNS）are presented.Q4-CNS can be regarded as an improved hybrid FE-meshless four-node quadrilateral element（FE-LSPIM QUAD4）, which is a hybrid FE-meshless method.Derivatives of Q4-CNS are continuous at nodes, so the continuous nodal stress can be obtained without any smoothing operation.It is found that,compared with the standard four-node quadrilateral element（QUAD4）,Q4- CNS can achieve significantly better accuracy and higher convergence rate.It is also found that Q4-CNS exhibits high tolerance to mesh distortion.Moreover,since derivatives of Q4-CNS shape functions are continuous at nodes,Q4-CNS is potentially useful for the problem of bending plate and shell models.
基金This project is supported by Provincial Natural Science foundation of Guangdong!(970516)
文摘How to automatically generate three-dimensional finite element Delaunay mesh by a peifected node connection method is introduced, where nodes are generated based on existing elements, instead of independence of node creation and elements generation in traditional node connection method. Therefore, Ihe the difficulty about how to automatically create nodes in the traditional method is overcome.
文摘A novel polygonal finite element method （PFEM） based on partition of unity is proposed, termed the virtual node method （VNM）. To test the performance of the present method, numerical examples are given for solid mechanics problems. With a polynomial form, the VNM achieves better results than those of traditional PFEMs, including the Wachspress method and the mean value method in standard patch tests. Compared with the standard triangular FEM, the VNM can achieve better accuracy. With the ability to construct shape functions on polygonal elements, the VNM provides greater flexibility in mesh generation. Therefore, several fracture problems are studied to demonstrate the potential implementation. With the advantage of the VNM, the convenient refinement and remeshing strategy are applied.
基金supported by the Natural Science Foundation of China China (Nos. 60533060, 10672032, and 10726067)the Science Foundation of Dalian University of Technology (No. SFDUT07001)
文摘Isopaxametric quadrilateral elements are widely used in the finite element method. However, they have a disadvantage of accuracy loss when elements are distorted. Spline functions have properties of simpleness and conformality. A 17onode quadrilateral element has been developed using the bivaxiate quaxtic spline interpolation basis and the triangular area coordinates, which can exactly model the quartic displacement fields. Some appropriate examples are employed to illustrate that the element possesses high precision and is insensitive to mesh distortions.
文摘A new method called node dynamic relaxation is proposed to simulate multilayer welding. A two dimensional plane strain model for multilayer welding is simulated and the results show that mesh distortion can be decreased, and it is also found that the node dynamic relaxation is a kind of method to calculate welding deformation accurately by comparing experiment results with simulation results.
文摘A locking-free rectangular Mindlin plate element with a new multi-resolution analysis (MRA) is proposed and a new finite element method is hence presented. The MRA framework is formulated out of a mutually nesting displacement subspace sequence whose basis functions are constructed of scaling and shifting on the element domain of basic full node shape function. The basic full node shape function is constructed by extending the split node shape function of a traditional Mindlin plate element to other three quadrants around the coordinate zero point. As a result, a new rational MRA concept together with the resolution level (RL) is constituted for the element. The traditional 4-node rectangular Mindlin plate element and method is a mono-resolution one and also a special case of the proposed element and method. The meshing for the monoresolution plate element model is based on the empiricism while the RL adjusting for the multiresolution is laid on the rigorous mathematical basis. The analysis clarity of a plate structure is actually determined by the RL, not by the mesh. Thus, the accuracy of a plate structural analysis is replaced by the clarity, the irrational MRA by the rational and the mesh model by the RL that is the discretized model by the integrated.
基金supported by the National Nature Science Foundation of China (Grant No 90510017)
文摘This paper proposes a new, simple and efficient method for nonlinear simulation of arch dam cracking from the construction period to the operation period, which takes into account the arch dam construction process and temperature loads. In the calculation mesh, the contact surface of pair nodes is located at places on the arch dam where cracking is possible. A new effective iterative method, the mixed finite element method for friction-contact problems, is improved and used for nonlinear simulation of the cracking process. The forces acting on the structure are divided into two parts： external forces and contact forces. The displacement of the structure is chosen as the basic variable and the nodal contact force in the possible contact region of the local coordinate system is chosen as the iterative variable, so that the nonlinear iterative process is only limited within the possible contact surface and is much more economical. This method was used to simulate the cracking process of the Shuanghe Arch Dam in Southwest China. In order to prove the validity and accuracy of this method and to study the effect of thermal stress on arch dam cracking, three schemes were designed for calculation. Numerical results agree with actual measured data, proving that it is feasible to use this method to simulate the entire process of nonlinear arch dam cracking.
基金financial support by the Open Foundation of Chongqing Key Laboratory of Geomechanics and Geoenvironment Protection(Logistical Engineering University)(No.GKLGGP 2013-02)
文摘A multi-resolution rectangular shell element with membrane-bending based on the Kirchhoff-Love theory is proposed. The multi-resolution analysis （MRA） framework is formulated out of a mutually nesting displacement subspace sequence, whose basis functions are constructed of scaling and shifting on the element domain of basic node shape functions. The basic node shape functions are constructed from shifting to other three quadrants around a specific node of a basic element in one quadrant and joining the corresponding node shape functions of four elements at the specific node. The MRA endows the proposed element with the resolution level （RL） to adjust the element node number, thus modulating structural analysis accuracy accordingly. The node shape functions of Kronecker delta property make the treatment of element boundary condition quite convenient and enable the stiffness matrix and the loading column vectors of the proposed element to be automatically acquired through quadraturing around nodes in RL adjusting. As a result, the traditional 4-node rectangular shell element is a mono-resolution one and also a special case of the proposed element. The accuracy of a structural analysis is actually determined by the RL, not by the mesh. The simplicity and clarity of node shape function formulation with the Kronecker delta property, and the rational MRA enable the proposed element method to be implemented more rationally, easily and efficiently than the conventional mono-resolution rectangular shell element method or other corresponding MRA methods.
