Assessment of New Adaptive Finite Elements Based on Carrera Unified Formulation for Meshes with Arbitrary Polygons

The meshing technique represents the capability to discretize the domain of interest, to fit the real physical continuum in the best possible way. The most used approach is the finite-element method (FEM), a numerical method to solve partial differential equations. To overcome the classical issues presented by FEM, other models are investigated. The goal is to allow the problem domain to be discretized by elements represented by arbitrary polygons, which can be concave and convex. Moreover, different polynomial consistency is sought within these methods with the possibility to handle non-conforming discretizations, mainly for local refinement and so on. This work aims to present the new adaptive elements, which are finite elements based on Carrera unified formulation, to demonstrate that all the previous capabilities can be done with these new elements, with easy implementation of the relative model. First, a classical patch test is done to investigate the mesh distortion sensitivity. Then, different study cases are presented with more complex meshes combining very distorted concave and convex elements.


Introduction
The meshing technique represents the capability to discretize the domain of interest, to fit the real physical continuum in the best possible way.In the engineering field, meshing is a crucial milestone in creating the model of the problem to solve, it is the starting point on which the simulation is performed and it is important that it can represent in the more accurate way the analyzed domain.The physical area investigated is not so regular, in general, and the meshing should account for it.
The most used technique is the finite-element method (FEM), which relies on a quite simple mathematical model capable of describing the field under investigation (structural and/or fluid dynamic) up to a certain level of accuracy.This method allows the construction of 2D and 3D regular elements with a prescribed number of nodes each [1].The main issues of this method are represented by: the constraint to regular presetting of elements; the difficulties managing different types of elements on the same domain, perhaps with a different number of nodes, therefore, a different approximation; the high computational cost due to the high number of elements used, to adapt when the real domain it is not geometrically regular or when enrichment of the mesh is required, because it is in the area of interest.
To alleviate meshing issues due to FEM, more advanced techniques can be used, resorting to numerical methods that are designed from the very beginning to provide arbitrary order of accuracy on more generally shaped elements.These techniques are based, for instance, on the virtual element method (VEM) or polytopal element method (PEM) [2].In fact, in contrast to FEM in which elements are typically triangular and quadrilateral in 2D or tetrahedral and hexahedral in 3D, the VEM permits arbitrary two-dimensional polygonal and three-dimensional polyhedral elements [3,4], having many aspects closely connected with FVM and in particular with mimetic finite differences (MFD) [5].This allows the problem domain to be discretized by elements represented by arbitrary polygons, which can be concave and convex.Moreover, different polynomial consistency is allowed within the method and non-conforming discretizations can be handled, mainly for local refinement and so on, representing the key aspect of this method [6].
The newly proposed adaptive elements, which are finite elements based on Carrera unified formulation (CUF), represent more agile and manageable elements, based on a well-known technique (FEM), easy to implement and can produce the same benefits obtained using VEM.As a matter of fact, the CUF, applied both to 1D and 2D finite elements, provides very accurate 3D-like solutions with much lower computational cost than the classical 3D elements, becoming an important tool for the analysis of structural problems.
Recently, in the framework of CUF, a novel approach called node-dependent kinematics (NDK) has been proposed to further increase the numerical efficiency of the models [7][8][9][10][11].By relating the kinematic assumptions to the chosen nodes, FE models with variable nodal kinematics can be built conveniently.Naturally, the critical zone with higher order assumptions can be bridged to the less-critical area modeled with adequate lower order theories [12], in such a way the numerical modeling procedure can further be simplified.
The NDK approach, along with the 3D modelling of non-orthogonal geometries, opens the possibility to create elements with a non-conventional number of nodes, with non-conventional geometrical shapes, and with different polynomial approximations along different directions, leading to a certain freedom in constructing proper elements.These elements are formulated considering that the integrals inside the matrices of the governing equations (stiffness and mass matrices), which involves the approximating functions, and the Jacobian matrix of the elements are computed in 3D form, although 1D and 2D CUF models are still employed.In conclusion, 3D elements can be recovered with a non-regular shape that can adapt to the edges of the domain considered or to adjacent portions which share one or more sides, when the geometrical constraints or the kinematic behaviour require that.
The paper is organized as follows: an introduction is presented in Sect.1; the 1D and 2D models based on CUF are briefly recalled in Sect.2; the derivation of governing equations from the Principle of Virtual Displacements by using both the finite-element method and the CUF is presented in Sect.3, where also the present approach is explained in the framework of NDK technique along with the 3D modelling; Sect. 4 contains the results obtained by the static and modal analysis of some useful study cases that demonstrate the capabilities of the present elements; finally, Sect. 5 summarises the conclusions of this work.

