A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue prob...A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue problem,the method is spectral-correct and spurious-free.Stability and error estimates are obtained,including the interpolation error estimates and the error estimates between the finite element solution and the exact solution.The method is suitable for singular solution as well as smooth solution,and consequently,the method is valid for nonconvex domains which may have a number of reentrant corners.Of course,the method is suitable for arbitrary quadrilaterals(under the usual shape-regular condition).展开更多
For k given graphs H_(1),...,H_(k) with k≥2,the k-color Ramsey number R(H_(1),...,H_(k)) represents the minimum integer N with the following property:if the edges of the complete graph K_(N) are colored with k colors...For k given graphs H_(1),...,H_(k) with k≥2,the k-color Ramsey number R(H_(1),...,H_(k)) represents the minimum integer N with the following property:if the edges of the complete graph K_(N) are colored with k colors,then there exists some i with 1≤i≤k such that K_(N) has a subgraph in color i isomorphic to H_(i).Let C_(m) be a cycle of length m and K_(1,n) a star of order n+1.In this paper,we systematically introduce the latest research progress on star-quadrilateral Ramsey numbers and provide an overview of Ramsey numbers concerning quadrilaterals,including multicolor cases.展开更多
In this paper, two new efficient fully nonconforming polynomial finite element methods on arbitrary convex quadrilaterals are constructed to solve the Brinkman problem. Both elements have 12 local degrees of freedom a...In this paper, two new efficient fully nonconforming polynomial finite element methods on arbitrary convex quadrilaterals are constructed to solve the Brinkman problem. Both elements have 12 local degrees of freedom and are uniformly first-order convergent. For the first element, we carefully design a shape function space of only degree four, which is convenient for practical computation. Meanwhile, the velocity solution obtained by our second element has O(h^(2)) convergence order in the case of Darcy limit. These new elements can be regarded as modifications of a known element to effectively improve the computational efficiency, only the shape function space is properly changed. We verify our theoretical findings by numerical examples.展开更多
A new quadrilateral finite element IQ4 is developed for the free vibration of carbon nanotube-reinforced composite(CNTRC)perforated plates with a central cutout.By enriching the membrane part and incorporating a proje...A new quadrilateral finite element IQ4 is developed for the free vibration of carbon nanotube-reinforced composite(CNTRC)perforated plates with a central cutout.By enriching the membrane part and incorporating a projected shear technique,the IQ4 element is proposed to address the known limitations of the standard Q4 element,such as shear locking and limited consistency in the coupling ofmembrane-bending components.The proposed element is formulated within the FSDT-based framework and assessed through benchmark tests to verify its convergence and accuracy.The governing equations are obtained via theweak formofHamilton’s principle.Particular attention is given to the influence of carbon nanotube volume fraction,distribution patterns,and boundary conditions on the fundamental frequency response of CNTRC plates with cutouts.In addition,a parametric study is conducted to assess the influence of cutout geometric configuration,shape,and size ratios on the vibrational response of the CNTRC plate.The numerical results demonstrate that the formulated IQ4 element provides stable and accurate estimations of natural frequencies,even in the presence of a cutout and the coupled effects of the non-uniform distribution of reinforcement through the plate thickness.The developed formulation is expected to contribute to the structural design and optimization of advanced lightweight systems,particularly in aerospace and mechanical engineering applications.展开更多
In the analysis of functionally graded materials （FGMs）, the uncoupled approach is used broadly, which is based on homogenized material property and ignores the effect Of local micro-structural interaction. The high...In the analysis of functionally graded materials （FGMs）, the uncoupled approach is used broadly, which is based on homogenized material property and ignores the effect Of local micro-structural interaction. The higher-order theory for FGMs （HOTFGM） is a coupled approach that explicitly takes the effect of micro-structural gradation and the local interaction of the spatially variable inclusion phase into account. Based on the HOTFGM, this article presents a quadrilateral element-based method for the calculation of multi-scale temperature field （QTF）. In this method, the discrete cells are quadrilateral including rectangular while the surface-averaged quantities are the primary variables which replace the coefficients employed in the temperature function. In contrast with the HOTFGM, this method improves the efficiency, eliminates the restriction of being rectangular cells and expands the solution scale. The presented results illustrate the efficiency of the QTF and its advantages in analyzing FGMs.展开更多
Based on the construction of reference e le ment and bilinear transformation, a quasi-Wilson element for arbitrary narrow q uadrilateral is presented. Using the interpolation Theorem for narrow quadrilate ral isoparam...Based on the construction of reference e le ment and bilinear transformation, a quasi-Wilson element for arbitrary narrow q uadrilateral is presented. Using the interpolation Theorem for narrow quadrilate ral isoparametric finite element and related methods, the bounds of interpolatio n error for arbitrary narrow quadrilateral quasi-Wilson element are obtained in case when the condition ρ K/h K≥σ 0】0 is not satisfied, where h K is the diameter of the element K and ρ K is the diameter of an ins cribed circle in K. The interpolation error is O(h2 K) in the L2( K)-norm and O(h K) in the H1(K) -norm provided that the in terpolated function belongs to H2(K).展开更多
Recently, some new quadrilateral finite elements were successfully developed by the Quadrilateral Area Coordinate （QAC） method. Compared with those traditional models using isoparametric coordinates, these new model...Recently, some new quadrilateral finite elements were successfully developed by the Quadrilateral Area Coordinate （QAC） method. Compared with those traditional models using isoparametric coordinates, these new models are less sensitive to mesh distortion. In this paper, a new displacement-based, 4-node 20-DOF （5-DOF per node） quadrilateral bending element based on the first-order shear deformation theory for analysis of arbitrary laminated composite plates is presented. Its bending part is based on the element AC-MQ4, a recent-developed high-performance Mindlin-Reissner plate element formulated by QAC method and the generalized conforming condition method; and its in-plane displacement fields are interpolated by bilinear shape functions in isoparametric coordinates. Furthermore, the hybrid post-rocessing procedure, which was firstly proposed by the authors, is employed again to improve the stress solutions, especially for the transverse shear stresses. The resulting element, denoted as AC-MQ4-LC, exhibits excellent performance in all linear static and dynamic numerical examples. It demonstrates again that the QAC method, the generalized conforming condition method, and the hybrid post-processing procedure are efficient tools for developing simple, effective and reliable finite element models.展开更多
A new 4-node quadrilateral flat shell element is developed for geometrically nonlinear analyses of thin and moderately thick laminated shell structures. The fiat shell element is constructed by combining a quadrilater...A new 4-node quadrilateral flat shell element is developed for geometrically nonlinear analyses of thin and moderately thick laminated shell structures. The fiat shell element is constructed by combining a quadrilateral area co- ordinate method （QAC） based membrane element AGQ6- II, and a Timoshenko beam function （TBF） method based shear deformable plate bending element ARS-Q12. In order to model folded plates and connect with beam elements, the drilling stiffness is added to the element stiffness matrix based on the mixed variational principle. The transverse shear rigidity matrix, based on the first-order shear deformation theory （FSDT）, for the laminated composite plate is evaluated using the transverse equilibrium conditions, while the shear correction factors are not needed. The conventional TBF methods are also modified to efficiently calculate the element stiffness for laminate. The new shell element is extended to large deflection and post-buckling analyses of isotropic and laminated composite shells based on the element independent corotational formulation. Numerical re- sults show that the present shell element has an excellent numerical performance for the test examples, and is applicable to stiffened plates.展开更多
Free vibration analysis of quadrilateral multilayered graphene sheets（MLGS） embedded in polymer matrix is carried out employing nonlocal continuum mechanics.The principle of virtual work is employed to derive the eq...Free vibration analysis of quadrilateral multilayered graphene sheets（MLGS） embedded in polymer matrix is carried out employing nonlocal continuum mechanics.The principle of virtual work is employed to derive the equations of motion.The Galerkin method in conjunction with the natural coordinates of the nanoplate is used as a basis for the analysis.The dependence of small scale effect on thickness,elastic modulus,polymer matrix stiffness and interaction coefficient between two adjacent sheets is illustrated.The non-dimensional natural frequencies of skew,rhombic,trapezoidal and rectangular MLGS are obtained with various geometrical parameters and mode numbers taken into account,and for each case the effects of the small length scale are investigated.展开更多
Two cases of the nested configurations in R3 consisting of two regular quadrilaterals are discussed. One case of them do not form central configuration, the other case can be central configuration. In the second case ...Two cases of the nested configurations in R3 consisting of two regular quadrilaterals are discussed. One case of them do not form central configuration, the other case can be central configuration. In the second case the existence and uniqueness of the central configuration are studied. If the configuration is a central configuration, then all masses of outside layer are equivalent, similar to the masses of inside layer. At the same time the following relation between r(the ratio of the sizes) and mass ratio b = m/m must be satisfied in which the masses at outside layer are not less than the masses at inside layer, and the solution of this kind of central configuration is unique for the given ratio (6) of masses.展开更多
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.展开更多
To simulate two-dimensional free-surface flows with complex boundaries directly and accurately, a novel VOF (Volume-of-fluid) method based on unstructured quadrilateral mesh is presented. Without introducing any compl...To simulate two-dimensional free-surface flows with complex boundaries directly and accurately, a novel VOF (Volume-of-fluid) method based on unstructured quadrilateral mesh is presented. Without introducing any complicated boundary treatment or artificial diffusion, this method treated curved boundaries directly by utilizing the inherent merit of unstructured mesh in fitting curves. The PLIC (Piecewise Linear Interface Calculation) method was adopted to obtain a second-order accurate linearized reconstruction approximation and the MLER (Modified Lagrangian-Eulerian Re-map) method was introduced to advect fluid volumes on unstructured mesh. Moreover, an analytical relation for the interface’s line constant vs. the volume clipped by the interface was developed so as to improve the method’s efficiency. To validate this method, a comprehensive series of large straining advection tests were performed. Numerical results provide convincing evidences for the method’s high volume conservative accuracy and second-order shape error convergence rate. Also, a dramatic improvement on computational accuracy over its unstructured triangular mesh counterpart is checked.展开更多
In this paper, a new method of topological cleanup for quadrilateral mesh is presented. The method first selects a patch of mesh around an irregular node. It then seeks the best connection of the selected patch accord...In this paper, a new method of topological cleanup for quadrilateral mesh is presented. The method first selects a patch of mesh around an irregular node. It then seeks the best connection of the selected patch according to its irregular valence using a new topological operation： small polygon reconnection （SPR）. By replacing the original patch with an optimal one that has less irregular valence, mesh quality can be improved. Three applications based on the proposed approach are enumerated： （1） improving the quality of a quadrilateral mesh, （2） converting a triangular mesh to a quadrilateral one, and （3） adapting a triangle generator to a quadrilateral one. The presented method is highly effective in all three applications.展开更多
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.展开更多
Quadrilateral finite elements for linear micropolar continuum theory are developed using linked interpolation.In order to satisfy convergence criteria,the newly presented finite elements are modified using the Petrov-...Quadrilateral finite elements for linear micropolar continuum theory are developed using linked interpolation.In order to satisfy convergence criteria,the newly presented finite elements are modified using the Petrov-Galerkin method in which different interpolation is used for the test and trial functions.The elements are tested through four numerical examples consisting of a set of patch tests,a cantilever beam in pure bending and a stress concentration problem and compared with the analytical solution and quadrilateral micropolar finite elements with standard Lagrangian interpolation.In the higher-order patch test,the performance of the first-order element is significantly improved.However,since the problems analysed are already describable with quadratic polynomials,the enhancement due to linked interpolation for higher-order elements could not be highlighted.All the presented elements also faithfully reproduce the micropolar effects in the stress concentration analysis,but the enhancement here is negligible with respect to standard Lagrangian elements,since the higher-order polynomials in this example are not needed.展开更多
Accuracy of a simulation strongly depends on the grid quality. Here, quality means orthogonality at the boundaries and quasi-orthogonality within the critical regions, smoothness, bounded aspect ratios and solution ad...Accuracy of a simulation strongly depends on the grid quality. Here, quality means orthogonality at the boundaries and quasi-orthogonality within the critical regions, smoothness, bounded aspect ratios and solution adaptive behaviour. It is not recommended to refine the parts of the domain where the solution shows little variation. It is desired to concentrate grid points and cells in the part of the domain where the solution shows strong gradients or variations. We present a simple, effective and com- putationally efficient approach for quadrilateral mesh adaptation. Several numerical examples are presented for supporting our claim.展开更多
To provide a iheoreiical basis for h-type .finire elemenl analysis with quadrilateralelements, in the present paper, the h-convergence of quadrilateral elements is established. whose related lemmas and theorems being ...To provide a iheoreiical basis for h-type .finire elemenl analysis with quadrilateralelements, in the present paper, the h-convergence of quadrilateral elements is established. whose related lemmas and theorems being presented, and therefore, theerror estimate problems are investigated.展开更多
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.展开更多
We propose a newmethod to generate surface quadrilateralmesh by calculating a globally defined parameterization with feature constraints.In the field of quadrilateral generation with features,the cross field methods a...We propose a newmethod to generate surface quadrilateralmesh by calculating a globally defined parameterization with feature constraints.In the field of quadrilateral generation with features,the cross field methods are wellknown because of their superior performance in feature preservation.The methods based on metrics are popular due to their sound theoretical basis,especially the Ricci flow algorithm.The cross field methods’major part,the Poisson equation,is challenging to solve in three dimensions directly.When it comes to cases with a large number of elements,the computational costs are expensive while the methods based on metrics are on the contrary.In addition,an appropriate initial value plays a positive role in the solution of the Poisson equation,and this initial value can be obtained from the Ricci flow algorithm.So we combine the methods based on metric with the cross field methods.We use the discrete dynamic Ricci flow algorithm to generate an initial value for the Poisson equation,which speeds up the solution of the equation and ensures the convergence of the computation.Numerical experiments show that our method is effective in generating a quadrilateral mesh for models with features,and the quality of the quadrilateral mesh is reliable.展开更多
With the development of advanced imaging technology, digital images are widely used. This paper proposes an automatic quadrilateral mesh generation algorithm for multi-colour imaged structures. It takes an original ar...With the development of advanced imaging technology, digital images are widely used. This paper proposes an automatic quadrilateral mesh generation algorithm for multi-colour imaged structures. It takes an original arbitrary digital image as an input for automatic quadrilateral mesh generation, this includes removing the noise, extracting and smoothing the boundary geometries between different colours, and automatic all-quad mesh generation with the above boundaries as constraints. An application example is provided to demonstrate the usefulness and effectiveness of the proposed approach.展开更多
基金supported by the National Natural Science Foundation of China(12401482)the second author was supported by the National Natural Science Foundation of China(12371371,12261160361,11971366)supported by the Open Research Fund of Hubei Key Laboratory of Computational Science,Wuhan University.
文摘A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue problem,the method is spectral-correct and spurious-free.Stability and error estimates are obtained,including the interpolation error estimates and the error estimates between the finite element solution and the exact solution.The method is suitable for singular solution as well as smooth solution,and consequently,the method is valid for nonconvex domains which may have a number of reentrant corners.Of course,the method is suitable for arbitrary quadrilaterals(under the usual shape-regular condition).
基金supported by NSFC(Nos.12161141003,11931006)supported by NSFC(Nos.11801520,12171436,12271489)supported by NSFC(No.11601527)。
文摘For k given graphs H_(1),...,H_(k) with k≥2,the k-color Ramsey number R(H_(1),...,H_(k)) represents the minimum integer N with the following property:if the edges of the complete graph K_(N) are colored with k colors,then there exists some i with 1≤i≤k such that K_(N) has a subgraph in color i isomorphic to H_(i).Let C_(m) be a cycle of length m and K_(1,n) a star of order n+1.In this paper,we systematically introduce the latest research progress on star-quadrilateral Ramsey numbers and provide an overview of Ramsey numbers concerning quadrilaterals,including multicolor cases.
基金Supported by the National Natural Science Foundation of China(Grant No.12201254)。
文摘In this paper, two new efficient fully nonconforming polynomial finite element methods on arbitrary convex quadrilaterals are constructed to solve the Brinkman problem. Both elements have 12 local degrees of freedom and are uniformly first-order convergent. For the first element, we carefully design a shape function space of only degree four, which is convenient for practical computation. Meanwhile, the velocity solution obtained by our second element has O(h^(2)) convergence order in the case of Darcy limit. These new elements can be regarded as modifications of a known element to effectively improve the computational efficiency, only the shape function space is properly changed. We verify our theoretical findings by numerical examples.
文摘A new quadrilateral finite element IQ4 is developed for the free vibration of carbon nanotube-reinforced composite(CNTRC)perforated plates with a central cutout.By enriching the membrane part and incorporating a projected shear technique,the IQ4 element is proposed to address the known limitations of the standard Q4 element,such as shear locking and limited consistency in the coupling ofmembrane-bending components.The proposed element is formulated within the FSDT-based framework and assessed through benchmark tests to verify its convergence and accuracy.The governing equations are obtained via theweak formofHamilton’s principle.Particular attention is given to the influence of carbon nanotube volume fraction,distribution patterns,and boundary conditions on the fundamental frequency response of CNTRC plates with cutouts.In addition,a parametric study is conducted to assess the influence of cutout geometric configuration,shape,and size ratios on the vibrational response of the CNTRC plate.The numerical results demonstrate that the formulated IQ4 element provides stable and accurate estimations of natural frequencies,even in the presence of a cutout and the coupled effects of the non-uniform distribution of reinforcement through the plate thickness.The developed formulation is expected to contribute to the structural design and optimization of advanced lightweight systems,particularly in aerospace and mechanical engineering applications.
基金National Natural Science Foundation of China (2009ZB52028,05C52013)Ph.D. Programs Foundation of Ministry of Education of China (20070287039)
文摘In the analysis of functionally graded materials （FGMs）, the uncoupled approach is used broadly, which is based on homogenized material property and ignores the effect Of local micro-structural interaction. The higher-order theory for FGMs （HOTFGM） is a coupled approach that explicitly takes the effect of micro-structural gradation and the local interaction of the spatially variable inclusion phase into account. Based on the HOTFGM, this article presents a quadrilateral element-based method for the calculation of multi-scale temperature field （QTF）. In this method, the discrete cells are quadrilateral including rectangular while the surface-averaged quantities are the primary variables which replace the coefficients employed in the temperature function. In contrast with the HOTFGM, this method improves the efficiency, eliminates the restriction of being rectangular cells and expands the solution scale. The presented results illustrate the efficiency of the QTF and its advantages in analyzing FGMs.
文摘Based on the construction of reference e le ment and bilinear transformation, a quasi-Wilson element for arbitrary narrow q uadrilateral is presented. Using the interpolation Theorem for narrow quadrilate ral isoparametric finite element and related methods, the bounds of interpolatio n error for arbitrary narrow quadrilateral quasi-Wilson element are obtained in case when the condition ρ K/h K≥σ 0】0 is not satisfied, where h K is the diameter of the element K and ρ K is the diameter of an ins cribed circle in K. The interpolation error is O(h2 K) in the L2( K)-norm and O(h K) in the H1(K) -norm provided that the in terpolated function belongs to H2(K).
基金The project is supported by the National Natural Science Foundation of China(10502028)the Special Foundation for the Authors of the Nationwide(China)Excellent Doctoral Dissertation(200242)the Science Research Foundation of China Agricultural University(2004016).
文摘Recently, some new quadrilateral finite elements were successfully developed by the Quadrilateral Area Coordinate （QAC） method. Compared with those traditional models using isoparametric coordinates, these new models are less sensitive to mesh distortion. In this paper, a new displacement-based, 4-node 20-DOF （5-DOF per node） quadrilateral bending element based on the first-order shear deformation theory for analysis of arbitrary laminated composite plates is presented. Its bending part is based on the element AC-MQ4, a recent-developed high-performance Mindlin-Reissner plate element formulated by QAC method and the generalized conforming condition method; and its in-plane displacement fields are interpolated by bilinear shape functions in isoparametric coordinates. Furthermore, the hybrid post-rocessing procedure, which was firstly proposed by the authors, is employed again to improve the stress solutions, especially for the transverse shear stresses. The resulting element, denoted as AC-MQ4-LC, exhibits excellent performance in all linear static and dynamic numerical examples. It demonstrates again that the QAC method, the generalized conforming condition method, and the hybrid post-processing procedure are efficient tools for developing simple, effective and reliable finite element models.
文摘A new 4-node quadrilateral flat shell element is developed for geometrically nonlinear analyses of thin and moderately thick laminated shell structures. The fiat shell element is constructed by combining a quadrilateral area co- ordinate method （QAC） based membrane element AGQ6- II, and a Timoshenko beam function （TBF） method based shear deformable plate bending element ARS-Q12. In order to model folded plates and connect with beam elements, the drilling stiffness is added to the element stiffness matrix based on the mixed variational principle. The transverse shear rigidity matrix, based on the first-order shear deformation theory （FSDT）, for the laminated composite plate is evaluated using the transverse equilibrium conditions, while the shear correction factors are not needed. The conventional TBF methods are also modified to efficiently calculate the element stiffness for laminate. The new shell element is extended to large deflection and post-buckling analyses of isotropic and laminated composite shells based on the element independent corotational formulation. Numerical re- sults show that the present shell element has an excellent numerical performance for the test examples, and is applicable to stiffened plates.
文摘Free vibration analysis of quadrilateral multilayered graphene sheets（MLGS） embedded in polymer matrix is carried out employing nonlocal continuum mechanics.The principle of virtual work is employed to derive the equations of motion.The Galerkin method in conjunction with the natural coordinates of the nanoplate is used as a basis for the analysis.The dependence of small scale effect on thickness,elastic modulus,polymer matrix stiffness and interaction coefficient between two adjacent sheets is illustrated.The non-dimensional natural frequencies of skew,rhombic,trapezoidal and rectangular MLGS are obtained with various geometrical parameters and mode numbers taken into account,and for each case the effects of the small length scale are investigated.
基金Supported by the NSF of China(10231010)Supported by the NSF of CQSXXY (20030104)
文摘Two cases of the nested configurations in R3 consisting of two regular quadrilaterals are discussed. One case of them do not form central configuration, the other case can be central configuration. In the second case the existence and uniqueness of the central configuration are studied. If the configuration is a central configuration, then all masses of outside layer are equivalent, similar to the masses of inside layer. At the same time the following relation between r(the ratio of the sizes) and mass ratio b = m/m must be satisfied in which the masses at outside layer are not less than the masses at inside layer, and the solution of this kind of central configuration is unique for the given ratio (6) of masses.
文摘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.
基金the National Natural Science Foundation ofChina under Grant No. 50779043, 50779045
文摘To simulate two-dimensional free-surface flows with complex boundaries directly and accurately, a novel VOF (Volume-of-fluid) method based on unstructured quadrilateral mesh is presented. Without introducing any complicated boundary treatment or artificial diffusion, this method treated curved boundaries directly by utilizing the inherent merit of unstructured mesh in fitting curves. The PLIC (Piecewise Linear Interface Calculation) method was adopted to obtain a second-order accurate linearized reconstruction approximation and the MLER (Modified Lagrangian-Eulerian Re-map) method was introduced to advect fluid volumes on unstructured mesh. Moreover, an analytical relation for the interface’s line constant vs. the volume clipped by the interface was developed so as to improve the method’s efficiency. To validate this method, a comprehensive series of large straining advection tests were performed. Numerical results provide convincing evidences for the method’s high volume conservative accuracy and second-order shape error convergence rate. Also, a dramatic improvement on computational accuracy over its unstructured triangular mesh counterpart is checked.
基金supported by the National Natural Science Foundation of China (10972006, 11172004)National Basic Research Program of China (2010CB832701)
文摘In this paper, a new method of topological cleanup for quadrilateral mesh is presented. The method first selects a patch of mesh around an irregular node. It then seeks the best connection of the selected patch according to its irregular valence using a new topological operation： small polygon reconnection （SPR）. By replacing the original patch with an optimal one that has less irregular valence, mesh quality can be improved. Three applications based on the proposed approach are enumerated： （1） improving the quality of a quadrilateral mesh, （2） converting a triangular mesh to a quadrilateral one, and （3） adapting a triangle generator to a quadrilateral one. The presented method is highly effective in all three applications.
基金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.
基金The research presented in this paper has been financially supported by the Croatian Science Foundation(Grants HRZZ-IP-11-2013-1631 and HRZZ-IP-2018-01-1732)Young Researchers'Career Development—Training of Doctoral Students,as well as a French Government Scholarship.
文摘Quadrilateral finite elements for linear micropolar continuum theory are developed using linked interpolation.In order to satisfy convergence criteria,the newly presented finite elements are modified using the Petrov-Galerkin method in which different interpolation is used for the test and trial functions.The elements are tested through four numerical examples consisting of a set of patch tests,a cantilever beam in pure bending and a stress concentration problem and compared with the analytical solution and quadrilateral micropolar finite elements with standard Lagrangian interpolation.In the higher-order patch test,the performance of the first-order element is significantly improved.However,since the problems analysed are already describable with quadratic polynomials,the enhancement due to linked interpolation for higher-order elements could not be highlighted.All the presented elements also faithfully reproduce the micropolar effects in the stress concentration analysis,but the enhancement here is negligible with respect to standard Lagrangian elements,since the higher-order polynomials in this example are not needed.
文摘Accuracy of a simulation strongly depends on the grid quality. Here, quality means orthogonality at the boundaries and quasi-orthogonality within the critical regions, smoothness, bounded aspect ratios and solution adaptive behaviour. It is not recommended to refine the parts of the domain where the solution shows little variation. It is desired to concentrate grid points and cells in the part of the domain where the solution shows strong gradients or variations. We present a simple, effective and com- putationally efficient approach for quadrilateral mesh adaptation. Several numerical examples are presented for supporting our claim.
文摘To provide a iheoreiical basis for h-type .finire elemenl analysis with quadrilateralelements, in the present paper, the h-convergence of quadrilateral elements is established. whose related lemmas and theorems being presented, and therefore, theerror estimate problems are investigated.
基金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.
基金supported by NSFC Nos.61907005,61720106005,61936002,62272080.
文摘We propose a newmethod to generate surface quadrilateralmesh by calculating a globally defined parameterization with feature constraints.In the field of quadrilateral generation with features,the cross field methods are wellknown because of their superior performance in feature preservation.The methods based on metrics are popular due to their sound theoretical basis,especially the Ricci flow algorithm.The cross field methods’major part,the Poisson equation,is challenging to solve in three dimensions directly.When it comes to cases with a large number of elements,the computational costs are expensive while the methods based on metrics are on the contrary.In addition,an appropriate initial value plays a positive role in the solution of the Poisson equation,and this initial value can be obtained from the Ricci flow algorithm.So we combine the methods based on metric with the cross field methods.We use the discrete dynamic Ricci flow algorithm to generate an initial value for the Poisson equation,which speeds up the solution of the equation and ensures the convergence of the computation.Numerical experiments show that our method is effective in generating a quadrilateral mesh for models with features,and the quality of the quadrilateral mesh is reliable.
基金supported by the Australian Research Council (ARC DP066620, LP0560932, and LX0989423)
文摘With the development of advanced imaging technology, digital images are widely used. This paper proposes an automatic quadrilateral mesh generation algorithm for multi-colour imaged structures. It takes an original arbitrary digital image as an input for automatic quadrilateral mesh generation, this includes removing the noise, extracting and smoothing the boundary geometries between different colours, and automatic all-quad mesh generation with the above boundaries as constraints. An application example is provided to demonstrate the usefulness and effectiveness of the proposed approach.
