The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES...The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.展开更多
We extend the differential quadrature element method (DQEM) to the buckling analysis of uniformly in-plane loaded functionally graded (FG) plates fully or partially resting on the Pasternak model of elastic support. M...We extend the differential quadrature element method (DQEM) to the buckling analysis of uniformly in-plane loaded functionally graded (FG) plates fully or partially resting on the Pasternak model of elastic support. Material properties of the FG plate are graded in the thickness direction and assumed to obey a power law distribution of the volume fraction of the constituents. To set up the global eigenvalue equation, the plate is divided into sub-domains or elements and the generalized differential quadrature procedure is applied to discretize the governing, boundary and compatibility equations. By assembling discrete equations at all nodal points, the weighting coefficient and force matrices are derived. To validate the accuracy of this method, the results are compared with those of the exact solution and the finite element method. At the end, the effects of different variables and local elastic support arrangements on the buckling load factor are investigated.展开更多
The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh disto...The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh distortions. In a previous work, we constructed an 8-node quadrilateral spline element L8 using the triangular area coordinates and the B- net method, which can be insensitive to mesh distortions and possess the second order completeness in the Cartesian co- ordinates. In this paper, a thin plate spline element is devel- oped based on the spline element L8 and the refined tech- nique. Numerical examples show that the present element indeed possesses higher accuracy than the DKQ element for distorted meshes.展开更多
A general method to construct locking free Reissner-Mindlin plate elements is presented. According to this method the shear strain is replaced by its proper interpolation polynomial, which corresponds to the Kirchoff ...A general method to construct locking free Reissner-Mindlin plate elements is presented. According to this method the shear strain is replaced by its proper interpolation polynomial, which corresponds to the Kirchoff conditions at the interpolation points as the thickness of plate tends to zero, so the element is locking free. We construct two triangular elements by this method - a 3-node element and a 6-node element. The numerical results are provided.展开更多
The present study aims to analyze free vibration of thin skew plates made of functionally graded material(FGM)by using the weak form quadrature element method.The material properties vary continuously through the thic...The present study aims to analyze free vibration of thin skew plates made of functionally graded material(FGM)by using the weak form quadrature element method.The material properties vary continuously through the thickness according to a power-law form.A novel FGM skew plate element is formulated according to the neutral surface based plate theory and with the help of the differential quadrature rule.For verifications,Numerical results are compared with available data in literature.Results reveal that the non-dimensional frequency parameters of the FGM skew plates are independent of the power-law exponent and always proportional to those of homogeneous isotropic ones when the coupling and rotary inertias are neglected.In addition,employing the physical neutral surface based plate theory is equivalent to using the middle plane based plate theory with the reduced flexural modulus matrix.展开更多
The draping behavior of fabric is simulated by using four node quadrilateral thin plate elements with finite rotation. The finite element formulation is based on the total Lagrangian approach. An exact representatio...The draping behavior of fabric is simulated by using four node quadrilateral thin plate elements with finite rotation. The finite element formulation is based on the total Lagrangian approach. An exact representation of finite rotation is introduced. The strain energy function accounting for the material symmetry is obtained by the tensor representation theory. To avoid shear locking, the assumed strain technique for transverse shear is adopted. The conjugate gradient method with a proposed line search algorithm is employed to minimize energy and reach the final shape of fabric. The draping behavior of a rectangular piece of fabric over a rectangular table is simulated. (Author abstract) 9 Refs.展开更多
Based on the strain formulation of the quasi-conforming finite element, displacement functions are constructed which have definite physical meaning, and a conclusion can be obtained that the coefficients of the consta...Based on the strain formulation of the quasi-conforming finite element, displacement functions are constructed which have definite physical meaning, and a conclusion can be obtained that the coefficients of the constant and the linear strain are uniquely determined, and the quasi-conforming finite element method is convergent to constant strain. There are different methods for constructing the rigid displacement items, and different methods correspond to different order node errors, and this is different from ordinary displacement method finite element.展开更多
Similar to the very vast prior literature on analyzing laminated composite structures,“higher-order”and“layer-wise higher-order”plate and shell theories for functionally-graded(FG)materials and structures are also...Similar to the very vast prior literature on analyzing laminated composite structures,“higher-order”and“layer-wise higher-order”plate and shell theories for functionally-graded(FG)materials and structures are also widely popularized in the literature of the past two decades.However,such higher-order theories involve(1)postulating very complex assumptions for plate/shell kinematics in the thickness direction,(2)defining generalized variables of displacements,strains,and stresses,and(3)developing very complex governing equilibrium,compatibility,and constitutive equations in terms of newly-defined generalized kinematic and generalized kinetic variables.Their industrial applications are thus hindered by their inherent complexity,and the fact that it is difficult for end-users(front-line structural engineers)to completely understand all the newly-defined generalized DOFs for FEM in the higher-order and layer-wise theories.In an entirely different way,very simple 20-node and 27-node 3-D continuum solid-shell elements are developed in this paper,based on the simple theory of 3D solid mechanics,for static and dynamic analyses of functionally-graded plates and shells.A simple Over-Integration(a 4-point Gauss integration in the thickness direction)is used to evaluate the stiffness matrices of each element,while only a single element is used in the thickness direction without increasing the number of degrees of freedom.A stress-recovery approach is used to compute the distribution of transverse stresses by considering the equations of 3D elasticity in Cartesian as well as cylindrical polar coordinates.Comprehensive numerical results are presented for static and dynamic analyses of FG plates and shells,which agree well,either with the existing solutions in the published literature,or with the computationally very expensive solutions obtained by using simple 3D isoparametric elements(with standard Gauss Quadrature)available in NASTRAN(wherein many 3D elements are used in the thickness direction to capture the varying material properties).The effects of the material gradient index,the span-to-thickness ratio,the aspect ratio and the boundary conditions are also studied in the solutions of FG structures.Because the proposed methodology merely involves:(2)standard displacement DOFs at each node,(2)involves a simple 4-point Gaussian over-integration in the thickness direction,(3)relies only on the simple theory of solid mechanics,and(4)is capable of accurately and efficiently predicting the static and dynamical behavior of FG structures in a very simple and cost-effective manner,it is thus believed by the authors that the painstaking and cumbersome development of“higher-order”or“layer-wise higher-order”theories is not entirely necessary for the analyses of FG plates and shells.展开更多
In this study, Artificial Neural Network has been employed for analysis of triangular plate with different geometrical and loading parameters. Plates, having different sizes of concentric holes are analyzed. Finite el...In this study, Artificial Neural Network has been employed for analysis of triangular plate with different geometrical and loading parameters. Plates, having different sizes of concentric holes are analyzed. Finite element analysis for 81 cases is carried out using ANSYS Workbench 15.0 software. Using these data of FEM analysis an Artificial Neural Network has been trained. The successfully trained network is further used for analysis of four new cases which are also validated by using ANSYS Workbench 15.0 software.展开更多
Current patch test for Mindlin plate element only satisfies the zero shear deformation condition.The patch test of non-zero constant shear for Mindlin plate problem cannot be performed.For shell element, the patch tes...Current patch test for Mindlin plate element only satisfies the zero shear deformation condition.The patch test of non-zero constant shear for Mindlin plate problem cannot be performed.For shell element, the patch test does not even exist.Based on the theory of enhanced patch test proposed by Chen W J (2006),the authors proposed the enhanced patch test function for Mindlin plate and thin cylindrical shell elements.This enhanced patch test function can be used to assess the convergence of the Mindlin plate and cylindrical thin shell elements.展开更多
Flange earrings of strong anisotropic sheet metals in deep-drawing process are numerically analyzed by the elastic-plastic large deformation finite element formulation based on a discrete Kirchhoff triangle plate shel...Flange earrings of strong anisotropic sheet metals in deep-drawing process are numerically analyzed by the elastic-plastic large deformation finite element formulation based on a discrete Kirchhoff triangle plate shell element model. A Barlat-Lian anisotropic yield function and a quasi-flow corner theory are used in the present formulation. The numerical results are compared with the experimental ones of cylindrical cup drawing process. The focus of the present researches is on the numerical analysis and the constraining scheme of the flange earring of circular sheets with strong anisotropy in square cup drawing process.展开更多
Some theorems of compactly supported non-tensor product form two-dimension Daubechies wavelet were analysed carefully. Compactly supported non-tensor product form two-dimension wavelet was constructed, then non-tensor...Some theorems of compactly supported non-tensor product form two-dimension Daubechies wavelet were analysed carefully. Compactly supported non-tensor product form two-dimension wavelet was constructed, then non-tensor product form two dimension wavelet finite element was used to solve the deflection problem of elastic thin plate. The error order was researched. A numerical example was given at last.展开更多
This article aims to propose a finite element formulation based on Quasi-3D theory for the static bending analysis of functionally graded porous(FGP)sandwich plates.The FGP sandwich plates consist of three layers incl...This article aims to propose a finite element formulation based on Quasi-3D theory for the static bending analysis of functionally graded porous(FGP)sandwich plates.The FGP sandwich plates consist of three layers including the bottom skin of homogeneous metal,the top skin of fully ceramic and the FGP core layer with uneven porosity distribution.A quadrilateral(Q4)element with nine degrees of freedom(DOFs)per node is derived and employed in analyzing the static bending response of the plate under uniform and/or sinusoidally distributed loads.The accuracy of the present finite element formulation is verified by comparing the obtained numerical results with the published results in the literature.Then,some numerical examples are performed to examine the effects of the parameters including power-law index k and porosity coefficient on the static bending response of rectangular FGP sandwich plates.In addition,a problem with a complicated L-shape model is conducted to illustrate the superiority of the proposed finite element method.展开更多
基金funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under Grant number 107.02-2019.330。
文摘The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.
文摘We extend the differential quadrature element method (DQEM) to the buckling analysis of uniformly in-plane loaded functionally graded (FG) plates fully or partially resting on the Pasternak model of elastic support. Material properties of the FG plate are graded in the thickness direction and assumed to obey a power law distribution of the volume fraction of the constituents. To set up the global eigenvalue equation, the plate is divided into sub-domains or elements and the generalized differential quadrature procedure is applied to discretize the governing, boundary and compatibility equations. By assembling discrete equations at all nodal points, the weighting coefficient and force matrices are derived. To validate the accuracy of this method, the results are compared with those of the exact solution and the finite element method. At the end, the effects of different variables and local elastic support arrangements on the buckling load factor are investigated.
基金supported by the National Natural Science Foundation of China(11001037,11102037,11290143)the Fundamental Research Funds for the Central Universities(DUT13LK07)
文摘The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh distortions. In a previous work, we constructed an 8-node quadrilateral spline element L8 using the triangular area coordinates and the B- net method, which can be insensitive to mesh distortions and possess the second order completeness in the Cartesian co- ordinates. In this paper, a thin plate spline element is devel- oped based on the spline element L8 and the refined tech- nique. Numerical examples show that the present element indeed possesses higher accuracy than the DKQ element for distorted meshes.
文摘A general method to construct locking free Reissner-Mindlin plate elements is presented. According to this method the shear strain is replaced by its proper interpolation polynomial, which corresponds to the Kirchoff conditions at the interpolation points as the thickness of plate tends to zero, so the element is locking free. We construct two triangular elements by this method - a 3-node element and a 6-node element. The numerical results are provided.
基金supported by the Priority Academic Program Development of Jiangsu Higher Education Institutions
文摘The present study aims to analyze free vibration of thin skew plates made of functionally graded material(FGM)by using the weak form quadrature element method.The material properties vary continuously through the thickness according to a power-law form.A novel FGM skew plate element is formulated according to the neutral surface based plate theory and with the help of the differential quadrature rule.For verifications,Numerical results are compared with available data in literature.Results reveal that the non-dimensional frequency parameters of the FGM skew plates are independent of the power-law exponent and always proportional to those of homogeneous isotropic ones when the coupling and rotary inertias are neglected.In addition,employing the physical neutral surface based plate theory is equivalent to using the middle plane based plate theory with the reduced flexural modulus matrix.
文摘The draping behavior of fabric is simulated by using four node quadrilateral thin plate elements with finite rotation. The finite element formulation is based on the total Lagrangian approach. An exact representation of finite rotation is introduced. The strain energy function accounting for the material symmetry is obtained by the tensor representation theory. To avoid shear locking, the assumed strain technique for transverse shear is adopted. The conjugate gradient method with a proposed line search algorithm is employed to minimize energy and reach the final shape of fabric. The draping behavior of a rectangular piece of fabric over a rectangular table is simulated. (Author abstract) 9 Refs.
文摘Based on the strain formulation of the quasi-conforming finite element, displacement functions are constructed which have definite physical meaning, and a conclusion can be obtained that the coefficients of the constant and the linear strain are uniquely determined, and the quasi-conforming finite element method is convergent to constant strain. There are different methods for constructing the rigid displacement items, and different methods correspond to different order node errors, and this is different from ordinary displacement method finite element.
基金This research is supported by the Mechanics Section,Vehicle Technology Division,of the US Army Research Labs.The support of National Natural Science Foundation of China(grant No.11502069)Natural Science Foundation of Jiangsu Province(grant No.BK20140838)is also thankfully acknowledged.
文摘Similar to the very vast prior literature on analyzing laminated composite structures,“higher-order”and“layer-wise higher-order”plate and shell theories for functionally-graded(FG)materials and structures are also widely popularized in the literature of the past two decades.However,such higher-order theories involve(1)postulating very complex assumptions for plate/shell kinematics in the thickness direction,(2)defining generalized variables of displacements,strains,and stresses,and(3)developing very complex governing equilibrium,compatibility,and constitutive equations in terms of newly-defined generalized kinematic and generalized kinetic variables.Their industrial applications are thus hindered by their inherent complexity,and the fact that it is difficult for end-users(front-line structural engineers)to completely understand all the newly-defined generalized DOFs for FEM in the higher-order and layer-wise theories.In an entirely different way,very simple 20-node and 27-node 3-D continuum solid-shell elements are developed in this paper,based on the simple theory of 3D solid mechanics,for static and dynamic analyses of functionally-graded plates and shells.A simple Over-Integration(a 4-point Gauss integration in the thickness direction)is used to evaluate the stiffness matrices of each element,while only a single element is used in the thickness direction without increasing the number of degrees of freedom.A stress-recovery approach is used to compute the distribution of transverse stresses by considering the equations of 3D elasticity in Cartesian as well as cylindrical polar coordinates.Comprehensive numerical results are presented for static and dynamic analyses of FG plates and shells,which agree well,either with the existing solutions in the published literature,or with the computationally very expensive solutions obtained by using simple 3D isoparametric elements(with standard Gauss Quadrature)available in NASTRAN(wherein many 3D elements are used in the thickness direction to capture the varying material properties).The effects of the material gradient index,the span-to-thickness ratio,the aspect ratio and the boundary conditions are also studied in the solutions of FG structures.Because the proposed methodology merely involves:(2)standard displacement DOFs at each node,(2)involves a simple 4-point Gaussian over-integration in the thickness direction,(3)relies only on the simple theory of solid mechanics,and(4)is capable of accurately and efficiently predicting the static and dynamical behavior of FG structures in a very simple and cost-effective manner,it is thus believed by the authors that the painstaking and cumbersome development of“higher-order”or“layer-wise higher-order”theories is not entirely necessary for the analyses of FG plates and shells.
文摘In this study, Artificial Neural Network has been employed for analysis of triangular plate with different geometrical and loading parameters. Plates, having different sizes of concentric holes are analyzed. Finite element analysis for 81 cases is carried out using ANSYS Workbench 15.0 software. Using these data of FEM analysis an Artificial Neural Network has been trained. The successfully trained network is further used for analysis of four new cases which are also validated by using ANSYS Workbench 15.0 software.
基金Supported by the National Natural Science Foundation of China(Grant Nos.50479058 and 10672032)
文摘Current patch test for Mindlin plate element only satisfies the zero shear deformation condition.The patch test of non-zero constant shear for Mindlin plate problem cannot be performed.For shell element, the patch test does not even exist.Based on the theory of enhanced patch test proposed by Chen W J (2006),the authors proposed the enhanced patch test function for Mindlin plate and thin cylindrical shell elements.This enhanced patch test function can be used to assess the convergence of the Mindlin plate and cylindrical thin shell elements.
基金The project supported by the National Natural Science Foundation of China (19832020)Provincial Natural Science Foundation of Jilin, China (200000519)
文摘Flange earrings of strong anisotropic sheet metals in deep-drawing process are numerically analyzed by the elastic-plastic large deformation finite element formulation based on a discrete Kirchhoff triangle plate shell element model. A Barlat-Lian anisotropic yield function and a quasi-flow corner theory are used in the present formulation. The numerical results are compared with the experimental ones of cylindrical cup drawing process. The focus of the present researches is on the numerical analysis and the constraining scheme of the flange earring of circular sheets with strong anisotropy in square cup drawing process.
文摘Some theorems of compactly supported non-tensor product form two-dimension Daubechies wavelet were analysed carefully. Compactly supported non-tensor product form two-dimension wavelet was constructed, then non-tensor product form two dimension wavelet finite element was used to solve the deflection problem of elastic thin plate. The error order was researched. A numerical example was given at last.
文摘This article aims to propose a finite element formulation based on Quasi-3D theory for the static bending analysis of functionally graded porous(FGP)sandwich plates.The FGP sandwich plates consist of three layers including the bottom skin of homogeneous metal,the top skin of fully ceramic and the FGP core layer with uneven porosity distribution.A quadrilateral(Q4)element with nine degrees of freedom(DOFs)per node is derived and employed in analyzing the static bending response of the plate under uniform and/or sinusoidally distributed loads.The accuracy of the present finite element formulation is verified by comparing the obtained numerical results with the published results in the literature.Then,some numerical examples are performed to examine the effects of the parameters including power-law index k and porosity coefficient on the static bending response of rectangular FGP sandwich plates.In addition,a problem with a complicated L-shape model is conducted to illustrate the superiority of the proposed finite element method.