Stretched grid method


The stretched grid method is a numerical technique for finding approximate solutions of various mathematical and engineering problems that can be related to an elastic grid behavior. In particular, meteorologists use the stretched grid method for weather prediction and engineers use the stretched grid method to design tents and other tensile structures.

FEM and BEM mesh refinement

In recent decades, the finite element and boundary element methods have become mainstays for industrial engineering design and analysis. Increasingly larger and more complex designs are being simulated using the FEM or BEM. However, some problems of FEM and BEM engineering analysis are still on the cutting edge. The first problem is that the reliability of the analysis strongly depends upon the quality of initial data generated at the pre-processing stage. It is known that automatic element mesh-generation techniques at this stage have become commonly used tools for the analysis of complex real-world models. With FEM and BEM increasing in popularity comes the incentive to improve automatic meshing algorithms. However, all of these algorithms can create distorted and even unusable grid elements. Several techniques exist that can take an existing mesh and improve its quality. For instance smoothing is one such method, which repositions nodal locations, so as to minimize element distortion. The stretched grid method allows one to obtain pseudo-regular meshes very easily and quickly in a one-step solution.
Assume that there is an arbitrary triangle grid embedded into a plane polygonal single-coherent contour and produced by an automeshing procedure. It may be assumed further that the grid considered as a physical nodal system is distorted by a number of distortions. It is supposed that the total potential energy of this system is proportional to the length of some -dimensional vector with all network segments as its components.
Thus, the potential energy takes the form
where is total number of segments in the network, is the length of segment number, and is an arbitrary constant.
The length of segment number may be expressed by two nodal coordinates as
It may also be supposed that coordinate vector of all nodes is associated with the non-distorted network, and the coordinate vector is associated with the distorted network. Vector may the be expressed as
The vector determination is related to minimization of the quadratic form by incremental vector ; that is,
where is the number of interior node of the area, and is the number of coordinate
After all transformations, we may write the following two independent systems of linear algebraic equations:
where is a symmetric matrix in the banded form similar to the global stiffness matrix of a FEM assemblage, and are incremental vectors of coordinates of all nodes at axes 1 and 2, and and are the right part vectors that are combined by coordinates of all nodes in axes 1 and 2.
The solution of both systems, keeping all boundary nodes conservative, obtains new interior node positions corresponding to a non-distorted mesh with pseudo-regular elements. For example, Fig. 2 presents the rectangular area covered by a triangular mesh. The initial auto mesh possesses some degenerative triangles. The final mesh produced by the SGM procedure is pseudo-regular without any distorted elements.
As the above systems are linear, the procedure reduces very quickly to a one-step solution. Moreover, each final interior node position meets the requirement of coordinate arithmetic mean of nodes surrounding it and meets the Delaunay criteria too. Therefore, the SGM has all the positive values peculiar to Laplacian smoothing and other kinds of smoothing approaches but is much easier and more reliable because of integer-valued final matrix representations. Finally, the SGM is applicable not only to 2D meshes but to 3D meshes consisting of any uniform cells as well as to mixed or transient meshes.

Minimum surface problem solution

Mathematically, a surface bounded by a closed curve is called minimal if its area is minimal amongst all the surfaces with the same boundary. The best-known minimal surface sample is a soap film bounded by a wire frame. Usually to create a minimum surface, a fictitious constitutive law, which maintains a constant prestress, independent of any changes in strain, is used. An alternative approximate approach is based on SGM.
The idea is to approximate a surface by a triangular grid. To converge such a grid to a grid with minimum area, one should solve the same two systems described above. Increments of the third nodal coordinates may be determined by a similar system at axis 3:
By solving all three systems simultaneously, one can obtain a new grid that will approximate desired minimal surface.
As an example, a catenoid calculated by the above approach is presented in Fig 3. The radii of the bounding rings and the height of the catenoid are equal to 1. The surface area of the catenoid determined by SGM is 2.9967189; the area of the ideal catenoid is 2.992.

Tensile fabric structures form finding

For structural analysis, the configuration of the structure is generally known beforehand. This is not the case for tensile structures such as tension fabric structures. Since the membrane in a tension structure possesses no flexural stiffness, its form or configuration depends upon initial prestressing and the loads that it is subjected to. Thus, the load-bearing behaviour and the shape of the membrane cannot be separated and cannot be generally described by simple geometric models only. The membrane shape, the loads on the structure, and the internal stresses interact in a nonlinear manner to satisfy the equilibrium equations.