Carrera Unified Formulation: 1D and 2D Models
The CUF (see [13,14]) describes the kinematic field in a unified manner that will be then exploited to derive the governing equations in a compact way, which are easy to implement unlike the VEM.This compact form represents an huge advantage with respect to classical theories, that have several limitations (as explained in many previous works related to CUF).The main of those limitations is that the kinematics should be theoretically enriched with an infinite number of terms (see Washizu [15]) to overcome physical inconsistencies and to predict high order effects.Studying real problems, this is not feasible.There is the necessity of truncating the expansion of the primary mechanical variables, along the smallest dimensions of the structure domain, to a given order, which is generally done when formulating structure's theories.Furthermore, the more the terms in the kinematics, the more the complexity of the formulation and resolution of the problem is.The CUF represents a solution for these issues, describing the governing equations in a compact and unified form, allowing the formulation of 1D and 2D models (beam and plate models, respectively).Let us consider a generic beam structure whose longitudinal axis, with respect to a Cartesian coordinate system, lays on the coordinate y, being its cross section defined in the xz-plane.The displacement field of one-dimensional models in CUF framework is described as a generic expansion of the generalized displacements (in the case of displacementbased theories) by arbitrary functions of the cross-sectional coordinates: where u = {u x , u y , u z } is the vector of 3D displacements and u = {u x , u y , u z } is the vector of general displacements, M is the number of terms in the expansion, denotes summation and the functions F (x, z) define the 1D model to be used.In fact, depending on the choice of F (x, z) func- tions, different classes of beam theories can be implemented, such as Euler-Bernoulli or Timoshenko beam models and so forth.
In the framework of plate theories, by considering the mid-plane of the plate laying in the xy-plane, CUF can be formulated in an analogous manner (see [16]): In the equation above, the generalized displacements are function of the mid-plane coordinates of the plate and the expansion is conducted in the thickness direction z.Again, depending on the choice of the F expanding functions, (1) various theories can be formulated, such as Kirchoff-Love or Reissner-Mindlin plate models and so forth.In this paper, only Lagrange expansions are considered, but other types of expansions based on the use of Taylor or Legendre polynomials can be conveniently employed as discussed in many previous works by Carrera et al. [8,10,14].

Lagrange Expansions
Lagrange expansions theories are based on the use Lagrange-type polynomials as generic functions over the cross section or along the thickness.
In the case of beam models, the cross section is divided into a number of local expansion sub-domains, whose polynomial degree depends on the type of Lagrange expansion employed.Four-nodes bilinear Q4, nine-nodes cubic Q9 and sixteen-nodes quartic Q16 polynomials can be used to formulate refined beam theories (see Carrera and Petrolo [17]).For example, the interpolation functions of a Q9 expansion are defined as where and vary over the cross section between −1 and +1 , and and represent the locations of the roots in the natural plane.The kinematic field of the single-Q9 beam theory is, therefore Refined beam models can be obtained by adopting higher order Lagrange polynomials or using a combination of Lagrange polynomials on multi-domain cross sections.More details about Lagrange-class beam models can be found in [17][18][19][20].
Plate theories based on Lagrange expansion are formulated by adopting 1D Lagrange polynomials as F thickness functions up to the desired order p: (3) where i denotes a summation and ranges from 2 to p + 1 .In this way two-nodes linear B2, three-nodes parabolic B3 and four-nodes cubic B4 polynomials can be used.Within the framework of 2D CUF theories, displacement fields like the one above have been used for the implementation of layerwise models for laminated structures (see, for example, [14,21,22]).
Lagrange polynomials have been extensively employed in the formulation of variable kinematics plate and shell theories in a unified framework by Kulikov and his co-workers.The readers are referred to the original papers for more details about 2D models based on Lagrange-type expansions, see, for example, [23,24].

Finite-Element Formulation
The main advantage of CUF is that it allows to write the governing equations and the related finite-element arrays in a compact and unified manner, which is formally an invariant with respect to the F functions.In the sections below, the math- ematical derivation of the fundamental nucleus (the invariant) of the stiffness matrix in the case of CUF 1D and 2D models is provided in detail.

Geometrical Relations and Constitutive Equations
In this section, the same notation and reference system as defined in Sect. 2 are adopted.Let the 3D displacement vector be defined as According to classical elasticity, stress and strain tensors can be organized in six-term vectors with no lack of generality.They read, respectively: Regarding to this expression, the geometrical relations between strains and displacements with the compact vectorial notation can be defined as ( 5) ( where, in the case of small deformations and angles of rotations, D is the following linear differential operator: On the other hand, for isotropic materials the relation between stresses and strains is obtained through the wellknown Hooke's law: where C is the isotropic stiffness matrix The coefficients of the stiffness matrix depend only on the Young's modulus, E, and the Poisson ratio, , and in the case of isotropic materials, they are:

Approximation of Displacements
In the case of 1D models, the discretization along the longitudinal axis of the beam is made by means of the finiteelement method.The generalized displacements are in this way described as functions of the unknown nodal vector, q i , and the 1D shape functions, N i ( ): where n elem is the number of nodes per element and the unknown nodal vector is defined as Similarly, the FEM discretization of generalized displacements on the midsurface of the plate can be written as follows: where 2D shape functions N i ( , ) are employed.Different sets of polynomials can be used to define FEM elements.Lagrange interpolating polynomials have been chosen in this work to generate both one-dimensional and two-dimensional elements.For the sake of brevity, their expression is not provided, but it can be found in the book by Carrera et al. [13], in which two-nodes (B2), three-nodes (B3) and four-nodes (B4) beam elements and four-nodes (Q4), nine-nodes (Q9) and 16-nodes (Q16) plate elements are described.
By combining the FEM approximation with the kinematic assumptions of the CUF, the 3D displacement field can be written as where the functions F and N i are defined according to the type of element (beam or plate).

Node-Dependent Kinematic Approach
In Eq. 16, the shape functions N i and the expanding func- tions F are independent.Recently, Carrera et al. [13, 25,  26] introduced a coupling by relating the expanding functions F to the shape functions N i by means of the following formalism: The difference of Eq. 17 from Eq. 16 is the additional superscript i of N i , which is now an index also of the function F .This definition introduces the dependency of the kinematic assumptions to the FE nodes, namely, the NDK.
As mentioned before, expanding functions with an increased order provide the chance to the better approximation of the structure responses.With NDK, the kinematic models can be refined locally on specific nodes which makes it easy to perform a local adaptable kinematic refinement with convenience.Different theories of structures will be blended naturally by the nodal shape functions within the element domain without any special coupling approaches.Meanwhile, no compatibility requirements for different nodal kinematics are needed.
The NDK technique can be applied to the efficient global-local modeling of structures.The kinematic model in the critical zone can be refined until the ideal accuracy is achieved while leaving the outlying area modeled with adequate lower order theories.Global-local FE model can be constructed conveniently without changing the meshes, and the same set of mesh grids can be re-used to build a family of models for concurrent global-local analyzes.This approach has been used in the efficient modeling of laminated structures in both 1D [25,26] and 2D [13,17,27] cases.

3D Modelling for Non-orthogonal Geometries
Classical beam and plate theories cannot be clearly employed in the modelling of beam with variable section or section not orthogonal to the axis, as well as plate with variable thickness or edges not orthogonal to the midsurface.In the framework of NDK approach, it is possible to extend the models of CUF to the modelling of this kind of geometries by incorporating the CUF kinematic assumption and the FEM discretization in a unique 3D approximation, as follows: where , , are the natural coordinates corresponding to x, y, z.In this expression, L i = (F i N i ) represents a non-con- ventional 3D shape function in which the order of expansion can be different along one of the spatial directions.Similarly, the virtual variation of displacements, that will be used in the derivation of governing equations below, can be approximated by where new summation indexes s and j have been introduced.
Considering this formalism, the volume integrals involved in the governing equations will be not split in 1D and 2D integrals as usual, but the functions (N i F ) and (N j F s ) will be handled as regular 3D shape functions; con- sequently, the Jacobian matrix relative to the transformation from natural coordinates , , to global coordinates x, y, z will be computed in 3D form.
The results will demonstrate that this particular approach allow us to save degrees of freedom with respect to the use of meshes based on classical 3D finite elements.Indeed, conventional 3D elements employ the same order of expansion in the three spatial directions and this usually implies some limitations on the choice of the aspect ratio of the element by leading to a detrimental increase of the number of elements used.

Governing Equations
The governing equations are obtained via the principle of virtual displacements.This variational statement sets as a necessary condition for the equilibrium of a structure that u = (F j s N j ) q sj = L sj ( , , ) q sj , the virtual variation of the internal work has to be equal to the virtual variation of external work, or The internal work is equivalent to the elastic strain energy: where V stands for the volume of the domain.By adopting the geometrical relation [Eq.( 8)], the constitutive law [Eq.(10)], the CUF kinematic field [Eq.( 1)] and the FEM discretization [Eq.( 13)], the internal work can be rewritten as where K sij is the 3 × 3 fundamental nucleus of the stiffness matrix.
The inertial work can be written as By adopting the CUF kinematic field and the FEM discretization (Eqs.( 18) and ( 19)), it becomes: where M sij is the 3 × 3 fundamental nucleus of the mass matrix.
The external work can be written as Following the same steps done for the internal work, it is possible to write also the external work in terms of nodal vector q i as where P i is the load vector in terms of shape functions N i and expanding functions F .
It should be noted that the formal expressions of the components of the fundamental nuclei are independent on the choice of the expanding functions F , which determine the theory of structure, and shape functions N i , which determine the numerical accuracy of the FEM approximation.This means that any beam or plate element can be automatically formulated by appropriately expanding the fundamental nuclei according to the indexes , s, i, and j.For more details about the expressions of fundamental nuclei of stiffness, mass and load matrices, the reader can refer to the work [28].
In this paper, linear, quadratic and cubic Lagrange-type shape functions are used for both F and N i functions. (20) Note that the acronyms B2 (linear), B3 (quadratic) and B4 (cubic) are used both for the 1D discretization of the axis in beam models and for the approximation through the thickness in plate models; as well as the acronyms Q4 (linear), Q9 (quadratic) and Q16 (cubic) are used both for the 2D approximation of the cross section in beam models and for the discretization of the midsurface in plate models.
According to these characteristics, the expansion order is set as a free input of the model that determines the number of unknowns to be solved.At the same time, the kinematic field can be enriched in particular zones of interest through the definition of multi-domain cross sections.

Numerical Results
In this section, some benchmark problems are considered for assessing the robustness of the proposed Adaptive Elements based on CUF.First, a set of problems involving a homogeneous and isotropic plate is used for assessing the element's performance, shown in Sect.4.1: the displacement sensitivity concerning distorted element shapes and the free-vibration analysis are considered.The influence of the boundary conditions is also addressed.Then, static and modal analyses are performed on a plate with distorted elements.The distorted mesh was created in a random way starting from the regular pattern: the new positions of nodes are obtained using a function which generates a random scalar drawn from the normal or Gaussian distribution having a mean equal to the previous position and different standard deviations.
For all the listed problems, both four-node and nine-node elements are considered.

Patch Test
This standard test is classically used to investigate the mesh distortion sensitivity in plate bending problems.The in-plane dimensions of the square plate considered are 100 × 100 m 2 , while the thickness is 5 m .The boundary conditions applied are: clamped-clamped (CC), all edges are fixed; simply supported (SS), all edges have no displacement in the thickness direction and the direction parallel to the edge.The load is expressed as a concentrated force at the centre of the top surface equal to F = 4 N .The material considered is isotropic with the following properties: E = 10.92Pa , = 0.30 and = 1 kg/m 3 .The starting meshes are a 4 × 4 Q4 elements and a 4 × 4 Q9 elements in the plane, with B3 approximation along the thickness in both cases.
Due to the symmetry of the problem, with symmetric geometry and symmetric boundary conditions, a quarter of the plate is considered to simplify the analysis: the full plate is divided into four equal sub-plates and the study is conducted only on one of them, as shown in Fig. 1.In particular, the bottom-left quarter is chosen and the following boundary conditions are applied, considering as u, v and w the displacements along x, y and z, respectively: bottom and left edges are fixed or simply-supported, on the top edge v = 0 , on the right edge u = 0 and on the top-right corner a force of F = 1 N is applied.

Mesh Distortion Sensitivity
The resulting mesh of the plate quarter is 2 × 2 based on Q4 or Q9 elements with B3 discretization along the thickness, see Fig. 2.
Note that from now on, some beam elements are introduced in the mesh to obtain distorted elements with different numbers of nodes along the edges.Then, three different configurations are analyzed: • 2xQ4 1xB2: in this configuration, the previous two elements at the top are substituted by a single 1D element with 2 nodes, having a Q9 approximation on the cross section; • 2xQ4 1xB3: as in the previous case, the two elements at the top are substituted by a single 1D element but having 3 nodes along the axis, where the cross section is expressed as a Q9 element; • 2xQ4 1xQ9: the number of nodes is identical with respect to the previous case, but the discretization of the two elements at the top is different, accounting for a single Q9 element having a B3 discretization along the thickness.
In Figs. 3 and 4, it is possible to visualize the described configurations: the black dots are the nodes used during the implementation, while the white dots represent the nodes added considering the thickness and the cross section of the plate and the beam elements, respectively.A color code is introduced for the different elements: the B2 element is colored heavenly, the B3 element is colored yellow, the Q4 element is colored light blue and the Q9 element is colored orange.
The emphasis now will be focused on the beam elements and on the possibility to freely move its nodes thanks to the 3D integration of its approximation function, as discussed in Sect.3.4.
The mesh distortion is characterized by the parameter s ∈ {−12, −10, −6, 0, 6, 10, 12} , which defines the coordi- nates of the central node of the plate quarter.Note that the parameter s can have negative and positive values both: this is not common in literature and, particularly in this work, it is used to control the property of the element: in

Static and Modal Analysis
The three different configurations are used to perform a static and modal analysis, as done in Ref. [29].Regarding the static analysis, the results are computed in terms of transverse displacement U z = u z (a, a, 0) (the center point on the midsurface of the entire plate) for the distorted mesh and it is normalized with respect to the reference value U 0 z calculated for the regular mesh ( s = 0 ).The values of U 0 z for the different meshes are compared with the 3D solution computed in Nastran for validation purposes.One can note that the 2x2Q9 solution is the most accurate with the highest number of degrees of freedom, the 2x2Q4 is the least accurate and the other solutions are in the middle according to their related degree of approximation.The meshes 2xQ4 1xQ9 and 2xQ4 1xB3 are very similar and they are compared to highlight that their solution are exactly the same, although the first uses a plate element (Q9) and the other one a beam element (B3).Moreover, one can see that the 2xQ4 1xB2 solution is more accurate than the 2x2Q4 one, although the number of degrees of freedom is the same: this is due to the use of quadratic approximation over the cross section of the B2 beam element.We can conclude that the present adaptive elements permit to increase the accuracy with respect to regular meshes with the same degrees of freedom (Table 1).
In addition, the results computed by varying the parameter s are reported in Figs. 6 and 7, where also different boundary conditions are considered (CC and SS, respectively).In both cases, these graphics demonstrate that the sensitivity to the distortion of the mixed meshes with beam elements is comparable to the regular meshes (2x2Q4 and 2x2Q9).The behaviors of 2xQ4 1xQ9 and 2xQ4 1xB3 meshes are again very similar; moreover, one can note that the use of B2 beam element in the mesh 2xQ4 1xB2 does not compromise the accuracy of the results when the mesh is distorted, although it permits to differentiate the number of nodes along the edges of the element (see Fig. 3b) by reducing the total number of the nodes in the element.
Finally, a free-vibration analysis is performed considering the different values of the parameter s.The results are shown in Figs. 8 and 9 in terms of natural frequencies associated with the first 16 modes: in this analysis, no boundary conditions are applied, hence only the non-null modes are reported.It is evident that the conclusions drawn before are confirmed: the solutions of the mixed meshes lie between the regular meshes, 2x2Q4 and 2x2Q9; the mixed meshes 2xQ4 1xQ9 and 2xQ4 1xB3 provide the same results, except in the case of s = −12 for SS boundary conditions; the solu- tion 2xQ4 1xB2 provide lower natural frequencies than the 2x2Q4 regular mesh, so one can suppose it is more accurate.Starting from these regular meshes, the elements are gradually distorted to demonstrate the capabilities of the presented adaptive elements.The nodes coordinates of some elements are varied in a random way, considering a Gaussian distribution having the mean equal to the original coordinate and the standard deviation increasing little by little.Stating that x 0 is the original coordinate of a node, dev the standard deviation and the mean value, the new coordinates are obtained as follow: where the Matlab Ⓡ function randn returns pseudorandom values drawn from the standard normal distribution, = 1 and the standard deviation varies in a range equal to dev ∈ {1, 2, 3, 5}.Moreover, in other examples, the nodes are moved intentionally to obtain specific shapes (for instance, a pentagon, which is a benchmark example for the virtual elements [3,5,30]).As a final example, a complete nonregular distorted mesh is considered.For all of these cases, a combination of plate and beam elements is exploited.
The meshes are reported in Fig. 11, where: • Case A: in Fig. 11a Mesh 9xQ9 1xQ4 6xB2 (DOFs = 387)-the mesh is obtained using a standard deviation equal to dev = 3; • Case B: in Fig. 11b Mesh 3xQ9 3xQ4 10xB2 (DOFs = 330)-the mesh is obtained creating some pentagon elements; • Case C: in Fig. 11c Mesh 3xQ9 3xQ4 10xB2 v2 (DOFs = 330)-an upgraded version of the previous mesh is obtained to highlight the presence of pentagon elements, leading to a more distorted mesh with the same number of DOFs; • Case D: in Fig. 11d Mesh 9xQ9 8xB3 (DOFs = 765)-a complete distorted mesh is obtained and modelled using high order accuracy in the approximation of the kinematics (quadratic interpolation); • Case E: in Fig. 11e Mesh 10xQ4 7xB2 (DOFs = 315)-a similar mesh to the previous one is considered, using a lower order accuracy of the kinematics approximation (linear interpolation), to investigate if it is possible to obtain comparable results with fewer degrees of freedom.

Static and Modal Analysis
The different meshes are investigated by performing static analysis and the results are reported in terms of displacements along the bottom free edge.In particular, in Fig. 12, the displacement magnitude is shown for the different meshes introduced above.Each mesh is analyzed and compared with the family of curves related to the Q4 and Q9 regular meshes.In particular, we consider: case A in which there is almost a balance between elements with high order accuracy (Q9) and elements with a lower order of accuracy (Q4, B2); case B in which the predominant elements are similar to the Q4 regular mesh; case C that is similar to the previous case, with the difference that now the pentagons are immersed in a more distorted mesh; cases D and E that, due to the severe distortion of the elements, tend to match the reference solutions related to Q9 distorted mesh with dev = 3 for the first case and related to Q4 distorted mesh with dev = 5 for the second case.In conclusion, it is evident that the results of the distorted meshes lie between the two reference curves of Q4 and Q9 regular meshes, with Fig. 12 Square plate-static analysis increasing accuracy according to the related number degrees of freedom.In a similar way to the static analysis, also for the modal investigation the frequencies are compared with respect to the solutions of the Q4 and Q9 regular meshes, as shown in Fig. 13.Note that the frequencies of the distorted meshes are perfectly bounded by those of the regular meshes, in particular the case D is very close to the full Q9 solution.

Conclusions
In this paper, the assessment of adaptive elements based on the CUF with respect to distorted mesh has been proposed.In the framework of CUF, which is an advanced tool that enables to formulate theories of structure in a unified manner, the NDK approach was exploited along with a 3D modelling approach.A kinematic refinement on some nodes is possible by allowing to choose different expansions through the same element.This aspect represents a key point to reduce the number of nodes inside the element and, consequentially, the computational cost by keeping the same accuracy.
When dealing with meshing technique, using the FEM, the main issue is the difficulty of approximate the physical domain with elements which must be triangular or quadrilateral in 2D and tetrahedral or hexahedral in 3D trying to contain, at the same time, the number of elements for computational cost reasons.The proposed adaptive elements can overcome these problems, offering a convenient solution.
The capabilities of these elements are already investigated in other works, see [28,31].This work aims to demonstrate that different polygons and polyhedra can be constructed using the adaptive elements, with a relatively low number of degrees of freedom, and to assess their performances when severe mesh distortions arise.Static and modal analyses are conducted, giving results in terms of displacements and natural frequencies.The outcomes are presented for distinct distorted meshes, combining elements with different kinematic approximations.The standard tests and the example problems analyzed prove that convenient distorted mesh (from a computational point of view) can be accomplished and the effect of distortion is comparable to classical elements.
Further studies shall concern the application of the adaptive elements to the analysis of shell structures, relying on the formulation of these elements in curvilinear coordinates, according to the recent author's work [32].
Funding Open access funding provided by Politecnico di Bari within the CRUI-CARE Agreement.

Fig. 13 Square plate-modal analysis
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http:// creat iveco mmons.org/ licen ses/ by/4.0/.
xz yz xy T = yy xx zz xz yz xy .

Fig. 1
Fig. 1 Isotropic plate with transversal concentrated load in the center point.Only a quarter of the plate is meshed with 2 × 2 elements

Fig. 5
Fig.5Quarter plate-mesh distortion defined by parameter s

Fig. 7 3 Fig. 8 Fig. 9
Fig. 7 Variation of the ratio U z ∕U 0 z for SS quarter plate

Table 1
Tip displacement Variation of the ratio U z ∕U 0 z for CC quarter plate