When designing a tension structure the shape is not known at the beginning of the process. Form-finding methods allow the designer to obtain an initial shape from given boundary conditions. Several form-finding methods for tension structures are already available in the technical literature; all of them posses certain limitations and drawbacks and no single method is optimal for all problems. The engineer may select the proper combination of methods best suited to the designer’s needs. In this paper it is proposed a combined method to achieve satisfactory equilibrium configurations for fabric tension structures. The force density method (FDM) implemented with topological mapping (TM) is used as a search engine for the preliminary design, and a procedure that employs nonlinear structural analysis is proposed for final refinement of the initial equilibrium configuration hence allowing the use of the same analysis tool for both refinement of the solution and analysis under loading.

Al diseñar una estructura tensada la forma inicial es normalmente desconocida. Los métodos de búsqueda de forma permiten al ingeniero obtener una geometría inicial dadas unas condiciones de contorno. Existen diferentes métodos de búsqueda de formas de equilibrio, pero todos tienen limitaciones y no existe uno único óptimo para cualquier tipo de problema. El ingeniero debe elegir la combinación de métodos que mejor se adapte a sus necesidades. En este artículo se propone un método combinado para generar configuraciones de equilibrio satisfactorias en estructuras tensadas. Como motor de búsqueda para el diseño preliminar se emplea el método de las densidades de fuerza (FDM) implementado con mallado en topología (TM), y se propone un procedimiento basado en análisis no lineal de estructuras para el refinamiento de la configuración inicial de equilibrio, permitiéndose así el empleo de las mismas herramientas tanto para el refinamiento de la solución inicial como para el cálculo posterior bajo condiciones de carga.

The design of tension structures is a relatively complex problem because of the relationship between the stress distribution and the final shape

Form-finding for tension structures is nowadays accomplished by numerical methods and computers. The parameters involved in the form-finding process are: i) type of element to model the fabric, ii) surface topology, iii) body forces, iv) surface tractions, v) nodal coordinates of the model, vi) geometry boundary conditions, and vii) internal stress distribution. Body forces and surface tractions are usually neglected during the shape-finding stage due to their insignificant influence in the geometrical shape at equilibrium. They can also complicate the initial equilibrium problem since the local direction and net magnitude of the loads may depend on an unknown initial surface shape. These loads are usually incorporated later into the stress analysis stage.

In most existing methods of equilibrium configuration, the shape is treated as an unknown. Then, it is necessary to introduce a set of parameters to ensure a unique solution. Ideally, one can use a nonlinear displacement analysis technique to solve the initial equilibrium problem using surface finite elements. In this case, the element stresses and the fixed nodes can be prescribed by the designer to guarantee a unique solution to the problem. The results of the equilibrium problem are the non-fixed nodes position. This technique has an important drawback: a geometrical initial mapping with a guess of the nodal coordinates is needed.

Several methods for solving the initial equilibrium problem (form-finding) have been developed to overcome some of the drawbacks associated with the nonlinear displacement method. In many of these methods a variety of limitations are imposed on the solution to transform the general nonlinear problem into a linear form. Siev and Eidelman

Haber and Abel

Other methods of form-finding are based on dynamic analysis. One of the main contributors to these methods was Barnes

Hernandez-Montes

Each of the methods mentioned above possesses certain limitations and drawbacks, therefore, no single solution method is optimal for all problems. The engineer may select the proper combination of methods best suited to the problem at hand and the designer’s needs. In this paper we take advantage of some of the existing methods to propose a combined method for computing satisfactory equilibrium configurations. The objective is to develop a straightforward methodology to achieve final equilibrium configurations for fabric tension structures, making use of stress control over surface elements, and taking advantage of the FDM with topological mapping as preliminary shape-finding.

The procedure presented in this paper is shown in

In order to make the paper self-contained an introduction to FDM and TM are included.

The starting point for the FDM

For a given pin-joint network with n nodes and m branches, the branch-node matrix C is a

The nodes Pi have coordinates (_{i}_{i}_{i}_{ix}_{iy}_{iz}_{j}_{j}

The great advantage of the FDM is the introduction of a parameter _{j}_{j}_{j})_{j}

Next, the _{x}, p_{y}_{z}