The preliminary design of tension structures involves the determination of an initial configuration, referred to as form finding. In addition to satisfying the equilibrium conditions, the initial configuration must accommodate both architectural and structural requirements. Further, the requirements of space and clearance should be met, the membrane principal stresses must be tensile to avoid wrinkling, and the radii of the double-curved surface should be small enough to resist out-of-plane loads and to insure structural stability. Several variations on form-finding approaches based on FEM have been developed to assist engineers in the design of tension fabric structures. All of them are based on the same assumption as that used for analysing the behaviour of tension structures under various loads. However, it might sometimes be preferable to use minimal surfaces in the design of tension structures.
The physical meaning of SGM consists in the convergence of the energy of an arbitrary grid structure with specified constraints to a grid that minimizes the sum of the distances between arbitrary pairs of grid nodes. It allows the minimum surface energy problem to be solved. The generalized formulation of SGM presupposes the ability to apply a set of external forces and rigid or elastic constraints to grid structure nodes, in order to model various external effects. We may obtain the following expression for such a formulation:
where
  • is the total number of grid segments,
  • is the total number of nodes,
  • is the length of segment,
  • is the stiffness of segment,
  • is the coordinate increment of node at axis,
  • is the stiffness of an elastic constraint in node at axis, and
  • is the external force on node at axis.

    Unfolding problem and cutting pattern generation

Once a satisfactory shape has been found, a cutting pattern may be generated. Tension structures are highly varied in their size, curvature and material stiffness. Cutting pattern approximation is strongly related to each of these factors. It is essential for a cutting-pattern generation method to minimize the approximation error and to produce reliable plane cloth data.
The objective is to develop the shapes described by these data, as close as possible to the ideal doubly curved strips. In general, cutting-pattern generation involves two steps. First, the global surface of a tension structure is divided into individual cloths. The corresponding cutting pattern at the second step can be found by simply taking each cloth strip and unfolding it onto a planar area. In the case of the ideal doubly curved membrane surface, the subsurface cannot be simply unfolded and must be flattened. For example, in, SGM has been used for the flattening problem solution.
The cutting-pattern generation problem is actually subdivided into two independent formulations. These are the generation of a distortion-free plane form unfolding each cloth strip and flattening double-curved surfaces that cannot be simply unfolded. From the position of differential geometry, both formulations are the same. We may consider it as an isometric mapping of a surface onto the plane area that will be conformal and equiareal simultaneously because of invariant angles between any curves and invariance of any pieces of area. In the case of a single-curved surface that can be unfolded precisely, equiareal mapping allows one to obtain a cutting pattern for fabric structure without any distortions. The second type of surfaces can be equiareally mapped only approximately, with some distortions of linear surface elements limited by the fabric properties. Assume that two surfaces are parameterized so that their first quadratic forms may be written as follows:
The condition of conformal mapping for two surfaces requires that
where is the ratio of the surface distortion due to conformal mapping.
It is known that the first quadratic form reflects the distance between two surface points and. When is close to 1, the above equation converges to the conditions of isometric mapping and equiareal mapping respectively because of invariant angles between any curves and invariance of any pieces of area. Remembering that the first stage of form-finding is based on a triangular mesh, and using the for the description of isometric and equiareal mapping of the minimum surface onto a plane area, we may write the following function, which is defined by the sum of integrals along segments of curved triangles:
where
  • is the total number of grid cells,
  • are the weight ratios,
  • is the total mapping residual, and
  • is a constant that does not influence the final result and may be used as a scale ratio.
Considering further weight ratios, we may transform the above equation into an approximate finite sum that is a combination of linear distances between nodes of the surface grid and write the basic condition of equiareal surface mapping as a minimum of following nonlinear function:
where
  • is the initial length of linear segment,
  • is the final length of segment number, and
  • is a distortion ratio close to 1 and may be different for each segment.
The initial and final lengths of segment number may be expressed as usual by two nodal coordinates as
where
  • are coordinates of nodes of the initial segment, and
  • are coordinates of nodes of the final segment.
According to the initial assumption, we can write for the plane surface mapping. The expression for vectors and may be written as

The definition of vector is as previously
After transformations, we may write the following two independent systems of nonlinear algebraic equations:
where all the parts of the system can be expressed as previously and and are vectors of pseudo-stresses at axes 1 and 2 that have the form
where
  • is the total number of nodes that surround node number, and
  • is the number of global axes.
The above approach is another form of SGM and allows one to obtain two independent systems of nonlinear algebraic equations that can be solved by any standard iteration procedure. The less the Gaussian curvature of the surface, the higher the accuracy of the plane mapping. As a rule, the plane mapping allows one to obtain a pattern with linear dimensions 1–2% less than corresponding spatial lines of a final surface. That is why it is necessary to provide the appropriate margins while patterning.
The typical sample of cutout—also called a gore or a patch—is presented in Figs. 9, 10, and 11.