Source: https://1library.co/document/zwvk2jgq-titulomodelaxe-sistema-fabricacion-maquetas-simulacion-dinamica-software-simulancla.html
Timestamp: 2020-07-13 20:51:55+00:00

Document:
TítuloModelaxe 3D do sistema de fondeo dun barco para a fabricación de maquetas e simulación dinámica co software SIMULANCLA
#Diego #Saya Carro
(1)AGRADECIMIENTOS Este trabajo realizado en la Escuela Politécnica Superior de Ferrol me ha permitido aprovechar una experiencia de aprendizaje muy gratificante, pudiendo conocer de cerca un importante sector en la comarca. La motivación de este proyecto presenta una importante oportunidad para la modernización en ciertas áreas de la industria naval y abriendo un apasionante camino en el desarrollo de estas herramientas. En primer lugar, quiero dar mi mayor agradecimiento a mi tutor Daniel Dopico Dopico por haber trabajado durante estos meses conmigo. Su disposición y su ayuda han sido una gran motivación para llevar a cabo este trabajo. Ha sido un placer poder trabajar durante este tiempo a su lado, por ello quiero agradecer de forma especial el apoyo que he recibido y el constante ánimo para mejorar. También quiero agradecer a Ana Álvarez García su importante. ayuda. Sus. explicaciones han sido vitales para poder llevar a cabo la fase de diseño y conocer en detalle los procedimientos y herramientas empleadas en esta disciplina. De manera muy especial, agradezco a mi familia todo el apoyo durante estos años. Especialmente mis padres, quienes me han ayudado en todo momento a seguir y esquivar las dificultades que se han ido presentando. Todo lo que pueda llegar a lograr es gracias a su esfuerzo. No quiero olvidarme de Alejandra, siendo su apoyo una importante fuente de ánimo. A todos los que me habéis ayudado, muchas gracias..
(3) Índice.
(7) Índice 1. Introducción ______________________________________ 1-1 1.1.. Antecedentes _____________________________________________________ 1-2. 1.2.. Motivación _______________________________________________________ 1-3. 2. Desarrollo del plano de formas del buque y sistema de fondeo _ 2-1 2.1.. Plano de formas del buque __________________________________________ 2-2. 2.2.. Sistema de fondeo _________________________________________________ 2-3. 2.2.1.. 2.3.. Parámetros elegidos para el sistema de fondeo ________________________________ 2-5. Modelo de CAD del buque __________________________________________ 2-7. 3. La dinámica de sistemas multicuerpo ____________________ 3-1 3.1.. Modelización en coordenadas mixtas _________________________________ 3-2. 3.1.1.. Las coordenadas naturales __________________________________________________ 3-2. 3.1.2.. Las coordenadas r elativas ___________________________________________________ 3-2. 3.1.3.. Las coordenadas mixtas ____________________________________________________ 3-3. 3.2.. El problema cinemático: problema de posición y velocidad inicial _________ 3-4. 3.2.1.. El problema de posición inicial _______________________________________________ 3-4. 3.2.2.. El problema de velocidad inicial _____________________________________________ 3-5. 3.3.. El problema dinámico: ecuaciones de la dinámica ______________________ 3-6. 3.3.1.. Introducción a la resolución de problemas dinámicos __________________________ 3-6. 3.3.2.. El primer problema dinámico: el problema de aceleración inicial ________________ 3-7. 3.3.3.. El problema dinámico en tiempo r eal: formulación A LI3-P ______________________ 3-9. 3.4.. Modelización elegida para el sistema de fondeo _______________________ 3-12. 3.4.1.. Definición del buque ______________________________________________________ 3-14. 3.4.2.. Definición de sólidos del sistema de fondeo __________________________________ 3-14. 3.4.3.. Grados de libertad ________________________________________________________ 3-18. 3.4.4.. El vector de variables ______________________________________________________ 3-19. 3.4.5.. Fuerzas __________________________________________________________________ 3-19. 4. SIMULANCLA: un software para la simulación del sistema de fondeo de buques ______________________________________ 4-1. V.
(8) 5. Simulación de la maniobra de anclas del buque diseñado ____ 5-1 5.1.. Configuración de la simulación _______________________________________5-2. 5.2.. Resultados de la simulación _________________________________________5-7. 6. Fabricación y simulación del modelo a escala_____________ 6-1 7. Conclusiones _____________________________________ 7-1 8. Bibliografía ______________________________________ 8-1 9. ANEJOS_________________________________________ 9-1 9.1.. VI. Anejo I: Plano de formas del buque ___________________________________9-1.
(9) Índice de figuras Figura 2-1. Referencias de un buque ...............................................................................................................2-2 Figura 2-2. Planos de un buque ........................................................................................................................2-2 Figura 2-3. Planos de formas.............................................................................................................................2-7 Figura 2-4. Planos de formas desplazados .....................................................................................................2-8 Figura 2-5. Modelado del casco ........................................................................................................................2-8 Figura 2-6. Modelado de proa ..........................................................................................................................2-9 Figura 2-7. Trazado de la brusca ................................................................................................................... 2-11 Figura 2-8. Divisiones para cálculo de b rusca ............................................................................................. 2-12 Figura 2-9. Modelado de la cubierta ............................................................................................................. 2-13 Figura 2-10. Modelado del escobén .............................................................................................................. 2-14 Figura 2-11. Grillete, caña y ancla ................................................................................................................. 2-14 Figura 3-1. Restricciones para los eslabones ............................................................................................... 3-14 Figura 3-2. Restricciones para el g rillete y la caña ..................................................................................... 3-16 Figura 3-3. Restricciones para la caña y las uñas del ancla ...................................................................... 3-17 Figura 3-4. Contacto no rmal .......................................................................................................................... 3-20 Figura 3-5. Contacto tangencial .................................................................................................................... 3-21 Figura 3-6. Ángulo del ancla con la caña ..................................................................................................... 3-24 Figura 4-1. Configuración de la simulación ....................................................................................................4-3 Figura 4-2. Configuración del buque................................................................................................................4-5 Figura 4-3. Cálculo de las características del elemento ................................................................................4-6 Figura 4-4. Configuración de la caña ...............................................................................................................4-7 Figura 4-5. Configuración del ancla .................................................................................................................4-7 Figura 4-6. Configuración de la cadena ...........................................................................................................4-9 Figura 4-7. Configuración del grillete y el molinete ......................................................................................4-9 Figura 4-8. Fuerzas representadas en la simulación .................................................................................. 4-10 Figura 4-9. Estructura de malla ..................................................................................................................... 4-11 Figura 5-1. Parámetro s de simulación .............................................................................................................5-2 Figura 5-2. Datos del buque ..............................................................................................................................5-3 Figura 5-3. Datos caña .......................................................................................................................................5-4 Figura 5-4. Datos ancla ......................................................................................................................................5-5 Figura 5-5. Datos de la cadena .........................................................................................................................5-6 Figura 5-6. Datos del grillete y el molinete .....................................................................................................5-7 Figura 5-7. Posición inicial .................................................................................................................................5-8 Figura 5-8. Posición de equilib rio .....................................................................................................................5-8. VII.
(10) Figura 5-9. Cuadro de control de la simulación ............................................................................................. 5-9 Figura 5-10. Fuerzas en la simula ción ............................................................................................................. 5-9 Figura 5-11. Evolución de la simulación ....................................................................................................... 5-10 Figura 5-12. Posición de estiba....................................................................................................................... 5-11 Figura 5-13. Posición ancla ............................................................................................................................. 5-12 Figura 5-14. Velo cidad de tiro ........................................................................................................................ 5-13 Figura 5-15. Recogida cadena ........................................................................................................................ 5-13 Figura 5-16. Fuerza del molinete ................................................................................................................... 5-14 Figura 6-1. Modelo a escala 1:30..................................................................................................................... 6-2. VIII.
(11) Índice de tablas. Tabla 2-1. Características del buque................................................................................................................2-1 Tabla 2-2. Tabla Numeral de Equipo................................................................................................................2-6 Tabla 3-1 Balance de grados de lib ertad. .................................................................................................... 3-18. IX.
(13) Capítulo 1 Introducción.
(15) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Introducción. 1. Introducción. El objetivo del presente proyecto es la validación del software SIMULANCLA, una novedosa tecnología que evita la necesidad de modelar el sistema de fondeo y todos los equipos que intervienen en esta maniobra, que permite el diseño de cascos de buques de grandes dimensiones. Las maniobras de fondeo en buques de gran tamaño entrañan muchos riesgos, tanto en el propio proceso de fondeo sobre el fondo marino como en la recogida de las cadenas y las anclas, que obligatoriamente deben quedar perfectamente estibadas en el casco sin que se produzcan situaciones de inestabilidad o atascos de las mismas antes de llegar a una estiba satisfactoria. El diseño de la zona del casco del buque sobre la que debe acomodarse el ancla en el proceso de estiba no es en absoluto trivial. Existen dos sistemas distintos en los buques objeto de esta investigación, el nicho y la regola, con formas geométricas complejas y especialmente diseñadas para este propósito. Hasta hace unos años, el estudio del sistema de fondeo de un buque se realizaba mediante el uso de maquetas, siendo costosas en tiempo y materiales. Además, limita n el estudio por los recursos necesarios y solamente ofrecen una aproximación general del comportamiento real. Con el uso de este software, se busca simplificar este anális is, permitiendo un estudio más cercano a la realidad, fiable y con unos tiempos de ejecución más rápidos. El software se ha desarrollado durante los últimos seis años y facilita a los ingenieros del astillero el diseño geométrico de la regola o nicho del buque, de tal forma que la maniobra para levar el ancla se lleve a cabo de manera satisfactoria, así como la elección de los eslabones que deben componer la cadena y todos los elementos mecánicos implicados en las maniobras.. 1-1.
(16) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Introducción. 1.1.. Antecedentes. Hasta hace unos años, se fabricaba, para cada diseño nuevo de buque que había que desarrollar, una maqueta de madera que representaba la maniobra de fondeo a escala 1:10. Estas maquetas son costosas en tiempo y no proporcionan la fiabilidad requerida, puesto que ni la escala ni los materiales en que están hechas son los reales, de modo que permiten sólo una evaluación cualitativa y aproximada de la bondad de la solución de diseño propuesta. La fabricación de una de estas maquetas supone aproximadamente 1500 horas de trabajo y la implicación de diversos departamentos. La nueva tecnología para el diseño de todos los elementos que intervienen en el proceso de fondeo supone un cambio de paradigma con respecto a la metodología empleada hasta la fecha en los astilleros, basada en las citadas maquetas de madera a escala del casco, cadena, ancla y de todos los sistemas que intervienen en el proceso de fondeo y leva de anclas. Hay diversos factores que se deben tomar en consideración para el diseño, como las distintas posiciones en las que se puede encontrar el ancla en el momento del izado y el posible balanceo y/o inclinación del buque durante el mismo. Asimismo, hay que asegurarse de que no se bloquea el mecanismo cuando el contacto entre el ancla y el casco se produce en una posición desfavorable. El diseño debe asegurar que no se produzca el bloqueo del ancla durante la maniobra y que ésta adopte la posición esperada en el casco del buque cuando esté completamente recogida. El software desarrollado permite simular, con gran detalle, la dinámica de las maniobras de leva de anclas, basándose en las técnicas de simulación de sistemas multicuerpo.. 1-2.
(17) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Introducción. 1.2.. Motivación. Con el fin de validar la maniobra de fondeo de un buque, el software SIMULANCLA va a permitir realizar una simulación completa y detallada sin la necesidad de disponer de una maqueta para comprobar la maniobra en cuestión. Gracias al software, será posible definir la maniobra con mayor similitud a la situación real. El empleo de una maqueta sólo permite una aproximación de la maniobra real, pues el empleo de otros materiales y una escala distintos a los reales, aportan una evaluació n cualitativa que solamente nos puede acercar a una idea del comportamiento real de los elementos que intervienen en el proceso. Sin embargo, el software SIMULANCLA permite realizar un estudio más completo, pudiendo analizar la dinámica de la maniobra, así como visualizar la maniobra como si de una maqueta se tratase, con la ventaja de poder conocer diferentes datos en tiempo de simulación. De este modo, será posible conocer con una mayor exactitud el sistema de fondeo de un buque, pudiendo comparar distintos modelos y comprobar que configuraciones se adaptan mejor a los requerimientos del equipo técnico. Todo esto sería muy difícil de materializar con el empleo de maquetas, pues requeriría de una inversión económica y horaria muy extensa. El desarrollo de este software además ha reducido el plazo necesario desde que se dispone de la estructura de la zona de proa hasta obtener las formas definitivas de la regola, nicho,. piezas fundidas. necesarias,. etc., así como la realización. de. modificaciones de manera fácil y rápida que permiten optimizar las maniobras.. 1-3.
(19) Capítulo 2 Desarrollo del plano de formas del buque y sistema de fondeo.
(21) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. 2. Desarrollo del plano de formas del buque y sistema de fondeo. Con el fin de realizar la simulación de la maniobra con el software SIMULANCLA, primero se desarrollará el modelo del buque, así como todos los componentes del sistema de fondeo necesarios para la simulación. Se propone para esta maniobra, un buque de apoyo a plataformas eólicas offshore. Las características principales del buque se pueden observar a continuación:. Tabla 2-1. Características del buque. El punto de partida para realizar el modelado, será el plano de formas que define la geometría del buque propuesto. El plano de formas es un plano fundamental que nos muestra las formas definitivas del buque. Se puede visualizar en el Anejo I: Plano de formas del buque. A partir de este plano de formas se generan las curvas que determinan el buque. Para ello habrá que ubicar las distintas curvas en el espacio, de tal manera que se ajusten a lo representado en el plano. Estas curvas serán la base para generar la superficie del buque. Una vez generado el casco se obtendrá el modelo CAD del buque, que posteriorme nte habrá que afinar para obtener el modelo final que será empleado en la simulación.. 2-1.
(22) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. Por último, se realizará el modelado de los distintos elementos que van a participar en la simulación.. 2.1.. Plano de formas del buque. Generalmente los buques se definen por las dimensiones del prisma regular en el que está contenido. La longitud del mismo se denominará eslora, el ancho será la manga y el alto se conoce como puntal. Al tratarse de un cuerpo móvil, considerando un observador en el centro del buque, se denomina proa al extremo situado ante él, popa al extremo tras él, babor al lado que está a su izquierda y estribor al lado que está a su derecha.. Figura 2-1. Referencias de un buque. Para definir la geometría del buque se definen el plano base, que representa la vista en planta; el plano transversal en la sección media, quedando a la izquierda las líneas de popa y a la derecha las líneas de proa; y el plano longitudinal por la línea de crujía, es decir, la línea que pasando por la popa y la proa divide el barco en dos bandas iguales.. Figura 2-2. Planos de un buque. 2-2.
(23) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. A partir de estos planos se determinan las siguientes secciones: . Las paralelas al plano de base se denominan líneas de agua, siendo la línea cero la línea de base o línea de agua.. . Las paralelas a los planos de los extremos de proa y popa, perpendiculares al plano base, se denominan cuadernas, empezando el cero por la perpendicular de popa. Desde el cero hasta popa se enumera con números negativos.. . Las paralelas al plano longitudinal y perpendiculares a los otros planos, se denominan secciones longitudinales, siendo uno la más próxima al plano longitudinal que pasa por crujía. El plano diametral, o sección media, será el que divide el buque en dos bandas iguales.. De este modo, el plano de formas está compuesto por las líneas de intersección de los planos de las secciones, tanto horizontales, transversales y longitudinales, descritas en el párrafo anterior; con la superficie exterior del buque, representadas en los tres planos de proyección. El procedimiento a llevar a cabo consiste en la translación de las formas del buque. Para ello se tomarán las curvas del plano de formas, se girarán y se colocarán en un punto de origen común. Este origen se encontrará en la intersección de la perpendicular de popa y de la línea de base. Posteriormente se acomodará cada una de estas curvas en su sitio, correspondiendo a la distancia obtenida al dividir la eslora con perpendiculares en 20 partes iguales, para construir una matriz de curvas en tres dimensiones. Con esta matriz se va a generar la superficie del buque. Una vez obtenida la superficie, se procederá a refinar la misma, en términos navales se denomina alisar las formas, para lograr que esta sea lo más aerodinámica posible.. 2.2.. Sistema de fondeo. El sistema de fondeo de un buque es el encargado de inmovilizar el mismo y de separarlo de la acción de las corrientes y el viento mediante elementos capaces de fijarse en el fondo del agua.. 2-3.
(24) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. Esta serie de elementos forman parte del buque y van unidos a él. Aunque puede incluir más elementos, las partes fundamentales de un sistema de fondeo son las siguientes: . Anclas. . Cadenas. . Molinete. . Escobén. . Estopor. . Tubos de las cadenas. . Caja de cadenas. El elemento principal es el ancla, encargado de fondear el buque por la acción de su propio peso y por sus uñas. Funciona de manera que, si es sometida a una fuerza paralela al fondo del agua, esta tiende a clavarse en el mismo. Para ello será necesario que el buque conserve una pequeña velocidad. Están fabricadas en hierro o acero y se unen al buque mediante un cabo o cadena. En cuanto a las cadenas, se ha ido generalizando el uso de cadenas de acero. En estas se pueden distinguir los siguientes elementos: grillete de ancla, eslabón final, eslabón grande, eslabón giratorio, eslabón ordinario y eslabón de unión. Haciendo referencia a los tipos de cadenas, pueden ser con contrete o sin contrete. El contrete es la pieza que se coloca en el interior del eslabón en sentido del eje menor para evitar su deformació n cuando está trabajando a tracción. Además, los eslabones deben ser normalizados y estandarizados para una mayor facilidad de sustitución cuando sea necesario. El molinete es el elemento que se encarga de elevar el ancla del buque. Éste gira sobre un eje horizontal provocando el paso de los eslabones hacia la caja de cadenas. El escobén es un conducto de sección circular o elíptica que se abre entre la cubierta y la armadura de un buque para permitir el paso de las cadenas y alojar la caña del ancla en la posición de estiba. El estopor es un elemento de acero colocado en la cubierta entre el molinete y el escobén. Este retiene la cadena evitando que la tensión actúe directamente sobre el 2-4.
(25) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. molinete. Todos los estopores tienen dos partes principales: una parte fija con una ranura que permite el paso del eslabón en posición vertical y una parte móvil que trata de retener al eslabón que se encuentra en la ranura. La caja de cadenas será el lugar del buque donde se estiban las cadenas de las anclas. Deben colocarse debajo del molinete y en las proximidades de su vertical. Las cadenas de cada ancla deben estibarse por separado, siendo necesaria la existencia de un mamparo que evite que se puedan mezclar. El fondo de la caja de cadenas debe hacerse de tal modo que permita el achique del agua y fango que se pueda ir depositando. 2.2.1. Parámetros elegidos para el sistema de fondeo Los medios de fondeo se calculan partiendo del numeral de equipo. El numeral de equipo es un número adimensional que sirve para determinar el equipo de anclas y cadenas en función de la reglamentación de las sociedades de clasificación. El número y el tamaño de las anclas apropiadas para un buque van a estar reguladas por el reglamento de dichas sociedades. Para el caso en cuestión, se hará uso del numeral de equipo de la sociedad Bureau Veritas. Se define como numeral de equipo al número ‘N’ resultante de la siguiente expresión: 2. NE  K  (length  breath  depth) 3 ,. (2.1). siendo length la eslora del buque, breath la manga y depth el calado. El coeficiente K podrá tener distinto valor. Para el caso de remolcadores sin restricción se toma K=1,3 y para el resto de buques se toma K=2,2. Consultando los datos de la Tabla 2-1, se obtiene el siguiente resultado para la ecuación (2.1): NE  2, 2  (176,1 29,8  9,9). 2. 3.  3063. (2.2). 2-5.
(26) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. Una vez obtenido el numeral de equipo, se obtendrán los distintos datos del sistema de fondeo consultando las tablas de la sociedad de clasificación correspondiente en función del numeral de equipo (NE) obtenido.. Tabla 2-2. Tabla Numeral de Equipo. De este modo se podrá conocer el número de anclas, su peso, el diámetro del eslabón de cadena, la longitud total de cadena, la carga de rotura, la longitud del cable de remolque y el número de amarras con su longitud y carga de rotura. Consultando los valores de la Tabla 2-2 se puede obtener la masa para el anclaje, fundamental para la simulación, que es de 9300 kg. De este modo, se hará uso de un ancla comercial que se ajuste a estas características. Cada línea horizontal de la tabla muestra el equipo de fondeo correspondiente a un buque cuyo numeral de equipo se encuentre entre los valores indicados. Una vez identificada la fila en la que se encuentra el numeral obtenido, los valores restantes de dicha fila muestran los valores del equipo de fondeo que le corresponde al buque en cuestión.. 2-6.
(27) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. 2.3.. Modelo de CAD del buque. Una vez esté completado el plano de formas del buque, este puede ser empleado para generar la superficie que da forma al casco del mismo. A partir de este sistema de curvas proyectadas en distintos planos, se podrá generar la superficie del casco empleando un programa de modelado 3D. Se hará uso del software Rhinoceros. El primer paso será importar los planos de formas y orientarlos de forma correcta en base a un punto común. Es importante que el plano transversal quede correctamente centrado respecto al plano longitudinal. De este modo, las líneas de proa y popa quedarán extendidas a un lado u otro del plano longitudinal para facilitar la obtención de la superficie. El plano base se representa en una mitad, pues este es simétrico. Este plano se importará a ambos lados del plano longitudinal para emplear como ayuda en la construcción del casco, aunque no será necesaria la utilización de sus curvas para generar la superficie. Sus curvas servirán de referencia para desplegar correctamente las curvas del plano base. En la Figura 2-3 se puede apreciar la correcta ubicación de los planos de formas en el espacio 3D:. Figura 2-3. Planos de formas. Una vez se hayan definido los planos correctamente, se tendrán que trasladar las curvas de proa y de popa a su posición real respecto al plano longitudinal, pues ahora están todas proyectadas en el plano transversal. Al finalizar este proceso, las curvas de proa 2-7.
(28) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. quedarán a un lado del plano longitudinal y las de popa al otro lado. Se tomará un conjunto de estas curvas, las de proa o las de popa, y se trasladarán al lado opuesto para que todas las curvas se encuentren alineadas.. Figura 2-4. Planos de formas desplazados. Tras haber ubicado cada curva en su lugar correspondiente, se podrá comenzar a dar forma a la superficie del casco. Se seleccionarán todas las curvas y, a través del comando ‘Transición’ en las opciones de superficie, se generará la superficie resultante del casco. Una vez obtenida dicha superficie, se realizará una simetría respecto al plano longitudinal para generar la otra mitad del casco.. Figura 2-5. Modelado del casco. El objetivo de la simulación se centra en la operación del sistema de fondeo, siendo la proa la parte de interés del casco para este estudio. Por ello, se han creado distintas. 2-8.
(29) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. capas que incluyen diferentes zonas del casco: la proa, la zona central y la popa. Con el fin de optimizar la simulación, únicamente se empleará la zona de proa para realizar la simulación, además de facilitar la visualización de la maniobra al focalizarla en esta zona. La superficie del casco de proa que será empleada en la simulación se puede apreciar en la Figura 2-6:. Figura 2-6. Modelado de proa. Con la superficie del casco definida, el siguiente paso será refinar la malla para obtener un diseño más limpio y, por otra parte, modelar la cubierta y el resto de elementos que intervienen en la simulación. Estos trabajos serán realizados con el software Solid Edge, que permitirá añadir estos elementos al casco creado con gran precisión y haciendo posible lograr las complicadas formas que toman algunos cuerpos. Es necesario exportar el modelo de Rhinoceros a un formato compatible con el nuevo software a emplear. Posteriormente se importará en Solid Edge para comenzar a refinar la superficie lograda. Será necesario recalcular la malla para que ésta presente una mayor uniformidad en su conjunto. De este modo, se evitarán posibles imperfeccio nes en aquellas zonas en las que el casco adopta formas más estilizadas. Es de vital importancia que la superficie final tenga un diseño limpio y óptimo para garantizar una buena interacción en la simulación de los distintos elementos que componen el sistema de fondeo. La existencia de imperfecciones en los cuerpos puede alterar los contactos entre estos, creando imprecisiones en los cálculos que no son deseables.. 2-9.
(30) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. En el diseño de la cubierta, es muy importante conocer los requerimientos a los que debe adaptarse para garantizar su correcta integración en el casco del buque. Durante la navegación, el buque está sometido a continuos movimientos ocasionados por el oleaje. Estos movimientos se deben a los seis grados de libertad, tres de traslación y tres de rotación: . Traslaciones o Traslación vertical. Este movimiento se debe a la flotación del barco, pudiendo ser de ascenso o descenso. o Traslación lateral. Este movimiento, conocido como ronza, se debe al desplazamiento lateral del barco de una banda a otra. o Traslación longitudinal. Este movimiento se debe al avance o retroceso a lo largo del plano longitudinal.. . Rotaciones o En eje vertical.. Esta rotación, conocida como virada, se debe. principalmente a la acción del timón, aunque puede ser ocasionada por otros fenómenos como el viento. o En eje transversal. Esta rotación, conocida como cabeceo cuando desciende la proa o arfada cuando esta asciende, se debe a la rotación en el eje transversal del buque, ocasionando desplazamientos verticales en proa o en popa. o En eje longitudinal. Esta rotación, conocida como escorado, se debe a la rotación del casco de una banda a otra. Durante la navegación el buque experimenta una combinación de estos movimientos, y el agua llega a alcanzar la cubierta en numerosas ocasiones, sobre todo en condiciones muy adversas. Si esta agua que se deposita en la cubierta no se desaloja, supone un riesgo muy importante, pues podría provocar el hundimiento del buque. Para facilitar la expulsión del agua que pueda alojarse en la cubierta, se dota a la misma de una curvatura que dirija el agua hacia los laterales del barco para su expulsión. Esta curvatura de la cubierta se denomina brusca y será necesario calcular su trazado en función de las medidas del buque. 2-10.
(31) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. Figura 2-7. Trazado de la brusca. Como se puede apreciar en la Figura 2-7, ‘b’ representa la altura en el centro de la cubierta y ‘l’ la mitad de la manga, pues la cubierta será simétrica respecto al plano longitudinal. El valor de ‘b’ estará situado entre 1/48 y 1/50 de la manga máxima, en función de la altura que se pueda emplear por cuestiones de diseño; en este caso es suficiente emplear un coeficiente de 1/50. El valor de ‘l’ será igual a la mitad de la manga máxima. De este modo, teniendo en cuenta los datos de la Tabla 2-1, los valores que se obtienen son los siguientes:. b  manga . l. 1  0,596m 50. 29,8  14,9m 2. (2.3). (2.4). El siguiente paso es establecer un número de divisiones para calcular la altura de la cubierta en dichas divisiones, para así generar la curva que dará forma a la superfic ie de la cubierta. Se toman al menos cuatro segmentos para obtener una precisión suficiente. Con este número, se crean las divisiones correspondientes en ‘l’ y en ‘b’. Las divisiones de ‘b’ se trasladan a la vertical en ‘0’, y cada una de ellas se une con la altura máxima en ‘4’, que representa la altura máxima de la cubierta. Cada línea generada va a determinar un punto de corte en cada división, dando la altura de la cubierta en dicho punto. Las divisiones resultantes se pueden apreciar en la Figura 2-8:. 2-11.
(32) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. Figura 2-8. Divisiones para cálculo de brusca. Una vez conocidos estos puntos se podrá calcular su valor, que representa la altura de la cubierta en cada uno. La expresión para obtener dicha altura viene dada por:. y1 . y2 . b 2  x1 l2. b 2  x2 l2. (2.5). (2.6). y3 . b 2  x3 l2. (2.7). y4 . b 2  x4 l2. (2.8). Para cada división, la altura correspondiente es la siguiente:. 2-12. x1  3725mm y1  37,25mm. (2.9). x2  7450mm y2  149mm. (2.10). x3  11175mm y3  335,25mm. (2.11). x4  14900mm y4  596mm. (2.12).
(33) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. Con estos datos se puede generar una curva en el software CAD para generar la superficie de la cubierta. Primero se realizará el boceto que represente la brusca con los datos calculados, para posteriormente generar la superficie a partir de dicho boceto.. Figura 2-9. Modelado de la cubierta. El siguiente elemento a definir es el escobén. Este elemento es de vital importancia en el sistema de fondeo. Forma la conexión entre la cubierta y el casco, atravesando el mismo para que la cadena pueda descender el ancla hasta el fondo marino. El primer parámetro fundamental para su diseño es el diámetro, pues debe permitir el paso de la cadena empleada y evitar que la caña del ancla pueda quedar atrapada cuando se procede al estibado. Además, debe permitir acomodar la caña cuando el ancla está completamente estibada. También requiere especial atención la inclinación que presentará, pues es fundamental para los contactos que los bordes sean suaves para facilitar el paso de los elementos. Unas uniones entre estos elementos muy remarcadas podrían propiciar rozamientos elevados que pueden dañar seriamente los equipos del sistema de fondeo. Por último, será necesario diseñar una prolongación del escobén. 2-13.
(34) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. que permita enfocar la cadena hacia el molinete, facilitando el traslado de la cadena a la propia caja de cadenas.. Figura 2-10. Modelado del escobén. Los modelos del grillete, la caña y las uñas del ancla han sido facilitados por Navantia, correspondiéndose a elementos de un buque real. Estos elementos cumplen con las características que se han ido definiendo anteriormente.. Figura 2-11. Grillete, caña y ancla. 2-14.
(35) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior Desarrollo del plano de formas del buque y sistema de fondeo. 2-15.
(37) Capítulo 3 La dinámica de sistemas multicuerpo.
(39) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. 3. La dinámica de sistemas multicuerpo. Un sistema multicuerpo o multibody system es la modelización física del conjunto de sólidos de un mecanismo y las relaciones y restricciones que existen entre ellos. A través de los sistemas multicuerpo podemos estudiar la cinemática, la dinámica inversa o la dinámica directa del mecanismo. El tipo de problema de interés para este proyecto es el de dinámica directa, que permite conocer la posición, velocidad y aceleración del sistema conocidas las fuerzas que actúan sobre los sólidos del mismo. El primer paso es la modelización del sistema multicuerpo. Para la aplicación de la modelización multicuerpo es precisa la parametrización de las posiciones de los sólidos mediante unas coordenadas apropiadas y la definición de las ecuaciones de restricción que ligan estas coordenadas. El siguiente paso es la resolución de la cinemática del sistema en la posición inic ia l, obteniendo unas posiciones y velocidades iniciales que cumplan las ecuaciones de restricción. A continuación se lleva a cabo un tercer paso es opcional: el planteamie nto de las ecuaciones de la estática del sistema y la resolución de la posición de equilib r io inicial del mismo. Finalmente, el último paso es el planteamiento y resolución de las ecuaciones del movimiento o ecuaciones dinámicas del sistema en cada paso de tiempo. Las entradas al problema dinámico son las fuerzas actuantes sobre el sistema y dichas fuerzas son tremendamente complejas de obtener, resultado de la evaluació n de complejos modelos de fuerzas, siendo su obtención un formidable problema en sí mismo. Es fundamental conocer la disciplina de dinámica de sistemas multicuerpo para comprender la base de cálculo del software SIMULANCLA. A partir de la solución multicuerpo, será posible la obtención de todos aquellos resultados que nos permitan analizar la dinámica de la maniobra de fondeo.. 3-1.
(40) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. 3.1.. Modelización en coordenadas mixtas. Las coordenadas son variables que definen completamente la posición de un mecanismo. Generalmente se trata de puntos, vectores, ángulos y distancias. Existen diferentes tipos de definición mediante coordenadas, y cada una de ellas va a dar lugar a distintos planteamientos del problema. En nuestro caso, se elegirá un sistema mixto de coordenadas, esto es, elegiremos simultáneamente dos tipos de definición para optimizar la solución del problema y aprovechar las ventajas de cada uno de ellos, minimizando así la posibilidad de que aparezcan problemas. Estos dos tipos de coordenadas serán las naturales y las relativas. Las primeras se emplearán para modelizar sólidos, y las segundas para actuadores (motores), definición de posiciones grados de libertad e introducción de ciertos tipos de fuerzas. 3.1.1. Las coordenadas naturales Las coordenadas naturales son básicamente las coordenadas cartesianas de puntos y vectores unitarios unidos a los sólidos del mecanismo. Esas coordenadas sitúan cada sólido del mecanismo Para cada sólido 3D se definen un mínimo de 12 coordenadas en total, así pues, dado que el número de grados de libertad de un sólido libre es 6, estas coordenadas son dependientes por propia definición. Aportan la ventaja de una definición sencilla y sistemática, con ecuaciones de restricción igualmente simples, ya que son lineales o cuadráticas, pero tienen el inconveniente de que, para obtener una modelización correcta, precisan de cierta experiencia. 3.1.2. Las coordenadas relativas Las coordenadas relativas se diferencian de las naturales en que, en lugar de referir los sólidos a un origen común a todos, se refieren cada uno de ellos a su antecesor en la cadena cinemática.. 3-2.
(41) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. Cuando se modelizan sistemas de cadena abierta exclusivamente en coordenadas relativas, se obtiene una cinemática sin restricciones ya que el número de coordenadas relativas necesarias es igual al número de grados de libertad del sistema. Cuando se modelizan sistemas de cadena cerrada exclusivamente con coordenadas relativas, éstas son redundantes, pero se obtiene un reducido número de ecuaciones de restricción; las necesarias para imponer las condiciones de cierre de lazo. Cuando se emplean las coordenadas relativas conjuntamente con otros sistemas de coordenadas, éstas son especialmente apropiadas para la definición de grados de libertad, actuadores o fuerzas asociadas a pares cinemáticos. El primer inconvenie nte de estas coordenadas es la necesidad de situar todos los elementos anteriores de la cadena para poder conocer la posición de un sólido, por las razones anteriorme nte descritas. El segundo inconveniente, es que las ecuaciones contienen términos trigonométr icos bastante complejos, lo cual aumenta considerablemente el coste computacional de la resolución del problema. 3.1.3. Las coordenadas mixtas Las coordenadas mixtas serán las empleadas en este proyecto. Utilizar dos sistemas de coordenadas diferentes nos permite aprovechar las ventajas de ambos. Por un lado, contaremos con la sencillez de las coordenadas naturales para la definición de las posiciones de los sólidos de forma independiente al resto de sólidos de la cadena cinemática, y por otro, tendremos una facilidad de definición para los grados de libertad, actuadores y fuerzas gracias a las coordenadas relativas. Tendremos que trabajar, entonces, con los puntos y vectores de las coordenadas naturales y con los ángulos y distancias de las coordenadas relativas.. 3-3.
(42) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. 3.2.. El problema cinemático: problema de posición y velocidad inicial. El problema cinemático consiste en la determinación de los posibles movimientos de un mecanismo concreto, independientemente de las fuerzas actuantes sobre él. Su resolución nos permite conocer la posición y velocidad iniciales del mecanismo, en función de las posiciones y velocidades de cada uno de los distintos grados de libertad. 3.2.1. El problema de posición inicial Antes de simular los movimientos posibles del mecanismo debemos conocer una posición que satisfaga todas las restricciones impuestas al mismo. Así, solucionar e mos este problema de posición inicial para evitar arrancar la simulación desde una posición que no cumpla las ecuaciones de restricción. Sea un sistema multicuerpo, cuya posición se define mediante el siguiente vector de coordenadas generalizadas:. q   q1 , q2 ,..., qm , z1, z2 ,..., zn . (3.1). Donde: qi: coordenadas dependientes zi: coordenadas independientes (variables que definen los GDL) Las coordenadas dependientes están relacionadas por un vector de m restriccio nes Φ(q) que consiste en una serie de ecuaciones no lineales. Este vector será desarrollado. en serie de Taylor entorno a la posición inicial. Una vez desechados los términos de orden superior a uno, tenemos: Φ(q)  0. (3.2). Φ(q)  Φ(q0 )  Φq (q0 )(q  q0 )  0. (3.3). 3-4.
(43) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. Para resolver esta última ecuación aplicaremos el método iterativo de NewtonRaphson: Φq (qi )(qi 1  qi )  Φ(qi ). (3.4). Φ q es la matriz jacobiana de las ecuaciones de restricción respecto a las coordenadas. generalizadas. Nótese que los valores de q z1 , z2 , zn correspondientes a los GDL son i 1 i conocidos de antemano z j  z j  0 , por tanto, las columnas correspondientes de la. matriz Jacobiana siempre van a ser multiplicadas por elementos nulos y no se calculan. Para arrancar las iteraciones del método de Newton-Raphson, se parte de un q0 estimado. Como criterio de convergencia se elige la norma-2 del vector de restricciones. Cuando en la iteración iésima sea menor que una cierta tolerancia emax se considera suficientemente aproximado el resultado. 1.  m 2 ei  Φ(qi )     2j   j 1 . (3.5) i. 3.2.2. El problema de velocidad inicial Lo que se pretende conseguir mediante la resolución de este problema son las velocidades de las coordenadas generalizadas para el instante inicial, que son las derivadas temporales primeras de dichas coordenadas. Derivando la ecuación (3.2) con respecto al tiempo: Φq (q)q  Φt  0. (3.6). Φq (q)q  Φt. (3.7). Del problema de posición inicial ya tenemos el jacobiano y las velocidades de los grados de libertad son conocidas de antemano, así que podemos resolver directamente la ecuación (3.6) para conocer las velocidades. 3-5.
(44) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. Este problema es lineal, por lo que no precisamos de un método iterativo.. 3.3.. El problema dinámico: ecuaciones de la dinámica. 3.3.1. Introducción a la resolución de problemas dinámicos Los dos conjuntos de ecuaciones más populares que nos permiten resolver el problema dinámico son las ecuaciones de Lagrange y las de Newton – Euler. Mientras que en estas últimas el número de ecuaciones es seis veces el número de sólidos, en las ecuaciones de Lagrange este número es proporcional al número de coordenadas generalizadas y generalmente muy inferior, por lo que se utilizarán éstas. Utilizando las coordenadas mixtas explicadas en el apartado 3.1.3 las ecuaciones de Lagrange son:. d  T  T  ΦqT λ  Q   dt  q  q. (3.8). ΦqT λ representa los esfuerzos requeridos para mantener las m restricciones entre las distintas variables dependientes y T es la energía cinética (función de la masa y la velocidad). Con las coordenadas elegidas se puede escribir de la forma siguiente :. 1 T  qTMq 2. (3.9). Sustituyendo la ecuación (3.9) en la ecuación (3.8)se obtiene:. Mq  ΦqT λ  Q Donde: M: matriz de masas q: vector de variables del problema Φ q: matriz jacobiana del vector de restricciones. 3-6. (3.10).
(45) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. λ: vector de multiplicadores de Lagrange Q: vector de fuerzas generalizadas El sistema (3.10) consiste en n ecuaciones, siendo n la cantidad de variables, sin embargo, existen m multiplicadores de Lagrange, tantos como ecuaciones de restricción. Por ello debemos añadir m ecuaciones de restricción. Como resultado, obtenemos un sistema compuesto por ecuaciones diferenciales y algebraicas o sistema DAE (Differential - Algebraic Equations) a ) Mq  ΦqT λ  Q  b)  Φ0. (3.11). Para resolver este sistema necesitaremos una formulación capaz de trabajar con él. Ya que la mayoría de los integradores están diseñados para sistemas de ecuaciones diferenciales ordinarias u ODE, así que previamente hemos de transformar el sistema DAE en un sistema ODE (Ordinary Differential Equations) o en un sistema que sea asimilable a un ODE. El sistema utilizado será una formulación de Lagrange aumentado para resolver el problema dinámico en tiempo real durante la simulación. Este sistema exige partir de un campo de aceleraciones inicial para arrancar la integración numérica. Dicho campo de aceleraciones se calculará con la formulación de los penalizadores (Bayo, García de Jalón & J. & Serna 1988). 3.3.2. El primer problema dinámico: el problema de aceleración inicial Este problema consiste en resolver el sistema (3.11) relacionando las variables del modelo con los multiplicadores de Lagrange. Así puede prescindirse de las ecuaciones (3.11) b) para obtener un sistema reducido con m ecuaciones menos en el que los multiplicadores de Lagrange no son incógnitas. Los multiplicadores están relacionados con las fuerzas que actúan en el mecanis mo para que las ecuaciones de restricción que ligan las variables se cumplan en todo instante. El método de los penalizadores (Bayo, García de Jalón & J. & Serna 1988) 3-7.
(46) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. impone que los multiplicadores sean proporcionales al incumplimiento de las restricciones y sus derivadas. Se sustituyen las restricciones por sistemas muelle-amortiguador con amortiguamie nto crítico, obteniendo la siguiente expresión para λ* :. λ  α  Φ  2Φ   2Φ . (3.12). Donde:.  : matriz de penalizadores, es decir la rigidez de los muelles correspondientes a cada restricción. Usualmente toma valores entre 106 y 1011 ..  : amortiguamiento relativo. Se le suele dar el valor ξ=1(amortiguamiento crítico).  : frecuencia natural del sistema. Se le suele dar el valor ω= 10. Sustituyendo la ecuación (3.12) en la ecuación (3.10) se obtiene:. Mq  ΦqTα  Φ  2Φ   2Φ   Q. (3.13). El vector de segundas derivadas de las restricciones Φ lo calculamos derivando dos veces el vector de restricciones con respecto al tiempo:.  Φ Φ   q dq  t dt  dΦ    Φ qΦ Φ  q t dt dt. (3.14). Donde Φ t es el vector de derivadas parciales temporales de las restricciones. Derivando una vez más en (3.14) se obtiene:. Φ  Φqq  Φqq  Φt. (3.15). Sustituyendo (3.15), en (3.13) y reagrupando términos se obtiene (3.16), las ecuaciones dinámicas del método de los penalizadores. Las incógnitas de este sistema son únicamente las aceleraciones q . 3-8.
(47) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo.  M  Φ αΦ  q  Q  Φ α Φ q  Φ  2Φ   Φ  T q. q. T q. 2. q. t. (3.16). 3.3.3. El problema dinámico en tiempo real: formulación ALI3-P Lo explicado anteriormente nos servirá para calcular el valor de los multiplicadores de Lagrange con suficiente precisión como para arrancar la formulación de Lagrange aumentado en index-3 con proyecciones ortogonales o ALI3-P (Index 3 Augmented Lagrange with Projections), que serán el conjunto de ecuaciones empleadas para resolver la dinámica en cada instante de tiempo. La formulación ALI3-P es tremendamente compleja y una descripción pormenorizada se sale del ámbito de este proyecto. La descripción completa de la última versión de la fomulación se puede consultar en (Dopico et al. 2014), aunque la descripción dada aquí corresponde a una versión previa de la formulación publicada en (Cuadrado et al. 2000). La expresión de la formulación ALI3-P es:. Mq  ΦqTαΦ  ΦqT λ*  Q. (3.17). Siendo: M: matriz de masas q: vector de variables del problema Φ q: matriz jacobiana del vector de restricciones λ* : vector de multiplicadores de Lagrange Q: vector de fuerzas generalizadas Los multiplicadores de Lagrange son obtenidos por un proceso iterativo λ*i 1  λ*i  αΦi 1 , i  0,1, 2.... (3.18). 3-9.
(48) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. λ *0 se inicializa con valor nulo en la primera iteración t=0. En los siguientes pasos de. tiempo t=nΔt, el valor de λ *0 será el valor de λ *n 1 . Cuando el proceso converja se obtendrán los valores exactos de λ *n . El integrador empleado es la regla trapezoidal implícita de paso simple. Las ecuaciones en velocidades y aceleraciones son las siguientes:. qn1 . 2 2  qn1  qˆ n con qˆ n    qn  qn  t  t . (3.19). qn1 . 4 ˆ con qˆ    4 q  4 q  q  q  q n  1 n n n n  2 n t 2 t  t . (3.20). Podemos comprobar que en esta última ecuación aparece el término q n perteneciente al paso de tiempo anterior. Esto justifica el cálculo del campo de aceleraciones inic ia l mediante el método de los penalizadores. En los siguientes pasos de tiempo. t   n  1 t se tomará el vector de aceleraciones correspondiente al estado t  ndt . Introduciendo las ecuaciones (3.19) y (3.20) en la ecuación (3.17) se obtiene el sistema:.  4  M  2 qn1  qˆ n   ΦTqn1  αΦn1  λ n1   Qn1  0  t . (3.21). 4 Mqn1  Mqˆ n  ΦTqn1 αΦn1  λ n1   Qn1  0 2 t. (3.22). Escalando por. Mqn1 . 4 se obtiene: t 2. t 2 Mqˆ n  ΦTqn1 αΦn1  ΦTqn1 λ n1  Qn1  0 4. . . (3.23). Para resolver el sistema no lineal (3.21) se usa el método iterativo de Newton-Raphson.. 3-10.
(49) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo.  f  q     qi   f (q)i  q  i. (3.24). Siendo:. f (q) . t 2  Mq  ΦTq αΦ  ΦTq λ*  Q 4. (3.25). La matriz tangente exacta es complicada de obtener, por tanto, se utiliza la matriz tangente aproximada siguiente:  f  q   t t 2  M  C  ΦTq αΦ q  K     2 4  q . (3.26). Donde:. K. Q : matriz de rigidez. q. C. Q : matriz de amortiguamiento. q. Este proceso conduce a una serie de posiciones qn+1 que satisfacen las ecuaciones del movimiento (3.21) y las ecuaciones de restricción. Sin embargo, es posible que no se cumpla Φ  0 y Φ  0 , por lo que para solventar este problema se emplearán las proyecciones ortogonales en velocidades y aceleraciones. * * Si q y q son las velocidades obtenidas después de la convergencia del método de. Newton-Raphson, las velocidades y aceleraciones depuradas se calculan según las siguientes expresiones:    t t 2 t t 2  * t 2 T T M  C  Φ αΦ  K q  M  C  Kq  Φ q αΦt   q q    2 4 2 4 4    . (3.27). 3-11.
(50) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo.   t t 2 T M  2 C  4  Φq αΦq  K  q      * t 2 T t t 2  M  C  K    q  Φ q α  Φ q q  Φt  2 4 4  . 3.4.. (3.28). Modelización elegida para el sistema de fondeo. Para modelizar cualquier sólido necesitamos cuatro entidades (12 coordenadas generalizadas), ya sean puntos o vectores, definiendo siempre al menos un punto y al menos un vector. Esto nos permitirá representar el sólido en cuestión en un espacio tridimensional y orientarlo, ya que construiremos, para cada sólido, un sistema de referencia formado por un punto y tres vectores linealmente independientes. Se puede demostrar que, definiendo los sólidos de esta manera, la matriz de masas y el vector de fuerzas generalizadas gravitatorias son constantes y no aparecerán en ellos términos de fuerzas de inercia dependientes de la velocidad o fuerzas dependientes de la velocidad respectivamente. Esto nos permite calcularla solamente en el instante inicial y prescindir de un cálculo continuo. El software SIMULANCLA emplea la biblioteca de sistemas multicuerpo MBSLIM (Multibody Systems en Laboratorio de Ingeniería Mecánica) para llevar a cabo la tarea de modelización del sistema de fondeo, planteamiento y resolución de las ecuaciones del movimiento. Simulancla emplea las siguientes restricciones entre las coordenadas naturales para definir un sólido rígido: 1. Distancia: para puntos se impone que la distancia entre dos puntos sea constante. En el caso de vectores, la norma del vector ha de ser la unidad. Es una sola ecuación.. p  p  p  p   l T. i. 3-12. j. i. j. 2 ij. 0. (3.29).
(51) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. viT vi  1  0. (3.30). 2. Producto escalar constante: el ángulo entre dos vectores se mantendrá constante. Si el sólido se define con más de un punto se tomarán como vectores los que unen esos puntos. Es una sola ecuación.. vi T v j  kk  0. (3.31). p  p . (3.32). T. i. j. vk  kl  0. p  p  p  p   k T. i. j. i. j. m. 0. (3.33). 3. Combinación lineal: los elementos de modelización del solido forman un subespacio. Los puntos y vectores definidos a mayores para facilitar la definición del modelo pertenecerán a ese subespacio (serán linealme nte dependientes) y además su posición en coordenadas locales permanecerá constante el tiempo, es decir los valores de los parámetros de combinac ió n lineal serán constantes. Con esta restricción ocurre lo mismo que con la de producto escalar, si se modeliza con más de un punto se tomará como vector el que los une. Se trata de una ecuación vectorial lo que equivale a tres ecuaciones. pi  p j   k 1  v l   k 2  v m   k 3  v n  0. (3.34). vi k1  vl k 2  vm k 3  vn  0. (3.35). Donde pi y vi son el punto o vector extra que no pertenece a la definición del sólido. Estas restricciones son incluidas automáticamente cuando se modeliza un sólido en MBSLIM. Ahora resta añadir las restricciones externas dependientes del mecanismo. Para definir los GDL utilizaremos coordenadas relativas, creando así las variables adicionales que nos permitirán el guiado. 3-13.
(52) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. 3.4.1. Definición del buque El suelo se define con el punto p1 y los vectores v1 , v2 y v3. Esto forma un sistema de coordenadas fijo al suelo del muelle, así que no es necesario incluir estas entidades en el vector de variables, pero sí precisaremos de ellas para la modelización gráfica. 3.4.2. Definición de sólidos del sistema de fondeo Eslabones de la cadena. Se describirán a continuación las variables y restricciones para dos eslabones consecutivos de la cadena.. Figura 3-1. Restricciones para los eslabones. 3-14.
(53) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. Variables de los elementos: pi. . Puntos:. . Vectores:. ,. p j pk ,. v i v j vk v m , , , .. Elementos modelización: Dos puntos y dos vectores por sólido. Hay un punto compartido entre eslabón y eslabón. Restricciones adicionales: vTi v k  0. vTj v m  0 Número de grados de libertad: El número de grados de libertad entre sólidos, teniendo en cuenta las restriccio nes automáticas y adicionales es de 2 entre eslabones de la cadena. Grillete y caña Se describirán a continuación las variables y restricciones para el grillete y la caña.. 3-15.
(54) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. Figura 3-2. Restricciones para el grillete y la caña. Variables elemento: . Puntos: pi, pj. . Vectores: vi, vj, vk , vm, vz. Elementos modelización: Un punto y dos vectores para el grillete y dos puntos y tres vectores para la caña. Restricciones: No son necesarias restricciones adicionales a las impuestas por la MBSLIM debido a que se comparte un punto y el vector del eje.. 3-16.
(55) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. Número de grados de libertad: Realizando el cómputo de las variables y restricciones impuestas, sale un grado de libertad entre el grillete y la caña. Caña y ancla Se describirán a continuación las variables y restricciones para la caña y el ancla.. Figura 3-3. Restricciones para la caña y las uñas del ancla. Variables elemento: . Puntos: pi, pj. . Vectores: vi, vj, vk , vm, vq, vz. . Ángulo  3-17.
(56) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. Elementos modelización: Dos puntos y tres vectores para la caña y un punto y tres vectores para las uñas del ancla. Restricciones: No son necesarias restricciones adicionales a las impuestas por la MBSLIM debido a que se comparte un punto y el vector del eje. Número de grados de libertad: Realizando el cómputo de las variables y restricciones impuestas, sale un grado de libertad entre la caña y las uñas del ancla. 3.4.3. Grados de libertad Haremos en primer lugar un balance de ecuaciones/variables para poder apreciar con mayor claridad los GDL.. Sólido. Grados de libertad relativos. Eslabón 1. 3. …. …. Eslabón “i”. 2. Caña. 1. Ancla. 1. Tabla 3-1 Balance de grados de libertad.. 3-18.
(57) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. 3.4.4. El vector de variables Ahora que conocemos las variables y restricciones que van a condicionar la simulación, podemos crear un vector, llamado vector de variables, que recogerá todas las variables que se introdujeron mediante la modelización de cada sólido más las variables correspondientes a la carga, más las que definen cada uno de los grados de libertad. Este vector q variará en cada paso integrable del tiempo, recogiendo el valor de cada variable. Tendrá la siguiente forma:. q  pi con i 2,3..,20 , vi con i 4,5..,16 , i con i 1,2...5 , s1 , s2 , s3 . (3.36). 3.4.5. Fuerzas Fuerzas normales de contacto El modelo de fuerzas de contacto normales empleado en este trabajo, es el modelo desarrollado por Flores en el 2011 en (Flores et al. 2011). El modelo es muy similar al modelo de Hunt-Crossley, desarrollado originalmente en (Hunt & Crossley 1975). La expresión de la fuerza normal del modelo de Flores es la siguiente:  8 1  Fn  kn n  1   5 .  n .  0 . . (3.37). donde k n es la rigidez equivalente del contacto y depende de la forma y el material de los sólidos que colisionan, n es el exponente de Hertz,   Rsph  pcenter  pcontact es la .. indentación,  es su derivada temporal,  0 es la velocidad normal relativa entre los sólidos en colisión cuando el contacto es detectado,. es el coeficiente de restitució n. y n es la dirección de la fuerza de contacto normal.. 3-19.
(58) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. Figura 3-4. Contacto normal. El valor de k n se puede calcular para dos sólidos en contacto con geometrías de esfera y plano respectivamente en los que la expresión para la rigidez puede ser expresada por. kn . 4 Rsph 3  sph   pln . (3.38). siendo Rsph el radio de la esfera en contacto con el plano y los parámetros de sus respectivos materiales  sph y  pln vienen dados por.  sph . 2 1  sph. Esph. ;  pln . 2 1  pln. E pln. (3.39). donde  y E son el cociente de Poisson y el módulo de Young o módulo elástico para cada uno de los dos materiales.. 3-20.
(59) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. Fuerzas de rozamiento de contacto. Figura 3-5. Contacto tangencial. El modelo de fuerzas de contacto tangenciales, o fuerzas de fricción, fue desarrollado en detalle en (Dopico et al. 2011) y está basado en la ley de Coulomb del rozamie nto seco incluyendo sticktion (adhesión) y un término de disipación viscosa. La forma general de la fuerza es la siguiente,. Ft  Fstick  (1 ) Fslide  visc vt. (3.40). En la expresión previa, los primeros dos términos constituyen el rozamiento seco, mientras que el tercer término representa el rozamiento viscoso. Para tener una transición suave entre sticking (adhesión) y slipping (deslizamiento), la fuerza de fricción se divide en dos componentes acopladas mediante una función suave, siguiendo las ideas propuestas en (Gonthier et al. 2004). El subíndice “t” indica aquí “tangencial”. En la ecuación (3.40), visc es el coeficiente de rozamiento viscoso, Fstick y Fslide son las componentes de sticktion o adhesión y slipping o deslizamiento, es una funció n suave de la velocidad tangencial en el punto de contacto, v t que está definida en. 3-21.
(60) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. función del punto central de la región de contacto, p contact y el vector normal al contacto, n, como sigue: .  .  vt  pcontact   nT pcontact  n  . (3.41). La citada función, tiene que cumplir las siguientes condiciones, 0;  1; . vt  vstick    vt  0  .  . (3.42). donde vstick es un parámetro del modelo que tiene en cuenta la velocidad de transició n entre adhesión y deslizamiento. Una buena elección para esta función de transición fue proporcionada por (Gonthier et al. 2004) y tiene la siguiente expresión:. . . 2  vtT vt / vstick.  e. (3.43). La ecuación (3.40) muestra que la fuerza total está compuesta de tres componentes: la fuerza de rozamiento seco de deslizamiento cuando las velocidades son significativas, la fuerza de rozamiento seco adhesivo a bajas velocidades y la fuerza de rozamie nto viscoso. La fuerza de rozamiento seco de deslizamiento viene dada por la expresión clásica de Coulomb, mientras que la fuerza de rozamiento seco adhesivo viene dada por las expresiones de los bristles, que son elementos viscoelásticos en forma de pelos o cerdas actuando entre los cuerpos en colisión. Las expresiones de las fuerzas de deslizamiento y adhesión vienen dadas por las ecuaciones (3.44) y (3.45):. Fslide. 3-22.  0;   vt din Fn v ; t . vt  0   vt  0 . (3.44).
(61) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. Fstick. 0;   m   f stick  I3  nnT   pcontact  p stick  ;   s. s  0   s  0 . (3.45). siendo  din el coeficiente de fricción dinámico, s  pcontact  pstick la deformación de los bristles, p stick el punto de adhesión, que inicialmente coincide con el centroide de la región de colisión en el instante en el que el contacto comenzó; I 3 es la matriz de m identidad de tamaño 3 x 3; f stick es la función que representa el comportamiento de los. bristles, m f stick  kstick s  cstick s. (3.46). siendo k stick y cstick los coeficientes de rigidez y amortiguamiento del modelo de adhesión. Sin embargo, hay un valor límite para la fuerza de adhesión:. Fstick  st Fn. (3.47). En (3.47),  st es el coeficiente de fricción estático que en general es más elevado que el dinámico. Si se excede este límite y (3.47) no se cumple, se producen dos consecuencias: en primer lugar, la ecuación (3.46) ya no es válida para el comportamiento de los bristles y tiene que ser sustituida por la ecuación (3.48); en segundo lugar, el punto de adhesión debe ser actualizado con la expresión (3.49): m f stick . I.   st Fn s. 3.  nn T   p contact  p stick .   F  v p stick  p contact   stick st n  t k stick   vt. (3.48). (3.49). El coeficiente  stick controla la deformación de los bristles cuando se alcanza la máxima fuerza. Físicamente el valor más razonable es stick  1, pero pequeñas variaciones con stick  1 pueden mejorar el comportamiento numérico del modelo. 3-23.
(62) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. Fuerza tope de ancla y caña La naturaleza de esta fuerza tiene una forma similar a la de un muelle-amortiguador que actúa por tramos. El parámetro fundamental es el ángulo que forman las uñas del ancla con la caña. Este ángulo tiene una rotación máxima, para este caso en particular, de 30º respecto a la vertical para evitar que las uñas puedan oscilar libremente y provocar problemas en el proceso de estibado.. Figura 3-6. Ángulo del ancla con la caña. La expresión de la fuerza queda:. f. lin m. si   max km  (  max )  cm   ;   km  (  max )  cm   ; si   max  0 si  max    max . (3.50). siendo km la constante elástica, cm la constante de amortiguamiento, φ el ángulo, φp su derivada, φ0 el ángulo natural y fm la fuerza de ancla y caña. Como se puede ver en (3.50), pueden presentarse tres casos. El primero es para.    máx , el segundo para   máx y el tercero para máx    máx . Para este último caso la fm será nula.. 3-24.
(63) UNIVERSIDADE DA CORUÑA Escola Politécnica Superior La dinámica de sistemas multicuerpo. 3-25.
(65) Capítulo 4 SIMULANCLA: un software para la simulación del sistema de fondeo de buques.
Modelado de Hueso Trabecular con Apoyo de Impresión de Réplicas 3D_Waldir L. Roque
Los resultados para un conjunto de muestras de imágenes µCT de hueso trabecular del radio distal han demostrado que una estructura con baja fracción de volumen, baja conectividad, alta

References: resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución