Comparison of results between SIMPACK® and the current method—Line 02.

## Abstract

This chapter is devoted to an application of a finite element method formulation to forecast the static and mechanical behavior of catenary-like structures subject to general force distributions, whose development was motivated by the need of installing assemblies of containment structures, called log boom lines, upstream a hydroelectric power plant to protect its integrity from the threats that logs carried through the river pose on it. Each log boom is modeled by a tridimensional truss element and the entire lines by assemblages of trusses. While the external forces, modeled with the aid of both simulations from computational fluid dynamics and experiments from a towing tank, originate from both the river stream and the logs that accumulate through the extension of the lines, the internal forces are calculated from classic expressions of solid mechanics; hence, the numerical method imposes equilibrium between them, which ultimately defines the geometry assembly. Verification and validation were performed at both model and prototype scales, and the results corroborated the accuracy of the tool for a series of flow conditions.

### Keywords

- finite element method
- fluid-structure interaction
- containment grids
- catenary-like
- numerical simulation

## 1. Introduction

The operation of a hydroelectric power plant at the Madeira River is compromised due to the presence of many logs carried through the water stream that can damage the power plant machinery and reduce its production efficiency. To prevent such damage, nine assemblies of containment structures, called log booms (Figures 1 and 2), were installed near the power plant and across the river to retain and deflect these logs. Nevertheless, as the influence of logs and stream can exert large loads over the assemblies, they are under the risk of structural failure; so a numerical tool based on the finite element method (FEM) was developed to predict the distributions of loads and equilibrium configurations of the lines, therefore assisting their surveillance and maintenance. The numerical tool was part of a research and development (R&D) project developed by the Institute for Technological Research (IPT), and it can be extended to similar problems in which the structures behave likewise, as is the case of fishing nets and cleaning structures [1]. Moreover, while commercial software could be used for the same purpose, the intrinsic drawbacks of large computation time, automation, and customization difficulties were generally not justified by the potential increase in accuracy, as the tool would eventually be integrated in a larger numerical simulator.

## 2. Numerical model

The mathematical model used in the present formulation is based on the works of Felippa and Gavin [2, 3, 4], whereas the numerical implementation was developed within the Matrix Laboratory (MATLAB) environment [5]. The material deformations are assumed elastic only (i.e., no plasticity or yielding) as each log boom is mainly composed of steel, making the unit respond linearly until approximately 80% of the yield stress [6], but the displacements are geometrically nonlinear, as the forces, the geometric stiffness, and the assembly equilibrium configuration depend on one another.

### 2.1 Mathematical model

Consider the three-dimensional truss element composed by six degrees of freedom (DoF) that represent all the possible translation directions its nodes have. Such element is used to model the structural behavior of each log boom (Figure 1).

With this numbering, the element constitutive stiffness matrix

while its geometric stiffness matrix

