HOME

TheInfoList



OR:

The stretched grid method (SGM) 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 structure A tensile structure is a construction of elements carrying only tension and no compression or bending. The term ''tensile'' should not be confused with tensegrity, which is a structural form with both tension and compression elements. Tensile ...
s.


FEM and BEM mesh refinement

In recent decades the
finite element The finite element method (FEM) is a popular method for numerically solving differential equations arising in engineering and mathematical modeling. Typical problem areas of interest include the traditional fields of structural analysis, heat ...
and
boundary element method The boundary element method (BEM) is a numerical computational method of solving linear partial differential equations which have been formulated as integral equations (i.e. in ''boundary integral'' form), including fluid mechanics, acoustics, ele ...
s (FEM and BEM) have become a mainstay 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 a reliability of engineering analysis that 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 which can take an existing mesh and improve its quality. For instance smoothing (also referred to as mesh refinement) is one such method, which repositions nodal locations, so as to minimize element distortion. The Stretched Grid Method (SGM) allows the obtaining of pseudo-regular meshes very easily and quickly in a one-step solution(see ). Let one assume that there is an arbitrary triangle grid embedded into plane polygonal single-coherent contour and produced by an automeshing procedure (see fig. 1) 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 \ n-dimensional vector with all network segments as its components. Thus, the potential energy takes the following form :\Pi = D\sum_^n ^2 where *\ n - total number of segments in the network, *\ R_j - The length of segment number \ j, *\ D - an arbitrary constant. The length of segment number \ j may be expressed by two nodal co-ordinates as :\ R = \sqrt It may also be supposed that co-ordinate vector \ of all nodes is associated with non-distorted network and co-ordinate vector \ is associated with the distorted network. The expression for vector \ may be written as :\ = \ + \ The vector \ determination is related to minimization of the quadratic form \ \Pi by incremental vector \, i.e. : \frac = 0 where *\ l - is the number of interior node of the area, *\ k - the number of co-ordinate After all transformations we may write the following two independent systems of linear algebraic equations : A\ = \ : A\ = \ where * A/math> - symmetrical matrix in the banded form similar to global stiffness matrix of FEM assemblage, *\ and \- incremental vectors of co-ordinates of all nodes at axes 1, 2, *\ and \ - the right part vectors that are combined by co-ordinates of all nodes in axes 1, 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 (left mesh). The final mesh (right mesh) produced by the SGM procedure is pseudo-regular without any distorted elements. As above systems are linear, the procedure elapses very quickly to a one-step solution. Moreover, each final interior node position meets the requirement of co-ordinate arithmetic mean of nodes surrounding it and meets the
Delaunay Delaunay is a French surname. Notable people with the surname include: People Arts * Catherine Delaunay (born 1969), French jazz clarinet player and composer * Charles Delaunay (1911–1988), French author and jazz expert * Joseph-Charles Dela ...
criteria too. Therefore, the SGM has all the positive values peculiar to Laplacian and other kinds of smoothing approaches but much easier and reliable because of integer-valued final matrices representation. Finally, the described above SGM is perfectly 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 the surface embedded into a non-plane closed curve is called minimum if its area is minimal amongst all the surfaces passing through this curve. The best-known minimum surface sample is a
soap film Soap films are thin layers of liquid (usually water-based) surrounded by air. For example, if two soap bubbles come into contact, they merge and a thin film is created in between. Thus, foams are composed of a network of films connected by Plate ...
bounded by 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. The alternative approximated approach to the minimum surface problem solution is based on SGM. This formulation allows one to minimize the surface embedded into non-plane and plane closed contours. The idea is to approximate a surface part embedded into 3D non-plane contour by an arbitrary triangle grid. To converge such triangle grid to grid with minimum area one should solve the same two systems described above. Increments of the third nodal co-ordinates may be determined additionally by similar system at axis 3 in the following way : A\ = \ Solving all three systems simultaneously one can obtain a new grid that will be the approximating minimal surface embedded into non-plane closed curve because of the minimum of the function \ \Pi where parameter \ j = 1,2,3. As an example the surface of
catenoid In geometry, a catenoid is a type of surface, arising by rotating a catenary curve about an axis (a surface of revolution). It is a minimal surface, meaning that it occupies the least area when bounded by a closed space. It was formally describe ...
which is calculated by the described above approach is presented in Fig 3. The radii of rings and the height of catenoid are equal to 1.0. The numerical area of catenoidal surface determined by SGM is equal to 2,9967189 (exact value is 2.992).


Tensile fabric structures form finding

For structural analysis, the configuration of the structure is generally known à priori. This is not the case for
tensile structure A tensile structure is a construction of elements carrying only tension and no compression or bending. The term ''tensile'' should not be confused with tensegrity, which is a structural form with both tension and compression elements. Tensile ...
s such as tension
fabric structure A fabric structure is a structure made of fabric, with or without a structural frame. The technology provides end users a variety of aesthetic free-form building designs. Custom-made structures are engineered and fabricated to meet worldwide stru ...
s. Since the membrane in a tension structure possesses no flexural stiffness, its form or configuration depends upon initial prestressing and the loads to which it is subjected. 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 non-linear 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 (aesthetics) and structural (strength and stability) 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 (work ). 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, as it is noted by some researchers it might sometimes be preferable to use the so-called ‘ minimal surfaces’ in the design of tension structures. The physical meaning of SGM consists in convergence of the energy of an arbitrary grid structure embedded into rigid (or elastic) 3D contour to minimum that is equivalent to minimum sum distances between arbitrary pairs of grid nodes. It allows the minimum surface energy problem solution substituting for finding grid structure sum energy minimum finding that provides much more plain final algebraic equation system than the usual FEM formulation. The generalized formulation of SGM presupposes a possibility to apply a set of outer forces and rigid or elastic constrains to grid structure nodes that allows the modelling of various outer effects. We may obtain the following expression for such SGM formulation :\Pi = \sum_^n D_j R_j^2 + \sum_^3 \left(\sum_^m C_ \Delta X_^2 - \sum_^m P_ \Delta X_\right) where *\ n - total number of grid segments, *\ m - total number of nodes, *\ R_j - length of segment number \ j, *\ D_j - stiffness of segment number \ j, *\ \Delta X_ - coordinate increment of node \ k at axis \ i, *\ C _ - stiffness of an elastic constrain in node \ k at axis \ i, *\ P _ - outer force in node \ k at axis \ i.


