Acta Polytechnica CTU Proceedings doi:10.14311/APP.2018.15.0036 Acta Polytechnica CTU Proceedings 15:36–40, 2018 © Czech Technical University in Prague, 2018 available online at http://ojs.cvut.cz/ojs/index.php/app EFFECT OF VISCOUS DAMPING IN NUMERICAL CALCULATION OF CABLE-MEMBRANE STRUCTURES BY THE METHOD OF DYNAMIC RELAXATION Miloš Huttner∗, Petr Fajman Czech Technical University in Prague, Faculty of Civil Engineering, Thakurova 7, 166 29 Prague 6, Czech Republic ∗ corresponding author: milos.huttner@fsv.cvut.cz Abstract. The aim of this article is to assess the speed of convergence of numerical calculation of cable-membrane structures using dynamic relaxation method in the process of finding critical damping pre-calculated on undamped system. The procedure is tested on four different constructions in six numerical cases. The variety of examples is as large as possible to demonstrate the greatest versatility of the test procedure. The efficiency of the procedure is evaluated based on the number of iterations. Keywords: Cable-membrane structures, dynamic relaxation, critical damping, viscous damping. 1. Introduction The design of cable-membrane structures is realized using numerical models of construction. Execution of the numerical analysis is not easy. The shape of the structure, unlike conventional structures, is not known in advance. The initial equilibrium shape of the structure must be found during the phase called the form-finding process. The second phase is the static response process. Both processes represent a geometrically non-linear case from a mathematical point of view, but each process has its own specifics. Different matrix methods have been developed to address both processes (for example force density method [1], update reference strategy [2], transient stiffness method, etc.). The method of dynamic re- laxation represents a very suitable alternative to ma- trix procedures [3]. In addition, the method can be easily used for both phases of the design of cable- membrane structures. Consuming CPU time with more unknowns takes approximately linearly, so the method is very well applicable to models with a large number of degrees of freedom [4]. The dynamic relaxation method builds on the princi- ple that the structure under free oscillation is directed to static equilibrium. The fictitious dynamic analysis is used to find static equilibrium. The artificial pa- rameters (masses, a damping factor, a time interval) influence the speed of convergence. The convergence is the fastest if the damping (a viscous damping is considered) is critical [4]. For fast convergence, it is necessary to find the fundamental oscillation fre- quency, which is not easy for geometrically non-linear cases. One possible way is to use a pre-calculation with an undamped oscillation mode. The aim of this article is to assess the speed of con- vergence of numerical calculation of cable-membrane structures in the process of finding critical damping pre-calculated on undamped system. The procedure is tested on four different constructions in six numerical cases. Both design processes (form-finding and static response) will be included. 2. Methodology Critical viscous damping coefficient can be found from [4]: C = 4πmf, (1) where C is viscous damping coefficient, m is mass (or vector of masses), f fundamental frequency of oscillation. The fundamental frequency of oscillation can be ob- tained pre-calculated on undamped system, with C set to zero. Based on the coordinate displacement record, the fundamental oscillation frequency is determined by the formula [4]: f = 1 N × ∆t , (2) where ∆t is time interval, N is the number of itera- tions required to complete one cycle of oscillations in a fundamental mode. As can be seen from the Equation 1, the critical damping factor is a multiple of the mass vector. This makes it easy to compare whether another damping coefficient is converging faster. This article investigates the determination of the damping coefficient based on the oscillation of the structure in the undamped mode (according to Equa- tion 2) with one cycle of oscillation. The procedure is a comparison with the otherwise set damping factor for a total of six computational cases (on four different structures). The efficiency of the process is evaluated based on the number of iterations. Other parameters are considered as follows: The time interval is selected one second. The mass vector is determined based on the stiffness of each node and 36 http://dx.doi.org/10.14311/APP.2018.15.0036 http://ojs.cvut.cz/ojs/index.php/app vol. 15/2018 Effect of viscous damping in numerical calculation... is set for each node separately (which is one of the most effective methods, for example, see [5]). 3. Testing There are a total of four reference examples, which are used together for a total of six cases. The variety of examples is as large as possible to demonstrate the greatest versatility of the test procedure. Structures are modeled in 3D space and include the form-finding (FF) and static response (SR) processes. Cables are numerically approximated using bar (FF and SR) and cable (SR) elements, see more [6]. Mem- branes are approximated using isoparametric triangu- lar elements (for FF and SR), see [3] or [7]. 3.1. Example 1 The structure is taken from [6]. A scheme of the model structure can be seen in Figure 1. The span L = 30.48 m. There have been many experiments on this real orthogonal cable net. Thus, numerical results can be compared with real-measured values. Cable elements are used to model cables. All cable segments have same cross-sectional ar- eas 1.4645 × 10−4 m2 and Young’s modulus 8.2737 × 1010 N/m2. The slack lengths of inner cables are 30.419 m, and the remaining cables are 31.76 m. Each of the four internal joints 4, 5, 8 and 9 carries a downward concentrated load 35.56 kN. A small self weigh of 1.46 × 10−3 kN/m for cables is considered. The structure has 12 degrees of freedom. Only the static response of the design is solved. Figure 1. Schema of example 1. Perspective view. 3.2. Example 2 The structure is taken from Lewis [4]. The results can be compared for both processes (form-finding and static response). A structure schema is shown in Figure 2. The symmetry of initial configuration about the x and y centre axes is maintained. The coordinates of node 29 are x = 12 m, y = 16 m and z = 3 m. Bar elements are used to approximate the cables. The structure has 45 degrees of freedom. For the form-finding process, initial geometry is gen- erated using a pre-tension force 90 kN in x direction, (except for the centre cable where it is 120 kN), a 30 kN in the y direction, plus a vertically point load 4.8 kN at all free nodes. For static response process the cross-sectional areas is assumed to be 350 mm2 in the x direction and 120 mm2 in the y direction. The elastic modulus is 160 kN/mm2. The structure is subjected to the further increase in load 2.0 kN at all free nodes. Figure 2. Schema of example 2. Perspective view. 3.3. Example 3 It is a fictitious membrane structure. A scheme of the model structure can be seen in Figure 3. The structure is supported throughout the circuit and the dimensions are L = 10 m and H = 4 m. To generate mesh, triangular elements with a maximum edge size of 1.25 meters are used. Thus, a model with 231 degrees of freedom is created. The model is tested for both design phases (FF and SR). For the FF process, a constant horizontal pre- tension σ = 2000 kN/m2 is used (the same for di- rections x and y). Thickness of membrane is 0.001 m. For a static response, a vertical flat load fz = 500 N/m2 is considered. The membrane pa- rameters are: the elastic modulus is 30 N/mm2, and the Poisson’s ratio is 0.2. Figure 3. Schema of example 3. Perspective view. 3.4. Example 4 It is a fictitious cable-membrane structure that has the parameters of a real construction. A scheme of the model structure can be seen in Figure 4. The structure is supported in nodes 1 to 6 and both top rings. The cables are around the perimeter of the structure. The dimensions of structure are L = 8 m and H = 3 m and the radius of the two top rings is 1 m. To generate mesh, triangular elements with a 37 Miloš Huttner, Petr Fajman Acta Polytechnica CTU Proceedings maximum edge size of 0.7 meters are used. Thus, a model with 1623 degrees of freedom is created. The model is tested only for form-finding process. For the FF process, a constant horizontal pre- tension σ = 1000 kN/m2 is used (the same for di- rections x and y). The thickness of membrane is 0.001 m. The horizontal pre-tension force 50 kN is considered in all cables. Figure 4. Schema of example 4. Perspective view. 4. Results All cases were tested in the author’s script for solving cable-membrane structures using the dynamic relax- ation method created in MATLAB®. The numerical accuracy of all calculations was controlled by a resid- ual limit force set at 0.0001 kN. 4.1. Example 1 The numerical finite element model (in final equi- librium state) for the static response is shown in Figure 5. Fundamental frequency was estimated at f = 0.0189 Hz, that corresponds to 0.237 multiple of the mass vector (α = 0.237). Figure 5. Example 1 (static response) in a final equilibrium state. The Figure 6 shows the number of iterations for damping coefficient values, when C = αm. It has been tested for α in the range from 0.05 to 0.6. It can be seen that the critical damping coefficient (red square) was set very well and that the number of iterations is very minimal. 4.2. Example 2 The numerical finite element model (in the initial and final state) for the form-finding process is shown in Figure 7. Fundamental frequency was estimated at Figure 6. Comparison of the number of iterations for the damping coefficient determined according to Equations 1 and 2 (red square) with different damping coefficients for Example 1 (SR). Figure 7. Example 2 (form-finding) in a final equi- librium state. f = 0.0926 Hz, that corresponds to 1.163 multiple of the mass vector (α = 1.163). The Figure 8 shows the number of iterations for damping coefficient values. It has been tested for α in the range from 0.8 to 1.6. It can be seen that the critical damping coefficient (red square) was set well and that the number of iterations is close to minimum. Figure 8. Comparison of the number of iterations for the damping coefficient determined according to Equations 1 and 2 (red square) with different damping coefficients for Example 2 (FF). Fundamental frequency for static response was es- timated at f = 0.0088 Hz, that corresponds to 0.111 multiple of the mass vector (α = 0.111). The Figure 9 shows the number of iterations for damping coefficient values. It has been tested for in the range from 0.04 to 0.22. It can be seen that the critical damping coefficient (red box) was set excellent 38 vol. 15/2018 Effect of viscous damping in numerical calculation... and that the number of iterations is minimal. Figure 9. Comparison of the number of iterations for the damping coefficient determined according to Equations 1 and 2 (red square) with different damping coefficients for Example 2 (SR). 4.3. Example 3 The numerical finite element model (in the initial and final state) for the form-finding process is shown in Figure 10. Fundamental frequency was estimated at f = 0.0151 Hz, that corresponds to 0.190 multiple of the mass vector (α = 0.190). Figure 10. Example 3 (form-finding) in a final equi- librium state. Figure 11. Comparison of the number of iterations for the damping coefficient determined according to Equations 1 and 2 (red square) with different damping coefficients for Example 3 (FF). The Figure 11 shows the number of iterations for damping coefficient values. It has been tested for in the range from 0.1 to 1.0. It can be seen that the critical damping coefficient (red square) was not set well and that the number of iterations is not close to minimum. Fundamental frequency for static response was es- timated at f = 0.0077 Hz, that corresponds to 0.096 multiple of the mass vector (α = 0.096). The Figure 12 shows the number of iterations for damping coefficient values. It has been tested for in the range from 0.04 to 0.22. It can be seen that the critical damping coefficient (red square) was set excellent and that the number of iterations is minimal. Figure 12. Comparison of the number of iterations for the damping coefficient determined according to 1 and 2 (red square) with different damping coefficients for Example 3 (SR). 4.4. Example 4 The numerical finite element model (in the initial and final state) for the form-finding process is shown in Figure 13. Fundamental frequency was estimated at f = 0.0417 Hz, that corresponds to 0.523 multiple of the mass vector (α = 0.523). Figure 13. Comparison of the number of iterations for the damping coefficient determined according to 1 and 2 (red square) with different damping coefficients for Example 4 (FF). The Figure 14 shows the number of iterations for damping coefficient values. It has been tested for α in the range from 0.2 to 0.9. It can be seen that the critical damping coefficient (red box was set well and that the number of iterations is close to minimum. 5. Conclusions This article investigated the determination of the damping coefficient based on oscillation of the struc- 39 Miloš Huttner, Petr Fajman Acta Polytechnica CTU Proceedings Figure 14. Comparison of the number of iterations for the damping coefficient determined according to 1 and 2 (red square) with different damping coefficients for Example 4 (FF). ture in the undamped mode with one cycle of oscil- lation. The procedure was a comparison with the otherwise set damping factor for a total of six com- putational cases (on four different structures). The efficiency of the process was evaluated based on the number of iterations. Base on the results presented in previous section, it can be concluded that the proposed method of deter- mining the dynamic relaxation damping coefficient is very suitable for cable-membrane structures. When setting a damping coefficient according to the pro- posed procedure, the minimum number of iterations (or near to minimum) was achieved in four of the six test cases. Excellent results are achieved by the static response analysis procedure with a construction close to the final state. It has been shown that only one oscillation cycle can be used to determine the fundamental frequency. Using multiple cycles would probably improve the frequency estimate, but the CPU time would increase. It could be tested in further studies. List of symbols C viscous damping coefficient [t s−1] m mass (or vector of masses) [t] f fundamental frequency of oscillation [Hz] ∆t time interval [s] N the number of iterations required to complete one cycle of oscillations in a fundamental mode [–] L the span of the structure [m] H the height of the structure [m] α multiple of the mass vector [–] Acknowledgements The results presented in this paper are outputs of the research project “SGS18/037/OHK1/1T/11- Development and application of advanced algorithms for numerical anal- ysis and modelling in mechanics of structures and ma- terials” supported by the Czech Technical University in Prague. References [1] K. Linkwitz, H. J. Schek. A new method of analysis of prestressed cable networks and its use on the roofs of the Olympic game facilities at Munich. IABSE, 1972. [2] E. R. Bletzinger. A general finite element approach to the form finding of tensile structures by the updated reference strategy. Int J Space Struct (14):131–146, 1999. [3] B. H. V. Topping, P. Iványi. Computer Aided Design of Cable Membrane Structures. Saxe-Coburg Publications, 2007. [4] W. J. Lewis. Tension structures: form and behaviour. Thomas Telford, 2003. [5] M. Huttner, J. Maca, P. Fajman. The efficiency of dynamic relaxation methods in static analysis of cable structures. Advances in Engineering Software (89):28–35, 2015. [6] H. Deng, Q. F. Jiang, A. S. K. Kwan. Shape finding of incomplete cable-strut assemblies containing slack and prestressed elements. Computers and Structures (83):21–22, 1767–1779, 2005. [7] M. R. Barnes. Form and stress engineering of tension structures. Structural Engineering Review (6 (3-4)):175–202, 1994. 40 Acta Polytechnica CTU Proceedings 15:36–40, 2018 1 Introduction 2 Methodology 3 Testing 3.1 Example 1 3.2 Example 2 3.3 Example 3 3.4 Example 4 4 Results 4.1 Example 1 4.2 Example 2 4.3 Example 3 4.4 Example 4 5 Conclusions List of symbols Acknowledgements References