A new kind of quadrilateral assumed stress hy- brid membrane element with drilling degrees of freedom and a traction-free inclined side has been developed based on an extended Hellinger-Reissner principle which is est...A new kind of quadrilateral assumed stress hy- brid membrane element with drilling degrees of freedom and a traction-free inclined side has been developed based on an extended Hellinger-Reissner principle which is established by expanding the essential terms of the assumed stress field as polynomials in the natural coordinates of the element. The homogeneous equilibrium equations are imposed in a variational sense through the internal displacements which are also expanded in the natural coordinates, while the tractionfree conditions along the inclined side are satisfied exactly. The use of such special element in the finite element solution is shown to be highly accurate when only a very coarse element mesh is used for plates with V-shaped rounded notches or inclined sides.展开更多
A set of basic deformation modes for hybrid stress finite elements are directly derived from the element displacement field. Subsequently, by employing the so-called united orthogonal conditions, a new orthogonalizati...A set of basic deformation modes for hybrid stress finite elements are directly derived from the element displacement field. Subsequently, by employing the so-called united orthogonal conditions, a new orthogonalization method is proposed. The result- ing orthogonal basic deformation modes exhibit simple and clear physical meanings. In addition, they do not involve any material parameters, and thus can be efficiently used to examine the element performance and serve as a unified tool to assess different hybrid elements. Thereafter, a convenient approach for the identification of spurious zero-energy modes is presented using the positive definiteness property of a flexibility matrix. More- over, based on the orthogonality relationship between the given initial stress modes and the orthogonal basic deformation modes, an alternative method of assumed stress modes to formulate a hybrid element free of spurious modes is discussed. It is found that the orthogonality of the basic deformation modes is the sufficient and necessary condition for the suppression of spurious zero-energy modes. Numerical examples of 2D 4-node quadrilateral elements and 3D 8-node hexahedral elements are illustrated in detail to demonstrate the efficiency of the proposed orthogonal basic deformation mode method.展开更多
The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identica...The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identical hybrid element. The Hilbert stress subspace of the assumed stress modes is established. So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt's method. Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency, is improved greatly. The numerical examples show that the method is effective.展开更多
The new methods to determine the zero-energy deformation modes in the hybrid elements and the zero-energy stress modes in their assumed stress fields are presented by the natural deformation modes of the elements. And...The new methods to determine the zero-energy deformation modes in the hybrid elements and the zero-energy stress modes in their assumed stress fields are presented by the natural deformation modes of the elements. And the formula of the additional element deformation rigidity due to additional mode into the assumed stress field is derived. Based on, it is concluded in theory that the zero-energy stress mode cannot suppress the zero-energy deformation modes but increase the extra rigidity to the nonzero-energy deformation modes of the element instead. So they should not be employed to assume the stress field. In addition, the parasitic stress modes will produce the spurious parasitic energy and result the element behaving over rigidity. Thus, they should not be used into the assumed stress field even though they can suppress the zero-energy deformation modes of the element. The numerical examples show the performance of the elements including the zero-energy stress modes or the parasitic stress modes.展开更多
In the process of shale gas exploitation,there exits two difficult problems:one is the real numerical simulation of a tremendous number of holes in actual shale;the other is the fluid–solid coupling problem involved ...In the process of shale gas exploitation,there exits two difficult problems:one is the real numerical simulation of a tremendous number of holes in actual shale;the other is the fluid–solid coupling problem involved in holes,where the difficulty of transition at the interface between the Eulerian grid and the Lagrangian grid becomes the most important.In response to these two problems,this paper establishes an element model with both fluid and solid.At the fluid–solid interface,the equilibrium condition of the surface force is introduced to obtain the modified complementary energy functional,and a new hybrid stress element with fluid is derived.The comparison of the simulation results with those of the ordinary commercial finite element software verifies the effectiveness and efficiency of this element,and proves its applicability in the problem of shale with numerous holes.Furthermore,this element can be extended to general problems of solid with fluid in.展开更多
A 3-dimensional hybrid stress element with a traction-free cylindrical surface based on amodified complementary energy principle has been derived for efficient and accurate analysis of stressconcentration around circu...A 3-dimensional hybrid stress element with a traction-free cylindrical surface based on amodified complementary energy principle has been derived for efficient and accurate analysis of stressconcentration around circular cutouts in thin to thick laminated composites. New expressions of sixstress components are developed by using three stress-functions in cylindrical co-ordinates, so that thehomogeneous equilibrium equations, the interlayer surface transverse-stresses and the traction-freeboundary condition on the cylindrical surface are satisfied exactly, while the interelement traction conti-nuity has been relaxed via the Lagrange multiplier method. Transverse-shear deformation effects areincorporated in each layer with displacement continuity enforced along interlayer surface. Selected ex-amples are used to demonstrate the efficiency and accuracy of the present special element.展开更多
A novel hybrid-stress finite element method is proposed for constructing simple 4-node quadrilateral plane elements, and the new element is denoted as HH4-3fl here. Firstly, the theoretical basis of the traditional hy...A novel hybrid-stress finite element method is proposed for constructing simple 4-node quadrilateral plane elements, and the new element is denoted as HH4-3fl here. Firstly, the theoretical basis of the traditional hybrid-stress elements, i.e., the Hellinger-Reissner variational principle, is replaced by the Hamilton variational principle, in which the number of the stress variables is reduced from 3 to 2. Secondly, three stress parameters and corresponding trial functions are introduced into the system equations. Thirdly, the displacement fields of the conventional bilinear isoparametric element are employed in the new models. Finally, from the stationary condition, the stress parameters can be expressed in terms of the displacement parameters, and thus the new element stiffness matrices can be obtained. Since the required number of stress variables in the Hamilton variational principle is less than that in the Hellinger-Reissner variational principle, and no additional incompatible displacement modes are considered, the new hybrid-stress element is simpler than the traditional ones. Furthermore, in order to improve the accuracy of the stress solutions, two enhanced post-processing schemes are also proposed for element HH4-3β. Numerical examples show that the proposed model exhibits great improvements in both displacement and stress solutions, implying that the proposed technique is an effective way for developing simple finite element models with high performance.展开更多
In this paper, on the basis of the incremental Reissner variational principle.a nonlinear finite element analysis has been accomplished and a formulation of hybrid stress element has been presented for incompressible ...In this paper, on the basis of the incremental Reissner variational principle.a nonlinear finite element analysis has been accomplished and a formulation of hybrid stress element has been presented for incompressible Mooney rubber-like materials. The corrected terms of the non-equilibrium force and the incompressibility deviation are considered in the formulation. The computed values of numerical example agree very closely with the exact solution.展开更多
Additive manufacturing(AM)of metals often results in parts with unfavorable mechanical properties.Laser peening(LP)is a high strain rate mechanical surface treatment that hammers a workpiece and induces favorable mech...Additive manufacturing(AM)of metals often results in parts with unfavorable mechanical properties.Laser peening(LP)is a high strain rate mechanical surface treatment that hammers a workpiece and induces favorable mechanical properties.Peening strain hardens a surface and imparts compressive residual stresses improving the mechanical properties of a material.This work investigates the role of LP on layer-by-layer processing of 3D printed metals using finite element analysis.The objective is to understand temporal and spatial residual stress development after thermal and mechanical cancellation caused by cyclically coupling printing and peening.Results indicate layer peening frequency is a critical process parameter affecting residual stress redistribution and highly interdependent on the heat generated by the printing process.Optimum hybrid process conditions were found to exists that favorably enhance mechanical properties.With this study,hybrid-AM has ushered in the next evolutionary step in AM and has the potential to profoundly change the way high value metal goods are manufactured.展开更多
In this paper, we discuss an adaptive hybrid stress finite element method on quadri- lateral meshes for linear elasticity problems. To deal with hanging nodes arising in the adaptive mesh refinement, we propose new tr...In this paper, we discuss an adaptive hybrid stress finite element method on quadri- lateral meshes for linear elasticity problems. To deal with hanging nodes arising in the adaptive mesh refinement, we propose new transition types of hybrid stress quadrilateral elements with 5 to 7 nodes. In particular, we derive a priori error estimation for the 5- node transition hybrid stress element to show that it is free from Poisson-locking, in the sense that the error bound in the a priori estimate is independent of the Lam~ constant A. We introduce~ for quadrilateral meshes, refinement/coarsening algorithms, which do not require storing the refinement tree explicitly, and give an adaptive algorithm. Finally, we provide some numerical results.展开更多
The sensitivity problem to mesh distortion and the low accuracy problem of the stress solutions are two inherent difficulties in the finite element method.By applying the fundamental analytical solutions (in global Ca...The sensitivity problem to mesh distortion and the low accuracy problem of the stress solutions are two inherent difficulties in the finite element method.By applying the fundamental analytical solutions (in global Cartesian coordinates) to the Airy stress function of the anisotropic materials,8-and 12-node plane quadrilateral hybrid stress-function (HS-F) elements are successfully developed based on the principle of the minimum complementary energy.Numerical results show that the present new elements exhibit much better and more robust performance in both displacement and stress solutions than those obtained from other models.They can still perform very well even when the element shapes degenerate into a triangle and a concave quadrangle.It is also demonstrated that the proposed construction procedure is an effective way for developing shape-free finite element models which can completely overcome the sensitivity problem to mesh distortion and can produce highly accurate stress solutions.展开更多
Superconvergence and recovery type a posteriori error estimators are analyzed for Pian and Sumihara's 4-node hybrid stress quadrilateral finite element method for linear elasticity problems. Superconvergence of or...Superconvergence and recovery type a posteriori error estimators are analyzed for Pian and Sumihara's 4-node hybrid stress quadrilateral finite element method for linear elasticity problems. Superconvergence of order O(h^(1+min){α,1}) is established for both the displacement approximation in H^1-norm and the stress approximation in L^2-norm under a mesh assumption, where α > 0 is a parameter characterizing the distortion of meshes from parallelograms to quadrilaterals. Recovery type approximations for the displacement gradients and the stress tensor are constructed, and a posteriori error estimators based on the recovered quantities are shown to be asymptotically exact. Numerical experiments confirm the theoretical results.展开更多
This paper presents eight-node solid-shell elements for geometric non-linear analyze of piezoelectric structures. To subdue shear, trapezoidal and thickness locking, the assumed natural strain method and an ad hoc mod...This paper presents eight-node solid-shell elements for geometric non-linear analyze of piezoelectric structures. To subdue shear, trapezoidal and thickness locking, the assumed natural strain method and an ad hoc modified generalized laminate stiffness matrix are employed. With the generalized stresses arising from the modified generalized laminate stiffness matrix assumed to be independent from the ones obtained from the displacement, an extended Hellinger-Reissner functional can be derived. By choosing the assumed generalized stresses similar to the assumed stresses of a previous solid ele- ment, a hybrid-stress solid-shell element is formulated. The presented finite shell element is able to model arbitrary curved shell structures. Non-linear numerical examples demonstrate the ability of the proposed model to analyze nonlinear piezoelectric devices.展开更多
文摘A new kind of quadrilateral assumed stress hy- brid membrane element with drilling degrees of freedom and a traction-free inclined side has been developed based on an extended Hellinger-Reissner principle which is established by expanding the essential terms of the assumed stress field as polynomials in the natural coordinates of the element. The homogeneous equilibrium equations are imposed in a variational sense through the internal displacements which are also expanded in the natural coordinates, while the tractionfree conditions along the inclined side are satisfied exactly. The use of such special element in the finite element solution is shown to be highly accurate when only a very coarse element mesh is used for plates with V-shaped rounded notches or inclined sides.
基金Project supported by the National Natural Science Foundation of China(No.10972188)the Fundamental Research Funds for the Central Universities of China(No.2010121073)the Scientific Program of Fujian Province of China(No.2007F3096)
文摘A set of basic deformation modes for hybrid stress finite elements are directly derived from the element displacement field. Subsequently, by employing the so-called united orthogonal conditions, a new orthogonalization method is proposed. The result- ing orthogonal basic deformation modes exhibit simple and clear physical meanings. In addition, they do not involve any material parameters, and thus can be efficiently used to examine the element performance and serve as a unified tool to assess different hybrid elements. Thereafter, a convenient approach for the identification of spurious zero-energy modes is presented using the positive definiteness property of a flexibility matrix. More- over, based on the orthogonality relationship between the given initial stress modes and the orthogonal basic deformation modes, an alternative method of assumed stress modes to formulate a hybrid element free of spurious modes is discussed. It is found that the orthogonality of the basic deformation modes is the sufficient and necessary condition for the suppression of spurious zero-energy modes. Numerical examples of 2D 4-node quadrilateral elements and 3D 8-node hexahedral elements are illustrated in detail to demonstrate the efficiency of the proposed orthogonal basic deformation mode method.
文摘The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identical hybrid element. The Hilbert stress subspace of the assumed stress modes is established. So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt's method. Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency, is improved greatly. The numerical examples show that the method is effective.
文摘The new methods to determine the zero-energy deformation modes in the hybrid elements and the zero-energy stress modes in their assumed stress fields are presented by the natural deformation modes of the elements. And the formula of the additional element deformation rigidity due to additional mode into the assumed stress field is derived. Based on, it is concluded in theory that the zero-energy stress mode cannot suppress the zero-energy deformation modes but increase the extra rigidity to the nonzero-energy deformation modes of the element instead. So they should not be employed to assume the stress field. In addition, the parasitic stress modes will produce the spurious parasitic energy and result the element behaving over rigidity. Thus, they should not be used into the assumed stress field even though they can suppress the zero-energy deformation modes of the element. The numerical examples show the performance of the elements including the zero-energy stress modes or the parasitic stress modes.
基金This work was supported by the National Natural Science Foundation of China(Grant Nos.11572142 and 12072135).
文摘In the process of shale gas exploitation,there exits two difficult problems:one is the real numerical simulation of a tremendous number of holes in actual shale;the other is the fluid–solid coupling problem involved in holes,where the difficulty of transition at the interface between the Eulerian grid and the Lagrangian grid becomes the most important.In response to these two problems,this paper establishes an element model with both fluid and solid.At the fluid–solid interface,the equilibrium condition of the surface force is introduced to obtain the modified complementary energy functional,and a new hybrid stress element with fluid is derived.The comparison of the simulation results with those of the ordinary commercial finite element software verifies the effectiveness and efficiency of this element,and proves its applicability in the problem of shale with numerous holes.Furthermore,this element can be extended to general problems of solid with fluid in.
基金The work was supported by the National Natural Science Foundation of China (Grant No. 100Tz064) .
文摘A 3-dimensional hybrid stress element with a traction-free cylindrical surface based on amodified complementary energy principle has been derived for efficient and accurate analysis of stressconcentration around circular cutouts in thin to thick laminated composites. New expressions of sixstress components are developed by using three stress-functions in cylindrical co-ordinates, so that thehomogeneous equilibrium equations, the interlayer surface transverse-stresses and the traction-freeboundary condition on the cylindrical surface are satisfied exactly, while the interelement traction conti-nuity has been relaxed via the Lagrange multiplier method. Transverse-shear deformation effects areincorporated in each layer with displacement continuity enforced along interlayer surface. Selected ex-amples are used to demonstrate the efficiency and accuracy of the present special element.
基金supported by the National Natural Science Foundation of China (10872108,10876100)the Program for New Century Excellent Talents in University (NCET-07-0477)the National Basic Research Program of China (2010CB832701)
文摘A novel hybrid-stress finite element method is proposed for constructing simple 4-node quadrilateral plane elements, and the new element is denoted as HH4-3fl here. Firstly, the theoretical basis of the traditional hybrid-stress elements, i.e., the Hellinger-Reissner variational principle, is replaced by the Hamilton variational principle, in which the number of the stress variables is reduced from 3 to 2. Secondly, three stress parameters and corresponding trial functions are introduced into the system equations. Thirdly, the displacement fields of the conventional bilinear isoparametric element are employed in the new models. Finally, from the stationary condition, the stress parameters can be expressed in terms of the displacement parameters, and thus the new element stiffness matrices can be obtained. Since the required number of stress variables in the Hamilton variational principle is less than that in the Hellinger-Reissner variational principle, and no additional incompatible displacement modes are considered, the new hybrid-stress element is simpler than the traditional ones. Furthermore, in order to improve the accuracy of the stress solutions, two enhanced post-processing schemes are also proposed for element HH4-3β. Numerical examples show that the proposed model exhibits great improvements in both displacement and stress solutions, implying that the proposed technique is an effective way for developing simple finite element models with high performance.
文摘In this paper, on the basis of the incremental Reissner variational principle.a nonlinear finite element analysis has been accomplished and a formulation of hybrid stress element has been presented for incompressible Mooney rubber-like materials. The corrected terms of the non-equilibrium force and the incompressibility deviation are considered in the formulation. The computed values of numerical example agree very closely with the exact solution.
基金supported by in part by the National Science Foundation through the awards CAREER #1846478 and STTR #1521188
文摘Additive manufacturing(AM)of metals often results in parts with unfavorable mechanical properties.Laser peening(LP)is a high strain rate mechanical surface treatment that hammers a workpiece and induces favorable mechanical properties.Peening strain hardens a surface and imparts compressive residual stresses improving the mechanical properties of a material.This work investigates the role of LP on layer-by-layer processing of 3D printed metals using finite element analysis.The objective is to understand temporal and spatial residual stress development after thermal and mechanical cancellation caused by cyclically coupling printing and peening.Results indicate layer peening frequency is a critical process parameter affecting residual stress redistribution and highly interdependent on the heat generated by the printing process.Optimum hybrid process conditions were found to exists that favorably enhance mechanical properties.With this study,hybrid-AM has ushered in the next evolutionary step in AM and has the potential to profoundly change the way high value metal goods are manufactured.
文摘In this paper, we discuss an adaptive hybrid stress finite element method on quadri- lateral meshes for linear elasticity problems. To deal with hanging nodes arising in the adaptive mesh refinement, we propose new transition types of hybrid stress quadrilateral elements with 5 to 7 nodes. In particular, we derive a priori error estimation for the 5- node transition hybrid stress element to show that it is free from Poisson-locking, in the sense that the error bound in the a priori estimate is independent of the Lam~ constant A. We introduce~ for quadrilateral meshes, refinement/coarsening algorithms, which do not require storing the refinement tree explicitly, and give an adaptive algorithm. Finally, we provide some numerical results.
基金supported by the National Natural Science Foundation of China(Grant No.10872108,10876100)the Program for New Century Excellent Talents in University(Grant No. NCET-07-0477)+1 种基金the National Basic Research Program of China(Grant No. 2010CB832701)ASFC
文摘The sensitivity problem to mesh distortion and the low accuracy problem of the stress solutions are two inherent difficulties in the finite element method.By applying the fundamental analytical solutions (in global Cartesian coordinates) to the Airy stress function of the anisotropic materials,8-and 12-node plane quadrilateral hybrid stress-function (HS-F) elements are successfully developed based on the principle of the minimum complementary energy.Numerical results show that the present new elements exhibit much better and more robust performance in both displacement and stress solutions than those obtained from other models.They can still perform very well even when the element shapes degenerate into a triangle and a concave quadrangle.It is also demonstrated that the proposed construction procedure is an effective way for developing shape-free finite element models which can completely overcome the sensitivity problem to mesh distortion and can produce highly accurate stress solutions.
基金supported by National Natural Science Foundation of China (Grant No. 11171239)Major Research Plan of National Natural Science Foundation of China (Grant No. 91430105)
文摘Superconvergence and recovery type a posteriori error estimators are analyzed for Pian and Sumihara's 4-node hybrid stress quadrilateral finite element method for linear elasticity problems. Superconvergence of order O(h^(1+min){α,1}) is established for both the displacement approximation in H^1-norm and the stress approximation in L^2-norm under a mesh assumption, where α > 0 is a parameter characterizing the distortion of meshes from parallelograms to quadrilaterals. Recovery type approximations for the displacement gradients and the stress tensor are constructed, and a posteriori error estimators based on the recovered quantities are shown to be asymptotically exact. Numerical experiments confirm the theoretical results.
基金Supported by the National Natural Science Foundation of China (Grant No. 10672111)the Major Project of the Natural Science Foundation of Jiangsu Province (Grant No. BK2006725)
文摘This paper presents eight-node solid-shell elements for geometric non-linear analyze of piezoelectric structures. To subdue shear, trapezoidal and thickness locking, the assumed natural strain method and an ad hoc modified generalized laminate stiffness matrix are employed. With the generalized stresses arising from the modified generalized laminate stiffness matrix assumed to be independent from the ones obtained from the displacement, an extended Hellinger-Reissner functional can be derived. By choosing the assumed generalized stresses similar to the assumed stresses of a previous solid ele- ment, a hybrid-stress solid-shell element is formulated. The presented finite shell element is able to model arbitrary curved shell structures. Non-linear numerical examples demonstrate the ability of the proposed model to analyze nonlinear piezoelectric devices.