where _{j}^{T}QC matrix is positive definite if we have a prestressed network (_{j}

An initial guess of the shape is needed in order to generate a network, for instance the approximate shape of the horizontal projection. With this horizontal projection it is easy to consider a network even when only the connections are picked as input for the FDM; this way of doing is called “geometrical mapping”, in order to differentiate it from the topological mapping method (TM).

The main feature of the topological mapping (TM) introduced by Hernández-Montes

The topology of the network is obtained with the following input data: the number of nodes in the second step, the number of steps and the type of relation between one step and the following one.

Three types of basic relations are defined between steps: A, B and C, see

In order to correctly close the topological mapping a classification of nodes and connections has been defined for these networks. For the basic relationships A, B and C (

The construction of the network begins with one node which is connected to the n nodes of the second step. However, this initial node may be suppressed whenever there is the need to create a new interior contour in the net. Moreover, the nodes and ring branches located in the last step constitute the contour of the equilibrium shape. The number of nodes on each step is double the amount on the previous one when relationship type A is used, and remains constant when using a relationship type B or C.

TM-FDM method assigns the fixed points to the nodes located at the last step through the calculation of the distance between consecutive fixed points and the perimeter formed by them in such a way that the distribution of nodes of the last ring is done proportionally to the real distance between the fixed points. This assignation is then performed automatically.

The number of steps used in TM depends on the desired degree of refinement. A large number of steps will result in a mesh formed by a higher density of triangles. Once the number of steps has been decided, it is desirable to automatically generate the topological relations. These topological relations can be analyzed so as to achieve one goal: to obtain a mesh as uniform as possible, i.e., a mesh having triangles whose areas are as similar as possible. We have conducted an empirical analysis in which, for a given shape of the contour and a certain topological relationship, we have calculated the deviation from the average of the area of the triangles forming the equilibrium configuration. This analysis was repeated for a sample consisting of:

Varying number of steps from 10-16.

All possible topological relationships (once the number of steps is fixed).

Different forms of the contour: regular and irregular.

From the obtained results, we have deduced a series of topological relations for the different number of steps. These relations are recommended in the generation of equilibrium configurations with the TM-FDM method and are shown in

NUMBER OF STEPS | Number of nodes for the second step | Topological relations from steps 3 to end | |||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|

3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | ||

10 | 6 | A | B | A | B | B | A | B | B | ||||||

11 | 6 | A | B | A | B | B | A | B | B | B | |||||

12 | 6 | A | B | A | B | B | A | B | B | B | B | ||||

13 | 6 | A | B | A | B | B | A | B | B | B | B | B | |||

6 | A | A | B | B | A | B | B | B | B | A | B | ||||

14 | 6 | A | B | A | B | B | A | B | B | B | B | B | B | ||

6 | A | A | B | B | A | B | B | B | B | A | B | B | |||

15 | 6 | A | B | A | B | B | A | B | B | B | B | B | B | B | |

6 | A | A | B | B | A | B | B | B | B | A | B | B | B | ||

16 | 6 | A | B | A | B | B | A | B | B | B | B | B | B | B | B |

6 | A | A | B | B | A | B | B | B | B | A | B | B | B | B |

The constant stress triangular element is used; it is described in

where _{x}_{y}_{xy}

A link is established between the stresses in a membrane model and the cable forces in its corresponding virtual triangulated network. First, a local coordinate system coplanar with the element is used. Afterwards, the components of the nodal forces (in a Cartesian system) are transformed into a skew system parallel to the edges of the element (

where ∆ is the area of the triangle. Let _{rs}_{rs}_{rs}_{rs}_{rs}

The shape of tension structures depends on their internal stress pattern. Nonlinear structural analysis is needed to obtain a new equilibrium configuration when loading (or a perturbation) is added to the Initial Equilibrium Surface Model. In this paper nonlinear analysis is employed with the hypothesis of large rotations and small strains, together with the constant stress triangular finite element previously described. New equilibrium configurations are obtained due to perturbations introduced in the system. Perturbations come from changes in the stress field, these changes generate non equilibrated nodal loads. New equilibrium configurations with subsequent geometric changes are obtained using the Stiffness Matrix Method

In case of the plane stress triangular finite element, the elastic stiffness matrix of the element (K^{e}_{E}) is

The components of the previous matrix are:

where E and ^{e}_{E})ij is the element in row i and column j of the matrix K^{e}_{E}.

The geometric stiffness matrix associated with the elastic stiffness matrix defined in ^{e}_{G}) is

where

Given an initial equilibrium configuration, usually a different equilibrium shape close to the initial shape might be obtained, changing both the geometry and the stress distribution of the initial structure. For example, it may be attractive to obtain an equilibrium configuration with a constant stress field, or other prescribed stress patterns, by changing the Initial Equilibrium Surface Model (derived from the initial cable model equilibrium configuration obtained by FDM and TM). It is possible to generate equilibrium configurations with an imposed stress distribution using the Stiffness Matrix Method. The procedure is based on the same displacement analysis techniques used to study the behaviour of tension structures under varying load condition. In this method geometrical nonlinear analysis techniques are used to compensate unbalanced forces. These unbalanced forces result from the imposition of a new stress field to the Initial Equilibrium Surface Model. The applied procedure contains the following steps, schematically described in

Step 1: Introduce the necessary data for the topology mapping: the coordinates of the fixed nodes, the topology of the pin-joint network and the force density values.

Step 2. Calculate an initial equilibrium configuration based on the force density method. The coordinates of the non-fixed nodes are obtained.

Step 3: Compute the membrane stresses using

Step 4: Modify the membrane stresses to the desired values. This introduces a perturbation in the system, and an unbalanced force vector (ΔF) appears.

Step 5: Assemble individual element stiffness matrices to form the elastic stiffness matrix K_{E} and the geometrical stiffness matrix K_{G} of the structure with the desired stresses.

Step 6: Equilibrium requires that the sum of all forces must be zero at each node, so a new equilibrium configuration of the structure is obtained by allowing the deformation of the structure under the unbalanced forces (ΔF).

where a is the vector of nodal displacements and ΔF is the vector of nodal forces resulting from the perturbation induced by the imposed stress field, which induces a non-equilibrium configuration. Each iteration requires the calculation of:

The new nodal positions, adding to the previous nodal position the nodal displacements a.

The vector of unbalanced forces ΔF, calculated with the stress field variation. Small variations are preferred for numerical convergence; the magnitude of these variations can be easily controlled with small values of the Young’s modulus E.

The error, which can be computed as the second norm of the vector ΔF (I have named it as error Type 1).

In case the error (Type 1) is bigger than an upper bound value, K_{E} and K_{G} are recalculated taking into account the new nodal position and the new stress field.

Step 7: Compute the vector of unbalanced forces ΔF with the final nodal position obtained in Step 6 and the desired stresses. Finally, compute the error as the second norm of the vector ΔF (Error Type 2). If this error is bigger than an upper bound value, go to Step 4, with the previous obtained nodal coordinates.

A tension structure with 8 fixed points is designed. The first step is the solution of the form-finding problem, TM and FDM is used. The necessary input data for this first step are the following:

Boundary conditions:

Number of fixed points: 8

Coordinates of the fixed points: (-6,-45,0), (6,-45,0), (27,-24,15), (27, 24,15), (6, 45,0), (-6, 45,0), (-27, 24,15), (-27, -24,15)

Topology:

Desired number of steps: 12

Type of network: obtained from

Force density values

For the inner branches: 1

For the external branches: 10

With these simple inputs the geometry corresponding to the initial equilibrium configuration is automatically obtained (

The cable model is transformed into a membrane model composed by triangular elements; the stress (units of force per length) field is obtained by

For simplicity, in this example an isotropic and uniform stress field has been chosen in the membrane. The fixed value of the principal stresses is 1 (force/length units). An imposed value of force of 48 (force units) has been given to the exterior cable. A Poisson’s ratio equal to zero and a Young’s modulus equal to 0.5 have been used. After few iterations (Newton-Raphson) from steps 4 to 7 a new equilibrium configuration is obtained. In order to fulfill with the requirements of error Type 2, the numerical procedure has gone four times through step 7.

1,15E+01 | 3,16E+00 | 4,47E-01 | 1,38E-02 | 2,18E-03 | 5,58E-04 | 1,55E-04 | 5,61E-05 | 1,84E-05 | 6,76E-06 | |

1,26E+00 | ||||||||||

1,26E+00 | 1,63E-02 | 4,01E-04 | 7,83E-06 | – | – | – | – | – | – | |

2,17E-01 | ||||||||||

2,17E-01 | 2,86E-03 | 4,65E-06 | – | – | – | – | – | – | – | |

8,71E-02 | ||||||||||

8,71E-02 | 9,67E-04 | 1,01E-06 | – | – | – | – | – | – | – | |

7,34E-02 | ||||||||||

7,35E-02 | 6,24E-04 | 7,78E-07 | – | – | – | – | – | – | – | |

6,89E-02 |

An extended application of topological mapping to obtain a final equilibrium configuration of a membrane structure under controlled stress conditions is presented. The procedure deals with a topology based pin-joint network that is easy to implement and the force density method (FDM). The procedure is versatile and it represents a simple and straightforward method to achieve optimal equilibrium configurations through the use of finite element analysis tools for both refinements of the solutions and for their analysis under loading conditions.