These matrices represent the stiffness of the truss due to its constitutiveness (as

An assembly is then modeled by the joint of **N** log booms, as in Figure 2; thus, to account for the assembly overall stiffness

Now, considering a function

The equilibrium configuration

Truncating Eq. (5) to the first order,

which, however, is a linear approximation of this nonlinear formulation. Given such nonlinearity, the problem must be solved numerically, and all the equations presented insofar are part of an iterative scheme, in which

As expected, by similarity,

### 2.2 Force calculations

The loads acting on the assembly originate from the hydrodynamic interaction and from the log accumulation and must be counterbalanced by the internal forces, which are a function of the mechanical properties of the material. Although the formulation is tridimensional, for the purposes by which the model was made, it is sufficient to obtain the equilibrium configuration on the

#### 2.2.1 External loads

To estimate the external forces, a database was created with the aid of computational fluid dynamics (CFD) (Figure 3). The software Siemens Star-CCM+, version 12.02.011, was used for these numerical simulations: both water and air were admitted incompressible, and the flow field was modeled using the unsteady Reynolds-averaged Navier-Stokes equations implicitly (at the first iterations, the movements were frozen to ensure fluid stabilization and then gradually released, still seeking stabilization); the interaction between the fluids was represented by the volume of fluid Eulerian multiphase model, and the multiphase conditions were defined by the flat wave model with a numerical damping, added in the momentum equation on the vertical direction, on both the inlet and outlet boundary conditions to minimize the effects of reflection and ensure better convergence; finally, the k-ω SST model was adopted.

A series of velocity magnitudes

Each part of the log boom that is free to rotate was assigned with a pair of coefficients (i.e., the grids and the chassis) such that the total hydrodynamic force acting on the trusses is retrieved based on them. So, for the

Once the hydrodynamic coefficients were obtained, the contribution from the logs that accumulate along the lines was assumed to be increments on the normal coefficients [10] and estimated with the aid of data from experiments at the towing tank of the IPT (Figure 6). Thus, considering already such increments, Eq. (9) accounts for the total external forces acting on the center of the element, which are equally distributed to its nodes and further transformed from local to global coordinates as well, according to the angle

So,

#### 2.2.2 Internal loads

The tension distribution throughout the trusses is calculated based on classic linear relations of solid mechanics, and it acts on the axial direction of the element. Therefore, for the same

This tension must be transformed to internal forces

### 2.3 Initial condition

Originally, the solution strategy was to adopt the line completely stretched along the

However, such quasi-static approach is computationally expensive because of the repetitiveness of the iterative procedure. Additionally,

in which

### 2.4 Static condensation

An intrinsic characteristic of structural problems is the presence of fixed degrees of freedom (F) in addition to the free ones (L). Thus, a common practice in FEM, known as static condensation, is to renumber the DoF seeking to reorder and split the linear system into two sub-systems. So, the displacement array can be written as:

and the unbalance array can be similarly split:

For the log boom assembly, the fixed DoF are related to the mooring points, so

Additionally, the matrix

Thus, from Eq. (4),

These sub-systems are solved recursively: the first of Eq. (18) is solved for

### 2.5 Sub-relaxation

Several times, the initial condition itself is not sufficient to ensure convergence because the average velocity used to calculate the catenary “weight” might not create a representative starting geometry. So, whenever the direct method fails, a sub-relaxation parameter is applied to the force unbalance as a way to prevent large deflections through the line, which is especially interesting during the first iterations; the restrained deflections are generally sufficient to adequately converge the simulation. The sub-relaxation

### 2.6 Post-processing

After numerical convergence, a series of variables can be outputted, such as the final geometry, the velocity magnitude and incidence angle acting over each log boom, the reaction forces, and the tension distribution along the length of the line

## 3. Results

In order to verify and validate the tool, thus assessing the numerical method developed, many simulations were compared to theoretical, experimental, and numerical data. As the ASME Verification & Validation Committee defines [11], validation determines the degree of accuracy that a numerical model has in representing the real world from its perspective, whereas verification evaluates its intrinsic errors and uncertainties. So, the committee states that verification precedes validation, as they must be performed to ensure the numerical reliability of the solution. In this work code and solution verifications were done by systematically refining the discretization and comparing the outputs with theoretical data from the problem of a catenary [12]. Following, the validation was done for two major cases: by comparing the solution with experimental results from a 1:10 scale model tested in the towing tank at the Naval and Ocean Engineering Laboratory of the IPT [13] and by comparing the solution with numerical simulations from SIMPACK® Multi-Body Simulation software in prototype scale [12, 14]; the later comparison had to be done against numerical data because it was unfeasible to instrument the in situ assemblies. Additionally, only hydrodynamic cases (i.e., without logs) were compared.

### 3.1 Verification from a theoretical catenary model

While evaluation of the error was performed through comparison with the catenary model, estimation of the error was done through classical Richardson extrapolation [12]. Code and solution verifications were performed by systematic grid refinement, with the coarsest mesh having 100 elements, the finest 800, and in between these two, successive meshes were created such that the number of trusses was consecutively increased by a value of approximately

Thus, according to Eq. (20),

The relative error is simply evaluated by calculating the difference between the theoretical and numerical values:

and the uncertainty

Finally, Figure 14 shows a comparison between numerical and theoretical geometries; as they practically overlap, the verification corroborates the precision of the proposed model. As a note, in order not to be much close to the correct catenary geometry, the verification was conducted with the line completely stretched along the

### 3.2 Model scale validation

Reduced model experiments were conducted at the facilities of the IPT in a 1:10 scale with the purpose of studying the hydrodynamic behavior of such structures [13]. To do so, the model, developed mainly in polycarbonate and polymer, was constructed with the aid of a laser cutting machine and a 3D printer; the assembly was ballasted by adding lead stripes along the structure but without changing its front area as to minimize drag interference; to test both symmetric and asymmetric setups, log booms were added from a five-unit symmetric configuration, while the left mooring point was offset by distances

Four S-Type S9M uniaxial load cells, produced by HBM, were used at each extremity, each of which with a nominal measure limit of 500 N, so the tensile forces acting on the leftmost and rightmost log booms could be acquired. Additionally, to understand how the stain behaves on each log boom and also how it is transferred throughout the structure, water strain gages, manufactured by Kyowa Electronic Instruments, were placed on one of the modules, as shown schematically in Figure 15. These locations were appropriately chosen with the aid of numerical simulations from finite element method: the criterion was to locate the stress concentration area that increased the measurement sensitivity while avoided the high strain gradients, since the latter can contribute to unreliable measurements.

Figures 16–18 compare experimental and numerical results for three different configurations as the function of the velocity. “LE” indicates “left extremity” and “RE” “right extremity,” while “Exp” and “Num” refer to experimental and numerical values. No numerical uncertainty was provided since the number of trusses was maintained equal the number of log booms.

The results agree satisfactorily: the differences are more evident as the free stream velocity increases; the numerical method reproduces the experimental observation that the tensile force is larger at the leftmost log boom, an expected trend since tensile loads are higher for elements that are parallel to the stream; the number of log booms seems not to affect the accuracy of the results, as the error remains practically the same regardless of the configuration; the larger differences are within

### 3.3 Prototype-scale validation

A prototype-scale validation is presented, and the results were compared to simulations from SIMPACK® Multi-Body Simulation software. The variables of comparison are again the reaction forces at the leftmost and rightmost anchor points, and the relative differences

#### 3.3.1 Line 02

The first prototype validation was performed for line 02 with a variable velocity profile, similar to the one presented in Figure 5. An example of converged solution for line 02 is presented in Figure 19, in which the equilibrium configuration, projected line, and velocity distribution are depicted.

The results from SIMPACK® and the current method are presented in [14]: the first simulation considers a 20% increase in the normal external force (

Table 1 shows that the maximum relative difference for the simulations that happen at the

SP (kN) | FEM (kN) | u (%) | |||
---|---|---|---|---|---|

LE | RE | LE | RE | LE | RE |

1025 | 814 | 1076 | 876 | 5.0 | 7.6 |

1239 | 1029 | 1305 | 1104 | 5.3 | 7.3 |

#### 3.3.2 Line 12

The second validation was performed for line 12 with several velocity profiles, also similar to the one presented in Figure 5. Analogously, Figures 19 and 20 depict an example of a converged solution.

Table 2 presents the comparison between SIMPACK® and the developed formulation. For this set of simulations, the maximum percentage difference happens at the

SP (kN) | FEM (kN) | u (%) | |||
---|---|---|---|---|---|

LE | RE | LE | RE | LE | RE |

1566 | 1234 | 1666 | 1295 | 6.4 | 4.9 |

1540 | 1204 | 1641 | 1267 | 6.6 | 5.2 |

1866 | 1479 | 1983 | 1543 | 6.3 | 4.3 |

2028 | 1578 | 2158 | 1659 | 6.4 | 5.2 |

1463 | 1124 | 1559 | 1181 | 6.6 | 5.1 |

1793 | 1387 | 1910 | 1459 | 6.5 | 5.2 |

1946 | 1487 | 2071 | 1564 | 6.4 | 5.2 |

1789 | 1382 | 1905 | 1454 | 6.5 | 5.2 |

1767 | 1359 | 1877 | 1424 | 6.2 | 4.8 |

1710 | 1306 | 1821 | 1374 | 6.5 | 5.2 |

2255 | 1748 | 2405 | 1840 | 6.6 | 5.3 |

1701 | 1314 | 1811 | 1381 | 6.4 | 5.1 |

1820 | 1414 | 1888 | 1425 | 3.7 | 0.8 |

#### 3.3.3 Lines 13

The last validation was conducted for the group of lines 13, i.e., for lines 13A, 13B, 13C, and 13D. Figure 21 is an example of a converged simulation for line 13C/13D.

Particularly for line 13, the magnitude and incidence angle remained constant throughout them, varying from simulation to simulation. In Table 3 these parameters are summarized: the first columns indicate the velocity range (VR) used in each simulation for the choice of the coefficients of Eq. (8), even though the actual velocity magnitude was sometimes not within that range.

Line 13A/13B | Line 13C/13D | |||
---|---|---|---|---|

VR (m/s) | V (m/s) | Angle (°) | V (m/s) | Angle (°) |

2.5 < V < 3.5 | 3.18 | 54.44 | 2.11 | 59.60 |

2.5 < V < 3.5 | 2.49 | 54.81 | 2.04 | 42.00 |

1.0 < V < 3.5 | 2.14 | 52.11 | 1.60 | 56.26 |

1.0 < V < 3.5 | 1.80 | 49.56 | 1.33 | 42.48 |

1.0 < V < 2.0 | 1.69 | 53.06 | 1.13 | 61.03 |

1.0 < V < 2.0 | 1.36 | 52.80 | 1.05 | 43.53 |

Table 4 presents the comparison between SIMPACK® and the developed formulation.

Line 13A/13B | Line 13C/13D | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|

SP (kN) | FEM (kN) | u (%) | SP (kN) | FEM (kN) | u (%) | ||||||

LE | RE | LE | RE | LE | RE | LE | RE | ||||

1550 | 1263 | 1655 | 1329 | 6.8 | 5.2 | 822 | 647 | 876 | 679 | 6.6 | 5.0 |

934 | 759 | 1000 | 801 | 7.0 | 5.5 | 1501 | 1318 | 1607 | 1401 | 7.1 | 6.3 |

745 | 618 | 793 | 650 | 6.5 | 5.2 | 537 | 436 | 572 | 458 | 6.5 | 5.1 |

577 | 486 | 615 | 512 | 6.7 | 5.5 | 600 | 525 | 639 | 556 | 6.6 | 5.9 |

429 | 358 | 455 | 376 | 6.0 | 5.2 | 209 | 165 | 221 | 173 | 5.8 | 4.7 |

280 | 234 | 297 | 247 | 6.2 | 5.4 | 341 | 301 | 361 | 318 | 6.0 | 5.5 |

The maximum percentage difference of the simulations happens again at the

## 4. Conclusions

A simple truss-based finite element method was proposed to simulate the load distributions and geometric configurations of assemblies that behave in a catenary-like manner, subject to external, variable loads. The formulation was tailored to the particular problem of log booms (structures that retain logs from reaching the machinery of hydropower plants) under the influence of river streams and the logs they convey. The nonlinear formulation imposes equilibrium between internal and external forces so that an iterative scheme must be numerically solved. The method was verified by comparison against analytical results from a theoretical catenary model: the relative error and uncertainty for the maximum and minimum forces were within 0.2%, while the mesh refinement order of convergence was close to 1. The tool was later validated against experimental model-scale data from the towing tank at the Institute for Technological Research [12] and prototype-scale numerical data from commercial software SIMPACK®, and all the results agree adequately [14]: for the experimental data validation, the numerical method was capable of reproducing the observations of the experiments, and the maximum relative discrepancy observed was about 6%. The differences are invariant to the increase in the assembly length, but they seem sensible to variations of the free stream. Likewise, the prototype-scale validations all show adequate agreement, regardless of the line configurations, and a maximum relative error, considering SIMPACK® as the reference, of less than 8%. These percentages corroborate the adequacy of the method for the purpose by which it was developed: to a have a fast, yet reliable, tool, to forecast the tension distributions throughout the lines such that it can be used as a simulator for the assessment of the structural safety of log boom assemblies.

## Acknowledgments

The authors would like to express their gratitude to the Santo Antonio Energy for funding the project (PD-06683-0116/2016), through the Research and Development Fund of Brazilian Electricity Regulatory Agency, and the Institute for Technological Research Foundation. They also would like to acknowledge the financial support and scholarship granted by the Coordination of Superior Level Staff Improvement, under projects 1784366 and 133846, and the Foundation for the institute for Technological Research.