Calculo de Formulas para Duque de Alba
Calculo de Formulas para Duque de Alba
Calculo de Formulas para Duque de Alba
RESUMEN
Los muelles son estructuras que sirven para el atraque de las embarcaciones. Una vez
que el barco llega al muelle es fijado a este con el fin de impedir que ocurran grandes
desplazamientos relativos entre ellos; para lograr este objetivo se hace uso de un sistema
de amarre formado por cabos y defensas. En esta situación los buques están expuestos a
diversas fuerzas externas producidas por la acción del viento, oleaje, corriente y demás
factores ambientales. En el presente trabajo se expone uno de los método que puede ser
utilizado para determinar las fuerzas estáticas de viento y corriente a las que eventualmente
puede estar sujeta una embarcación, así como también se presentan algunos conceptos
dinámicos de interés como las ecuaciones de movimiento, y se discuten las dificultades de
resolución, por vía estrictamente teórica, de la magnitud, el carácter, y la importancia
relativa de las acciones de oleaje.
Uno de los aspectos a considerar para el correcto diseño de las estructuras de amarre y
atraque son las fuerzas que le transmiten los cabos, cadenas de anclaje y defensas que
sujetan la embarcación. El cálculo de estas fuerzas puede resultar largo y complejo; para
facilitar el cálculo de estas fuerzas se ha desarrollado un programa denominado Amarre que
funciona en Visual Basic aplicado a Excel y que permite obtener las fuerzas de amarre, en
los cabos y defensas, de una manera rápida, sencilla y confiable.
Amarre utiliza el método de la rigidez aplicando el análisis matricial para obtener las
fuerzas y deformaciones del sistema de amare; para ello fue necesario desarrollar y deducir
las matrices y ecuaciones para el sistema específico (el buque y sus elementos de amarre);
considerando el pretensado de los cabos, la no linealidad del comportamiento de los
materiales que constituyen las defensas y la no linealidad de las cadenas de anclaje
(catenarias).
El usuario debe introducir las coordenadas globales del centro de gravedad del buque y
las distintas fuerzas externas aplicadas en su centro de gravedad. Para los cabos de amarre
es necesario introducir las coordenadas globales de cada bita, las coordenadas locales de
cada escoben, el pre-tensado inicial de cada cabo, y la curva de comportamiento Fuerza-
Deformación. La curva consiste en doce pares de puntos, donde las fuerzas deben ser
dadas en orden ascendente. Los datos de entrada de las defensas son: posición X (en
coordenadas globales) de cada defensa, longitud, curva de comportamiento Fuerza-
Deformación que sigue las mismas especificaciones que la curva de los cabos. Para las
cadenas de anclaje es necesario introducir: las coordenadas globales de cada bita, las
coordenadas locales de cada escoben. Dado que las cadenas tienen un comportamiento de
catenaria, para definir su curva será preciso introducir: la profundidad del agua, la longitud
Denise Da Costa
–i–
Resumen
de la cadena, peso sumergido, fuerza horizontal máxima, el peso del ancla y su factor de
seguridad, y en caso de existir, el peso del sinker.
Denise Da Costa
– ii –
Resumen
RESUMO
Os cais são estruturas que servem para o atraque de embarcações. O navio que chega
ao cais é ligado a este com o fim de impedir que ocorram grandes deslocamentos relativos
entre eles; para alcançar este objectivo utiliza-se um sistema de amarre formado por cabos
e defensas. Nesta situação as embarcações estão expostas a diversas forças produzidas
pela acção do vento, ondas, corrente e demais factores ambientais; na presente dissertação
expõe-se um dos métodos que pode ser utilizado para determinar as forças estáticas de
vento e corrente às quais eventualmente pode estar sujeita uma embarcação, e também se
apresentam alguns conceitos dinâmicos de interesse como as equações de movimento, e
discutem-se as dificuldades de resolução, por via estritamente teórica, da magnitude,
carácter e importância relativa das acções das ondas.
Amarre utiliza o método da rigidez aplicando a análise matricial para obter as forças e
deformações do sistema de amarração. É necessário desenvolver e deduzir as matrizes e
equações para o sistema específico (o navio e os seus elementos de amarração);
considerando o pré-esforço dos cabos, a não linearidade do comportamento dos materiais
que constituem as defensas e a não linearidade das cadeias de ancoragem (catenárias).
Denise Da Costa
– iii –
Resumen
Denise Da Costa
– iv –
Resumen
ABSTRACT
Docks are structures for berthing and mooring ships. Once the ship arrives at the quay it
is fixed with the purpose of diminishing the relative displacement between them both; in
order to achieve this objective, a mooring system is used. In this situation the ships are
exposed to diverse external forces produced by the action of the wind, waves, current and
other environmental factors; in the present these is exposed one of the method that can be
used to determine the static wind and current forces to which a vessel possibly can be
subject, as well as some dynamic concepts of interest like the equations of motion, and the
resolution difficulties are discussed, via strictly theoretical, of the magnitude, the character,
and the relative importance of wave loads.
For the correct design of a mooring structure it is necessary to know the forces
transmitted by the hawsers, chains and fenders. The calculation of these forces can be long
and complex; to facilitate the calculation, a software called “Amarre” (which works in Visual
Basic with Excel) has been developed. This software allows obtaining the mooring and
fender forces in a simple and reliable way.
Amarre uses the matrix analysis to obtain the forces and deformations of the complex
mooring and fender system. For that it was necessary to develop and to deduce the matrices
and equations for the specific system: the ship and its mooring’s elements; considering the
prestressed of the moorings lines, the material nonlinearity of the behavior of the fenders
and the nonlinearity of the anchorage chains (catenary).
The user has to introduce the gravity center coordinates of the ship and the external
forces applied at the gravity center. Hawser data requested are: mooring point coordinates in
global system, chock coordinates in local system, preload and load-extension curve. The
curve consists of twelve pairs of load-deflection points, that must be given in ascending
order. Fender input data are: the X coordinate of each fender, initial length and the load-
extension curve which follow the same rules at the curve for hawsers. Data for catenary
chains are: chock coordinate in local system, buoy coordinates (X, Y) in global system. For
computing the anchor chain load-extension catenary curve are necessary the following
characteristics of the chain system: water depth, lengths and unit weights of upper and lower
chain sections, the maximum horizontal load, weight of anchor and safety factor, and sinker
weight.
This data is processed in the program. The preliminary step is to compute the catenary
curve of the behavior for each chain according to an iterative process; then local coordinates
of chocks are transformed into global coordinates to obtain the direction cosines of each
element and with that in hand, their stiffness matrix, assembling therefore the stiffness matrix
Denise Da Costa
–v–
Resumen
of the structure. If the structure is unstable, the program will give a message alerting of that
fact. If this is not the case, the local forces of each member are determined. Due to the non
linear behavior of the curves it is necessary to use an iterative process, in which the value of
the elements stiffness will be updated at each iteration, until obtaining convergence. In all
this iterative process is verified the structure stability, as well as that no hawsers are working
in compression and that no fenders are working in traction. After obtaining the convergence
of the nonlinear problem, the forces in the mooring points are determined. Finally the results
are presented: displacement of the ship, forces in each element and forces in the mooring
points. This values are important for designing the mooring structures.
To improve the engineer criteria when designing a mooring structure, a parametric study
was realized comparing the most important characteristics of the mooring systems – material
nonlinearity of the behavior of the fenders, moorings lines configurations, vessel type and
anchorage chains contribution –
Denise Da Costa
– vi –
Agradecimientos
AGRADECIMIENTOS
Gracias a Dios....
....por la maravillosa familia que tengo y por el amor recibido.
Gracias a mis amigos; Alvaro Alfonzo-Larrain, Giomar Iglesias, Julia Estévez, Kelvyn
García, Néstor García, Luís Felipe París y Vasco De Freitas....
…por el apoyo constante y por la fuerza para seguir adelante. Gracias por superar la barrera
de la distancia y continuar aquí conmigo.
Gracias a mis amigos; Ana Isabel Cruz, Pedro Miguel Martins y Sara Cruz…
...por hacerme reír en los momentos malos y por compartir conmigo mis éxitos y alegrías.
Gracias por la confianza. Gracias por las meriendas con discusiones filosóficas.
Denise Da Costa
– vii –
Agradecimientos
Denise Da Costa
– viii –
Índice
INDICE
RESUMEN i
AGRADECIMIENTOS vii
INDICE ix
INDICE DE TABLAS Y FIGURAS xvii
GLOSARIOS DE TÉRMINOS xxi
1 INTRODUCCIÓN 01
1.1 ASPECTOS GENERALES 01
1.2 DESCRIPCIÓN DEL PROBLEMA 01
1.3 OBJETIVOS 02
1.4 HIPÓTESIS 03
1.4.1 Problema estático o quase-estático 03
1.5 DESCRIPCIÓN DEL CONTENIDO DE LA TESIS 03
– ix –
Índice
4 SISTEMAS DE DEFENSAS 61
4.1 TIPOS 61
4.1.1 Defensas de madera 61
4.1.2 Sistemas de pilares 62
4.1.2.1 Pilares de madera 62
4.1.2.2 Madera colgada 63
4.1.2.3 Pilares de acero 63
4.1.2.4 Pilares de concreto 63
4.1.3 Sistemas de goma 64
4.1.3.1 Defensas en compresión 64
4.1.3.2 Defensa en corte 65
4.1.3.3 Flexibles 66
4.1.3.4 Defensas de neumático 66
4.1.4 Sistemas hidráulicos e hidroneumáticos 67
4.1.4.1 Amortiguador hidráulico 67
4.1.4.2 Defensa flotante hidroneumática 67
4.1.5 Resortes 68
4.2 SELECCIÓN DEL TIPO DE DEFENSA 68
4.2.1 Condiciones de exposición del puerto 68
4.2.2 Atraque versus Amarre 68
4.2.3 Distancia máxima permitida entre el buque y el muelle 69
4.2.4 Tipo de muelle 69
4.3 CRITERIOS BÁSICOS DE DISEÑO 70
4.3.1 Procedimiento general de diseño 70
4.3.2 Pilares de defensa 71
4.3.3 Sistemas de defensas resistentes 72
Denise Da Costa
–x–
Índice
5 SISTEMAS DE AMARRE 73
5.1 DEFINICIÓN 73
5.2 ESTRUCTURAS DE ATRAQUE Y AMARRE 74
5.2.1 Plataformas 74
5.2.2 Duques de atraque 75
5.3 CABOS DE AMARRE 76
5.3.1 Tipos 76
5.3.2 Materiales 76
5.3.3 Características de los cabos de amarre 77
5.4 CADENAS DE ANCLAJE 78
5.4.1 Comportamiento 78
5.4.2 Materiales 78
5.4.2.1 Cadenas 78
5.4.2.2 Cuerdas de alambre 79
5.5 CRITERIOS BÁSICOS DE DISEÑO 80
5.5.1 Estructuras de amarre 80
5.5.2 estructuras de atraque 80
5.5.3 Cabos 81
5.5.3.1 Cabos de popa, proa y laterales 81
5.5.3.2 Cabos “resorte longitudinal” 81
6 METODOLOGÍA DE ANÁLISIS 83
6.1 METODO DE LA RIGIDEZ 83
6.1.1 Matriz de rigidez de un miembro 83
6.1.2 Matriz de transformación de desplazamientos y fuerzas 86
6.1.3 Matriz de transformación de desplazamiento 87
6.1.4 Matriz de transformación de fuerzas 94
6.1.5 Matriz de rigidez global de un miembro 95
6.1.6 Matriz de rigidez de la estructura 97
6.1.7 Matriz K22 100
6.1.8 Matriz K12 100
6.2 CATENARIAS 101
6.2.1 Ecuación general de la catenaria 101
6.2.2 Catenaria completamente desarrollada 104
6.2.2.1 Características geométricas y mecánicas 104
6.2.3 Catenaria no completamente desarrollada 107
6.2.4 Fuerza horizontal vs Distancia horizontal 108
6.2.4.1 Caso 1 109
6.2.4.2 Caso 2 109
6.2.4.3 Caso 3 109
Denise Da Costa
– xi –
Índice
– xii –
Índice
Denise Da Costa
– xiii –
Índice
Denise Da Costa
– xiv –
Índice
10 REFERENCIAS 205
ANEXOS 209
ANEXO A – MANUAL DEL USUARIO 209
Denise Da Costa
– xv –
Índice
Denise Da Costa
– xvi –
Índice de tablas y figuras
Denise Da Costa
– xvii –
Índice de tablas y figuras
4 SISTEMAS DE DEFENSAS
Figura 4.1 – Sistema de pilares de madera 62
Figura 4.2 – Sistema de defensa de Madera Colgante 63
Figura 4.3 – Defensa Cilíndrica 64
Figura 4.4 – Defensa en corte tipo Raykin 65
Figura 4.5 – Defensa de Neumático 67
5 SISTEMAS DE AMARRE
Figura 5.1 – Sistema de Amarre 73
Figura 5.2 – Distribución de la fuerza de impacto en los pilares inclinados de una 75
plataforma rígida
76
Figura 5.3 – Arreglo típico de los cabos en un buque
Figura 5.4 – Representación de la pérdida de rigidez de los cabos con tres cordones 77
por formación de lazos
80
Figura 5.5 – Cabo de anclaje compuesto por cadena y cuerda de alambre
Figura 5.6 – Ángulos horizontales y verticales para los cabos de amarre 82
6 METODOLOGÍA DE ANÁLISIS
Figura 6.1 – Estado general de cargas y deformaciones para una barra 84
Figura 6.2 – Desplazamiento permitido en el nodo i e impedido en el nodo j 84
Figura 6.3 – Desplazamiento permitido en el nodo j e impedido en el nodo i 84
Figura 6.4 – Posición relativa del sistema de coordenadas locales respecto al sistema 86
global
87
Figura 6.5 – Definición de los grados de libertad del sistema
88
Figura 6.6 – Representación del modelo en análisis con sus grados de libertad,
coordenadas globales y locales
89
Figura 6.7 – Sistema con todos los grados de libertad impedidos, excepto el
desplazamiento en 1
90
Figura 6.8 – Sistema con todos los grados de libertad impedidos, excepto el
desplazamiento en 2
91
Figura 6.9 – Sistema con todos los grados de libertad impedidos, excepto la rotación
en 3
92
Figura 6.10 – Sistema con todos los grados de libertad impedidos, excepto el
desplazamiento en 4
93
Figura 6.11 – Sistema con todos los grados de libertad impedidos, excepto el
desplazamiento en 5
94
Denise Da Costa
– xviii –
Índice de tablas y figuras
7 PROGRAMA AMARRE
Tabla 7.1 – Mensajes de error del programa Amarre 146
8 ESTUDIO PARAMÉTRICO
Tabla 8.1 – Cálculo de las fuerzas ambientales para el tanquero de 85.000 DWT 150
Tabla 8.2 – Cálculo de las fuerzas ambientales para el tanquero de 130.000 DWT 152
Tabla 8.3 – Comportamiento de la defensa tipo Foam expresado en % de fuerza vs % 154
de deformación
Tabla 8.4 – Fuerza máxima resistida por cada defensa según el tipo y tamaño (D1) 154
Tabla 8.5 – Comportamiento de la defensa tipo D & Square expresado en % de 155
fuerza vs % de deformación
Tabla 8.6 – Fuerza máxima resistida por cada defensa según el tipo y tamaño (D2) 156
Tabla 8.7 – Comportamiento de la defensa tipo Super Cone expresado en % de 157
fuerza vs % de deformación
Tabla 8.8 – Fuerza máxima resistida por cada defensa según el tipo y tamaño (D3) 157
Tabla 8.9 – Características de los cabos de nylon 159
Tabla 8.10 – Datos Fuerza vs %deformación para los 4 cabos a considerar 160
Denise Da Costa
– xix –
Índice de tablas y figuras
Denise Da Costa
– xx –
Glosario de términos
GLOSARIO DE TERMINNOS
Amarra: Hacer firme un cabo. También afirmar el barco a un muelle o boya (3)
Bita: Pieza metálica de uno o más brazos que sirve para amarrar un cabo o cable (3)
Codera: Cabo que se hace firme en la popa y que sirve para amarrar la embarcación o para
cambiar su orientación (3)
Defensas: Elementos que se colocan para preservar el casco de los choques y/o roces con
otros barcos o el muelle (3)
Escoben: En los buques se denomina así al orificio practicado en la proa, por el cual salen
las cadenas de las anclas y amarras (3)
Eslora: Longitud máxima del casco del buque de proa a popa (MOPU, 1990)
Manga: Anchura máxima de la sección transversal del casco del buque (MOPU, 1990)
(1) http://www.elmundo.es/diccionarios/
(2) http://www.librerianautica.com/diccionario_nautico.html
Denise Da Costa
– xxi –
Glosario de términos
Denise Da Costa
– xxii –
Introducción
INTRODUCCION
Las estructuras portuarias son una de las estructuras más importantes de un país
costero, dado que permiten la comunicación nacional e internacional de personas y
mercancía; contribuyendo al desarrollo y crecimiento del país. Por lo tanto, el correcto
diseño y análisis de este tipo de estructuras es de primordial importancia.
Por muchos que sean los cuidados que se tengan para realizar las maniobras de
atraque y amarre, debido a las condiciones extremas donde estas se encuentran, las
estructuras portuarias están expuestas a sufrir grandes daños si no se toman las
precauciones necesarias; cabe destacar que, al quedar fuera de servicio una estructura de
este tipo, las pérdidas serán mucho mayores que cualquier gasto adicional para protegerla y
asegurar su funcionamiento.
Por todo lo expuesto anteriormente surge la idea de crear un programa que permita la
determinación de las fuerzas generadas durante el amarre y defensa de un buque a un
muelle. Amarre es un programa que determina las fuerzas que ocurren en los cabos y
defensas bajo distintas condiciones de carga, permitiendo mejorar el criterio del ingeniero en
la toma de decisiones al momento de diseñar una estructura de atraque y amarre. El
programa Amarre utiliza el método de la rigidez aplicando el análisis matricial para obtener
las fuerzas y deformaciones del sistema de amarre incluyendo las defensas; para ello fue
necesario desarrollar y deducir las matrices y ecuaciones para el sistema específico: el
buque y sus elementos de amarre.
Los muelles son estructuras que sirven para el atraque de las embarcaciones. Una vez
que el barco llega al muelle es fijado a este con el fin de impedir que ocurran grandes
desplazamientos relativos entre ellos, para lograr esto se hace uso de un sistema de amarre
formado por cabos y defensas.
Denise Da Costa
–1–
Introducción
Uno de los aspectos a considerar para el correcto diseño de las estructuras de amarre y
atraque son las fuerzas que le transmiten los cabos y defensas que sujetan la embarcación.
El cálculo de estas fuerzas puede resultar largo y complejo. Primeramente se asume que el
sitio del muelle, el buque, y la configuración de amarre están dados antes del comienzo del
diseño detallado. Generalmente es necesario utilizar varias disposiciones de amarre porque
las cargas en los cabos y en las estructuras de atraque y amarre no son conocidas en la
primera etapa del diseño. El ingeniero proyectista tendrá que analizar varios tipos de buques
para una configuración de amarre dada, de modo a optimizar el diseño de las estructuras de
atraque y amarre.
Para facilitar el cálculo y el criterio del ingeniero al momento de diseñar una estructura
de amarre y atraque, se ha desarrollado en este Trabajo un programa denominado Amarre
que funciona en Visual Basic aplicado a Excel y permite obtener las fuerzas de amarre de
una manera rápida, sencilla y confiable.
1.3 OBJETIVOS
Denise Da Costa
–2–
Introducción
1.4 HIPÓTESIS
Las cargas inducidas por las olas pueden ser predominantes comparadas con las
fuerzas de viento y corriente en los casos donde el buque se encuentre amarrado en
ambientes desprotegidos. Mientras que en locales de amarre en áreas protegidas, como
puertos y muelles, estas fuerzas disminuyen, y el viento y la corriente comienzan a dominar.
En el presente trabajo, se asume que el buque se encuentra amarrado en un puerto o
muelle y por consiguiente las fuerzas predominantes en el análisis serán las fuerzas
estáticas de viento y corriente; en el capítulo 2 se presentan las ecuaciones para determinar
dichas fuerzas.
El análisis cuantitativo de fuerzas inducidas por las olas está más allá del alcance de
esta tesis; sin embargo, en el capítulo 3 se presentan algunos conceptos dinámicos de
interés como las ecuaciones de movimiento en varias situaciones, también se discuten las
dificultades de resolución, por vía estrictamente teórica, así como la magnitud, el carácter, y
la importancia relativa de cargas de oleaje.
Denise Da Costa
–3–
Introducción
En el capítulo 7 – Programa amarre – están los detalles del programa desarrollado para
realizar los cálculos en las fuerzas de amarre para buque amarrados en puertos; en este
capítulo se explican los principios del programa, cuales son los datos de entrada y de salida,
y el procedimiento de las diversas subrutinas que componen el algoritmo del programa.
Denise Da Costa
–4–
Fuerzas estáticas de viento y corriente
Las principales propiedades del viento y del agua son presentadas a continuación.
Denise Da Costa
–5–
Fuerzas estáticas de viento y corriente
Surge – En la dirección X
Sway – En la dirección Y
Heave – En la dirección Z
Roll – Rotación alrededor del eje X
Pitch – Rotación alrededor del eje Y
Yaw – Rotación alrededor del eje Z
– Figura 2.1 –
Perspectiva del buque con sus correspondientes movimientos (UFC, 2005)
Denise Da Costa
–6–
Fuerzas estáticas de viento y corriente
– Figura 2.2 –
Nomenclatura para las fuerzas de viento y corriente (NFEC, 1986)
Denise Da Costa
–7–
Fuerzas estáticas de viento y corriente
2.3 VIENTO
Existe una extensa variedad de métodos para determinar la fuerza del viento sobre un
buque amarrado. El procedimiento adoptado es tomado de Gregory P. Tsinker (1995)
conjuntamente con Naval Facilities Engineering Command (NFEC, 1986)
– Figura 2.3 –
Definición de términos (UFC, 2005)
Denise Da Costa
–8–
Fuerzas estáticas de viento y corriente
⋅ ρ a ⋅ Vw2 ⋅ Ay ⋅ C yw ⋅ f yw (θ w )
1
Fyw = (2.1)
2
donde,
El coeficiente de arrastre de la fuerza del viento depende del área del casco y de la
superestructura del buque, y es calculado utilizando la siguiente expresión
2 2
VS V
⋅ AS + H ⋅ AH
V VR
C yw =C⋅ R (2.2)
Ay
donde,
C yw = Coeficiente de arrastre lateral de la fuerza del viento.
C = Coeficiente empírico. (Ver tabla 2.2)
VS = Velocidad promedio del viento sobre la superestructura.
VR = Velocidad promedio del viento a 10 m sobre el nivel del mar (velocidad de referencia).
AS = Proyección longitudinal del área de la superestructura.
AH = Proyección longitudinal del área del casco (Hull; index: H).
Denise Da Costa
–9–
Fuerzas estáticas de viento y corriente
VS V
Los valores de y H son determinados usando las siguientes ecuaciones:
VR VR
1/ 7
V S hS
= (2.3)
V R hR
1/ 7
V H hH
= (2.4)
V R hR
donde,
VS
= Promedio normalizado de la velocidad del viento sobre la superestructura.
VR
hS = Altura promedio de la superestructura.
hR = Altura de referencia (10 m)
VH
= Promedio normalizado de la velocidad del viento sobre el casco.
VR
hH = Altura promedio del casco.
Buque C Ejemplos
De casco grande 0.82 Portaaviones.
Típico 0.92 Barcos con una superestructura de área moderada.
Con superestructura extensa 1.02 Cruceros, destroyers.
– Tabla 2.2 –
Coeficiente de viento para distintos tipos de buques (UFC, 2005)
sin (5 ⋅ θ w )
sin (θ w ) −
f yw (θ w ) = 20 (2.5)
1
1−
20
La ecuación (2.5) es positiva para ángulos comprendidos entre 0º < θ w < 180º y es
negativa para ángulos de viento entre 180º < θ w < 360º . En la figura 2.4 se muestra el
gráfico correspondiente a distintos valores de la función para distintos ángulos de viento.
Denise Da Costa
– 10 –
Fuerzas estáticas de viento y corriente
0 0
1,20
20 0,31 1,00
40 0,69 0,80
0,60
60 0,96
0,40
80 1,00
0,20
100 1,00
fyw
0,00
120 0,96 -0,20 0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360
-0,40
140 0,69
-0,60
160 0,31
-0,80
180 0 -1,00
260 -1,00
280 -1,00 Los buques amarrados en puertos y estructuras similares están
300 -0,96 resguardados (ie, locales donde hay muy poca agitación marítima) bien
sea por los propios cabos de amarre, bien sea por los edificios en el
320 -0,69
área de servicio del puerto y/o por estructuras portuarias similares. El
340 -0,31
ingeniero proyectista debe considerar las posibilidades de
360 0
abrigo/resguardo en la determinación del área proyectada del buque.
La fuerza longitudinal del viento es definida como la componente de la fuerza del viento
paralela a la línea central del buque y es determinada usando la siguiente ecuación
⋅ ρ a ⋅ Vw2 ⋅ Ax ⋅ C xw ⋅ f xw (θ w )
1
Fxw = (2.6)
2
donde:
Fxw = Fuerza longitudinal del viento.
ρa = Masa volumétrica del aire.
– 11 –
Fuerzas estáticas de viento y corriente
El coeficiente de arrastre longitudinal de la fuerza del viento varia dependiendo del tipo
y características de buque. Adicionalmente este coeficiente depende del ángulo del viento;
Cuando el ángulo θw = 0º se denomina condición de viento de cabeza (proveniente de la
C xwS
El ajuste de estos coeficientes para todos los buques de carga y petroleros con las
cubiertas obstaculizadas consiste en adicionar 0.08 para C xwB y C xwS
Denise Da Costa
– 12 –
Fuerzas estáticas de viento y corriente
f wx (θ w ) = − cos(φ ) (2.7)
donde,
90º
φ (− ) = ⋅ θ w Para θ w < θ wz
θ wz
⋅ (θ w − θ wz ) + 90º
90º
θ (+ ) = Para θ w > θ wz
180º − θ wz
forma.
sin 5 ⋅ γ
− sin γ −
f xw (θ w ) = 10 (2.8)
1
1−
10
donde,
90º
γ (− ) = ⋅ θ w + 90º Para θ w < θ wz
θw
90º 90º⋅θ wz
γ (+ ) = ⋅ θ w + 180º − Para θ w > θ wz
180º −θ w 180º −θ wz
⋅ ρ a ⋅ Vw2 ⋅ Ay ⋅ L ⋅ C xyw (θ w )
1
M xyw = (2.9)
2
Denise Da Costa
– 13 –
Fuerzas estáticas de viento y corriente
donde,
– Figura 2.5 –
Coeficientes de momento (NFEC, 1986)
θ ⋅ 180
C xyw (θ w ) = −a1 ⋅ sin w Para 0º < θ w < θ z (2.10a)
θz
C xyw (θ w ) = a 2 ⋅ sin (λ ⋅ (θ w − θ z )) Para θ z < θ w < 180º (2.10b)
Denise Da Costa
– 14 –
Fuerzas estáticas de viento y corriente
donde,
180º
λ=
180º − θ z
Tipo de buque θz a1 a2
Liner 80 0,075 0,14
Portador 90 0,068 0,072
Tanque 95 0,077 0,07
Crucero 90 0,064 0,05
Destroyer 68 0,02 0,12
– Tabla 2.5 –
Variables para determinar el momento debido al viento (UFC, 2005)
2.4 CORRIENTE
Denise Da Costa
– 15 –
Fuerzas estáticas de viento y corriente
donde,
Vc = Velocidad de la corriente.
LwL = Longitud de la línea de flotación del buque.
T = Calado promedio del buque.
C yc = Coeficiente transversal de fricción para corriente.
θc = Angulo de la corriente.
– Figura 2.6 –
Relación entre el calado del buque (T) y la profundidad del agua (d)
K
T
C yc = C 0 + (C1 − C 0 ) ⋅ (2.12)
d
Denise Da Costa
– 16 –
Fuerzas estáticas de viento y corriente
donde,
C 0 = Coeficiente de fuerza de resistencia para grandes profundidades. Cuando T/d ≈ 0,
este coeficiente es calculado según:
C 0 = 0,22 ⋅ x
L2wL ⋅ Am
x=
B ⋅V
Am = C m ⋅ B ⋅ T
donde,
LwL = Longitud del buque en la línea de agua.
Am = Área de la sección transversal sumergida en la sección media del buque.
B = Acho máximo del buque en la línea de flotación.
V = Volumen sumergido del buque.
C1 = Coeficiente de fuerza de resistencia para pequeñas profundidades. Para T/d = 1 y
corrientes de 1,5 m/s o menores, se recomienda C1 = 3,2
T = Calado promedio del buque.
d = Profundidad del agua.
k = Exponente adimensional.
k = 3 Para buques de carga, portaviones y submarinos
k = 2 Para los casos restantes.
Los valores del coeficiente C m son dados en las tablas características de cada buque.
El procedimiento descrito se refiere apenas a condiciones normales de diseño con
corrientes moderadas cuya velocidad no supere los 1,5 m/s (2,9 nudos) y buques atracados
en canales anchos o puertos.
donde,
Fxc = Fuerza estática longitudinal debida a la corriente.
Fx FORM = Fuerza estática longitudinal debido a la resistencia de forma.
Denise Da Costa
– 17 –
Fuerzas estáticas de viento y corriente
donde,
ρw = Masa volumétrica del agua.
Vc = Velocidad de la corriente.
B = Ancho máximo del buque en la línea de flotación.
T = Calado promedio del buque.
C xcb = Coeficiente longitudinal. = 0,1
θc = Angulo de la corriente.
donde,
S = Área de la superficie mojada.
C xca = Coeficiente de fricción lateral en la dirección longitudinal.
D
S = 1,7 ⋅ T ⋅ LwL +
T ⋅γ w
LwL = Línea de flotación.
D = Desplazamiento del buque.
γw = Densidad por peso del agua.
0,075
C xca =
((log10 RN ) − 2)2
R N = Número de Reynolds.
Denise Da Costa
– 18 –
Fuerzas estáticas de viento y corriente
Vc ⋅ LwL ⋅ cos(θ c )
RN =
ν
ν = Viscosidad cinemática del agua.
donde,
AP = Área de las láminas de la hélice.
C PROP = Coeficiente de resistencia de la hélice. C PROP = 1
ATPP
AP =
0,838
ATPP = Área total proyectada de la hélice.
LwL ⋅ B
ATPP =
AR
AR = Área adimensional para hélices. En la tabla 2.6 se presentan algunos valores típicos.
Buque AR
Destroyer 100
Crucero 160
Portador 125
Carga 240
Tanque 270
Submarino 125
– Tabla 2.6 –
Valores de AR para distintos tipos de buques (UFC, 2005)
Denise Da Costa
– 19 –
Fuerzas estáticas de viento y corriente
e
M xyc = Fyc ⋅ LwL ⋅ c (2.17)
LwL
donde,
ec
= Cociente entre la excentricidad y la longitud de la línea de flotación.
LwL
ec = excentricidad de Fyc
ec
= a + b ⋅θ c Para 0º < θ c < 180º (2.18a)
LwL
ec
= a + (b ⋅ (360º −θ c )) Para 180º < θ c < 360º (2.18b)
LwL
donde,
a = Intercepción con el eje Y (ver tabla 2.7)
b = Pendiente por grado (ver tabla 2.7)
Buque a b Notas
Series 60 -0,291 0,00353 Buques de carga.
FFG -0,201 0,00221
CVE-55 -0,168 0,00189
SS-212 -0,244 0,00255 Submarinos.
– Tabla 2.7 –
Variables para determinar el cociente de excentricidad (UFC, 2005)
Denise Da Costa
– 20 –
Conceptos dinámicos importantes
El siguiente apartado fue desarrollado con base en la teoría presentada por Wilson,
1984.
Se asume que la estructura o parte de ella es rígida, y posee una masa m la cual es
significativamente mayor que la masa de los demás restantes elementos flexibles que
limitan su movimiento. El movimiento será restringido al plano donde el desplazamiento
∑F v = m ⋅ v&& (3.1)
∑M O = J O ⋅ θ&& (3.2)
J O = ∫ r 2 dm (3.3)
m
Denise Da Costa
– 21 –
Conceptos dinámicos importantes
tiene que:
J O = J G + m ⋅ hG2 (3.4)
El tipo de fuerzas que constituyen las sumatorias en el lado izquierdo de las ecuaciones
(3.1) y (3.2) son clasificadas en tres categorías. La primera categoría son las fuerzas
aplicadas; incluyendo peso propio, las fuerzas de viento, oleaje y corriente como el resto de
las fuerzas ambientales; las fuerzas externas aplicadas dependientes del tiempo en la
dirección v son denotadas como P1 (t ) y el momento total dependiente del tiempo alrededor
de un punto fijo O, debido a las fuerzas externas aplicadas en el sistema descrito por la
coordenada θ , es denotado como M O (t ) . La segunda categoría son las reacciones, estas
fuerzas de resorte o restauradoras ocurren cuando existe un apoyo flexible en el modelo de
cuerpo rígido; si el movimiento es únicamente descrito por la coordenada v , entonces la
reacción es denotada como una función q (v) en esta coordenada; si el movimiento es
únicamente descrito por la coordenada θ, entonces el momento de reacción actuante
alrededor del punto fijo O es denotado como q (θ ) . En la tercera categoría está la fuerza de
fricción y/o de amortecimiento que es asumida como dependiente solamente de la
velocidad; análogamente con la fuerza de resorte, las respectivas fuerzas de
Denise Da Costa
– 22 –
Conceptos dinámicos importantes
3.2 OLAS
Esta acción lleva al encrespamiento de la superficie del mar. El viento ejerce, de esta
forma, una fuerza de empuje en la parte posterior de la ola y una fuerza de succión en la
parte frontal. Estas fuerzas son mayores como mayor sea la ola, dado que son
proporcionales al área de contacto. En el ciclo descrito, más viento implica una ola mayor y
una ola mayor genera más superficie aumentando la acción del viento; únicamente termina
cuando por efecto de la gravedad de la masa de agua, las fuerzas se equilibran; asumiendo
que la velocidad del viento es constante durante un período de tiempo suficientemente
grande.
Por lo expuesto anteriormente se puede concluir que las características de las olas
dependen esencialmente de la velocidad del viento que las genera, pudiendo establecer una
relación entre los dos fenómenos.
Las partículas de agua que forman una ola no experimentan movimiento de translación,
pero si de oscilación describiendo una trayectoria aproximadamente circular. La posición
relativa de una partícula que se encuentra en la superficie, en relación con las adyacentes,
varía con el tiempo, ocupando una posición nueva en el círculo que describe. La forma
creada por la superficie del agua se llama perfil de la ola; este perfil varia y son estas
variaciones que dan la sensación de traslación del agua, es decir, que esta avanza en la
misma dirección que la ola. En la realidad lo que se mueve es la forma de la ola sin que las
Denise Da Costa
– 23 –
Conceptos dinámicos importantes
– Figura 3.1 –
Representación esquemática de olas.
gT 2
forma: λ=
2π
Esta descripción es válida para aguas profundas, donde la interacción con el fondo
marino no afecta el movimiento del agua. Cuando la profundidad, d , es menor que la mitad
λ
de la longitud de la ola, d< , la fricción con el fondo marino causa una fuerza que se
2
Denise Da Costa
– 24 –
Conceptos dinámicos importantes
opone a la fuerza que general a ola; por esto, la altura y longitud de la ola disminuyen y la
trayectoria que describen las partículas es distorsionada, pasando de circular a elíptica.
Disminuyendo la profundidad del agua, las fuerzas de fricción con el fondo generan un
desequilibrio de manera que las partículas que ocupan una posición más elevada trasponen
a las que se encuentran más enfrente dando lugar al rebatimiento de la ola y a la liberación
de una gran cantidad de energía.
Para determinar la acción de las olas sobre una estructura es necesario determinar los
parámetros de cálculo de la ola y después calcular el efecto de la ola en la estructura.
Las diversas medidas de olas en varios locales permiten deducir las distribuciones de
frecuencias de los parámetros de la ola y con estos determinar los valores que tienen el
período de ocurrencia deseado.
La siguiente tabla relaciona alturas y longitudes de olas con la velocidad del viento que
las genera. Los valores presentados son válidos para aguas profundas (Minikin, 1963).
Denise Da Costa
– 25 –
Conceptos dinámicos importantes
Se pueden distinguir dos casos extremos del efecto de las olas en la mayoría de las
estructuras marítimas, siendo los siguientes:
Choque de la ola en una pared vertical rígida provocando el quiebre de la ola, como
en los rompe olas y muros costeros.
Efecto de la ola en elementos estructurales relativamente pequeños y que no
alteran el movimiento de la ola, como son las estacas.
La figura 3.2 representa una perspectiva de los métodos que describen las olas
superficiales gravíticas, resaltando que existen dos abordajes para el problema.
El abordaje determinístico revela ser más útil en la descripción de los efectos de corta
duración y consiste en métodos que pueden ser descritos como analíticos. Los métodos
analíticos se basan, básicamente, en las teorías clásicas de las olas y pueden ser divididos
en lineales y no lineales. Por otro lado, el abordaje probabilístico, al cual no se hará
referencia, es más útil en la previsión del comportamiento de la ola a largo plazo.
– Figura 3.2 –
Esquema de las teorías de ola
Denise Da Costa
– 26 –
Conceptos dinámicos importantes
Las olas planas, o olas con gran longitud, son olas independientes en la coordenada
perpendicular a la dirección de propagación de la ola. De esta forma, las crestas y valles
pueden ser representados con líneas rectas en el plano horizontal.
En la figura 3.3 se representa una ola con forma armónica simple, tal que para cada
instante de tiempo, t , tiene una forma sinusoidal referente al nivel medio del agua.
– Figura 3.3 –
Representación esquemática de una ola.
ηt = 0 = A cos(kx ) (3.7)
donde,
H
A= , es la amplitud de la ola.
2
k es el número de ola que transforma el término kx a radianes.
2π
k= (3.8)
λ
Denise Da Costa
– 27 –
Conceptos dinámicos importantes
La expresión general de η para una ola sinusoidal con movimiento progresivo con
η ( x, t ) = A cos k ( x − ct ) (3.9)
– Figura 3.4 –
Variación con el tiempo de la elevación instantánea de la superficie
donde T es el período de la ola y representa el valor particular del tiempo que hace al
coseno igual a 1, dando el tiempo que demora la próxima cresta a llegar
2π
kc = = ω [rad / s ] (3.12)
T
donde ω es la frecuencia angular de la ola. Con las ecuaciones (3.8) y (3.11) se obtiene:
λ = cT (3.13)
Denise Da Costa
– 28 –
Conceptos dinámicos importantes
η ( x, t ) = A cos(kx − ωt ) (3.14)
La teoría de las olas lineal es también conocida por la teoría de Airy o teoría de primer
orden; fue desarrollada por Airy (1845) y Laplace (1816), es la más importante de las teorías
clásicas porque es basada en los espectros de olas de las teorías probabilísticas.
Considerando una ola plana, simple y armónica, definida en la figura 3.3 que se
propaga en el sentido positivo de x. Se presentan a continuación las nueve suposiciones de
la teoría linear.
∂u ∂w
− =0 (3.15a)
∂z ∂x
Denise Da Costa
– 29 –
Conceptos dinámicos importantes
∂u ∂w
+ =0 (3.15b)
∂x ∂z
∂u 1 ∂p
=− (3.15c)
∂t ρ ∂x
∂w 1 ∂p
=− −g (3.15d)
∂t ρ ∂z
∂η
w= , para z = 0 (3.16)
∂t
w = 0 , para z = − d (3.17)
p = pa , para z = 0 (3.18)
Tomando:
η = A cos(kx − ωt ) (3.19)
Quedaría
cosh k ( z + d )
u = Aω cos(kx − ωt )
sinh (kd )
(3.20)
sinh k ( z + d )
w = Aω sin (kx − ωt )
sinh (kd )
(3.21)
cosh k ( z + d )
p = pa − ρgz + ρgA cos(kx − ωt )
cosh (kd )
(3.22)
Denise Da Costa
– 30 –
Conceptos dinámicos importantes
Las ecuaciones (3.19) a (3.22) satisfacen las ecuaciones (3.15) a (3.17). Estas
ecuaciones también satisfacen las condiciones (3.16) y (3.18), llevando a que el número de
ola y frecuencia satisfagan la relación:
ω λ
c= = (3.24)
k T
se tiene que:
tanh (kd )
g
c2 = (3.25)
k
Parámetro Fórmula
4π 2 A cosh k ( z + d )
Aceleración horizontal de la partícula u& = sin (kx − ωt )
T2 sinh (kd )
4π 2 A sinh k ( z + d )
Aceleración vertical de la partícula w& = cos(kx − ωt )
T2 sinh (kd )
cosh k ( z + d )
p = ρgA cos(kx − ωt )
cosh (kd )
Presión dinámica
1
g 2
Celeridad de la ola c = tanh (kd )
k
1 2kd
cg = c ⋅ 1 +
2 sinh (2kd )
Velocidad de grupo
– Tabla 3.2 –
Cuadro resumen de la teoría lineal.
Denise Da Costa
– 31 –
Conceptos dinámicos importantes
La explicación teórica de esta teoría pierde interés en los objetivos de este trabajo;
sirviendo apenas como una introducción histórica de otras teorías no lineales.
De este modo, la teoría Trocoidal sirve de ligación con la teoría de olas oscilatorias de
amplitud finita desarrollada por Stokes (1845), Levi-Cevite (1926), Struik (1926), Havelock
(1914) entre otros (Myers et al., 1969; Gaythwaite, 1981; Sarpkaya & Isaacson, 1981).
De esta forma, las ecuaciones que describen el movimiento pueden ser escritas de la
misma forma que para la teoría lineal.
∂w ∂u
− =0 (3.26a)
∂x ∂z
∂u ∂u
+ =0 (3.26b)
∂x ∂z
∂u ∂u ∂u 1 ∂p
+u +ω =− (3.26c)
∂t ∂x ∂z ρ ∂x
∂w ∂w ∂w 1 ∂p
+u +ω =− −g (3.26d)
∂t ∂x ∂z ρ ∂z
Denise Da Costa
– 32 –
Conceptos dinámicos importantes
∂p ∂p ∂p
+u +w =0 (3.28)
∂t ∂x ∂z
Denise Da Costa
– 33 –
Conceptos dinámicos importantes
Parámetro Fórmula
H 2 π cosh kd
cos(kx − ωt ) +
H
η= ×
Perfil de la superficie 2 8 λ sinh 3 kd
× [2 + cosh (2kd ) ⋅ cos[2(kx − ωt )]]
π cosh k ( z + d )
u=H cos(kx − ωt ) +
Velocidad horizontal de
T sinh (kd )
la partícula 3 2 π 2 cosh k ( z + d )
+ H cos[2(kx − ωt )]
4 Tλ sinh 4 (kd )
π sinh k ( z + d )
w = H sin (kx − ωt ) +
Velocidad vertical de la T sinh (kd )
partícula 3 2 π 2 sinh 2k ( z + d )
+ H sin [2(kx − ωt )]
4 Tλ sinh 4 (kd )
2π 2 H cosh k ( z + d )
u& = sin (kx − ωt ) +
Aceleración horizontal de T2 sinh (kd )
la partícula π 3 cosh 2k ( z + d )
+ 3H 2 2 sin[2(kx − ωt )]
T λ sinh (kd )
4
2π 2 H sinh k ( z + d )
w=
& sin (kx − ωt ) +
Aceleración vertical de la T2 sinh (kd )
partícula π 3 sinh 2k ( z + d )
+ 3H 2 2 cos[2(kx − ωt )]
T λ sinh (kd )
4
H cosh k ( z + d )
p = ρg cos(kx − ωt ) +
2 cosh (kd )
Presión dinámica 3 2 ρ gπ 1 cosh 2k ( z + d ) 1
+ H × − ×
λ sinh (2kd ) sinh kd
2
4 3
× cos[2(kx − ωt )]
1
g 2
Celeridad de la ola c = tanh (kd )
k
– Tabla 3.3 –
Cuadro resumen. Teoría de amplitud finita de Stokes.
Las expresiones de orden superior son simplemente aquellas en las cuales las
aproximaciones de los efectos de los términos correctivos son desarrolladas hasta el
término de orden correspondiente. En principio la teoría de Stokes es usada en un orden
suficientemente elevado, debe ser adecuada en la descripción de olas para cualquier valor
de profundidad. En la práctica esto es solamente posible para olas en aguas profundas. En
aguas poco profundas los términos correctivos son extensos, provocando así que la
Denise Da Costa
– 34 –
Conceptos dinámicos importantes
convergencia sea lenta; es necesario un gran número de términos para obtener un grado de
precisión uniforme.
η0 λ
3
Ur = (3.29)
λ d
Existen otras formulaciones clásicas que necesitan menos términos para obtener la
precisión deseada, pero se escapan de los objetivos de este trabajo.
El problema que muchas veces se presenta es cual es la mejor teoría para una dada
situación. El gráfico propuesto por Le Mehaute (1976) es bastante práctico y cubre toda la
gama de olas de forma permanentemente posibles, a pesar de que es basado en
consideraciones subjetivas. El gráfico de Le Mehaute es bastante útil en los cálculos
preliminares dado que apunta inmediatamente para la aplicabilidad de determinada teoría
(Sarpkaya & Isaacson, 1981).
En la figura 3.6 se presenta un gráfico que resulta de la sobre posición del gráfico de Le
Mehaute con el de Dean (1974) y es actualmente uno de los más aceptados para la
verificación de la validad de determinada teoría de olas. Puede fácilmente realizarse una
comparación con la figura 3.5, teniendo en cuenta que se encuentran en sistemas de
unidades diferentes.
Denise Da Costa
– 35 –
Conceptos dinámicos importantes
– Figura 3.5 –
Límites de validad según Le Mehaute.
– Figura 3.6 –
Límite de validad de las teorías de olas de Le Mehaute sobrepuesta a la propuesta por Dean
Denise Da Costa
– 36 –
Conceptos dinámicos importantes
Una de las formas utilizadas para estimar las fuerzas máximas de olas sobre las
estructuras es hacer los cálculos basados en una única ola de diseño. Para una teoría de
ola en particular, con una altura y período de ola seleccionados de acuerdo con la
localización de la estructura, se determina el campo de presiones y la componente
horizontal de la velocidad y aceleración de la partícula. Con estas informaciones de
movimiento, se determina la distribución de los dos parámetros que gobiernan el
movimiento, Re (número de Reynolds) y KC (número de Keulegan-Carpenter) para las
varias componentes estructurales; el régimen de escurrimiento es determinado; los
coeficientes de arrastre, inercia y difracción son convenientemente escogidos (Sarpkaya &
Isaacson, 1981; Gaythwaite, 1981).
πD 2 γ 1 γ
f = CM ⋅ ⋅ ⋅ u& + CD ⋅ ⋅ D ⋅ u ⋅ u (3.30)
4 g 2 g
En el caso en que la ola pasa alrededor de un obstáculo, el desvío local del flujo de
agua produce efectos muy semejantes a los del viento; es decir, el flujo se separa del
obstáculo creando áreas de arrastre y de succión dando lugar a una fuerza de arrastre
global que es función del área expuesta, de la densidad del fluido y del cuadrado de la
velocidad, así como de un coeficiente que depende, a su vez, de la forma del objeto (Barros,
1999-b).
1 γ
fD = ⋅ CD ⋅ ⋅ D ⋅ u ⋅ u (3.30a)
2 g
donde,
C D = coeficiente de arrastre.
γ = peso volumétrico del agua.
g = aceleración de la gravedad.
D = diámetro del tubo o eje de la sección expuesta al flujo
u = velocidad horizontal para la cual tendería la partícula de fluido en el punto en estudio si
no existiese el obstáculo.
Denise Da Costa
– 37 –
Conceptos dinámicos importantes
γ πD 2
fM = C M ⋅ ⋅ u& (3.30b)
g 4
donde,
γ πD 2 γ πD 2 γ πD 2
fM = Cm ⋅ ⋅ u& + ⋅ u& = (Cm + 1) ⋅ u&
g 4 g 4 123 g 4 (3.31)
14 4244 3 14243 CM
Fuerza de masa hidrodinám ica Fuerza de Froude - Krylov
∂
∑F X =
∂t ∫ V
ρ u dV + ∫ ρ u ⋅ u dA0
A0
(3.32)
Denise Da Costa
– 38 –
Conceptos dinámicos importantes
Definiendo q como la fuerza de reacción por unidad de longitud del cilindro, se tiene
que ∑F X = − q ⋅ ∆z . Cabe destacar que el signo negativo representa que la fuerza es de
reacción (Sarpkaya & Isaacson, 1981).
– Figura 3.7 –
Esquema de aplicación de las fuerzas.
∂
∫ ρ u dV ≅ − ρD 2 ∆z u& (3.33)
∂t V
El signo negativo es utilizado en este caso dado que el fluido acelera negativamente
dentro del volumen que rodea al disco. La segunda integral representa el flujo de momento
en x, o la diferencia entre el momento exterior y el momento interior a través del área
A0 = D ⋅ ∆z . Así:
∫ ρ u ⋅ u dA = ∫
0
ext
ρ u ⋅ u dA0 − ∫ ρ u ⋅ u dA0 = 0 − ρ u ⋅ u ⋅ D ⋅ ∆z
int
(3.34)
− q ⋅ ∆z ≅ − ρ ⋅ D 2 ⋅ ∆z ⋅ u& − ρ ⋅ D ⋅ u ⋅ u ⋅ ∆z
⇔ (3.35)
q ≅ ρ ⋅ D ⋅ u& + ρ ⋅ D ⋅ u ⋅ u
2
Denise Da Costa
– 39 –
Conceptos dinámicos importantes
3.3 SEICHE
La superficie del agua exhibe una serie de nodos y de antinodos con respecto a la
columna del agua. Los antinodos son regiones donde el movimiento vertical es máximo y las
velocidades horizontales son mínimas. Donde la longitud de la onda es suficientemente
mayor que la longitud del buque, el buque atracado en el antinodo experimentará una
subida apacible y caerá con el período de onda estacionaria. En el nodo, el buque estará
conforme a una acción de movimiento longitudinal (surge) horizontal periódica debido a las
corrientes. La figura 3.8 ilustra un sistema unidimensional de onda estacionaria típico
(NFEC, 1981).
– Figura 3.8 –
Sistema unidimensional de onda estacionaria.
Denise Da Costa
– 40 –
Conceptos dinámicos importantes
El problema general consiste en solucionar los seis movimientos acoplados del buque
amarrado bajo la carga del oleaje y comportamiento no lineales de los cabos y de defensas
que lo sujetan. Cuando un buque está sometido a la acción de seiche, aunque la amplitud
de seiche es generalmente baja, el buque experimenta un mayor movimiento en la dirección
longitudinal (surge). Así en la práctica solamente los movimientos surge (longitudinal al eje)
y sway (transversal al eje) son normalmente importante para los propósitos del diseño de la
amarradura (NFEC, 1981).
FW + FD + Fm + FI = 0 (3.36)
donde,
FW = Fuerza de la ola (por diferencia de presión)
FD = Fuerza de fricción.
Fm = Fuerza de amarre.
FI = Fuerza de inercia.
Denise Da Costa
– 41 –
Conceptos dinámicos importantes
Los resultados de tal análisis se pueden presentan en una figura similar a la figura 3.9; en
esta figura la amplitud de movimiento de la oleada del buque se traza contra el período de
seiche para una amplitud dada. Según lo demostrado en la figura, el buque puede resonar
en los períodos críticos del seiche, causando movimientos grandes en oleada. Después de
la determinación de la magnitud de movimiento de la oleada, se pueden cuantificar los
daños posibles del sistema, buque y cabos.
– Figura 3.9 –
Amplitud del movimiento Surge de un buque sujeto a la acción seiche
(NFEC, 1981)
Una aproximación para la ecuación básica del movimiento desacoplado surge fue
presentada por Barros (1997) y proviene de la segunda ley de Newton
donde,
&x& = Aceleración del centro de masa del buque.
m = masa del buque.
f s = Fuerza longitudinal no lineal (rigidez longitudinal) para cabos de amarre de popa y
proa.
fs = c ⋅ xn
donde,
c = Coeficiente de rigidez no lineal.
x = Traslación en la dirección surge del centro de masa.
n = Exponente positivo.
f d = Fuerza amortiguadora en la dirección longitudinal.
f d = b ⋅ (u − x& )
Denise Da Costa
– 42 –
Conceptos dinámicos importantes
donde,
b = Coeficiente del amortiguamiento por viscosidad.
(u − x& ) = Velocidad relativa.
c m = Coeficiente de masa añadida del casco sumergido.
0,26 ⋅ π ⋅ B
cm =
L
donde,
B es el ancho del buque.
L es la longitud del buque entre perpendiculares.
donde, a es la masa virtual total (masa del buque más masa añadida), expresada como
a = (1 + c m ) m
De acuerdo con Barros (1997) la ecuación no lineal (3.38) puede ser expresada con el
U ⋅ cm
E [x(t )] = &x& + 2 ⋅ D ⋅ ω ⋅ x& + ω 2 ⋅ x n + ⋅ ω S ⋅ sin (ω S t ) − 2 ⋅ D ⋅ ω ⋅ U ⋅ cos(ω S t ) = 0
1 + cm
(3.39)
donde,
b
D=
2 ⋅ a ⋅ω
ω 2 = c a : Frecuencia no lineal pseudo natural.
Denise Da Costa
– 43 –
Conceptos dinámicos importantes
donde, X V es la amplitud del buque provocada por el movimiento surge inducido por
seiche y ε es el ángulo de fase.
U
2
cm 2
2
(X n −1
φn − η ) + (2 ⋅ D ⋅η )
2 2 2
= (2 ⋅ D ) +
2
η (3.41)
ω ⋅ XV 1 + cm
V
Para un determinado sistema de amarre, el análisis del movimiento en surge puede ser
aproximado mediante la siguiente metodología seguida por Barros (1997, 2001-b):
el valor absoluto de las amplitudes X V es ploteado como una función del parámetro de
frecuencia η . El parámetro de la curva de respuesta es la amplitud de seiche S .
Denise Da Costa
– 44 –
Conceptos dinámicos importantes
– Figura 3.10 –
Curva de respuesta para un tanque de aceite en situación standard (Barros, 1997)
Es sabido que para periodos pequeños del buque, corresponden altas velocidades en
las partículas de la onda debido a las grandes amplitudes de seiche. Para pequeñas
amplitudes de seiche, el periodo del buque es afectado por la velocidad de las partículas de
onda larga, aunque la resonancia en surge ocurra, sus consecuencias son menos
importantes porque las fuerzas en los cabos de amarre son pequeñas (Barros, 1997; 2001-
b).
La figura 3.11 ilustra la curva de respuesta para el movimiento surge de un buque para
las amplitudes de seiche correspondiente a los periodos T = 20 y 40 segundos. En el gráfico
se incluye también el valor límite para la amplitud en surge del buque (2,40 m) en donde
probablemente ocurre la ruptura de los cabos de amarre. Una característica evidenciada en
los diagramas es la respuesta cíclica debido al aumento y disminución de las amplitudes de
seiche, posibles para cada periodo de seiche.
– Figura 3.11 –
Respuesta cíclica para periodos peligrosos de sieche (Barros, 1997)
Denise Da Costa
– 45 –
Conceptos dinámicos importantes
– Figura 3.12 –
Seis grados de libertad de un cuerpo flotante.
Un análisis completo del movimiento del buque requiere la solución simultánea de seis
ecuaciones. Dado que el buque no posee doble simetría, el plano de flotación es asimétrico
transversalmente y el buque tiende a rotar surgiendo así la necesidad de estudiar de forma
acoplada ambos movimientos, heave y pitch. Estos últimos son los únicos movimientos
acoplados del sistema porque para los movimiento sway, surge y yaw no hay
descompensación hidrostática, es decir, para estos movimientos no hay variación del centro
de masa del buque.
El método de análisis utilizado en esta sección está basado en la teoría presentada por
Mc. Cormick conocido como la teoría de tiras o elemento diferencial (strip theory).
La suposición básica de esta teoría es que el flujo dependiente del tiempo adyacente a
la superficie flotante de las secciones verticales de un cuerpo flotante es bidimensional.
Entonces se asume que las fuerzas y momentos resultantes en estas secciones o
elementos diferenciales se pueden sumar a lo largo de la longitud total del cuerpo para
obtener las fuerzas y momentos totales. Además, estas dos suposiciones se juntan con la
hipótesis de que los movimientos de un cuerpo no alteran los movimientos de la partícula en
la ola superficial, aunque los movimientos de la partícula influencian los movimientos del
cuerpo.
Denise Da Costa
– 46 –
Conceptos dinámicos importantes
Un cuerpo flotante, como el ilustrado en la figura 3.12, está sujeto a las siguientes
fuerzas: F1 , reacción de inercia del fluido; F2 , fuerza de restauración hidrostática; F3 ,
fuerza de amortiguamiento; y F4 , fuerza inducida por la onda superficial. La fuerza F4 es de
hecho un caso especial de las fuerzas de inercia e hidrostáticas y actúa como la función de
fuerzas del sistema.
d 2z
m = F1 + F2 + F3 + F4 (3.42)
dt 2
d 2θ
IY = M1 + M 2 + M 3 + M 4 (3.43)
dt 2
donde,
I Y = Momento de inercia
M i = Momento en pitch debido a la fuerza Fi
Puesto que el líquido adyacente al buque está acelerado, la reacción de inercia del
fluido actúa como una fuerza externa de la siguiente manera:
= − [m w (ξ )wb ]
dF1 D
dξ Dt
∂ ∂ ∂ ∂
= − [m w (ξ )wb ] − (− Vi + v j + wk ) ⋅ i + j + k [mw (z& + ξθ − Vθ )]
dF1
dξ ∂t ∂ξ ∂y ∂z
dF1
dξ
( )
= − m w &z& + ξθ&& − 2Vθ& + V
dm w
∂ξ
(
z& + ξθ& − Vθ ) (3.44)
donde,
m w (ξ ) = Masa añadida por unidad de longitud para la sección transversal del cuerpo.
ξ ≈ ξ&& ≈ 0
Denise Da Costa
– 47 –
Conceptos dinámicos importantes
El momento alrededor del centro de gravedad del cuerpo flotante resultado de la fuerza
de movimiento inducido es:
dM 1 dF
=ξ 1 (3.45)
dξ dξ
L/2
dF1
F1 = ∫ dξ
−L / 2
dξ
L/2 L/2 L/2
dmw
F1 = − ∫ m w dξ &z& − ∫ ξ m w dξ θ&& + V ∫ ξ + 2m w dξ θ& (3.46)
−L / 2 −L / 2 −L / 2
dξ
Y,
L/2
dF1
M1 = ∫ ξ dξ
−L / 2
dξ
L/2 L/2 L/2
dm w
M 1 = − ∫ ξ m w dξ &z& − ∫ ξ mw dξ θ&& + V
2
∫ ξ dξ z& +
−L / 2 −L / 2 −L / 2
dξ
(3.47)
dmw
L/2 L/2
dmw
+ V ∫ ξ 2 + 2ξ mw dξ θ& − V 2 ∫ ξ dξ θ
−L / 2 d ξ −L / 2
d ξ
Para ilustrar el uso de las ecuaciones (3.46) y (3.47), considere una esfera medio
sumergida como se muestra en la figura 3.13. La tira (o el elemento diferencial) tiene la
forma de un disco semicilíndrico, cuya expresión para masa añadida es:
B′
2
m w (ξ ) = ρ ⋅ π
1
2 2
Denise Da Costa
– 48 –
Conceptos dinámicos importantes
2
L
donde, según la figura 3.13(c), B ′(ξ ) = 2 − ξ 2
2
Así,
1 L 2
mw = ρ ⋅ π − ξ 2
2 2
Y,
dm w
= −ρ ⋅ π ⋅ ξ
dξ
– Figura 3.13 –
Esfera semi sumergida.
1 1
F1 = − ⋅ π ⋅ ρ ⋅ L3 ⋅ &z& + ⋅ V ⋅ π ⋅ ρ ⋅ L3 ⋅ θ&
12 12
Denise Da Costa
– 49 –
Conceptos dinámicos importantes
1 L2
M1 = − ⋅ ρ ⋅ π ⋅ L3 θ&& + V z& − V 2 ⋅ θ
12 20
dF2
= − ρ ⋅ g ⋅ B ′(ξ )ζ (3.48)
dξ
dF2
M2 = ξ (3.49)
dξ
L/2 L/2
F2 = − ρ ⋅ g ∫ B ′(ξ ) dξ z − ρ ⋅ g ∫ ξ B′(ξ )dξ θ (3.50)
−L / 2 −L / 2
L/2
dF2
M2 = ∫ ξ dξ
−L / 2
dξ
L/2 L/2
M 2 = −ρ ⋅ g ∫ ξ B ′(ξ ) dξ z − ρ ⋅ g ∫ξ B ′(ξ ) dξ θ
2
(3.51)
−L / 2 −L / 2
Denise Da Costa
– 50 –
Conceptos dinámicos importantes
π ⋅ ρ ⋅ g ⋅ L2 ⋅ z
F2 = −
4
1
M 2 = − π ⋅ ρ ⋅ g ⋅ L4 ⋅ θ
64
3.4.3 Amortiguamiento
– Figura 3.14 –
Onda de amortiguamiento.
ρ ⋅ g ⋅ a d2 ⋅ c g
E& dξ = dξ (3.52)
2
Denise Da Costa
– 51 –
Conceptos dinámicos importantes
c ωn ω n2
tiene que c g = = en que, k→
2 2⋅k g
ρ ⋅ g 2 ⋅ a d2
E& dξ = dξ (3.53)
4 ⋅ ωn
4 ⋅ Ab 2 ⋅ Ab ⋅ ω n
dFd wb = dFd = dFd (3.54)
T π
donde,
T = Período del cuerpo y de la ola.
wb = Velocidad promedio.
Las amplitudes de la ola y del cuerpo no son necesariamente iguales.
Dado que la energía perdida por el cuerpo debe ser igual a la energía ganada por el
agua, las expresiones (3.53) y (3.54) pueden ser utilizadas para expresar la igualdad
donde el factor 2 es debido a los dos conjuntos de olas generadas (incidente y transmitida).
= = (3.56)
dξ wb 4 ⋅ ω n2
Denise Da Costa
– 52 –
Conceptos dinámicos importantes
dN π 2 ⋅ ρ ⋅ g 2 (a d Ab )
2
dFd 1
= = (3.57)
dξ wb dξ 8 ⋅ ω n3
dF3 dN
=− wb (3.58)
dξ dξ
dM 3 dF dN
= ξ 3 = −ξ wb (3.59)
dξ dξ dξ
L/2
dN
F3 = − ∫ wb dξ
−L / 2 dξ
π 2 ⋅ρ ⋅ g2 L/2
F3 = − ∫ A (ξ ) ⋅ w dξ
2
(3.60)
8 ⋅ ω n3
b
−L / 2
L/2
dN
M3 = − ∫ ξ ⋅ wb dξ
−L / 2 dξ
π 2 ⋅ρ ⋅ g2 L/2
M3 = − ∫ ξ ⋅ A (ξ ) ⋅ w dξ
2
(3.61)
8 ⋅ ω n3
b
−L / 2
donde,
ad
A (ξ ) = : Razón entre la amplitud de la onda y la amplitud de movimiento del elemento
Ab
diferencial.
Denise Da Costa
– 53 –
Conceptos dinámicos importantes
π 2 ⋅ρ ⋅ g2 2
F3 = − A ⋅ L( z& − Vθ ) (3.62)
8 ⋅ ω n3
π 2 ⋅ρ ⋅ g2 2 3 &
M3 = − A ⋅Lθ (3.63)
96 ⋅ ω n3
Las fuerzas actuantes en un cuerpo flotante que son inducidas por las olas son de
naturaleza inercial e hidrostática, sin embargo, difieren de las descritas en los apartados
3.4.1 y 3.4.2 porque en sus expresiones está explicita de dependencia del tiempo (son
independientes del movimiento del cuerpo)
– Figura 3.15 –
Cuerpo flotante.
∂φ
p w ds dξ = − ρ ds dξ (3.64)
∂t
Denise Da Costa
– 54 –
Conceptos dinámicos importantes
′
dF4 a B / 2
= ∫ p w dy
dξ − B′ / 2
′
B /2
dF4 a
= ρ ⋅ a w ⋅ g ⋅ cos(kξ − wt ) ∫ e kz dy (3.67a)
dξ − B′ / 2
La segunda fuerza inducida por la ola es el cambio hidrostático debido al paso de la ola,
siendo expresada por la forma
dF4b
= ρ ⋅ g ⋅ B ′ ⋅η
dξ
dF4b
= ρ ⋅ g ⋅ B ′ ⋅ a ⋅ cos(kξ − ωt ) (3.67b)
dξ
La fuerza total inducida por la ola en el elemento diferencial es obtenida por la suma de
las ecuaciones (3.67a) y (3.67b):
dF4
=
d
(F4 a + F4b )
dξ dξ
B′ / 2 kz
= ρ ⋅ g ⋅ a w ∫ e + B ′ [cos(kξ ) ⋅ cos(ωt ) + sin (kξ ) ⋅ sin (ωt )]
dF4
(3.68)
dξ − B′ / 2
Denise Da Costa
– 55 –
Conceptos dinámicos importantes
dM 4 dF
=ξ 4 (3.69)
dξ dξ
L/2
B′ / 2 kz
F4 = ρ ⋅ g ⋅ a w ∫ ∫ e dy + B ′(ξ ) ⋅ [cos(kξ ) ⋅ cos(ωt ) + sin (kξ ) ⋅ sin (ωt )]dξ
− L / 2 − B′ / 2
(3.70)
L/2
B′ / 2 kz
M 4 = ρ ⋅ g ⋅ a w ∫ ξ ∫ e dz + B ′(ξ ) ⋅ [cos(kξ ) ⋅ cos(ωt ) + sin (kξ ) ⋅ sin (ωt )] dξ
− L / 2 − B′ / 2
(3.71)
B′
z = − sin (α ) = y ⋅ tan (α )
2
donde, B′ = 2 ( L 2 )2 − ξ 2
π ⋅ ρ ⋅ g ⋅ a w ⋅ L2 k ⋅L
F4 = 1 − cos(ωt )
2 6
π ⋅ ρ ⋅ g ⋅ a w ⋅ k ⋅ L4
M4 = sin (ωt )
32
Denise Da Costa
– 56 –
Conceptos dinámicos importantes
Las ecuaciones (3.46), (3.50), (3.60) y (3.70) son combinadas con la ecuación (3.42)
para obtener la ecuación del movimiento heave
L/2
a= ∫m w dξ
−L / 2
L/2
dN π 2 ⋅ρ ⋅ g2 L/2
b= ∫ d ξ = ∫ A 2 dξ
−L / 2
dξ 8 ⋅ ωn −L / 2
3
L/2
c = ρ⋅g ∫ B ′ dξ
−L / 2
L/2
d= ∫ξ ⋅ m
−L / 2
w dξ
L/2
dmw π 2 ⋅ρ ⋅ g2 L/2
e = −V ∫ dξ ξ + 2 ⋅ m
w d ξ + ∫ξ ⋅ A
2
dξ
−L / 2 8 ⋅ ω n3 −L / 2
V ⋅π 2 ⋅ ρ ⋅ g 2
L/2 L/2
h = ρ ⋅ g ∫ ξ ⋅ B ′ dξ − ∫A
2
dξ
−L / 2 8 ⋅ ω n3 −L / 2
L/2
B′ / 2 kz
FA = ρ ⋅ g ⋅ a w ∫ ′ cos(kξ ) dξ
∫ e dy + B
− L / 2 − B′ / 2
L/2
B′ / 2 kz
FB = ρ ⋅ g ⋅ a w ∫ ′ sin (kξ ) dξ
∫ e dy + B
− L / 2 − B′ / 2
F0 = FA2 + FB2
− FB
γ = tan −1
FA
Las ecuaciones (3.47), (3.51), (3.61) y (3.71) son combinadas con la ecuación (3.43)
para obtener la ecuación del movimiento pitch
Denise Da Costa
– 57 –
Conceptos dinámicos importantes
L/2
A= ∫ξ ⋅ m w dξ
2
−L / 2
dmw
L/2
π 2 ⋅ρ ⋅ g2 L/2
B = −V ∫ ξ ξ + 2 ⋅ m w dξ + ∫ξ
2
⋅ A 2 dξ
−L / 2
dξ 8 ⋅ ω n3 −L / 2
V ⋅ π 2 ⋅ ρ ⋅g 2
L/2 L/2 L/2
dmw
C =V 2 ∫ ξ dξ + ρ ⋅ g 4 ∫ ξ 2 ⋅ B ′ dξ − ∫ξ ⋅ A
2
dξ
−L / 2
dξ −L / 2 8 ⋅ ω n3 −L / 2
D=d
L/2
dm w π 2 ⋅ρ ⋅ g2 L/2
E = −V ∫ ξ dξ + ∫ξ ⋅ A
2
dξ
−L / 2
dξ 8 ⋅ ω n3 −L / 2
L/2
H = ρ⋅g ∫ ξ ⋅ B ′ dξ
−L / 2
B′ / 2 kz
L/2
M A = ρ ⋅ g ⋅ a w ∫ ξ ∫ e dy + B ′ cos(kξ ) dξ
− L / 2 − B′ / 2
Movimiento heave:
π ⋅ ρ ⋅ g ⋅ a w ⋅ L2 k⋅L π ⋅ ρ ⋅ L3 π 2 ⋅ ρ ⋅ g2 ⋅ A2 ⋅ L
1 − cos(ωt ) = m + &z& + z& +
2 6 12 8 ⋅ ω n3
π ⋅ ρ ⋅ g ⋅ L2 V ⋅ π ⋅ ρ ⋅ L3 & V ⋅ π 2 ⋅ ρ ⋅ g 2 ⋅ A 2 ⋅ L
+ z− θ− θ
4 12 8 ⋅ ω n3
Movimiento pitch:
π ⋅ ρ ⋅ g ⋅ a w ⋅ k ⋅ L4 π ⋅ ρ ⋅ L5 && π 2 ⋅ ρ ⋅ g 2 ⋅ A 2 ⋅ L3 &
1 − k ⋅ L sin (ωt ) = I y +
2
θ + θ+
32 15 240 96 ⋅ ω n3
π ⋅ ρ ⋅ g ⋅ L4 V 2 ⋅ π ⋅ ρ ⋅ L3 V ⋅ π ⋅ ρ ⋅ L3
+ − θ + z&
64 12 12
– 58 –
Conceptos dinámicos importantes
CASO 1
Para cuerpos con simetría en el plano y – z se obtienen los siguientes resultados
CASO 2
L/2 L/2
∫ ξ ⋅ A dξ = A
2 2
∫ ξ dξ = 0
−L / 2 −L / 2
CASO 3
Cuando el cuerpo es simétrico y no existe corriente o no hay movimiento adelante del
cuerpo, entonces e=E=0
CASO 4
Cuando el cuerpo está completamente sumergido, c = H = 0 y también cualquier
término que tenga B ′ en la integral es cero excepto en los límites de las integrales para
FA ; FB ; M A ; M B
CASO 5
Para cuerpos simétricos con A constante y V = 0 los movimientos heave y pitch son
independientes desde que d = e = h = D = E = H = 0 . En este caso las ecuaciones de
movimiento son desacopladas.
Las cargas en los elementos de amarre pueden ser evaluadas mediante métodos de
análisis estático o dinámico. El análisis estático es apropiado cuando el sistema (buque y
sus elementos de amarre) no experimentará movimientos dinámicos significativos para las
acciones consideradas en el diseño. Por otro lado, el análisis dinámico es necesario cuando
se prevé que el sistema estará sujeto a grandes movimientos. Desafortunadamente, es
difícil determinar, a priori, si el sistema experimentará cargas dinámicas importantes; por lo
Denise Da Costa
– 59 –
Conceptos dinámicos importantes
que esta determinación debe ser hecha para cada problema en particular, considerando las
condiciones ambientales del lugar, el tipo de buque y el arreglo de amarre (Tsinker, 1995).
En la tabla 3.4 se resume los casos más comunes de problemas dinámicos y estáticos
de modo a asistir al ingeniero proyectista en identificar cuando será necesario un análisis
dinámico; recordando que esta tabla es simplemente una guía.
Sin embargo, el presente capítulo es de importancia en el sentido que enriquece la tesis con
miras a posterior utilización en desarrollos futuros.
Denise Da Costa
– 60 –
Sistemas de defensas
SISTEMAS DE DEFENSAS
Sin embargo, las defensas no son solamente para proteger al sistema del impacto,
también sirven como protección tanto para el buque como para el muelle mientras el buque
está amarrado a este.
4.1 TIPOS
Defensas de madera.
Sistemas de pilares.
Defensas de goma.
Sistemas hidráulicos e hidroneumáticos.
Resortes.
Las defensas de madera son el sistema de defensas más simple que existe; consisten
en troncos de madera con fibra de alta dureza y resistencia al impacto. La energía de
absorción se traduce en el aplaste de las fibras y la flexión del grupo de troncos. La pérdida
de fuerza puede ser causada por una mala construcción provocando que la madera se
astille o por consideraciones de altas energías (Tsinker, 1986).
Este sistema es frecuentemente combinado con otro tipo de sistema para proveer
mayor área de contacto incrementando así la eficiencia del sistema. Deben ser utilizadas
superficies de goma para proteger al buque de posibles daños ocasionados por el muelle
durante el proceso de atraque o de posibles daños ocasionados por el movimiento relativo
con otros barcos mientras se encuentran amarrados (Tsinker, 1986).
Denise Da Costa
– 61 –
Sistemas de defensas
Este sistema emplea pilares de madera dirigidos a lo largo del muelle; el tope de los
pilares puede estar sin apoyo lateral o apoyado en varios niveles de fijación por medio de
apoyos y de cuñas. El uso de uno o varios apoyos depende de la longitud del pilar y de la
variación de la marea (Derucher & Heins, 1979).
– Figura 4.1 –
Sistema de pilares de madera.
Denise Da Costa
– 62 –
Sistemas de defensas
– Figura 4.2 –
Sistema de defensa de Madera Colgante.
Los pilares de acero son ocasionalmente utilizados cuando la profundidad del agua es
superior a 12 metros o cuando las solicitaciones debido a las fuerzas son elevadas. Estos
sistemas tienen como ventaja una mayor resistencia ante grandes fuerzas y cuando las
condiciones del suelo son adversas, este tipo de pilares suele presentarse como solución.
Las desventajas son, vulnerabilidad a la corrosión y alto costo.
Los pilares de concreto reforzado no son una practica común debido a su limitada
capacidad para resistir las fuerzas de tracción; sin embargo, se han utilizado pilares de
concreto pretensado con refuerzos intermedios de goma en el nivel superior. La ventaja es
que estos pilares resisten la deterioración natural y biológica. Las desventajas son su baja
resistencia a tracción y la corrosión limitada del acero de refuerzo a través de las grietas.
Denise Da Costa
– 63 –
Sistemas de defensas
– Figura 4.3 –
Defensa Cilíndrica (Fentek, catalogo)
Denise Da Costa
– 64 –
Sistemas de defensas
Los tubos de goma cubiertos que cuelgan de los tabiques herméticos sólidos del buque
se pueden utilizar como sistema de defensa en compresión. La capacidad de absorción de
energía de tal sistema puede ser variada usando los tubos en capas solas o dobles, o
variando el tamaño del tubo. La absorción de la energía de un tubo cilíndrico es casi
directamente proporcional a la fuerza del buque hasta que la deflexión iguala
aproximadamente mitad del diámetro externo; después de esto, la fuerza aumenta mucho
más rápidamente que la absorción de energía. Por lo tanto, debe ser utilizada una defensa
suficientemente grande de modo que la energía de atraque sea absorbida sin requerir una
desviación de tal magnitud que dé lugar a un aumento desproporcionado de la fuerza.
Consisten en una serie de bloques de goma entre láminas de acero formando una
especie de sándwich que será montado en un soporte de acero que se une a las placas
mediante agujeros y otra placa de acero en el tope que sirve de amortiguador de la energía
de impacto. Las defensas en corte tienen una capacidad de absorción de energía de 100%,
el único problema es que son muy rígidas para pequeñas embarcaciones y las placas de
acero tienen a corroerse; sin embargo, dada su elevada capacidad de absorción de energía
son ideales para grandes buques (Derucher & Heins, 1979), (Tsinker, 1986).
– Figura 4.4 –
Defensa en corte tipo Raykin
Denise Da Costa
– 65 –
Sistemas de defensas
Estos tipos de defensas, especialmente las denominadas Raykin, poseen una rigidez
muy elevada lo que trae como desventaja no poder ser utilizadas para embarcaciones
pequeñas; otras desventajas son la posible corrosión de las placas de acero y problemas
con el enlace entre la placa y la goma.
4.1.3.3 Flexibles
Este sistema consiste en un bloque de goma arqueado enlazado entre dos placas de
acero. Pueden ser instalados en muelles abiertos (más flexibles) o muelles reforzados (más
rígidos), duques de atraque, o incorporado en un sistema de pilares. La energía de impacto
es absorbida doblando la columna de compresión arqueada; cuando se aplica una fuerza de
impacto, se acumula una carga relativamente alta con una desviación pequeña; las hebillas
todavía con desviaciones más pequeñas, mantienen una carga virtualmente constante
sobre el rango de las deformaciones por el doblez
Las ventajas de este tipo de defensa son, su alta absorción de energía y que las
fuerzas transmitidas al muelle son relativamente pequeñas. Estas defensas tienen como
desventajas, posible destrucción de la ligación de contacto entre las placas de acero y la
goma y problemas de fatiga del material.
Estas defensas son diseñadas para absorber la energía de impacto, por compresión del
aire contenido en ellas; y pueden ser utilizadas tanto en estructuras de atraque como
defensas incorporadas en el buque. Consisten en un grupo de neumáticos de caucho y
ruedas capaces de rodar libremente alrededor de un eje fijo o flotante; donde la parte fija es
diseñada para la colocación de tabiques herméticos en concreto, la parte flotante consiste
en dos a cinco neumáticos. La capacidad de absorción de energía y la resistencia a las
fuerzas de impacto dependen del tamaño y número de neumáticos y de la presión de aire
inicial en ellos instalada cuando son inflados.
Las defensas de neumáticos pueden ser aseguradas al muelle o duque de alba por
medio de cabos o cadenas que son atadas a los anillos que se encuentran en los extremos,
otra forma de asegurarlas al muelle es introduciendo la defensa en una red y colgarla por
varios puntos a lo largo de su longitud. Cuando estas defensas son usadas en duques de
Denise Da Costa
– 66 –
Sistemas de defensas
Las ventajas de este sistema es que es conveniente para ambos casos, atraque y
amarre de buques; además de ser utilizados para la protección de las esquinas de los
muelles. En contra partida presentan un alto costo de mantenimiento.
– Figura 4.5 –
Defensa de Neumático (Fentek, catalogo)
Este sistema consiste en un cilindro lleno de aceite u otro fluido, dispuesto de tal forma
que cuando es presionado el émbolo por el impacto, el líquido es desplazado a través de un
orificio, no variable o variable, en un depósito situado a una cota superior. Cuando es
disipada la energía de impacto, la alta presión dentro del cilindro crea una fuerza que coloca
al émbolo de nuevo en su posición original y el líquido fluye nuevamente dentro del cilindro
por gravedad. Este sistema se utiliza comúnmente donde existen viento severo, olas y
corriente.
En este sistema, una defensa de goma es llenada de agua o agua y aire; absorbiendo
así la energía de impacto, por resistencia viscosa o por la compresión del aire. Esta defensa
cumple los requisitos de la defensa ideal pero es económicamente poco viable por su alto
Denise Da Costa
– 67 –
Sistemas de defensas
costo inicial y de mantenimiento. Al igual que los amortiguadores hidráulicos, este tipo de
defensa presenta características de absorción de energía favorables para los casos de
atraque y amarre de buques.
4.1.5 Resortes
Existe una extensa variedad de factores que afectan la selección apropiada del sistema
de defensa. Estos factores incluyen: el ambiente local de la marina, la exposición del puerto,
clase y configuración del buque, velocidad y dirección del buque en el momento de atraque,
ayuda disponible en el muelle, tipo de estructura de atraque, e incluso las habilidades del
piloto o capitán del buque.
– 68 –
Sistemas de defensas
energía, tales como las defensas tipo Raykin o los sistemas de defensas de goma en
compresión axial. Este es el caso cuando el atraque ocurre bajo la acción de vientos, de
corrientes y olas, sin la ayuda de remolcadores.
Para localizaciones donde el comportamiento del buque amarrado es el principal factor, son
utilizadas defensas suaves combinadas con cabos de amarre obteniendo así una reducción
en las fuerzas de amarre y minimización de los movimientos del buque. Un tipo de defensas
suaves (flexible), como por ejemplo las defensas de goma en compresión radial, tiende a
aumentar el período natural de oscilación del sistema de amarre, por eso la resonancia
asociada a las olas de largo período o “seiches” puede ser evitado. Esto es aplicable en
puertos en los cuales el atraque es auxiliado por remolcadores; pero la oscilación del agua
en la zona protegida del puerto debido a la acción de los “seiches” es un factor importante
en la selección de la defensa.
La distancia requerida por el sistema de defensa debe ser limitada para evitar
inconvenientes durante el proceso de carga y descarga; generalmente el límite máximo es
de 1 a 1,5 metros. Sin existir ningún problema si el sistema de defensa es para un tanque
petrolero que implique solamente una fuente de combustible.
Para un muelle abierto, cualquier tipo de sistemas de defensa puede ser aplicable. Para
un muelle sólido, se debe considerar el uso de defensas resistentes de modo a reducir al
mínimo daños en el buque.
Denise Da Costa
– 69 –
Sistemas de defensas
Denise Da Costa
– 70 –
Sistemas de defensas
Coherentemente con los requisitos de fuerzas, el espaciamiento entre pilares debe ser
de 3 a 3,5 metros en casos de servicio ligero; para cruceros y auxiliares el espacio debe
estar entre 2 y 2,7 metros y para servicios pesados es recomendable espacios de 1,5 a 2
metros entre pilares.
Las esquinas externas y las expuestas deben ser protegidas con grupos de pilares de
defensa. Para pequeñas embarcaciones, incluyendo destroyers hasta un máximo de 3.000
toneladas, estas esquinas son protegidas con grupos de 7 a 9 pilares dispuestos en dos
filas. Para muelles donde atracan buques mayores, las esquinas expuestas se deben
recubrir con pilares de defensa dispuestos en dos filas escalonadas en arco circular con
radio de 4 metros; excepto en los casos donde los buques son mucho mayores y/o la
estructura este localizada de forma tal que este sumamente expuesta, son necesarias tres o
cuatro filas de pilares.
Las defensas tipo resorte o de goma pueden ser reforzadas con pilares de madera
dispuestos en dos filas atornillados a otros elementos horizontales de refuerzo (wale) que
soportan el sistema de absorción de energía. Si se utilizan defensas de goma tubulares, los
pilares de defensa son dispuestos en dos filas separadas y ligados por cuerdas de acero en
varias vueltas. El número de pilares varia entre 20 a 40, ligados por cadenas y cabos para
restringir movimientos longitudinales o laterales del conjunto.
– 71 –
Sistemas de defensas
realizado un dragado mayor, debe realizarse una penetración mayor. Si existe una capa de
material blando, con espesura inferior a 3 metros, por encima de un suelo firme, los pilares
de defensa deben penetrar el suelo firme por lo menos 2,5 metros de profundidad. Para
suelos firmes a gran profundidad, debajo de una capa espesa de material blando, los pilares
deben alcanzar una penetración por lo menos igual a la determinada por el peso del martillo
de penetración y debe asegurar una capacidad resistente de 2 a 3 toneladas (conforme a la
fórmula de estacas). La experiencia indica que la capacidad resistente aumenta después de
la conclusión de la penetración, de modo a asegurar la resistencia necesaria de los pilares.
La inclinación de los pilares no debe exceder de 1/6 hacia el lado del muelle.
Para defensas tipo resorte o defensas de goma, consistentemente con los requisitos de
fuerza, el espaciamiento vertical de estas defensas puede ser incrementado por encima de
los límites de espaciamiento previamente nombrados para cada tipo de buque.
Para defensas tipo Raykin y tipo amortiguador, el espaciamiento debe ser conforme a la
fuerza y requisitos de absorción de energía. Para defensas elásticas sin incluir dispositivos
de goma, son necesarias superficies de metal o de madera (ambas con gran resistencia por
fricción) para resistir las fuerzas transversales.
Denise Da Costa
– 72 –
Sistemas de amarre
SISTEMAS DE AMARRE
Una vez que el buque llega al muelle es amarrado a este mediante un sistema de
amarre para evitar que ocurran grandes desplazamientos relativos entre ambos. La elección
de un sistema de amarre depende del tipo de rompeolas que se utilizarán y de las
condiciones del sitio. Esta investigación está principalmente referida a los cabos usados en
el sistema y no a las anclas. Se discuten los tipos de estructuras de amarre y atraque, la
configuración de los cabos de amarre y las cadenas de anclaje.
5.1 DEFINICIÓN
– Figura 5.1 –
Sistema de Amarre
Denise Da Costa
– 73 –
Sistemas de amarre
Las estructuras de amarre están diseñadas para soportar las cargas de amarre;
mientras que las estructuras de atraque están diseñadas para resistir ambas solicitaciones,
las de amarre y las de impacto. Las estructuras más comunes en estos sistemas son:
plataformas y duques de atraque. Las estacas aisladas de amarre son frecuentemente
usadas como estructuras de amarre en muelles flotantes.
5.2.1 Plataformas
Las plataformas pueden presentar dos tipos de comportamiento, rígido o flexible. Las
plataformas rígidas consisten en una placa de concreto soportada por un conjunto de pilares
de concreto verticales e inclinados; esta estructura está diseñada para resistir grandes
cargas laterales que se transfieren, en mayor parte, a los pilares inclinados. Cuando un
buque atraca en estas plataformas, la carga de impacto es distribuida en los pilares
inclinados de forma tal que los pilares que se encuentran en la parte posterior resisten
cargas de compresión y los pilares frontales resisten cargas de tracción, provocando un
efecto de levantamiento de estos (Ver figura 5.2). Debido a que la capacidad de los pilares
para resistir el efecto de levantamiento es generalmente pequeña comparada con la
capacidad del suelo para resistir las compresiones, es recomendable que la placa de
concreto tenga una espesura suficiente de manera a resistir parte de la carga que tiende a
levantar a los pilares; en consecuencia los pilares verticales son diseñados para resistir el
peso de la placa. El sistema de defensas de estas plataformas consiste en pilares de
madera dispuestos en fila en la parte frontal y en grupo en las esquinas; cabe recordar que
en las plataformas rígidas toda la energía de impacto es absorbida por el sistema de
defensas. En suelos poco resistentes es usual el uso de plataformas rígidas en vez de
duques de atraque.
Denise Da Costa
– 74 –
Sistemas de amarre
– Figura 5.2 –
Distribución de la fuerza de impacto en los pilares inclinados de una plataforma rígida.
Las plataformas flexibles están formadas por pilares tubulares de acero y por una placa
que puede ser de acero o de concreto. Dado el carácter flexible de estas estructuras, los
pilares de acero tienden a deformarse, absorbiendo así parte de la energía cinética; de está
forma las solicitaciones en las defensas son menores que en las plataformas rígidas, puesto
que ahora las fuerzas de impacto son resistidas por el sistema de defensas y por la
estructura (NFEC, 1986).
Los duques de atraque o duques de alba son estructuras aisladas, de amarre y/o de
atraque, que se encuentran fuera del muelle. Como su nombre lo indica estas estructuras
son comúnmente más usadas como estructuras de defensa, ya que son capaces de resistir
grandes cargas de atraque, producidas al momento del impacto del buque con estas.
Los duques de atraque pueden ser flexibles o rígidos. Los duques de alba rígidos
consisten en un único pilar tubular de acero o de concreto, relleno con material granular;
que sirve de soporte a la placa de concreto provista con defensas de goma. Estos tipos de
duques por ser rígidos requieren de un sistema de defensas capaz de absorber la energía
de impacto. Y son diseñados para el atraque de grandes buques.
Los duques de atraque flexibles son formados por una serie de pilares de madera
constructiva (tres, siete, diecinueve o más) y asegurados en el extremo superior con
cuerdas. Uno de los pilares centrales posee una longitud mayor de manera que sirva de
punto de amarre del cabo. Debido al carácter flexible de estas estructuras, no es necesario
Denise Da Costa
– 75 –
Sistemas de amarre
un sistema de defensas, sin embargo, la perdida de uno de los pilares, por un impacto
brusco del buque, causar una redistribución de las fuerzas en los demás pilares y la
consecuente debilitación del sistema. Por tanto, no es recomendable el uso de este tipo de
duques como estructura de atraque. Surge entonces otra modalidad de los duques flexibles,
que consiste en una serie de pilares de acero, con un sistema de defensas de madera o de
goma, logrando así una mayor resistencia a las cargas laterales. Estos duques son
típicamente usados para el atraque de buques de menores proporciones.
5.3.1 Tipos
Los buques son asegurados al muelle mediante cabos de amarre; estos cabos son
clasificados según el tipo de carga que resisten (longitudinal o lateral) y según su posición
respecto al buque. Los cabos de proa, de popa y laterales resisten cargas laterales y
momento. Los cabos “resortes longitudinales” (spring) resisten cargas longitudinales. Las
cargas predominantes en un sistema de amarre son las laterales, en consecuencia existe un
mayor número de cabos de proa, de popa y laterales que cabos “resorte longitudinal”
(spring).
– Figura 5.3 –
Arreglo típico de los cabos en un buque.
5.3.2 Materiales
Existen cuatro tipos de cabos para asegurar un boque a un muelle; estos son de: Fibra
natural, fibra sintética, alambre y cadena. Los más usados son los cabos de fibra natural y
fibra sintética. Las cadenas son usadas para amarrar buques inactivos.
Denise Da Costa
– 76 –
Sistemas de amarre
5.3.3 Características
Las características más importantes para los cabos son su resistencia, elasticidad y
construcción.
Resistencia: Para un mismo diámetro con distintos materiales, se tiene que los
cabos de alambre son los más resistentes seguidos de los cabos de fibra sintética
y luego los de fibra natural.
Elasticidad: La elasticidad de los cabos de determina mediante un gráfico carga vs
porcentaje de elongación. De los distintos tipos de cabos se tiene que los de
alambre son considerablemente menos elásticos que los demás.
Construcción: Los cabos sintéticos varían de acuerdo a su construcción
En general, los cabos de ocho cordones y doble trenza son preferibles (Flory et al.,
1977). Los cabos de tres cordones no son recomendables, pues tienen tendencia
a formar lazos, lo cual reduce gradualmente la fuerza del cabo, resultado de torcer
un cabo puesto en la dirección opuesta a la dirección de configuración. Este
efecto, es semejante a tener nodos en los cabos.
– Figura 5.4 –
Representación de la pérdida de rigidez de los cabos con tres cordones por formación de
lazos.
Denise Da Costa
– 77 –
Sistemas de amarre
5.4.1 Comportamiento
5.4.2 Materiales
Las cadenas son utilizadas en la mayoría de los casos de anclaje en comparación con
otros tipos de cabos como fibra sintética, fibra natural, etc; porque las cadenas tienen una
alta resistencia a la abrasión y presentan características favorables de absorción de
energía.
5.4.2.1 Cadenas
En el cálculo de las catenarias se puede observar que las cadenas al tener un peso por
unidad de longitud mayor que el de las cuerdas de alambre, poseen mayor capacidad de
absorción de energía. Mientras que la cadena comportándose como catenaria tiene esta
característica muy deseable de brindar buena capacidad de absorción de energía, presenta
también la desventaja de ser tan fuerte como su eslabón lo permita.
Existen tres tipos principales de cadenas de anclaje; fundidas, soldadas (flash butt-
welded), y dilok. Estos tipos difieren unos de otros en el método de manufactura y en su
resistencia. Tanto las cadenas fundidas como las soldadas (flash butt-welded) son utilizadas
en sistemas de amarre; mientras que las cadenas dilok son utilizadas principalmente como
cadenas del buque.
Denise Da Costa
– 78 –
Sistemas de amarre
Dilok: Este es un tipo de cadena forjada que no requiere soldadura o adhesión del
metal durante la fabricación. Cada eslabón consiste en dos partes, un macho y
una hembra, la fabricación consiste en perforar y calentar la parte hembra del
eslabón; la parte macho es enlazada a través de la próxima conexión y es inserida
fría en la hembra, que es martillada sobre la parte macho. De este proceso resulta
una ligación de resistencia prácticamente uniforme siendo generalmente más
resistente que las ligaciones por fundición o soldadura.
Es evidente que las cadenas tipo dilok son más susceptibles al fracaso que las cadenas
fundidas o soldadas. La naturaleza de la construcción de este tipo de cadenas permite la
posibilidad que el agua penetre a través de la junta de cierre, causando problemas de
corrosión que no son detectados con una simple inspección visual.
Las cuerdas de alambre son una alternativa. Estos cabos de anclaje son más ligeros
que las cadenas con la misma resistencia a tensión y en consecuencia también poseen una
menor capacidad de absorción de energía, durante el comportamiento como catenarias. Sin
embargo, las cuerdas de alambre son mucho más elásticas que las cadenas teniendo así
una respuesta tan buena como la de una cadena de la misma fuerza.
Denise Da Costa
– 79 –
Sistemas de amarre
– Figura 5.5 –
Cabo de anclaje compuesto por cadena y cuerda de alambre.
Denise Da Costa
– 80 –
Sistemas de amarre
Los buques llevan normalmente sus propios cabos de amarre; por lo tanto, el diseño de
un muelle raramente envuelve estos elementos. Sin embargo, el tamaño, la longitud, la
composición, y el arreglo de los cabos afectan el diseño del muelle. Por lo tanto, es
necesario determinar el tipo y el arreglo de los cabos del buque(s) para el cual se diseña el
muelle.
La configuración de los cabos depende sobre todo del uso del muelle, espacio
disponible para los cabos, tamaño de la embarcación, cargas de viento y oleaje, fuerzas de
amarre y de impacto; tipo, número y tamaño de los cabos disponibles en el buque (Tsinker,
1994).
Para estos cabos, el ángulo vertical con la horizontal debe ser menor a 25º; este ángulo
depende de la altura de la estructura de amarre y de la longitud de los cabos. Los cabos de
popa, de proa y laterales resisten mejor las cargas laterales cuando el ángulo vertical es
más pequeño (Ver figura 5.6).
Para los cabos de popa y de proa, el ángulo horizontal, con la línea central transversal
del buque debe ser menor que 45º. Para los cabos laterales, este ángulo debe ser menos
que 15º (Argeschou, 1984; OCIMF, 1997) Como se indica en la figura 5.6.
Para estos cabos, el ángulo vertical con la horizontal debe ser menor que 25º. A
medida que este ángulo sea menor, los cabos “resorte longitudinal” son más eficientes
porque resisten menor las cargas longitudinales (Argeschou, 1984; OCIMF, 1997) (Ver
figura 5.6).
El ángulo horizontal, con el eje longitudinal del buque debe ser menor que 10º (Ver
figura 5.6)
Denise Da Costa
– 81 –
Sistemas de amarre
– Figura 5.6 –
Ángulos horizontales y verticales para los cabos de amarre (OCIMF, 1997).
Denise Da Costa
– 82 –
Metodología de análisis
METODOLOGÍA DE ANÁLISIS
Inicialmente, para obtener las fórmulas necesarias para el cálculo, se aceptará que los
elementos de sujeción (cabos y defensas) pueden resistir tanto tracción como compresión;
más adelante, se explicará con detalle como esta falsa aceptación es remediada por un
procedimiento en el análisis.
Para hallar la matriz de rigidez del miembro, que representa las relaciones fuerza-
desplazamiento entre los extremos del mismo, se establece un sistema de coordenadas
locales x’ - y’ ubicando su origen en el extremo del miembro y alineando el eje x’ con el eje
del miembro, el estado general de cargas y deformaciones para una barra se muestra en la
figura 6.1
Denise Da Costa
– 83 –
Metodología de análisis
– Figura 6.1 –
Estado general de cargas y deformaciones para una barra.
La matriz de rigidez establece las relaciones que existen entre las cargas y los
desplazamientos, para ello, el caso general de carga será descompuestos en dos casos,
cada uno asociado al desplazamiento de un extremo (Hibbeler, 1997).
Caso 1
– Figura 6.2 –
Desplazamiento permitido en el nodo i e impedido en el nodo j
Caso 2
– Figura 6.3 –
Desplazamiento permitido en el nodo j e impedido en el nodo i
Las ecuaciones de equilibrio obtenidas, para cada caso, son las siguientes
A⋅ E
f i´ = di (6.1a)
L
Denise Da Costa
– 84 –
Metodología de análisis
− A⋅ E
f j´ = di (6.1b)
L
− A⋅ E
f i´´ = dj (6.2a)
L
A⋅ E
f j´´ = dj (6.2b)
L
A⋅ E A⋅ E
f i = f i´ + f i´´ = di − dj (6.3a)
L L
A⋅ E A⋅ E
f j = f j´ + f j´´ = − di + dj (6.3b)
L L
fi EA 1 − 1 d i
= ×
fj L −1 1 d j
f = k′ ⋅ d (6.4)
EA 1 −1
k′ = =
L −1 1
q B EA 1 − 1 d B
= ×
qE L −1 1 dE
q = k′⋅ d (6.4a)
Denise Da Costa
– 85 –
Metodología de análisis
– Figura 6.4 –
Posición relativa del sistema de coordenadas locales respecto al sistema global.
XE − XB
λx = Cosθx = (6.5a)
( X E − X B ) 2 + (YE − YB ) 2
YE − YB
λy = Cosθy = (6.5b)
( X E − X B ) 2 + (YE − YB ) 2
donde:
XE = X escoben
YE = Y escoben
XB = X bita
YB = Y bita
Denise Da Costa
– 86 –
Metodología de análisis
– Figura 6.5 –
Definición de los grados de libertad del sistema
Denise Da Costa
– 87 –
Metodología de análisis
– Figura 6.6 –
Representación del modelo en análisis con sus grados de libertad, coordenadas globales y
locales
Considerando por separado cada uno de los desplazamientos globales se analizan los
desplazamientos generados según las coordenadas locales.
Denise Da Costa
– 88 –
Metodología de análisis
dE = D1 * Cos θx
dE = D1 * λX
– Figura 6.7 –
Sistema con todos los grados de libertad impedidos, excepto el desplazamiento en 1
Denise Da Costa
– 89 –
Metodología de análisis
dE = D2 * Cos θy
dE = D2 * λY
– Figura 6.8 –
Sistema con todos los grados de libertad impedidos, excepto el desplazamiento en 2
Denise Da Costa
– 90 –
Metodología de análisis
R senD3 ∼ R*D3
dE = R*D3 * sen ( θx - αx )
dE = R*D3 ( sen θx ⋅ cos αx - cos θx ⋅ sen αx )
R cos αx = x´E
R sen αx = y´E
– Figura 6.9 –
Sistema con todos los grados de libertad impedidos, excepto la rotación en 3
Denise Da Costa
– 91 –
Metodología de análisis
dB = D4 * Cos θx
dB = D4 * λx
– Figura 6.10 –
Sistema con todos los grados de libertad impedidos, excepto el desplazamiento en 4
Denise Da Costa
– 92 –
Metodología de análisis
dB = D5 * Cos θy
dB = D5 * λY
– Figura 6.11 –
Sistema con todos los grados de libertad impedidos, excepto el desplazamiento en 5
d B = D4 ⋅ λ X + D5 ⋅ λY
(6.6a)
d E = D1 ⋅ λ X + D2 ⋅ λY + D3 ⋅ a
(6.6b)
Denise Da Costa
– 93 –
Metodología de análisis
D4
D5
d B λX λY 0 0 0
= × D1
dE 0 0 λX λY a
D2
D3
d =T ⋅D (6.7)
donde:
λX λY 0 0 0
T=
0 0 λX λY a
– Figura 6.12 –
Representación de las fuerzas en coordenadas locales del sistema.
Denise Da Costa
– 94 –
Metodología de análisis
Análogamente, considerando por separado cada una de las fuerzas aplicadas en los
ejes locales y analizando sus componentes en coordenadas globales, se pueden obtener
las siguientes ecuaciones:
Q4 = q B ⋅ Cosθ X (6.8a)
Q5 = q B ⋅ Cosθ Y (6.8b)
Q1 = FX = q E ⋅ Cosθ X = q E ⋅ λ X (6.8c)
Q2 = FY = q E ⋅ Cosθ Y = q E ⋅ λY (6.8d)
Q3 = q E ⋅ (λY ⋅ x ′e − λ X ⋅ y ′e ) = q E ⋅ a (6.8e)
Q4 λ X 0
Q5 λY 0
q
Q1 = 0 λX × B
qE
Q2 0 λY
Q3 0 a
Q =TT ⋅q (6.9)
Las defensas son consideradas como miembros que solo impiden el movimiento
vertical en la coordenada X en que se encuentren, por lo tanto, sus cosenos directores
embarcación.
q = k′⋅ d (6.4a)
Expresión matricial que relaciona las fuerzas f aplicadas al miembro con los
desplazamientos d que ocurren como consecuencia estas; ambas en coordenadas locales.
d =T ⋅D (6.7)
Denise Da Costa
– 95 –
Metodología de análisis
Q =TT ⋅q (6.9)
Expresión matricial que relaciona las fuerzas Q en coordenadas globales con sus
respectivas componentes en coordenadas locales.
Con estas expresiones podemos obtener la relación entre las fuerzas y los
desplazamientos, ambos en coordenadas globales para cada elemento.
Q = T T ⋅ k′⋅ d (6.10a)
Q = T T ⋅ k′⋅T ⋅ D (6.10b)
λX 0
λY 0
EA 1 −1 λX λY 0 0 0
k= 0 λX × ×
L −1 1 0 0 λX λY a
0 λY
0 a
λ2X λ X ⋅ λY − λ2X − λ X ⋅ λY − a ⋅ λX
EA
K= λ X ⋅ λY λY2 − λ X ⋅ λY −λ 2
Y − a ⋅ λY
L
− λ2X − λ X ⋅ λY λ 2
X λ X ⋅ λY a ⋅ λX
Denise Da Costa
– 96 –
Metodología de análisis
Sustituyendo k
Q4 λ2X λ X ⋅ λY − λ2X − λ X ⋅ λY − a ⋅ λ X D4
Q5 λ X ⋅ λY λY2 − λ X ⋅ λY − λY2 − a ⋅ λY D5
EA
Q1 = − λ2X − λ X ⋅ λY λ2X λ X ⋅ λY a ⋅ λ X × D1
L
Q2 − λ X ⋅ λY − λY2 λ X ⋅ λY λY2 a ⋅ λY D2
Q3 − a ⋅ λX − a ⋅ λY a ⋅ λX a ⋅ λY a 2
D3
Q = k ⋅D (6.11)
QB EA K 11 K 12 DB
= ×
QE L K 21 K 22 DE
Una vez obtenida la matriz de rigidez global de todos los elementos es necesario
ensamblarlas en el orden apropiado para poder encontrar la matriz K de rigidez de la
estructura entera. El proceso de combinar las matrices de miembro depende de una
cuidadosa identificación de los elementos de cada matriz. La matriz de rigidez global de la
estructura tendrá un orden que será igual al número de grados de libertad total que esta
posea. Cuando se ensamblan las matrices k , cada uno de sus elementos se colocará en su
misma posición de fila y columna en la matriz K de rigidez de la estructura (Hibbeler, 1993,
1997). En particular, cuando dos o más miembros están conectados al mismo nodo,
entonces algunos de los elementos de cada una de las matrices k se asignaran a la misma
posición en la matriz K , tal como en el problema en estudio donde todos los miembros
convergen a un mismo cuerpo rígido que es la embarcación. Cuando esto ocurre, los
elementos asignados a la posición común deben sumarse entre si algebraicamente.
Denise Da Costa
– 97 –
Metodología de análisis
Q4 x x 0 0 0 0 0 0 0 0 x x x D4
Q5 x x 0 0 0 0 0 0 0 0 x x x D5
Q6 0 0 x x 0 0 0 0 0 0 x x x D6
Q7 0 0 x x 0 0 0 0 0 0 x x x D7
Q8 0 0 0 0 x x 0 0 0 0 x x x D8
Q9 0 0 0 0 x x 0 0 0 0 x x x ⋅
⋅ = 0 0 0 0 0 0 ⋅ ⋅ 0 0 x x x × ⋅
⋅ 0 0 0 0 0 0 ⋅ ⋅ 0 0 x x x ⋅
⋅ 0 0 0 0 0 0 0 0 x x x x x ⋅
Qn 0 0 0 0 0 0 0 0 x x x x x Dn
Q1 x x x x x x x x x x xxx xxx xxx D1
Q2 x x x x x x x x x x xxx xxx xxx D2
Q3 x x x x x x x x x x xxx xxx xxx D3
Q = K ⋅D (6.12)
Como hemos identificado los grados de libertad del buque seguidamente, podremos
subdividir la expresión anterior Q = K ⋅ D de la siguiente manera
QB K K 12 DB
= 11 ×
QE K 21 K 22 DE
donde:
QB = K 11 DB + K 12 DE (6.13)
QE = K 21 DB + K 22 DE (6.14)
Denise Da Costa
– 98 –
Metodología de análisis
Dado que el muelle es una estructura fija, y los desplazamientos en las bitas son cero
( DB = 0 ), las ecuaciones (6.13) y (6.14) quedan de la siguiente manera
QB = K 12 DE (6.13a)
QE = K 22 DE (6.14a)
DE = K 22−1 QE (6.14b)
Una vez conocidos los desplazamientos del buque, es fácil obtener las fuerzas
generadas en el extremo de la bita a través de la ecuación (6.13a)
Para determinar las fuerzas en los elementos de amarre, se debe sustituir la ecuación
(6.7) en la ecuación (6.4a) para cada miembro, donde D4 y D5 serán las componentes en
X e Y en el extremo cercano a la bita; de esta manera se obtiene la expresión matricial
D4
D5
q B EA 1 − 1 λ X λY 0 0 0
= × × D1
qE L −1 1 0 0 λX λY a
D2
D3
q = k′ ⋅T ⋅ D (6.15)
De los resultados obtenidos con esta expresión, debido a la convención de eje local que
hemos adoptado (Bita – Escoben), un valor positivo de q E indica que el elemento está a
tensión.
La expresión general (6.15) se aplica a todos los elementos por igual; donde D4 y D5
se refiere a los desplazamientos en los extremos cercanos a la bita para X e Y
respectivamente. Como las bitas se encuentran fijas al muelle y no sufren ningún
desplazamiento, la expresión (6.15) para el problema en estudio queda de la siguiente
manera
Denise Da Costa
– 99 –
Metodología de análisis
0
0
q B EA 1 − 1 λ X λY 0 0 0
= × × D1
qE L −1 1 0 0 λX λY a
D2
D3
6.1.7 Matriz K22
Como la ecuación (6.14a) lo indica, es la matriz que relaciona las cargas aplicadas a la
embarcación con los desplazamientos de la misma. Para hallar esta matriz es necesario
calcular la matriz K22 de cada miembro y luego sumarlas: K22 = Σ k22
λ2X λ X ⋅ λY a ⋅ λX
EA
K 22 =∑ λ X ⋅ λY λY2 a ⋅ λY
L
a ⋅ λX a ⋅ λY a2
La matriz K12 esta compuesta por todas las matrices k12 de cada miembro escritas en
forma de columna, es decir
K 12 (miembro"1")
K 12 (miembro"2")
K 12 (miembro"3")
K 12 = .
.
.
K 12 (miembro" n")
EA − λ2X − λ X ⋅ λY − a ⋅ λX
k12 =
L − λ X ⋅ λY − λY2 − a ⋅ λY
Denise Da Costa
– 100 –
Metodología de análisis
Es importante señalar que tanto la matriz K22 como la K12, son función únicamente de
la geometría inicial del problema y de las características de los elementos de amarre: área
(A) y módulo de elasticidad (E).
6.2 CATENARIAS
Una catenaria es la curva formada por un cable flexible de densidad uniforme que
cuelga a partir de dos puntos sobre su peso propio. Ejemplos de estas curvas podemos
observar en los cables de los puentes de suspensión y en los cables unidos a los postes del
teléfono (http://math.fullerton.edu/)
Los cabos de anclaje asumen la forma geométrica de una catenaria entre el fondo
marino y el escoben del buque (punto de amarre).
La figura 6.13 representa la definición de las partes de una catenaria usadas para el
análisis de estas curvas. Para cualquier punto (x, z) se cumple lo siguiente:
V = T ⋅ Senθ (6.16a)
H = T ⋅ Cosθ (6.16b)
T = w⋅ z (6.16c)
donde,
V = Fuerza vertical en el punto (x, z)
H = Fuerza horizontal en el punto (x, z)
T = Tracción en el punto (x, z)
θ = Angulo del cabo de anclaje con la horizontal.
w = Peso sumergido por unidad de longitud.
Denise Da Costa
– 101 –
Metodología de análisis
– Figura 6.13 –
Definición de las partes de una catenaria, para el análisis (U.S. Navy, 1985).
– Figura 6.14 –
Porción de un cabo de anclaje (catenaria)
Denise Da Costa
– 102 –
Metodología de análisis
dH
dx = 0 (6.17)
dx
de lo que se concluye que la fuerza horizontal (H) es constante a lo largo de todo el eje X.
dV
dV = dx = γ ⋅ ds (6.18)
dx
dz
V =H (6.19)
dx
2
dV d 2z ds dz
= H 2 =γ = γ 1+ (6.19a)
dx dx dx dx
dz
Designando t= y sustituyendo en la ecuación (6.19a), queda
dx
dt
H = γ 1+ t2 (6.19b)
dx
dt γ γ ⋅x
= dx → arc sh t = + C1
1+ t 2 H H
dz γ ⋅x
t= = sh + C1 (6.19c)
dx H
Denise Da Costa
– 103 –
Metodología de análisis
H γ ⋅x
z ( x) = ∫ t dx = ⋅ch + C1 + C 2 (6.20)
γ H
donde las constantes C1 y C 2 son determinadas con las condiciones de borde en el punto
A del cabo de anclaje. Dependiendo de estas condiciones de borde, existen dos tipos de
catenarias; catenaria completamente desarrollada y catenaria no completamente
desarrollada. A continuación se explica con detalles cada una de estas catenarias (Barros,
1990).
Esta curva catenaria es tangente con el eje horizontal, es decir, el ángulo del cabo de
anclaje con la horizontal es cero (θ = 0º). Por ende la componente vertical, de la fuerza de
tracción, en el origen es nula (V0 = 0). Sustituyendo estas condiciones de borde en la
ecuación general de la catenaria (ec. 6.20) se tiene que:
H
C2 = − (6.20a)
γ
H γ ⋅ x
z ( x) = ⋅ ch −1
γ H
(6.21)
La figura 6.15 representa una curva catenaria entre los puntos A(0, 0) situado en el
fondo marino y el punto B(x, z). Donde, la línea roja representa la catenaria límite y la línea
azul representa una catenaria con determinado peso y tracción.
Denise Da Costa
– 104 –
Metodología de análisis
– Figura 6.15 –
Catenaria entre los puntos A y B.
γ ⋅x V
tg α ( x ) =
dz
= sh = (6.21)
dx H H
La longitud de la catenaria entre el punto de origen A (0, 0) y el punto genérico B (x, z),
es dada por:
γ ⋅x γ ⋅x
x x x 2 x
dz H
S X = ∫ ds = ∫ dx + dz = ∫ 2 2
1 + dx = ∫ ch dx = ⋅ sh (6.22)
0 0 0 dx 0
H γ H
2
γ ⋅x γ ⋅x H γ ⋅x
x x x
H
∫ x ds ∫ x ch H
dx x ⋅ ⋅ sh
γ
− ⋅ ch
H γ H
dx
0
XG = 0
= 0
=
SX SX SX
γ ⋅x γ ⋅z
ch −1
H H H H = x− H ⋅ z
XG = x − ⋅ = x− ⋅
γ γ ⋅x γ γ ⋅ SX γ SX
sh
H H
Denise Da Costa
– 105 –
Metodología de análisis
H z
XG = x − ⋅ (6.23a)
γ SX
H γ ⋅ x γ ⋅x γ ⋅x γ ⋅x
x x x x
H H
∫ z ds ∫0 γ ⋅ ch H − 1 ⋅ ch H dx ∫γ ⋅ ch 2
H
dx − ∫
γ
⋅ ch
H
dx
ZG = 0
= = 0 0
SX SX SX
1 H 1 H γ ⋅x H 1 H 1 H H
⋅ ⋅ x + ⋅ ⋅ S X ⋅ ch − ⋅ SX ⋅ ⋅ x + ⋅ z + ⋅ S X − ⋅ S X
2 γ 2 γ H γ 2 γ 2 γ γ
ZG = =
SX SX
1 H 1 H H
⋅ ⋅ x + ⋅ z + ⋅ S X − ⋅ S X
2 γ 2 γ γ
ZG = (6.23b)
SX
donde,
γ ⋅x V γ ⋅ SX
α = arc tg sh = arc tg = arc tg
H H H
γ ⋅x 2 γ ⋅x γ ⋅z
2
T 2 = H 2 + V 2 = H 2 + (γ ⋅ S X ) = H 2 ⋅ 1 + sh 2 = H ⋅ ch = H 2 ⋅ + 1
2 2
H H H
por lo tanto,
γ ⋅ z
T = H ⋅ + 1 = H + γ ⋅ z (6.24b)
H
Denise Da Costa
– 106 –
Metodología de análisis
γ ⋅x γ ⋅x
SX SX SX
ds H H d H
E P = ∫ z ⋅ γ ds = γ ∫ z dx = γ ∫ γ ⋅ ch − ⋅ sh dx (6.25)
0 0
dx 0
H γ dx γ H
⋅ γ ⋅ z ⋅ S X − ⋅ H ⋅ (S X − x )
1 1
EP = (6.25a)
2 2
Esta curva forma un ángulo con el eje horizontal diferente de cero (θ ≠ 0º).
dz V0
Para x = 0; z = 0; = . Sustituyendo en la ecuación (6.19c)
dx H
V0 V
= sh C1 → C1 = arc sh 0 (6.26a)
H H
⋅ ch (C1 ) + C 2
H
0=
γ
2
H V H H V
C 2 = − ⋅ ch arc sh 0 = − ⋅ 1 + sh 2 C1 = − ⋅ 1 + 0
γ H γ γ H
2
H V
C2 = − ⋅ 1+ 0 (6.26b)
γ H
Denise Da Costa
– 107 –
Metodología de análisis
H γ ⋅ x V0
2
V0
z ( x ) = ⋅ ch
+ arc sh − 1 + (6.27)
γ H H H
– Figura 6.16 –
Definición de las variables de un cabo de anclaje, como catenaria (NFEC, 1985).
Denise Da Costa
– 108 –
Metodología de análisis
Las ecuaciones para determinar la distancia horizontal (R) en cada caso son
presentadas a continuación
6.2.4.1 Caso 1
2
S S2
R = L1 + L2 − S 2 + C 2 ⋅ log 2 + + 1 (6.28a)
C2 C2
donde,
S 2 = D ⋅ (D + 2 ⋅ C 2 )
6.2.4.2 Caso 2
a + a2 +1
R = L1 + C 2 ⋅ log 1 1 (6.28b)
v + v2 +1
donde,
2
1 D 4 D S2
v = ⋅ ⋅ − +1 −
2 C2 S2
2
C2 C 2
C2
L
a1 = v + 2
C2
6.2.4.3 Caso 3
S1
2 a + a2 +1
S1
R = L1 − S1 + C1 ⋅ log + + 1 + C 2 ⋅ log 2 2 (6.28c)
C1
C1 a1 + a1 + 1
2
Denise Da Costa
– 109 –
Metodología de análisis
donde,
S1 w
a1 = +
C1 H
L2
a 2 = a1 +
C2
(a )
2
S1
D = C1 ⋅ + 1 − 1 + C 2 ⋅ 2
2 + 1 − a12 + 1
C1
6.2.4.4 Caso 4
a + a2 +1 a + a2 +1
R = C1 ⋅ log 1 1 + C 2 ⋅ log
3 3 (6.28d)
v + v2 +1 a + a2 +1
2 2
donde,
L1
a1 = v +
C1
w
a 2 = a1 +
H
L
a3 = a 2 + 2
C2
D = C1 ⋅ (a 2
1 )
+ 1 − v 2 + 1 + C2 ⋅ (a 2
3 + 1 − a 22 + 1 )
Mediante el método de Newton-Raphson se determina el valor de v
donde,
a1 ; a 2 ; a 3 = Variables intermedias en la ecuación de la catenaria.
C1 ; C 2 = Constantes de la catenaria.
D = Profundidad del agua.
H = Fuerza horizontal de la cadena de anclaje.
L1 = Longitud de la sección inferior de la cadena.
L2 = Longitud de la sección superior de la cadena.
R = Distancia horizontal entre el ancla y el escoben.
S1 ; S 2 = Longitud de las secciones levantadas.
w = Peso sumergido del sinker.
Denise Da Costa
– 110 –
Metodología de análisis
El método de mínimos cuadrados asume que la mejor aproximación a una curva dada
es otra que tiene la suma mínima de las desviaciones ajustadas (menor error cuadrado) de
un sistema de datos.
d 1 = Y1 − f ( X 1 )
d 2 = Y2 − f ( X 2 )
....
d n = Yn − f ( X n )
n n 2
Dada la función f ( X ) = a o + a1 X
n 2 n 2
∂R n
= 2∑ [Yi − (a 0 + a1 X i )] = 0 (6.31a)
∂a 0 i =1
Denise Da Costa
– 111 –
Metodología de análisis
∂R n
= 2∑ [Yi − (a 0 + a1 X i )]X i = 0 (6.32a)
∂a1 i =1
a0 =
∑Y ⋅ ∑ X − ∑ X ⋅ ∑ X
i i
2
i i ⋅ Yi
n∑ X − (∑ X )
2
(6.31b)
2
i i
n∑ X i ⋅ Yi − ∑ X i ⋅ ∑ Yi
a1 =
n∑ X i2 − (∑ X i )
2
(6.32b)
Dada la función f ( X ) = a o + a1 X + a 2 X 2
[ ( )]
n 2 n 2
∂R
La condición para que el residuo sea mínimo es: =0
∂a i
∂R
[ ( )]
n
= −2∑ Yi − a0 + a1 X + a 2 X 2 = 0 (6.34)
∂a 0 i =1
∂R
[ ( )]
n
= −2∑ Yi − a0 + a1 X + a 2 X 2 X = 0 (6.35)
∂a1 i =1
∂R
[ ( )]
n
= −2∑ Yi − a0 + a1 X + a 2 X 2 X 2 = 0 (6.36)
∂a 2 i =1
n ∑X i ∑X i
2
a0 ∑Y i
∑ Xi ∑X i
2
∑X i
3
× a1 = ∑ X ⋅Y i i
∑ X i2 ∑X i
3
∑X i
4
a2 ∑ X ⋅Y i
2
i
Denise Da Costa
– 112 –
Metodología de análisis
f ( x ) = b ⋅ e a⋅ x (6.37)
f ( x) = b ⋅ x a (6.38)
( )
n 2
( )
n 2
R = ∑ Yi − b ⋅ X i
a
Para el caso de la ecuación (6.38)
i =1
f ( x ) = b ⋅ e a⋅ x (6.39)
Denise Da Costa
– 113 –
Metodología de análisis
n∑ X i ⋅ Ln(Yi ) − ∑ X i ⋅ ∑ Ln(Yi )
a=
n∑ X i2 − (∑ X i )
2
(6.40a)
Ln(b ) =
∑ Ln(Y ) ⋅ ∑ X − ∑ X ⋅ ∑ X
i i
2
i i ⋅ Ln(Yi )
n∑ X − (∑ X )
2
(6.40b)
2
i i
f ( x) = b ⋅ x a (6.41)
Ln(b ) =
i i i i i
n∑ Ln( X ) − [∑ Ln( X )]
2
(6.42b)
2
i i
El tipo de aproximación por trozos utilizando polinomios cúbicos entre cada par de
puntos sucesivos se llama interpolación “cubic spline” y es una de las técnicas más
populares actualmente usada. Debido a que en la forma general de un polinomio cúbico
están involucradas cuatro constantes, se tiene suficiente flexibilidad en el procedimiento de
interpolación por “cubic spline” para asegurar no sólo que la interpolación sea continua y
derivable en el intervalo sino que también posea una derivada segunda continua. Vale la
Denise Da Costa
– 114 –
Metodología de análisis
pena recordar, sin embargo, que no se asume que las derivadas de la interpolación
coincidan con las de la función (que se está aproximando) ni siquiera en los puntos
sucesivos (Burden, et al., 1978).
Dada la función f definida en [a, b] y una serie de puntos a = X0 < X1 < X2 < … < Xn = b,
la interpolación por “cubic spline” (S) para la función f satisface las siguientes condiciones:
S es un polinomio cúbico, denotado por Sj, en el intervalo [Xj ; Xj+1] para cada j = 0, 1, 2,
…, n-1
i) S ( x j ) = f ( x j ) para cada j = 0, 1, 2, …., n
b. S ′( x0 ) = f ′( x0 ) y S ′( x n ) = f ′( x n ) → Clamped condition
S j ( x) = a j + b j ( x − x j ) + c j ( x − x j ) 2 + d j ( x − x j ) 3 (6.43)
Denise Da Costa
– 115 –
Metodología de análisis
50
40
30
% Fuerza
20
10
0
0,0 3,0 6,0 9,0 12,0 15,0
% deformación
– Figura 6.17 –
Curva Fuerza versus Deformación para distintos tipos de cabos de amarre.
Denise Da Costa
– 116 –
Metodología de análisis
6.5.2 Defensas
δ δ
2
F = a 0 + a1 ⋅ + a 2 ⋅ (6.46)
L L
∂F δ
= EA = a1 + 2 ⋅ a 2 ⋅ (6.46a)
∂ (δ / L ) L
δ
en el instante inicial =0 → EA0 = a1
L
Para que esta aproximación tenga validez física, a1 debe ser positivo.
Denise Da Costa
– 117 –
Metodología de análisis
Dados los datos (δ/L ; F) de la curva de comportamiento, estos son aproximados a una
ecuación exponencial del tipo:
δ
a⋅
F = b⋅e L
(6.47)
δ
Ln(F ) = Ln(b ) + a ⋅ (6.47a)
L
δ
∂F a ⋅
= EA = b ⋅ a ⋅ e L (6.47b)
∂ (δ / L )
δ
para el instante inicial =0 → EA0 = b ⋅ a
L
Esta aproximación sólo tiene validez física cuando el producto b ⋅ a es positivo.
Dados los datos (δ/L ; F) de la curva de comportamiento, estos son aproximados a una
ecuación exponencial del tipo
δ
a
F = b ⋅ (6.48)
L
Denise Da Costa
– 118 –
Metodología de análisis
δ
Ln(F ) = Ln(b ) + a ⋅ Ln (6.48a)
L
( a −1)
∂F δ
= EA = b ⋅ a ⋅ (6.48b)
∂ (δ / L ) L
δ
para el instante inicial =0 → EA0 = 0
L
Derivando:
S ′( x) = bi − 2 ⋅ ci ⋅ xi + 3 ⋅ d i ⋅ xi2 + 2 ⋅ (ci − 3 ⋅ d i ⋅ xi ) ⋅ x + 3 ⋅ d i ⋅ x 2
Denise Da Costa
– 119 –
Metodología de análisis
δ EA0 = b0 − 2 ⋅ c0 ⋅ x0 + 3 ⋅ d 0 ⋅ x02
=0 →
L EA0 = b0
y1 − y 0 ( x1 − x 0 ) ⋅ (c1 + 2 ⋅ c0 )
b0 = −
x1 − x0 3
y1 x1 ⋅ (c1 + 2 ⋅ c0 )
b0 = −
x1 3
Para que la aproximación tenga validez físicamente, se tiene que cumplir que:
b0 > 0
y1 x1 ⋅ (c1 + 2 ⋅ c0 )
− >0
x1 3
3 ⋅ y1 − x12 ⋅ (c1 + 2 ⋅ c0 ) > 3 ⋅ x1 > 0
3 ⋅ y1 > x12 ⋅ (c1 + 2 ⋅ c0 )
3 ⋅ y1
(c1 + 2 ⋅ c0 ) <
x12
Denise Da Costa
– 120 –
Metodología de análisis
6.5.2.5 Conclusión:
En la figura 6.20 se presenta una curva para defensas la cual fue aproximada a una
parábola por mínimos cuadrados y también aproximada exponencialmente según la
aproximación semi-logarítmica. Pudiéndose observar las limitaciones de dichos métodos de
aproximación.
– Figura 6.18 –
Curva de comportamiento para defensas aproximada por varios métodos.
Por lo tanto, en el programa Amarre los datos de las curvas de comportamiento de los
cabos y defensas son aproximados por trozos rectos, para cada par de puntos,
determinando así la rigidez axial instantánea de cada elemento. Cabe destacar, que la
aproximación por rectas, a pesar de ser algo menos exacta, ofrece valores de EA inicial
siempre positivos; lo cual no sucede siempre con las aproximaciones por cubic spline.
Denise Da Costa
– 121 –
Metodología de análisis
Evidentemente que los cabos son elementos que si bien no resisten corte ni momento,
tampoco son capaces de resistir compresión. Por otro lado las defensas sólo trabajan a
compresión.
La incapacidad de los cabos para resistir compresión y de las defensas para resistir
tracción hace necesario un procedimiento particular de cálculo para adaptar el método de la
rigidez a este tipo de problemas.
Este problema consiste en restringir los movimientos del barco con la ayuda de vínculos
ficticios ubicados en el centro de masa y colocar los cabos con sus respectivas fuerzas de
tensión inicial.
Para obtener las reacciones en el vínculo ficticio, se deben descomponer las fuerzas
del tensado inicial y llevarse al centro de masa de la embarcación.
Denise Da Costa
– 122 –
Metodología de análisis
Las siguientes ecuaciones expresan las reacciones de los vínculos ficticios en el centro
de masa debido al tensado inicial:
n n
FX′ = ∑ Ti ⋅ λ X = ∑ FXi′ (6.49)
i =1 i =1
n n
FY′ = ∑ Ti ⋅ λY = ∑ FYi′ (6.50)
i =1 i =1
n
M ′ = ∑ (Fy i′ ⋅ xi′ − Fxi′ ⋅ y i′ ) (6.51)
i =1
Para obtener las reacciones en el extremo cercano a la bita, para este caso primario,
basta con descomponer la tensión inicial de los cabos en coordenadas globales
Q4 = −T ⋅ λ X (6.52)
Q5 = −T ⋅ λY (6.53)
Q1 = FX − FX´ (6.54)
Q2 = FY − FY´ (6.55)
Q3 = M − M ´ (6.56)
Denise Da Costa
– 123 –
Metodología de análisis
La solución del problema real viene dado por la suma de las soluciones del problema
primario más el problema complementario.
En definitiva el problema es tratado como si la embarcación, una vez alineada con las
defensas, se inmovilizara con vínculos artificiales, luego se tensaran las cuerdas, luego se
aplicarán las cargas de viento y oleaje, y por último se eliminarán los vínculos ficticios que
inicialmente lo inmovilizaban. Evidentemente que al desaparecer estos vínculos, si la
embarcación se encuentra bien sujeta, esta se desplazará hasta alcanzar una nueva
posición de equilibrio. La situación de fuerzas para esta última posición de equilibrio es la
que se determina en este trabajo.
Denise Da Costa
– 124 –
Metodología de análisis
– Figura 6.21 –
Esquema representativo de la consideración del pretensado de los cabos de amarre.
Denise Da Costa
– 125 –
Metodología de análisis
6.7.4 Ejemplo
Solución:
Denise Da Costa
– 126 –
Metodología de análisis
El caso complementario consiste en la estructura original sin las cargas internas. Aquí
se deben resolver el problema únicamente con las cargas externas originales (las cuales
son nulas en este ejemplo) junto con las reacciones obtenidas en el caso primario pero
aplicadas con signo contrario.
Denise Da Costa
– 127 –
Metodología de análisis
Del álgebra matricial se sabe que al haber dos o más filas linealmente dependientes en
la matriz, su determinante es nulo. Y del análisis estructural se sabe que cuando se está en
presencia de un problema inestable no es posible obtener un número de ecuaciones
linealmente independientes igual al número de incógnitas involucradas.
Denise Da Costa
– 128 –
Metodología de análisis
Denise Da Costa
– 129 –
Metodología de análisis
Denise Da Costa
– 130 –
Programa Amarre
PROGRAMA AMARRE
7.1 DESCRIPCIÓN
7.2.1 Buque
Denise Da Costa
– 131 –
Programa Amarre
En la figura 7.1 se representan de forma gráfica los datos de entrada del buque.
– Figura 7.1 –
Datos de entrada del Buque.
7.2.2 Cabos
Número de cabos.
X Bita: Coordenada en X global de la bita donde se amarra el cabo [m].
Y Bita: Coordenada en Y global de la bita donde se amarra el cabo [m].
x´ escoben: Coordenada en X local, referida al centro de masa del buque, del
escoben a donde llega el cabo [m].
y´ escoben: Coordenada en Y local, referida al centro de masa del buque, del
escoben a donde llega el cabo [m].
Curva Fuerza vs Deformación: Depende del material del cabo; y es introducida por
el usuario mediante doce pares de puntos (δ/L ; F), donde las fuerzas deben ser
dadas en orden ascendente.
Tensado Inicial: Fuerza de pretensado al momento de amarre [T].
En la figura 7.2 se representan de forma gráfica los datos de entrada del cabo.
Denise Da Costa
– 132 –
Programa Amarre
– Figura 7.2 –
Datos de entrada del Cabo.
7.2.3 Defensas
– Figura 7.3 –
Datos de entrada de la defensa.
Denise Da Costa
– 133 –
Programa Amarre
Número de cadenas.
X Boya: Coordenada en X global de la boya [m].
Y Boya: Coordenada en Y global de la boya [m].
x´ escoben: Coordenada en X local, referida al centro de masa del buque, del
escoben a donde llega la cadena [m].
y´ escoben: Coordenada en Y local, referida al centro de masa del buque, del
escoben a donde llega la cadena [m].
Curva de comportamiento: Para generar la curva de comportamiento de las
cadenas de anclaje y determinar la fuerza horizontal máxima del sistema de
anclaje, el usuario debe introducir: profundidad del agua (D), longitud de la parte
superior de la cadena (L2), longitud de la parte inferior de la cadena (L1), peso por
unidad de longitud de cada una de las partes de la cadena y fuerza máxima que la
cadena puede resistir dependiendo del material, peso del sinker (w), peso del
ancla y ángulo horizontal en el ancla.
En la figura 7.4 se representa de forma gráfica, los datos de entrada para el sistema de
anclaje con cadenas.
– Figura 7.4 –
Datos de entrada de las cadenas.
7.3 PROCEDIMIENTO
Denise Da Costa
– 134 –
Programa Amarre
A continuación se explica con detalle las subrutinas más importantes del programa.
El programa considera el tensado inicial de los cabos; para lo cual se resuelve el caso
complementario sumándole a las fuerzas externas aplicadas al buque el tensado de los
cabos con signo contrario y trasladado al centro de masa de la embarcación. Una vez
obtenido los resultados de este caso, se suma con el caso primario, el cual sólo contiene las
fuerzas de tensado inicial.
Denise Da Costa
– 135 –
Programa Amarre
El programa Amarre aproxima las curvas de comportamiento por trozos de recta; esto
es, asume que cada par de puntos forman una recta y calcula la pendiente de esta para así
tener el valor de EA para cada par puntos.
Una vez obtenidos los valores de EA para cada par de puntos, el programa asigna el
valor de EA inicial para cada elemento y calcula la matriz de rigidez de la estructura así
como las fuerzas en cada cabo y defensa. Después de obtener el equilibrio del sistema,
Amarre verifica si el valor de la fuerza está o no dentro del rango asumido; caso no esté, el
programa verifica en que rango se encuentra la fuerza y determina el nuevo valor de EA; si
el nuevo valor de EA tiene un error mayor del 2% respecto al valor anterior, el programa
sustituye este valor por el nuevo y vuelve a determinar la matriz de rigidez de la estructura y
todas las fuerzas en cada elemento. Este proceso iterativo es realizado hasta conseguir la
convergencia del problema no lineal.
Dado que los cabos y las defensas son elementos que no trabajan a compresión ni a
tracción, respectivamente, se deben eliminar del cálculo aquellos cabos cuya fuerza no sea
de tracción y aquellas defensas cuya fuerza no sea de compresión. Para esto se creó una
subrutina que elimina todos los términos en la matriz de rigidez de la estructura
pertenecientes al cabo más comprimido o a la defensa más traccionada, y recalcula todas
las fuerzas en los demás elementos; esta subrutina se repite tantas veces como cabos a
compresión y defensas a tracción existan.
Las fuerzas en las bitas son calculadas en coordenadas globales. En el caso que más
de un cabo este amarrado a la misma bita el programa lo detectará ya que contiene una
subrutina para comparar las coordenadas de todas las bitas y contabilizar una sola vez
aquellas cuyas coordenadas X e Y sean iguales, de este modo los resultados son
presentados de forma total, es decir, se determina la sumatoria de fuerzas parciales de cada
cabo que llega a la bita para obtener la fuerza total de la bita. Esta fuerza, como bien es
sabido, es la misma fuerza del cabo (en el extremo cercano a la bita) en sentido contrario y
descompuesta en los ejes X e Y globales del sistema.
Denise Da Costa
– 136 –
Programa Amarre
Arroja las fuerzas de amarre en las bitas, en coordenadas X e Y globales. Estas son las
fuerzas que debe resistir la estructura de amarre
El mismo ejemplo será estudiado con rigor en el próximo capítulo (Estudio Paramétrico)
presentando el análisis de resultados.
Denise Da Costa
– 137 –
Programa Amarre
Como quiera que el sistema temporal de amarre estará en uso sólo por un período
limitado de tiempo, la opción de amarre con ancla implicará el uso de un cabo largo de proa.
La solución definitiva, representada por la reconstrucción del duque de amarre dañado,
retomará el esquema de amarre originalmente diseñado.
– Figura 7.5 –
Datos de entrada del buque.
Denise Da Costa
– 138 –
Programa Amarre
– Figura 7.6 –
Datos de entrada de los cabos.
Denise Da Costa
– 139 –
Programa Amarre
– Figura 7.7 –
Datos de entrada de la curva de comportamiento de los cabos.
– Figura 7.8 –
Datos de entrada de las defensas.
Denise Da Costa
– 140 –
Programa Amarre
– Figura 7.9 –
Datos de entrada de la curva de comportamiento de las defensas.
Denise Da Costa
– 141 –
Programa Amarre
– Figura 7.10 –
Datos de entrada del cabo + cadena de anclaje.
Denise Da Costa
– 142 –
Programa Amarre
– Figura 7.11 –
Datos de entrada del sistema cadena y ancla.
El programa Amarre posee dos versiones de salida de datos, una en formato de Visual
Basic y otra en formato de Excel para mayor comodidad de impresión y tratamiento de
datos.
– Figura 7.12 –
Salida de datos en Visual Basic.
Denise Da Costa
– 143 –
Programa Amarre
– Figura 7.13 –
Salida de datos en Excel.
Denise Da Costa
– 144 –
Programa Amarre
Inicio
Datos de Entrada
Geométricos.
Fuerzas
aplicadas al buque en
su centro de masa.
Curvas
"Fuerza vs
Deformación"
Faltan Datos.
¿Los datos están NO Introduzca los Datos del
completos ? problema
SI
Advertencia:
SI
Determinante = 0 La estructura es inestable.
Se suspende la ejecución del programa
NO
SI Defensas a Tracción.
Eliminar
Cabos a Compresión.
NO
NO Cálculo de las
Error > 2% fuerzas globales
en las Bitas.
SI
Resultados:
Sustituir EA (i) por Desplazamiento
EA(i+1) globales del buque.
Versión
Fuerzas locales
Imprimible de los
en cada miembro.
resultados.
Fuerzas totales en
las bitas.
Fin
Denise Da Costa
– 145 –
Programa Amarre
Denise Da Costa
– 146 –
Estudio paramétrico
ESTUDIO PARAMETRICO
Algunos de los datos utilizados en los problemas fueron ofrecidos por la empresa
NOUEL (Caracas, Venezuela), dedicada a la ingeniería de obras marítimas y portuarias;
permitiendo así la realización de este capítulo con datos de proyectos verdaderos y
agregando más valor al trabajo aquí realizado.
Denise Da Costa
– 147 –
Estudio paramétrico
– Figura 8.1 –
Configuración de amarre 1.
– Figura 8.2 –
Configuración de amarre 2.
– Figura 8.3 –
Configuración de amarre 3.
Denise Da Costa
– 148 –
Estudio paramétrico
– Figura 8.4 –
Configuración de amarre 4.
8.1.2 Buques
Los buques que serán utilizados para realizar el presente estudio paramétrico son de
85.000 y 130.000 toneladas de peso muerto; tanto los buques como las condiciones
ambientales fueron obtenidos de dos proyectos reales realizados por la empresa NOUEL, la
cual facilitó esta información a fin de dar sentido práctico a los ejemplos aquí realizados.
– Figura 8.5 –
Dimensiones de buques.
Denise Da Costa
– 149 –
Estudio paramétrico
Coeficientes
Areas
– Tabla 8.1 –
Cálculo de las fuerzas ambientales para el tanquero de 85.000 DWT
Denise Da Costa
– 150 –
Estudio paramétrico
⋅ ρ a ⋅ Vw2 ⋅ Ax ⋅ C xw ⋅ f xw (θ w )
Viento: 1
Fxw =
2
⋅ ρ a ⋅ Vw2 ⋅ Ay ⋅ C yw ⋅ f yw (θ w )
1
Fyw =
2
⋅ ρ a ⋅ Vw2 ⋅ Ay ⋅ L ⋅ C xyw (θ w )
1
M xyw =
2
e
M xyc = F yc ⋅ L wL ⋅ c
L wL
Combinaciones
Denise Da Costa
– 151 –
Estudio paramétrico
Coeficientes
Areas
– Tabla 8.2 –
Cálculos de las fuerzas ambientales sobre el buque de 130.000 DWT.
Denise Da Costa
– 152 –
Estudio paramétrico
⋅ ρ a ⋅ Vw2 ⋅ Ax ⋅ C xw ⋅ f xw (θ w )
Viento: 1
Fxw =
2
⋅ ρ a ⋅ Vw2 ⋅ Ay ⋅ C yw ⋅ f yw (θ w )
1
Fyw =
2
⋅ ρ a ⋅ Vw2 ⋅ Ay ⋅ L ⋅ C xyw (θ w )
1
M xyw =
2
e
M xyc = F yc ⋅ L wL ⋅ c
L wL
Combinaciones
Denise Da Costa
– 153 –
Estudio paramétrico
8.1.3 Defensas
Fuerza Deformación
(%) (%)
0 0
5 10
10 15
20 24
30 30
40 36
50 41
60 45
70 49
80 53
90 57
100 60
– Tabla 8.3 –
Comportamiento de la defensa tipo Foam expresado en
% de Fuerza vs % de deformación (Fentek, 2002)
– Tabla 8.4 –
Fuerza máxima resistida por cada defensa según el tipo y tamaño (Fentek, 2002)
Denise Da Costa
– 154 –
Estudio paramétrico
CURVA 1
Curva de Comportamiento
Defensa tipo Foam
Deformación
Fuerza (t) relativa δ/L
(%) 25,0
0,000 0,000
20,0
1,015 10,000
2,030 15,000
Fuerza (t)
15,0
4,060 24,000
6,090 30,000
10,0
8,120 36,000
10,150 41,000 5,0
12,180 45,000
14,210 49,000 0,0
16,240 53,000 0 20 40 60 80
18,270 57,000
% deformación
20,300 60,000
– Figura 8.6 –
Curva de comportamiento de la defensa 1 (tipo Foam).
Fuerza Deformación
(%) (%)
0 0
10 2,5
18 5
20 8
30 18
40 25
50 31
60 37
70 41
80 45
90 48
100 50
– Tabla 8.5 –
Comportamiento de la defensa tipo D & square expresado en
% de Fuerza vs % de deformación (Fentek, 2002)
Denise Da Costa
– 155 –
Estudio paramétrico
– Tabla 8.6 –
Fuerza máxima resistida por cada defensa según el tipo y tamaño (Fentek, 2002)
CURVA 2
Curva de Comportamiento
Defensa D & Square
Deformación
Fuerza (t) relativa δ/L
(%) 25,00
0,000 0,000
2,100 2,500 20,00
3,780 5,000
Fuerza (t)
15,00
4,200 8,000
6,300 18,000
10,00
8,400 25,000
10,500 31,000
5,00
12,600 37,000
14,700 41,000 0,00
16,800 45,000 0 10 20 30 40 50 60
18,900 48,000
% deformación
21,000 50,000
– Figura 8.7 –
Curva de comportamiento de la defensa 2 (tipo Square).
Denise Da Costa
– 156 –
Estudio paramétrico
Fuerza Deformación
(%) (%)
0 0
75 20
89 25
97 30
100 35
98 40
92 45
77 55
73 60
77 65
91 70
100 72
– Tabla 8.7 –
Comportamiento de la defensa tipo Super Cone expresado en
% de Fuerza vs % de deformación (Fentek, 2002)
SCN 500
R (kN) R (t) E (kNm)
E 0.9 164 16,7 36,5
E 1.0 182 18,6 40,5
E 1.1 187 19,1 41,9
E 1.2 191 19,5 43,2
E 1.3 196 20,0 44,6
E 1.4 200 20,4 45,9
E 1.5 205 20,9 47,3
E 1.6 209 21,3 48,6
E 1.7 214 21,8 50,0
E 1.8 218 22,2 51,3
E 1.9 223 22,7 52,7
E 2.0 227 23,1 54,0
– Tabla 8.8 –
Fuerza máxima resistida por cada defensa según el tipo y tamaño (Fentek, 2002)
Denise Da Costa
– 157 –
Estudio paramétrico
CURVA 3
Curva de Comportamiento
Defensa Super Cone
Deformación
Fuerza (t) relativa δ/L
(%) 25,00
0,000 0,000
15,300 20,000 20,00
18,156 25,000
Fuerza (t)
15,00
19,788 30,000
20,400 35,000
10,00
19,992 40,000
18,768 45,000
5,00
15,708 55,000
14,892 60,000 0,00
15,708 65,000 0 20 40 60 80
18,564 70,000
% deformación
20,400 72,000
– Figura 8.8 –
Curva de comportamiento de la defensa 3 (tipo Super Cone).
8.1.4 Cabos
Curva de Comportamiento
Fuerza Deformación Cabos de nylon
(%) (%)
0 0,0 120
5 4,3
100
10 7,5
15 9,6 80
Fuerza (%)
20 11,5
60
25 13,2
30 14,5 40
35 15,7
20
40 17,0
60 22 0
80 26 0,0 5,0 10,0 15,0 20,0 25,0 30,0 35,0
– Figura 8.9 –
Comportamiento de los cabos de nylon expresado en
% de Fuerza vs % de deformación (Puget Sound Rope, 2003)
Denise Da Costa
– 158 –
Estudio paramétrico
– Tabla 8.9 –
Características de los cabos de nylon (Puget Sound Rope, 2003)
Denise Da Costa
– 159 –
Estudio paramétrico
CURVA 1 CURVA 2
Deformación Deformación
Fuerza (t) relativa δ/L Fuerza (t) relativa δ/L
(%) (%)
CURVA 3 CURVA 4
Deformación Deformación
Fuerza (t) relativa δ/L Fuerza (t) relativa δ/L
(%) (%)
– Tabla 8.10 –
Datos Fuerza vs % deformación para los 4 cabos a considerar.
Denise Da Costa
– 160 –
Estudio paramétrico
De esta forma cada caso en estudio queda definido por un código donde se incorporan
las 4 letras que definen la variable a utilizar; por ejemplo, Buque de 85.000 DWT amarrado
con una configuración como la representada en la figura 8.2, utilizando cabos de nylon de Ø
= 56 mm y defensas tipo square; el código seria: B1c2.C3D2
Como variables fijas fueron utilizadas, buque de 85.000 DWT amarrado según la
configuración presentada en la figura 8.3 con cabos de nylon de Ø = 104 mm.
Denise Da Costa
– 161 –
Estudio paramétrico
Configuración 3
Buque
X CM 118,5 m
Y CM 51,0 m
Defensa X Li
1 78,00 1,5
2 111,00 1,5
3 146,00 1,5
8.2.2 Resultados
18,00
16,00
14,00
Compresión (t)
12,00
10,00
8,00
6,00
4,00
2,00
0,00
1 2 3
Defensa
– Figura 8.10 –
Fuerzas en las defensas variando el tipo de material.
Denise Da Costa
– 162 –
Estudio paramétrico
Fuerza en el cabo 1
3,00
2,50
2,00
Tensión (t)
1,50
1,00
0,50
0,00
1
Cabo
– Figura 8.11 –
Fuerza en el cabo 1 para el análisis del comportamiento de las defensas.
50,80
50,60
50,40
Y CM
50,20
50,00
49,80
49,60
118,30 118,35 118,40 118,45 118,50 118,55
X CM
– Figura 8.12 –
Posición del centro de masa del buque para el análisis del comportamiento de las defensas.
Denise Da Costa
– 163 –
Estudio paramétrico
A continuación se presentan los resultados obtenidos para los seis casos analizados.
8.2.2.1 Defensa 1
B1c3.C4D1 – Lineal
B1c3.C4D1 – No lineal
Denise Da Costa
– 164 –
Estudio paramétrico
8.2.2.2 Defensa 2
B1c3.C4D2 – Lineal
Denise Da Costa
– 165 –
Estudio paramétrico
B1c3.C4D2 – No lineal
Denise Da Costa
– 166 –
Estudio paramétrico
8.2.2.3 Defensa 3
B1c3.C4D3 – Lineal
B1c3.C4D3 – No lineal
Denise Da Costa
– 167 –
Estudio paramétrico
En la figura 8.10 podemos observar que la compresión en las defensas varia para el
mismo material dependiendo si el análisis es lineal o no lineal. De las tablas de resultados
se puede observar que para la defensa 1 (foam) existen variaciones hasta de 18% entre el
análisis lineal y no lineal; para la defensa 2 (square) existen diferencias entre 11% y 19% en
las compresiones de las defensas cuando se varia de análisis lineal a no lineal y por último
podemos observar que la compresión en las defensas tipo 3 (super cone) varia en un
máximo de 31% debido a la alta no linealidad que presentan provocando que el análisis
lineal introduzca un alto grado de error.
La no linealidad de las defensas también influye en las fuerzas resistidas por los cabos
que a su vez serán transmitidas a las bitas; analizando el cabo #1, por ser este uno de los
más representativos y estar siempre trabajando, se puede concluir que a medida que la
curva de comportamiento de las defensas presenta una alta no linealidad las fuerzas en los
cabos presentan mayores variaciones de un análisis a otro, es decir, para el tipo de defensa
1 (foam) se obtienen diferencias en las fuerzas de tensión de los cabos de 4% mientras que
para el caso de defensa 3 (super cone) la diferencia entre las tensiones de los cabos es del
orden de 31%
Por todo lo expuesto se puede concluir que un análisis más riguroso incluyendo la no
linealidad en el comportamiento de las defensas es necesario ya que se obtienen resultados
de fuerzas significativamente diferentes y que pueden condicionar el diseño del muelle.
Denise Da Costa
– 168 –
Estudio paramétrico
En este caso las fuerzas externas aplicadas fueron colocadas con signo positivo en la
dirección Y global de forma que las defensas no trabajan por estar en tracción, pero en
contra partida se puede realizar un mejor análisis del comportamiento de los cabos por estar
todos trabajando simultáneamente.
Debido a que los cabos de amarre son unos de los factores primordiales en la
configuración de amarre, también se analizaron tres tipos de configuraciones.
Buque
X CM 126,6 m
Y CM 53,8 m
Denise Da Costa
– 169 –
Estudio paramétrico
Configuración 2
Buque
X CM 133,39 m
Y CM 56,50 m
Configuración 3
Buque
X CM 118,5 m
Y CM 51,0 m
Denise Da Costa
– 170 –
Estudio paramétrico
8.3.2 Resultados
56,00
55,50
55,00
Y CM
54,50
54,00
53,50
125,00 125,50 126,00 126,50 127,00
X CM
– Figura 8.13 –
Posición del CM del buque para el análisis de los cabos (Configuración 1)
57,00
56,50
56,00
Y CM
55,50
55,00
54,50
54,00
130,00 130,50 131,00 131,50 132,00 132,50 133,00 133,50 134,00
X CM
– Figura 8.14 –
Posición del CM del buque para el análisis de los cabos (Configuración 3)
Denise Da Costa
– 171 –
Estudio paramétrico
59,50
59,00
58,50
58,00
57,50
57,00
56,50
56,00
Y CM
55,50
55,00
54,50
54,00
53,50
53,00
52,50
52,00
51,50
51,00
50,50
114,50 115,00 115,50 116,00 116,50 117,00 117,50 118,00 118,50 119,00
X CM
– Figura 8.15 –
Posición del CM del buque para el análisis de los cabos (Configuración 3)
8.3.2.1 Configuración 1
B1c1.C1
Desplazamientos globales del buque
Desplazamiento en X: -1,08 m
Desplazamiento en Y: 1,87 m
Rotación: 1,06 º
Posición de equilibrio
X CM: 125,52
Y CM: 55,67
B1c1.C2
Desplazamientos globales del buque
Desplazamiento en X: -0,31 m
Desplazamiento en Y: 0,53 m
Rotación: 0,30 º
Posición de equilibrio
X CM: 126,29
Y CM: 54,33
Denise Da Costa
– 172 –
Estudio paramétrico
B1c1.C3
Desplazamientos globales del buque
Desplazamiento en X: -0,09 m
Desplazamiento en Y: 0,15 m
Rotación: 0,00 º
Posición de equilibrio
X CM: 126,51
Y CM: 53,95
B1c1.C4
Desplazamientos globales del buque
Desplazamiento en X: -0,03 m
Desplazamiento en Y: 0,05 m
Rotación: 0,03 º
Posición de equilibrio
X CM: 126,57
Y CM: 53,85
Cabo Tensión
1 4,71
2 5,05
3 5,21
4 8,45
5 0,65
6 5,93
7 6,88
8 7,05
Denise Da Costa
– 173 –
Estudio paramétrico
8.3.2.2 Configuración 2
B1c2.C1
Desplazamientos globales del buque
Desplazamiento en X: -2,93 m
Desplazamiento en Y: -2,02 m
Rotación: 0,49 º
Posición de equilibrio
X CM: 130,46
Y CM: 54,48
B1c2.C2
Desplazamientos globales del buque
Desplazamiento en X: -0,83 m
Desplazamiento en Y: -0,57 m
Rotación: 0,00 º
Posición de equilibrio
X CM: 132,56
Y CM: 55,93
B1c2.C3
Desplazamientos globales del buque
Desplazamiento en X: -0,24 m
Desplazamiento en Y: -0,17 m
Rotación: 0,00 º
Posición de equilibrio
X CM: 133,15
Y CM: 56,33
B1c2.C4
Denise Da Costa
– 174 –
Estudio paramétrico
Cabo Tensión
1 3,44
2 3,03
3 4,79
4 3,05
5 6,30
6 3,56
7 4,73
8 3,89
9 5,12
10 4,73
8.3.2.3 Configuración 3
B1c3.C1
Desplazamientos globales del buque
Desplazamiento en X: -3,58 m
Desplazamiento en Y: 7,63 m
Rotación: 1,83 º
Posición de equilibrio
X CM: 114,92
Y CM: 58,63
Denise Da Costa
– 175 –
Estudio paramétrico
B1c3.C2
Desplazamientos globales del buque
Desplazamiento en X: -1,01 m
Desplazamiento en Y: 2,16 m
Rotación: 0,52 º
Posición de equilibrio
X CM: 117,49
Y CM: 53,16
B1c3.C3
Desplazamientos globales del buque
Desplazamiento en X: -0,29 m
Desplazamiento en Y: 0,62 m
Rotación: 0,00 º
Posición de equilibrio
X CM: 118,21
Y CM: 51,62
B1c3.C4
Desplazamientos globales del buque
Desplazamiento en X: -0,10 m
Desplazamiento en Y: 0,21 m
Rotación: 0,00 º
Posición de equilibrio
X CM: 118,40
Y CM: 51,21
Denise Da Costa
– 176 –
Estudio paramétrico
Cabo Tensión
1 5,89
2 7,41
3 6,40
4 4,10
5 10,57
6 8,85
Bien para cualquiera de las tres configuraciones, se obtuvieron iguales valores para las
fuerzas en los cabos y las correspondientes fuerzas transmitidas a las bitas
independientemente del tipo de cabos. Estos resultados no son de extrañar porque las
fuerzas en las reacciones dependen de las fuerzas aplicadas y la posición en los apoyos.
Para sistemas hiperestáticos considerando el buque como elemento rígido (como es el
caso) las fuerzas son siempre transmitidas de la misma forma a los elementos flexibles
(cabos) y la distribución de esfuerzos depende de la configuración de amarre; si
consideramos el buque como un elemento flexible (en el concepto de ingeniería civil)
obtendríamos pequeñas diferencias en las fuerzas resistidas por los elementos al variar la
rigidez de estés.
Analizando las fuerzas en las bitas que amarran los cabos 1 y 2 y la última bita de cada
sistema, por estar estas ubicadas en una posición semejante, podemos observar que las
fuerzas resistidas por las bitas en el caso de la configuración 3 son 10% y hasta 80%
mayores que en el caso de las configuraciones 1 y 2.
Denise Da Costa
– 177 –
Estudio paramétrico
Comparando las fuerzas máximas de los cabos de amarre en cada caso, se puede
destacar que la configuración 3 tiene fuerzas 40% mayores que la configuración 2 y 20%
mayores a las fuerzas máximas de los cabos en la configuración 1. A su vez, los cabos en la
configuración 1 poseen fuerzas máximas 25% mayores a las fuerzas máximas de los cabos
en la configuración 2. De este modo quedaría la configuración 2 con las menores fuerzas de
tensión en los cabos.
En conclusión, los sistemas de amarre con mayor número de cabos tienen una
distribución de esfuerzos más favorable provocando, en general, menores fuerzas en las
bitas y menor desplazamiento del buque. La rigidez de los cabos apenas contribuye a
obtener posiciones finales del buque más próximas de la posición inicial a medida que esta
aumenta.
En este estudio se comparan los dos buques para la misma configuración y los mismos
materiales constituyentes de los elementos de amarre y defensas. Con el objetivo de
verificar las fuerzas en las bitas para cada caso y estudiar el efecto que podría tener atracar
y amarrar buques de diferentes categorías en el mismo puerto; se utilizó la configuración de
amarre 1 por ser clásica, con cabos de naylon Ø = 104 mm por ser los más rígidos y dos
tipos de defensas, curva 2 y 3 que presentan curvas de comportamiento con diferentes
grados de no linealidad y por último las fuerzas aplicadas en cada buque corresponden a
sus respectivas fuerzas estáticas equivalentes para el caso de buque cargado por ser este
el caso más desfavorable.
Denise Da Costa
– 178 –
Estudio paramétrico
Configuración 1
Defensa X Li
1 93,9 1,5
2 126,6 1,5
3 159,3 1,5
Defensa X Li
1 93,9 1,5
2 126,6 1,5
3 159,3 1,5
Denise Da Costa
– 179 –
Estudio paramétrico
8.4.2 Resultados
54,00
53,50
Y CM
53,00
52,50
126,50 126,55 126,60 126,65 126,70
X CM
– Figura 8.16 –
Posición del CM del buque para el análisis de comparación de buques.
– Figura 8.17 –
Fuerza en los cabos “spring” para el análisis de comparación de buques.
Denise Da Costa
– 180 –
Estudio paramétrico
8.4.2.1 Buque 1
B1c1.C4D2
Desplazamientos globales del buque
Desplazamiento en X: 0,00 m
Desplazamiento en Y: -0,26 m
Rotación: -0,02 º
Posición de equilibrio
X CM: 126,60
Y CM: 53,54
B1c1.C4D3
Desplazamientos globales del buque
Desplazamiento en X: 0,04 m
Desplazamiento en Y: -0,43 m
Rotación: -0,13 º
Posición de equilibrio
X CM: 126,64
Y CM: 53,37
Denise Da Costa
– 181 –
Estudio paramétrico
8.4.2.2 Buque 2
B2c1.C4D2
Desplazamientos globales del buque
Desplazamiento en X: -0,05 m
Desplazamiento en Y: -1,02 m
Rotación: -0,61 º
Posición de equilibrio
X CM: 126,55
Y CM: 52,78
Denise Da Costa
– 182 –
Estudio paramétrico
B2c1.C4D3
Denise Da Costa
– 183 –
Estudio paramétrico
Tanto en la figura 8.17 como en las tablas de resultados se puede comprobar que hay
una gran diferencia en las fuerzas de reacciones en los elementos y en las bitas, cuando
son comparados los dos tipos de buques en estudio. Obteniendo fuerzas de más del triple
para el caso del buque 2; nótese que la fuerza máxima en los cabos no es superada en
ninguno de los casos en estudio porque los cabos a ser utilizados son los de mayor
resistencia; cabe destacar que el uso de cabos con diámetro inferior a 52 mm sería riesgoso
pues para el caso de amarrar el buque de 130.000 DWT estos romperían, poniendo en
causa el sistema de amarre.
Para esto se agregaron dos cadenas de popa y proa (sentido longitudinal) y dos
cadenas de anclaje en sentido lateral. Las cuatro cadenas utilizadas tienen una longitud de
150,8 m ancladas al fondo marino con un ancla tipo Stockless de 6.000 lb de peso.
Denise Da Costa
– 184 –
Estudio paramétrico
Defensa X Li
1 93,9 1,5
2 126,6 1,5
3 159,3 1,5
Denise Da Costa
– 185 –
Estudio paramétrico
8.5.3 Resultados
60
50
40
Tensión (t)
30
20
10
0
1 2 3 4 5
Cabo
– Figura 8.18 –
Tensión en los cabos para el análisis de comparación de buque
e incluyendo cadenas en el buque 2.
Denise Da Costa
– 186 –
Estudio paramétrico
Denise Da Costa
– 187 –
Estudio paramétrico
La opción de amarre alterno por boya, cadena y ancla consiste en una boya, un tramo
de 110 m de cadena de acero de 2,5´´ de diámetro y un ancla tipo Danfor de 15.000 lb de
peso. Para ubicar correctamente el ancla se debe considerar que el cabo y la cadena se
ubicarán frente a los duques de amarre existentes, tal como se indica en las figuras 8.19 y
8.20.
Con el fin de realizar también una comparación con las opciones de amarre estudiadas
en el apartado 8.3 Comportamiento de los cabos, fue utilizado el buque 1 (85.000 DWT) y
cabos de amarre tipo 4.
Denise Da Costa
– 188 –
Estudio paramétrico
– Figura 8.19 –
Amarre alterno 1.
– Figura 8.20 –
Amarre alterno 2.
Denise Da Costa
– 189 –
Estudio paramétrico
Buque
XCM YCM
133,39 56,50
Buque
XCM YCM
126,6 53,8
Denise Da Costa
– 190 –
Estudio paramétrico
Denise Da Costa
– 191 –
Estudio paramétrico
8.6.3 Resultados
56,9
56,8
56,7
Y CM
56,6
56,5
56,4
133,2 133,3 133,4
X CM
– Figura 8.21 –
Posición del centro de masa del buque para el análisis entre
amarre sin cadenas de anclaje y amarre alterno
(Configuración 2).
Denise Da Costa
– 192 –
Estudio paramétrico
54,2
54,1
54
Y CM
53,9
53,8
53,7
126,5 126,6 126,7
X CM
– Figura 8.22 –
Posición del centro de masa del buque para el análisis entre
amarre sin cadenas de anclaje y amarre alterno
(Configuración 1).
Denise Da Costa
– 193 –
Estudio paramétrico
Denise Da Costa
– 194 –
Estudio paramétrico
Denise Da Costa
– 195 –
Estudio paramétrico
Del estudio realizado se puede observar que el sistema alterno compuesto por boya,
cadena y ancla, para ambos ejemplos, resiste tanto las solicitaciones en la condición del
buque vacío como las solicitaciones en la condición del buque cargado, obteniendo un
factor de seguridad de 1,5 en la condición más desfavorable, lo que permite asegurar el uso
de este tipo de sistemas en situaciones de emergencia o alternas como las presentadas en
ambos casos.
Comparando las soluciones presentadas en el apartado 8.3 con las aquí desarrolladas
como amarres alternos, podemos notar que la posición final del buque varia
sustancialmente entre la solución de amarre sin ancla con la solución de amarre alterno,
existiendo en esta última desplazamientos de hasta 6 veces mayores que en la primera;
valores debidos a que el sistema con boya y cadena de anclaje es un sistema muy elástico
comparado al sistema amarrado en una estructura rígida (ver figuras 8.21 y 8.22).
Por otra parte se puede constatar que las fuerzas resistidas por los cabos en ambas
situaciones de amarre varían entre 8% y 60% comprobando así que la distribución de
esfuerzos varía de forma poco influyente, ya que los cabos no superan de ninguna forma
sus valores de resistencia máximos. También se puede observar que el sistema de boya y
cadena de anclaje coparticipa en la distribución de esfuerzos, resistiendo parte de las
fuerzas externas a ser resistidas por el sistema y proporcionando estabilidad, la cual fue
puesta en riesgo en el momento que uno o más de los duques de amarre quedó fuera de
utilización.
Denise Da Costa
– 196 –
Estudio paramétrico
Los datos del sistema flexible, ancla, tipo de cadena y características, así como parte
de la configuración de amarre fueron obtenidos de ejemplos desarrollados por NFEC (1985).
Para el sistema flexible se utilizó en el sentido longitudinal un tramo de 178,3 m de cadena
de 1¼´´ de diámetro con ancla tipo Stato de 3.307 lb de peso; en el sentido lateral se utilizó
un tramo de 164,6 m de cadena de 2¼´´ de diámetro con ancla tipo Stato de 12.894 lb de
peso.
Fuerzas
Fx -5,06 t
Fy 30,16 t
Mxy 298,92 t.m
Buque
XCM YCM
186,70 71,00
Buque
XCM YCM
186,70 71,00
Denise Da Costa
– 197 –
Estudio paramétrico
Denise Da Costa
– 198 –
Estudio paramétrico
Denise Da Costa
– 199 –
Estudio paramétrico
8.7.4 Resultados
72,00
71,80
71,60
Y CM
71,40
71,20
71,00
70,80
186,40 186,50 186,60 186,70 186,80
X CM
– Figura 8.23 –
Posición final del buque al comparar sistema de amarre rígido y flexible.
B1c4.C4
Denise Da Costa
– 200 –
Estudio paramétrico
Cabo Tensión
1 0,00
2 13,79
3 16,37
4 5,06
B1c4.C4
Al comparar la posición final del buque para ambos sistemas, se puede verificar que el
sistema flexible presenta desplazamientos en la dirección X de 66,7% mayores que los
desplazamientos para el sistema rígido, de igual forma el sistema flexible tiene
desplazamientos en la dirección Y de 16,5% mayores que los desplazamientos obtenidos en
el sistema rígido. Lo que lleva a concluir que los sistemas flexibles presentan el
inconveniente de permitir mayores desplazamientos del buque, factor de gran importancia
debido a que el objetivo del sistema de amarre es impedir grandes desplazamientos
relativos y absolutos del buque.
Denise Da Costa
– 201 –
Estudio paramétrico
Denise Da Costa
– 202 –
Conclusiones y Recomendaciones
CONCLUSIONES Y RECOMENDACIONES
9.1 CONCLUSIONES
Conocer el tipo de buques que irán a atracar en una determinada estructura de amarre
es importante en el estudio previo de la misma porque, según resultados obtenidos en el
capítulo 8, la diferencia de toneladas entre buques provoca grandes diferencias en las
solicitaciones de las bitas poniendo en causa la seguridad del sistema. Es importante
recordar que la incorporación de cadenas de anclaje al sistema no es garantía para
asegurar el amarre de buques mayores a muelles que no han sido diseñados para este fin;
y en algunos casos el amarre a boya presenta el inconveniente de ser inoperante debido a
que es un sistema muy elástico que debería trabajar en conjunto con el resto de los cabos
que, al estar sujetos a estructuras, representa un sistema muy rígido. Además, del punto de
vista operativo no es tan seguro como el anterior ya que los operadores deberán acceder a
la boya para amarrar y soltar los cabos respectivos y esta operación, en condiciones de
oleaje medio a fuerte, es riesgosa.
Cabe destacar que, según la OCIMF y otras publicaciones que manejan temas
marítimos, la tendencia moderna para el amarre de buques es la de evitar el uso de los
cabos largos tanto de proa como de popa. El sistema de amarre más eficiente, en forma
esquemática, se ilustra en la Figura 5.6 (capítulo 5 – Sistemas de amarre).
Denise Da Costa
– 203 –
Conclusiones y Recomendaciones
9.2 RECOMENDACIONES
Denise Da Costa
– 204 –
Referencias
REFERENCIAS
ARGESCHOU, H, et al (1984) – Planning and design of ports and marine terminals, John
Wiley & Sons, USA.
BARROS, Rui (1997) – Approximate longitudinal motions of moored vessels under seiches,
XVIII Congreso Ibero latinoamericano de métodos computacionales en ingeniería, XVIII
CILAMCE, Brasilia, Brasil.
BARROS, Rui (1999-a) – First order heave response of a free floating spar buoy in a random
th
seaway, 7 EPMESC, Universidad de Macau, Macau, China.
BARROS, Rui (1999-b) – Acções hidrodinâmicas sobre cilindros verticais e pipelines sob a
acção da corrente e da ondulação, Congresso Luso-Moçambicano de Engenharia, Maputo,
Moçambique.
BARROS, Rui (2001-a) - Random response of a spar buoy in heave, Artigo publicado na
Revista Engenharia: Estudo e Pesquisa, Rio de Janeiro, Brasil.
BURDEN et al. (1978) – Numerical Analysis, Prindle Weber & Schmidt, Boston,
Massachusetts, USA
CHAPRA & CANALE (1998) – Numerical Methods for Engineers: with programming and
software applications, 3ra edición, Mac Graw Hill.
Denise Da Costa
– 205 –
Referencias
DERUCHER & HEINS (1979) – Bridge and pier protective systems and devices, Marcel
Dekker, New York, USA.
GAYTHWAITE, John (1981) – The marine environment and structural design, Van Nostrand
Reinhold Company, USA.
HIBBELER, Russell C (1993) – Statics and Mechanics of Materials, Prentice Hall College.
HIBBELER, Russell C (1997) – Análisis Estructural, 3era edición, Prentice Hall Collage,
Mexico.
MYERS, Jonh et al. (1969) – Handbook of ocean and underwater engineering, Mc. Graw
Hill, USA.
Naval Facilities Engineering Command (NFEC) (1984) – Harbors. Design Manual 26.1,
Alexandria, Virginia.
Naval Facilities Engineering Command (NFEC) (1985) – Fleet Mooring. Design Manual 26.5,
Alexandria, Virginia.
Naval Facilities Engineering Command (NFEC) (1986) – Fixed Mooring. Design Manual
26.4, Alexandria, Virginia.
Oil Companies International Marine Forum (OCIMF) (1997) – Mooring equipment guidelines
(Second Edition), USA.
Denise Da Costa
– 206 –
Referencias
R.C. Minikin (1963) – Wind, Waves and Maritime Structures. Studies in Harbor Making and
in the Protection of Coasts, 2 ed. Londres.
RODRIGUES, C. et al. (2005) – Torre eólica offshore, Proyecto de estructuras, FEUP, Porto,
Portugal.
ROSARIO DEL, Z. & PEÑALOZA, S (2003) – Guía para la elaboración formal de reportes de
investigación, Universidad Católica Andrés Bello, Caracas, Venezuela.
TSINKER, Gregory (1986) – Floating Ports, Gulf publishing company. Houston, Texas, USA.
TSINKER, Gregory (1995) – Marine structures engineering, Chapman & Hall, USA.
http://www.efunda.com
http://mathworld.wolfram.com
http://www.fishing-catalog.com/docks/fend_off.htm
http://www.schuylerrubber.com/fenders1/sr3d-2_soft.html
http://www.vulcanhammer.net/marine/
Denise Da Costa
– 207 –
Referencias
Denise Da Costa
– 208 –
Anexos
ANEXO A
Denise Da Costa
– 209 –
Anexos
Denise Da Costa
– 210 –
Anexos
4.- Para introducir los datos del problema oprima el botón Entrar.
Denise Da Costa
– 211 –
Anexos
6.- Dentro de “Cabos” y “Defensas”, usted tiene la opción de crear la curva “Fuerza
vs Deformación” de cada elemento. Para lograr esto usted debe:
6.1.- Oprimir el botón “Curva” el cual dará acceso a una hoja de cálculo donde se
podrán introducir los doce puntos característicos de dichas curvas, como se ilustra en la
siguiente figura:
Denise Da Costa
– 212 –
Anexos
6.2.- Introducir los doce puntos característicos de cada curva que desea utilizar.
6.3.- Oprimir el botón “Aceptar” para regresar a la introducción de datos.
6.4.- Para indicar cual curva corresponde a cada cabo o defensa, simplemente
seleccione el número de la curva en la introducción de datos.
7.- Dentro de la pestaña “Cabo + Cadena” podrá introducir los datos necesarios
para definir sistemas que incluyan boya, cadena y ancla.
7.1.- Para el correcto análisis de las cadenas debe oprimir el botón “Catenaria” el
cual dará acceso a una hoja de cálculo donde podrá introducir los datos específicos que
permiten determinar la fuerza máxima de estos sistemas.
7.2.- Después de introducir los datos es preciso oprimir el botó “Analizar” y luego el
botón “Aceptar”
8.- Después de introducir todos los datos completos, oprima el botón “Calcular” para
obtener los resultados.
NOTA: Tenga presente las unidades de los datos de entrada, en cada caso se indican
cuales son las unidades en las que trabaja Amarre
Denise Da Costa
– 213 –
Anexos
4.- Si desea imprimir los resultados, oprima el botón: “Ver versión Imprimible”
Denise Da Costa
– 214 –
Anexos
Denise Da Costa
– 215 –