Unfolding problem and cutting pattern generation

Once a satisfactory shape has been found, a
cutting pattern Cutting is the separation or opening of a physical object, into two or more portions, through the application of an acutely directed force. Implements commonly used for cutting are the knife and saw, or in medicine and science the scalpel and mi ...
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 possible approximation 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 on a planar area. In the case of the ideal doubly curved membrane surface the subsurface cannot be simply unfolded and they must be flattened. For example, in,Popov, E.V
Cutting pattern generation for tent type structures represented by minimum surfaces
The Transactions of the Canadian Society for Mechanical Engineering, Univ. of Alberta, vol. 22, N 4A, 1999, pp. 369–377.
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. Studying the problem carefully one can notice that from the position of differential geometry both formulations are the same. We may consider it as an
isometric mapping In mathematics, an isometry (or congruence, or congruent transformation) is a distance-preserving transformation between metric spaces, usually assumed to be bijective. The word isometry is derived from the Ancient Greek: ἴσος ''isos'' me ...
of a surface onto the plane area that will be
conformal mapping In mathematics, a conformal map is a function that locally preserves angles, but not necessarily lengths. More formally, let U and V be open subsets of \mathbb^n. A function f:U\to V is called conformal (or angle-preserving) at a point u_0\i ...
and equiareal mapping simultaneously because of invariant angles between any curves and invariance of any pieces of area. In the case of single-curved surface that can be unfolded precisely equi-areal mapping allows one to obtain a cutting pattern for fabric structure without any distortions. The second type of surfaces can be equi-areal mapped only approximately with some distortions of linear surface elements limited by the fabric properties. Let us assume that two surfaces are parameterized so that their first quadratic forms may be written as follows :I_1 = E_1(u, v) \operatornameu^2 + 2F_1(u, v) \operatornameu \operatornamev + G_1(u, v) \operatornamev^2 :I_2 = E_2(u, v) \operatornameu^2 + 2F_2(u, v) \operatornameu \operatornamev + G_2(u, v) \operatornamev^2 The condition of
conformal mapping In mathematics, a conformal map is a function that locally preserves angles, but not necessarily lengths. More formally, let U and V be open subsets of \mathbb^n. A function f:U\to V is called conformal (or angle-preserving) at a point u_0\i ...
for two surfaces as is formulated in differential geometry requires that :\sqrt = \lambda \sqrt where \ \lambda 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 \ (u,v) and \ (u + \operatornameu, v + \operatornamev). When \ \lambda-ratio is close to 1 the above eqn converges to condition of isometric mapping and to equi-areal 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 triangular mesh of a surface and using th
method of weighted residuals
for the description of isometric and equi-areal 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 :\Pi = D\sum_^n \oint_ w_j \left(\lambda \sqrt - \sqrt\right)^2 \operatornames where *\ n - total number of grid cells, *\ w_j - weight ratios, *\ \Pi - the total mapping residual, *\ D - the constant that does not influence the final result and may be used as a scale ratio. Considering further weight ratios \ w_j = 1 we may transform eqn. into approximate finite sum that is a combination of linear distances between nodes of the surface grid and write the basic condition of equi-areal surface mapping as a minimum of following non-linear function :\Pi = D\sum_^n \oint_ w_j \left(\lambda R_j - L_j\right)^2 \operatornames where *\ R_j - initial length of linear segment number \ j, *\ L_j - final length of segment number \ j, *\ \lambda - distortion ratio close to 1 and may be different for each segment. The initial and final lengths of segment number \ j may be expressed as usual by two nodal co-ordinates as : R = \sqrt : L = \sqrt where *\ X_ - co-ordinates of nodes of the initial segment, *\ x_ - co-ordinates of nodes of the final segment. According to the initial assumption we can write \ x_ = x_ = 0 for the plane surface mapping. The expression for vectors \ and \ with co-ordinate increments term use may be written as :\ = \ + \ The vector \ definition is made as previously : \frac = 0 After transformations we may write the following two independent systems of non-linear algebraic equations : A\ = \ + \ : A\ = \ + \ where all the parts of the system can be expressed as previously and \ and \ are vectors of pseudo-stresses at axes 1, 2 that has the following form :\ = -\left\ where *\ N - total number of nodes that surround node number \ t, *\ l - the number of global axes. The above approach is another form of SGM and allows the obtaining of two independent systems of non-linear algebraic equations that can be solved by any standard iteration procedure. The less Gaussian curvature of the surface is the higher the accuracy of the plane mapping is. As a rule, the plane mapping allows 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 cut out — also called a cutout, a
gore (segment) A gore is a sector of a curved surface or the curved surface that lies between two close lines of longitude on a globe and may be flattened to a plane surface with little distortion. The term has been extended to include similarly shaped piec ...
, or a patch — is presented in Figs. 9, 10, 11.


See also

* Mesh generation * Types of mesh


References


External links


K3-Tent system for tensile fabric structures formfinding and cutting patterning

Kubantent corp
Mesh generation Numerical differential equations {{Mesh generation, state=autocollapse