Based on the first-order shear deformation theory,a 3-node co-rotational triangular finite element formulation is developed for large deformation modeling of non-smooth,folded and multi-shell laminated composite struc...Based on the first-order shear deformation theory,a 3-node co-rotational triangular finite element formulation is developed for large deformation modeling of non-smooth,folded and multi-shell laminated composite structures.The two smaller components of the mid-surface normal vector of shell at a node are defined as nodal rotational variables in the co-rotational local coordinate system.In the global coordinate system,two smaller components of one vector,together with the smallest or second smallest component of another vector,of an orthogonal triad at a node on a non-smooth intersection of plates and/or shells are defined as rotational variables,whereas the two smaller components of the mid-surface normal vector at a node on the smooth part of the plate or shell(away from non-smooth intersections)are defined as rotational variables.All these vectorial rotational variables can be updated in an additive manner during an incremental solution procedure,and thus improve the computational efficiency in the nonlinear solution of these composite shell structures.Due to the commutativity of all nodal variables in calculating of the second derivatives of the local nodal variables with respect to global nodal variables,and the second derivatives of the strain energy functional with respect to local nodal variables,symmetric tangent stiffness matrices in local and global coordinate systems are obtained.To overcome shear locking,the assumed transverse shear strains obtained from the line-integration approach are employed.The reliability and computational accuracy of the present 3-node triangular shell finite element are verified through modeling two patch tests,several smooth and non-smooth laminated composite shells undergoing large displacements and large rotations.展开更多
The current research of supporting and transmission system in flywheel energy storage system(FESS) focuses on the low consumption design. However, friction loss is a non-negligible factor in the high-speed but lightwe...The current research of supporting and transmission system in flywheel energy storage system(FESS) focuses on the low consumption design. However, friction loss is a non-negligible factor in the high-speed but lightweight FESS energy and momentum storage with mechanical-type supporting system. In order to realize the support system without mechanical loss and to maximize the e ciency of the flywheel battery, a permanent magnet biased magnetic bearings(PMBMB) is applied to the FESS with the advantages of low loss, high critical speed, flexible controllability and compact structure. In this frame, the relevant research of three degrees of freedom(3-DOF) PMBMB for a new type FESS is carried out around the working principle, structural composition, coupling characteristics analysis, mathematical model, and structural design. In order to verify the performance of the 3-DOF PMBMB, the radial force mathematical model and the coupling determination equations of radial two DOF are calculated according to an equivalent magnetic circuit, and radial–axial coupling is analyzed through finite element analysis. Moreover, a control system is presented to solve the control problems in practical applications. The rotor returns to the balanced position in 0.05 s and maintains stable suspension. The displacement fluctuation is approximately 40 μm in the y direction and 30 μm in the x direction. Test results indicate that the dynamic rotor of the proposed flywheel energy storage system with PMBMB has excellent characteristics, such as good start-of-suspension performance and stable suspension characteristics. The proposed research provides the instruction to design and control a low loss support system for FESS.展开更多
By considering the characteristics of deformation of rotationally periodic structures under rotationally periodic loads, the periodic structure is divided into some identical substructures in this study. The degrees-o...By considering the characteristics of deformation of rotationally periodic structures under rotationally periodic loads, the periodic structure is divided into some identical substructures in this study. The degrees-of-freedom (DOFs) of joint nodes between the neighboring substructures are classified as master and slave ones. The stress and strain conditions of the whole structure are obtained by solving the elastic static equations for only one substructure by introducing the displacement constraints between master and slave DOFs. The complex constraint method is used to get the bifurcation buckling load and mode for the whole rotationally periodic structure by solving the eigenvalue problem for only one substructure without introducing any additional approximation. The finite element (FE) formulation of shell element of relative degrees of freedom (SERDF) in the buckling analysis is derived. Different measures of tackling internal degrees of freedom for different kinds of buckling problems and different stages of numerical analysis are presented. Some numerical examples are given to illustrate the high efficiency and validity of this method.展开更多
A new kind of super parametric finite elements for geometric nonlinear analysis of plates and shells is presented.Besides the nodes on the middle surface, additional virtual nodes are used to determine the normal to ...A new kind of super parametric finite elements for geometric nonlinear analysis of plates and shells is presented.Besides the nodes on the middle surface, additional virtual nodes are used to determine the normal to the middle surface.There are three displacement d.o.f.for each node of the element, and two transverse shear strains are taken as additional independent d.o.f.for each node on the middle surface.Therefore, the element is suitable for large rotation analysis of plates and shells. It can also be easily applied to the analysis of laminated and sandwich shells and plates.Numerical examples are given to show the accuracy and efficiency of the element.展开更多
In this article,we propose a new finite element spaceΛh for the expanded mixed finite element method(EMFEM)for second-order elliptic problems to guarantee its computing capability and reduce the computation cost.The ...In this article,we propose a new finite element spaceΛh for the expanded mixed finite element method(EMFEM)for second-order elliptic problems to guarantee its computing capability and reduce the computation cost.The new finite element spaceΛh is designed in such a way that the strong requirement V h⊂Λh in[9]is weakened to{v h∈V h;d i v v h=0}⊂Λh so that it needs fewer degrees of freedom than its classical counterpart.Furthermore,the newΛh coupled with the Raviart-Thomas space satisfies the inf-sup condition,which is crucial to the computation of mixed methods for its close relation to the behavior of the smallest nonzero eigenvalue of the stiff matrix,and thus the existence,uniqueness and optimal approximate capability of the EMFEM solution are proved for rectangular partitions in R d,d=2,3 and for triangular partitions in R 2.Also,the solvability of the EMFEM for triangular partition in R 3 can be directly proved without the inf-sup condition.Numerical experiments are conducted to confirm these theoretical findings.展开更多
Analyzing static and dynamic problems including composite structures has been of high significance in research efforts and industrial applications.In this article,equivalent single layer approach is utilized for dynam...Analyzing static and dynamic problems including composite structures has been of high significance in research efforts and industrial applications.In this article,equivalent single layer approach is utilized for dynamic finite element procedures of 3D composite beam as the building block of numerous composite structures.In this model,both displacement and strain fields are decomposed into cross-sectional and longitudinal components,called consistent geometric decomposition theorem.Then,the model is discretized using finite clement procedures.Two local coordinate systems and a global one are defined to decouple mechanical degrees of freedom.Furthermore,from the viewpoint of consistent geometric decomposition theorem,the transformation and element mass matrices for those systems are introduced here for the first time.The same decomposition idea can be used for developing element stiffiness matrix.Finally,comprehensive validations are conducted for the theory against experimental and numerical results in two case studies and for various conditions.展开更多
A finite element model for the supercavitating underwater vehicle was developed by employing 16-node shell elements of relative degrees of freedom.The nonlinear structural dynamic response was performed by introducing...A finite element model for the supercavitating underwater vehicle was developed by employing 16-node shell elements of relative degrees of freedom.The nonlinear structural dynamic response was performed by introducing the updated Lagrangian formulation.The numerical results indicate that there exists a critical thickness for the supercavitating plain shell for the considered velocity of the vehicle.The structure fails more easily because of instability with the thickness less than the critical value,while the structure maintains dynamic stability with the thickness greater than the critical value.As the velocity of the vehicle increases,the critical thickness for the plain shell increases accordingly.For the considered structural configuration,the critical thicknesses of plain shells are 5 and 7 mm for the velocities of 300 and 400 m/s,respectively.The structural stability is enhanced by using the stiffened configuration.With the shell configuration of nine ring stiffeners,the maximal displacement and von Mises stress of the supercavitating structure decrease by 25% and 17% for the velocity of 300 m/s,respectively.Compared with ring stiffeners,longitudinal stiffeners are more significant to improve structural dynamic performance and decrease the critical value of thickness of the shell for the supercavitating vehicle.展开更多
In this paper, a computation method has been developed so as to compare the finite element method (FEM) with the test results directly. The structure is divided into the 'master' and 'slave' degrees of...In this paper, a computation method has been developed so as to compare the finite element method (FEM) with the test results directly. The structure is divided into the 'master' and 'slave' degrees of freedom. The simplified model can be obtained with modal reduction. Then the design sensitivity analysis of the eigenvalues and eigenvectors has been carried out using the modal frequency and modal shape of the test. A two-story frame structure and a jacket model structure have been calculated. Meanwhile, the modified coefficient, the FEM computational and experimental values have been given. It has been shown that the FEM model modified using the test modal value is efficient.展开更多
基金This work was supported by National Natural Science Foundation of China under Grant 11672266.
文摘Based on the first-order shear deformation theory,a 3-node co-rotational triangular finite element formulation is developed for large deformation modeling of non-smooth,folded and multi-shell laminated composite structures.The two smaller components of the mid-surface normal vector of shell at a node are defined as nodal rotational variables in the co-rotational local coordinate system.In the global coordinate system,two smaller components of one vector,together with the smallest or second smallest component of another vector,of an orthogonal triad at a node on a non-smooth intersection of plates and/or shells are defined as rotational variables,whereas the two smaller components of the mid-surface normal vector at a node on the smooth part of the plate or shell(away from non-smooth intersections)are defined as rotational variables.All these vectorial rotational variables can be updated in an additive manner during an incremental solution procedure,and thus improve the computational efficiency in the nonlinear solution of these composite shell structures.Due to the commutativity of all nodal variables in calculating of the second derivatives of the local nodal variables with respect to global nodal variables,and the second derivatives of the strain energy functional with respect to local nodal variables,symmetric tangent stiffness matrices in local and global coordinate systems are obtained.To overcome shear locking,the assumed transverse shear strains obtained from the line-integration approach are employed.The reliability and computational accuracy of the present 3-node triangular shell finite element are verified through modeling two patch tests,several smooth and non-smooth laminated composite shells undergoing large displacements and large rotations.
基金Supported by National Natural Science Foundation of China(Grant Nos.51707082,51877101,51607080)Jiangsu Provincial Natural Science Foundation of China(Grant Nos.BK20170546,BK20150510)+1 种基金China Postdoctoral Science Foundation(Grant No.2017M620192)Priority Academic Program Development of Jiangsu Higher Education Institutions
文摘The current research of supporting and transmission system in flywheel energy storage system(FESS) focuses on the low consumption design. However, friction loss is a non-negligible factor in the high-speed but lightweight FESS energy and momentum storage with mechanical-type supporting system. In order to realize the support system without mechanical loss and to maximize the e ciency of the flywheel battery, a permanent magnet biased magnetic bearings(PMBMB) is applied to the FESS with the advantages of low loss, high critical speed, flexible controllability and compact structure. In this frame, the relevant research of three degrees of freedom(3-DOF) PMBMB for a new type FESS is carried out around the working principle, structural composition, coupling characteristics analysis, mathematical model, and structural design. In order to verify the performance of the 3-DOF PMBMB, the radial force mathematical model and the coupling determination equations of radial two DOF are calculated according to an equivalent magnetic circuit, and radial–axial coupling is analyzed through finite element analysis. Moreover, a control system is presented to solve the control problems in practical applications. The rotor returns to the balanced position in 0.05 s and maintains stable suspension. The displacement fluctuation is approximately 40 μm in the y direction and 30 μm in the x direction. Test results indicate that the dynamic rotor of the proposed flywheel energy storage system with PMBMB has excellent characteristics, such as good start-of-suspension performance and stable suspension characteristics. The proposed research provides the instruction to design and control a low loss support system for FESS.
文摘By considering the characteristics of deformation of rotationally periodic structures under rotationally periodic loads, the periodic structure is divided into some identical substructures in this study. The degrees-of-freedom (DOFs) of joint nodes between the neighboring substructures are classified as master and slave ones. The stress and strain conditions of the whole structure are obtained by solving the elastic static equations for only one substructure by introducing the displacement constraints between master and slave DOFs. The complex constraint method is used to get the bifurcation buckling load and mode for the whole rotationally periodic structure by solving the eigenvalue problem for only one substructure without introducing any additional approximation. The finite element (FE) formulation of shell element of relative degrees of freedom (SERDF) in the buckling analysis is derived. Different measures of tackling internal degrees of freedom for different kinds of buckling problems and different stages of numerical analysis are presented. Some numerical examples are given to illustrate the high efficiency and validity of this method.
文摘A new kind of super parametric finite elements for geometric nonlinear analysis of plates and shells is presented.Besides the nodes on the middle surface, additional virtual nodes are used to determine the normal to the middle surface.There are three displacement d.o.f.for each node of the element, and two transverse shear strains are taken as additional independent d.o.f.for each node on the middle surface.Therefore, the element is suitable for large rotation analysis of plates and shells. It can also be easily applied to the analysis of laminated and sandwich shells and plates.Numerical examples are given to show the accuracy and efficiency of the element.
基金supported by NSF of China grant 11971276H.Chen was supported by NSF of China grants 12171287,10971254 and 11471196+1 种基金H.Wang was supported by the ARO MURI Grant W911NF-15-1-0562by the National Science Foundation under Grant DMS-2012291.
文摘In this article,we propose a new finite element spaceΛh for the expanded mixed finite element method(EMFEM)for second-order elliptic problems to guarantee its computing capability and reduce the computation cost.The new finite element spaceΛh is designed in such a way that the strong requirement V h⊂Λh in[9]is weakened to{v h∈V h;d i v v h=0}⊂Λh so that it needs fewer degrees of freedom than its classical counterpart.Furthermore,the newΛh coupled with the Raviart-Thomas space satisfies the inf-sup condition,which is crucial to the computation of mixed methods for its close relation to the behavior of the smallest nonzero eigenvalue of the stiff matrix,and thus the existence,uniqueness and optimal approximate capability of the EMFEM solution are proved for rectangular partitions in R d,d=2,3 and for triangular partitions in R 2.Also,the solvability of the EMFEM for triangular partition in R 3 can be directly proved without the inf-sup condition.Numerical experiments are conducted to confirm these theoretical findings.
文摘Analyzing static and dynamic problems including composite structures has been of high significance in research efforts and industrial applications.In this article,equivalent single layer approach is utilized for dynamic finite element procedures of 3D composite beam as the building block of numerous composite structures.In this model,both displacement and strain fields are decomposed into cross-sectional and longitudinal components,called consistent geometric decomposition theorem.Then,the model is discretized using finite clement procedures.Two local coordinate systems and a global one are defined to decouple mechanical degrees of freedom.Furthermore,from the viewpoint of consistent geometric decomposition theorem,the transformation and element mass matrices for those systems are introduced here for the first time.The same decomposition idea can be used for developing element stiffiness matrix.Finally,comprehensive validations are conducted for the theory against experimental and numerical results in two case studies and for various conditions.
文摘A finite element model for the supercavitating underwater vehicle was developed by employing 16-node shell elements of relative degrees of freedom.The nonlinear structural dynamic response was performed by introducing the updated Lagrangian formulation.The numerical results indicate that there exists a critical thickness for the supercavitating plain shell for the considered velocity of the vehicle.The structure fails more easily because of instability with the thickness less than the critical value,while the structure maintains dynamic stability with the thickness greater than the critical value.As the velocity of the vehicle increases,the critical thickness for the plain shell increases accordingly.For the considered structural configuration,the critical thicknesses of plain shells are 5 and 7 mm for the velocities of 300 and 400 m/s,respectively.The structural stability is enhanced by using the stiffened configuration.With the shell configuration of nine ring stiffeners,the maximal displacement and von Mises stress of the supercavitating structure decrease by 25% and 17% for the velocity of 300 m/s,respectively.Compared with ring stiffeners,longitudinal stiffeners are more significant to improve structural dynamic performance and decrease the critical value of thickness of the shell for the supercavitating vehicle.
文摘In this paper, a computation method has been developed so as to compare the finite element method (FEM) with the test results directly. The structure is divided into the 'master' and 'slave' degrees of freedom. The simplified model can be obtained with modal reduction. Then the design sensitivity analysis of the eigenvalues and eigenvectors has been carried out using the modal frequency and modal shape of the test. A two-story frame structure and a jacket model structure have been calculated. Meanwhile, the modified coefficient, the FEM computational and experimental values have been given. It has been shown that the FEM model modified using the test modal value is efficient.