Depositos Sismo
Depositos Sismo
Depositos Sismo
Agradecimientos
INDICE GENERAL.
1. OBJETIVOS 11
2. INTRODUCCIN 13
2.1 QU ES UN DEPSITO? 13
2.2 TIPOS DE DEPSITOS 13
2.2.1 Segn forma 13
2.2.2 Segn presin de trabajo 14
2.2.3 Segn altura 15
2.2.4 Abiertos o cerrados 16
2.2.5 Segn material 17
2.2.6 Segn orientacin. 17
2.3 COMPORTAMIENTO DE DEPSITOS NO ELEVADOS
DURANTE UN SISMO 18
2.4. EVOLUCION DEL ANLISIS SISMICO DE DEPSITOS 22
2.4.1 Estado Del Arte. 22
3. METODOLOGA 25
3.1 MODELO DE ANALTICO. 25
3.1.1 Modelo Housner para un depsito circular 27
3.2 MODELO NUMRICO 29
3.2.1 Generalidades 29
3.2.2 Aproximacin del modelo con elementos finitos. 30
3.2.3 Seleccin del elemento fluido 31
6. CONCLUSIONES 112
BIBLIOGRAFA 117
ANEXOS 119
INDICE DE FIGURAS.
Figura 39. Segundo modo del sistema depsito fluido (modo 305 del anlisis) 63
Figura 40. Tercer modo del sistema depsito fluido (modo 337 del anlisis) 63
Figura 41. Tensin de Von Mises en el contorno del depsito 64
Figura 42. Tensin de Traccin Anular en el contorno del depsito 65
Figura 43.Tensin de Compresin Vertical en el contorno del depsito 65
Figura 44. Depsito de relacin H/D = 0.2 67
Figura 45. Modelo ANSYS del depsito relacin H/D = 0.2 68
Figura 46. Tensin de traccin anular para el depsito relacin H/D = 0.2 69
Figura 47. Tensin de compresin vertical para el depsito relacin H/D = 0.2 69
Figura 48. Tensin de Von Mises para el depsito relacin H/D = 0.2 70
Figura 49. Depsito de relacin H/D = 1 71
Figura 50. Modelo ANSYS del depsito relacin H/D = 1 72
Figura 51. Tensin de traccin anular para el depsito relacin H/D = 1 73
Figura 52. Tensin de compresin vertical para el depsito relacin H/D = 1 73
Figura 53. Tensin de Von Mises para el depsito relacin H/D = 1 74
Figura 54. Depsito de relacin H/D = 2 75
Figura 55. Modelo ANSYS del depsito relacin H/D = 2 76
Figura 56. Tensin de traccin anular para el depsito relacin H/D = 2 77
Figura 57. Tensin de compresin vertical para el depsito relacin H/D = 2 77
Figura 58. Tensin de Von Mises para el depsito relacin H/D = 2 78
Figura 59. Estructura para el depsito elevado 81
Figura 60. Representacin del depsito de relacin H/D = 0.4 y elevado 1m. 82
Figura 61. Modo 1 del fluido(1 del anlisis) para el depsito relacin
H/D = 0.4 elevado 1m 84
Figura 62. Modo 2 del fluido(5 del anlisis) para el depsito relacin
H/D = 0.4 elevado 1m 84
Figura 63. Modo 3 del fluido(13 del anlisis) para el depsito relacin
H/D = 0.4 elevado 1m 85
Figura 64. Modo 1 del sistema fluido-estructura para el depsito relacin
H/D = 0.4 elevado 1m 85
Figura 65. Modo 2 del sistema fluido-estructura para el depsito relacin
H/D = 0.4 elevado 1m 86
Figura 66. Tensin de traccin anular para el depsito de relacin
H/D = 0.4 y elevado 1m. 87
Figura 67. Tensin de compresin vertical para el depsito de relacin
H/D = 0.4 y elevado 1m. 87
Figura 68. Tensin de Von Mises para el depsito de relacin
H/D = 0.4 y elevado 1m. 88
Figura 69. Representacin del depsito de relacin H/D = 0.4 y elevado 5m. 89
Figura 90. Modo 3 del fluido(13 del anlisis) para el depsito relacin
H/D = 1 elevado 5m 106
Figura 91. Modo 1 del sistema fluido-estructura para el depsito relacin
H/D = 1 elevado 5m 106
Figura 92. Modo 2 del sistema fluido-estructura para el depsito relacin
H/D = 1 elevado 5m 107
Figura 93. Tensin de traccin anular para el depsito de relacin
H/D = 1 y elevado 5m. 108
Figura 94. Tensin de compresin vertical para el depsito de relacin
H/D = 1 y elevado 5m. 108
Figura 95. Tensin de Von Mises para el depsito de relacin
H/D = 1 y elevado 5m. 109
INDICE DE TABLAS.
1
Objetivos
1. OBJETIVOS.
2
Introduccin
2. INTRODUCCIN.
2.1 QU ES UN DEPSITO?.
Existen varios criterios para clasificar los depsitos. Una de las maneras
ms bsicas de clasificarlos, es considerando la forma del mismo. Algunas de
las formas ms comunes son, rectangular, cuadrada y cilndrica.
- Daos en el techo del depsito causado por los modos de chapoteo del
lquido.
3
Metodologa
3. METODOLOGA.
Efectos impulsivos:
1.7 * R
tanh
M0 = H
*M
1 .7 * R
H
M
H 0 = 0.38 * H 1 + 1
M0
Efectos convectivos:
1.8 * H
0.71 * tanh
M1 = R
*M
1 .8 * H
R
2
R*M
2
M R
1
R
H 1 = H 1 0.21 + 0.55 0.15
M1 H H H * M1
y
2
4.75 * g * M 1 * H
K=
M * R2
Donde:
R = Radio del depsito.
H = altura del lquido contenido en el depsito.
g = aceleracin de la gravedad (9,81 m/s2).
M = masa total del lquido contenido en el depsito.
M0 = masa impulsiva del lquido.
H0 = altura a la cual est aplicada la masa impulsiva (M0).
M1 = masa convectiva del lquido.
H1 = altura a la cual est aplicada la masa convectiva (M1).
K = rigidez requerida para la oscilacin de la masa convectiva (M1).
3.2.1 Generalidades.
Los nodos del fluido estn acoplados con los nodos de la pared del
depsito en la direccin radial, mientras que los movimientos relativos en la
direccin tangencial y vertical pueden ocurrir libremente.
Este elemento tridimensional est definido por ocho nodos, teniendo tres
grados de libertad en cada nodo: traslaciones nodales, en x, y, z
respectivamente. Adems este elemento tiene resortes especiales, que tienden
a retardar los movimientos hidrodinmicos ya que en la superficie libre del
En general, los cambios de presin del fluido, asociados con los modos
convectivos, son muy pequeos y, a menudo, son despreciables en los
problemas de interaccin fluido-estructura, resueltos por el mtodo de
elementos finitos. Sin embargo, dentro del modelo de elementos finitos, el
comportamiento convectivo puede ser fcilmente incluido, al agregar resortes
verticales bajo la primera capa de elementos finitos a partir de la superficie
libre. Este efecto ya lo tiene incorporado el software utilizado.
4
Estudio analtico del depsito
apoyado sobre el terreno
Geometra Valor
Dimetro 9m
Altura 3.5 m
Nivel del agua 3m
Espesor paredes 5 mm
Espesor techo 4 mm
Espesor del suelo 5 mm
Tabla 1. Caractersticas geomtricas del depsito de relacin H/D = 0.4.
Propiedad Valor
Composicin H2O
Densidad 1000 Kg/m3
Modulo elasticidad ( 20C ) 2200 MPa
Tabla 4. Propiedades del fluido a almacenar.
Donde:
Masa total lquido:
K = 320539,43 N
m
X 1 = 1.14 m
W LIQ R 2 R * W LIQ
2
R
1
X 2 = H 1 0.21 + 0.55 0.15
W2 H H H * W2
X 2 = 1.67 m
a c = S * * ab
Donde:
- ab: Valor caracterstico de la aceleracin horizontal de la
superficie del terreno, que se define en el mapa de peligrosidad
ssmica (figura 22).
- Para
* ab 0.4 * g
S = 1
siendo:
C: Coeficiente de terreno.
I 1,0
II 1,3
III 1,6
IV 2,0
Tabla 5. Coeficiente de terreno.
ab = 0.04 * g Madrid
Si T T A (T ) = 1 + 1.5 T T
A
Si T A T TB (T ) = 2.5
Si T TB (T ) = K * C T
A
Donde:
(T) : Valor del espectro normalizado de respuesta elstica.
K = 1 (Madrid. Figura 22).
C = 1.3 (Coeficiente de suelo).
TA y TB son los perodos caractersticos del espectro de respuesta, de valores:
TA = K * C = 0.13
10
TB = K * C = 0.52
2.5
Q = Z * I (C1 * WS + C1 * W R + C1 * W1 + C 2 * W2 )
M = Z * I (C1 * WS * X S + C1 * W R * H T + C1 * W1 * X 1 + C 2 * W2 * X 2 )
Donde:
Q : Corte Basal.
M : Momento Volcante.
Z : Factor de zona ssmica. (factor de la aceleracin ssmica Horizontal).
I: Factor de Importancia que puede ser tomado como 1,0 a no ser que se
especifique otra cosa. En todo caso, el valor de I, no debe ser superior a
1,25 y este valor slo debe aplicarse a depsitos estratgicos en caso de
terremotos.
C1, C2 : Coeficientes ssmicos de masas impulsivas y convectivas
respectivamente.
Ws : Peso total de las paredes del depsito.
Xs : Distancia del fondo del depsito al centro de gravedad del depsito.
WR : Peso total del techo (Planchas, vigas, etc.).
WT : Peso del fluido.
HT : Distancia del fondo del depsito al centro de gravedad del techo.
W1 : Peso impulsivo.
X1 : Distancia desde el fondo del depsito al centro de accin de la fuerza
ssmica aplicada sobre W1.
W2 : Peso Convectivo (Modo fundamental).
X2 : Distancia desde el fondo del depsito al centro de accin de la fuerza
ssmica aplicada sobre W2.
D : Dimetro nominal del depsito.
H : Altura del Fluido.
WLIQ = WT = 190851,75 Kg
W1 = 73936.75 Kg
W2 = 94136.79 Kg
X 1 = 1.14 m
X 2 = 1.67 m
WS = 3938.61 Kg
WR = 2968.77 Kg
Para obtener el coeficiente ssmico C2 la norma nos dice que su valor es:
0.4
2.75 * A0 T 0.05
n
C2 =
g *R T*
Donde:
A0 es la aceleracin efectiva mxima A0 = 0.04g
R es el factor de modificacin de respuesta R = 4
1.33 0.4
2.75 * 0.04 g 0.35 0.05
C2 = = 4.91e 3
g *4 3.42 0.02
Q = 1049.87 Kg
3 .5
M = 0.04 * 1 0.32 * 3938.6 * + 0.32 * 2968.7 * 3.5 + 0.32 * 73936 .7 * 1.14 + 0.004 * 94136 .7 * 1.67
2
M = 1325.2 Kg * m
Segn la API-650:
wL = 99 * t b Fby * G * H N
m
Donde:
wL = peso mximo del contenido del depsito que puede ser usado para
resistir el momento de vuelco, en N / m.
tB = espesor de la placa de fondo en el depsito, en mm.
FBY = lmite de fluencia del acero, MPa.
G= gravedad de diseo especfico del lquido almacenado.
H= nivel de diseo mximo de lquido, en m.
D= dimetro nominal del depsito, en m.
wL = 99 * 5 206 *1 * 3 = 12305.5 N
m
196 * G * H * D = 5292 N
m
Entonces:
wL = 5292 N
m
M 1325.2
F= = 2 = 0.022
D (wL + wt ) 9 (540 + 205.45)
2
b 226.28
= =
12 * t 12 * 5
= 3.77 Kg
cm 2
83 * t G * H * D2
Fa = Si 44
D t2
83 * t G * H * D2
Fa = + 7.5 G * H Si < 44
2 .5 D t2
Donde:
G= gravedad de diseo especfico del lquido almacenado.
H= nivel de diseo mximo de lquido, m.
D= dimetro nominal del depsito.
G * H * D 2 1* 3 * 9 2
= = 9.72 < 44 entonces:
t2 52
83 * t 83 * 5
Fa = + 7.5 G * H = + 7.5 1 * 9
2 .5 D 2 .5 * 9
Fa = 31.43MPa
Para la modelizacion del depsito con ANSYS, se han utilizado los tipos
de elementos y las hiptesis descritas en el capitulo 3. A continuacin se
muestran una serie de imgenes capturadas de la interfaz del programa
ANSYS que representan el depsito y las fases de mallado del mismo.
Parmetro Valor
Tensin lmite 235
elstico fy (N/mm2) (t < 16 mm)
H 220 mm
B 110 mm
e1 9.2 mm
E 5.9 mm
Figura 27. Perfil viga IPE. Ix 2772e4 mm4
Iy 204.9e4 mm4
Tabla 6. Parmetros viga IPE 220.
z VM r
Traccin anular Compresin vertical Tensin Von Mises radial
Valor
252 258 275 0
(Kg/cm2)
Tabla 7. Tensiones obtenidas en el anlisis esttico.
Segn los valores para los periodos obtenidos para cada modo de
vibracin y comparando estos con los periodos pico en el espectro de
respuesta dado por la norma, se puede concluir que el modo predominante de
vibracin ser un modo del fluido solo, por lo que se vera durante el sismo el
agua vibrar segn este modo.
Donde:
j son las races de la funcin J1, donde J1 es la funcin de Bessel de primera
especie. Los primeros cuatro valores de son:
j 1 2 3 4
j 1.84119 5.33145 8.53633 11.706
Tabla 9. Races de la funcin J1.
Figura 33. Primer modo convectivo del fluido (modo 1 del anlisis).
Figura 34. Segundo modo convectivo del fluido (modo 6 del anlisis).
Figura 35. Tercer modo convectivo del fluido (modo 13 del anlisis).
Figura 36. Cuarto modo convectivo del fluido (modo 23 del anlisis).
Los resultados para las formas modales del sistema depsito-fluido, que
se obtienen a continuacin, se resumen en la tabla 11.
Kh 1 E
f0 =
2 H
Donde:
Kh es una constante, que se determina a travs de la relacin H/R y de t /
(1000*R), mediante la grafica siguiente:
Segn lo anterior:
f 0 = 18.44 s 1
1 1 1 1
= + +
(f )
'
0
2
f0
2
(f ) (f )
'
F
2 '
S
2
1 KF
f F' =
2 Wte
1 KS
f S' =
2 Wte
K F = 3e10 Kg
s2
R * E *t 4.5 * 190e9 * 5e 3
KS = =
2(1 + ) H 2(1 + 0.3) 3
K F = 1.72e9 Kg
s2
Wte es la masa del techo:
Wte = R 2 * t techo * acero = * 4.5 2 * 7960 = 2025.57 kg
Entonces :
1 KF
f F' = = 505.9 s 1
2 Wte
1 KS
f S' = = 121.14s 1
2 Wte
f 0' = 18.21s 1
Figura 38. Primer modo del sistema depsito fluido (modo 277 del anlisis).
Figura 39. Segundo modo del sistema depsito fluido (modo 305 del anlisis).
Figura 40. Tercer modo del sistema depsito fluido (modo 337 del anlisis).
Tensin de
Corte basal Momento de vuelco
METODO compresin
(Kg) (Kg*m)
(Kg/cm2)
ANSYS 930.2 1382 3.34
API 650 1049.94 1325.2 3.77
Tabla 12. Comparacin de los resultados entre ANSYS y norma API-650.
Geometra Valor
Dimetro 12 m
Altura 2.4 m
Nivel del agua 1.7 m
Espesor paredes 5 mm
Espesor techo 4 mm
Espesor suelo 5 mm
Tabla 14. Geometra depsito relacin H/D = 0.2.
Figura 46. Tensin de traccin anular para el depsito relacin H/D = 0.2.
Figura 47. Tensin de compresin vertical para el depsito relacin H/D = 0.2.
Figura 48. Tensin de Von Mises para el depsito relacin H/D = 0.2.
Compresin
Momento de
Mtodo Corte basal (Kg.) vertical mxima
vuelco (Kg*m)
(Kg/cm2)
ANSYS 262.53 2339.7 3.8
Tabla 17. Resultados anlisis espectral para el depsito relacin H/D = 0.2.
Geometra Valor
Dimetro 6.5 m
Altura 6.5 m
Nivel del agua 5.75 m
Espesor paredes 5 mm
Espesor techo 4 mm
Espesor suelo 5 mm
Tabla 18. Geometra depsito relacin H/D = 1.
Compresin
Momento de
METODO Corte basal (Kg) vertical mxima
vuelco (Kg*m)
(Kg/cm2)
ANSYS 3370.5 938.7 10.1
Tabla 21. Resultados anlisis espectral para el depsito relacin H/D = 1.
Geometra Valor
Dimetro 5m
Altura 10 m
Nivel del agua 9.7 m
Espesor paredes 5 mm
Espesor techo 4 mm
Espesor suelo 5 mm
Tabla 22. Geometra depsito relacin H/D = 2.
Compresin
Momento de
Mtodo Corte basal (Kg) vertical mxima
vuelco (Kg*m)
(Kg/cm2)
ANSYS 7913.2 2280.8 26.4
Tabla 25. Resultados anlisis espectral para el depsito relacin H/D = 2.
Parmetro
Tensin Tensin
Corte Momento Tensin de
Depsito compresin traccin
Basal Volcante Von Mises
vertical anular
(Kg) (Kg*m) (MPa)
(MPa) (MPa)
H/D = 0.2 262.52 2339 0.38 0.91 0.81
H/D = 0.4 930 1382 0.32 1.3 1.23
H/D = 1 3370 938.7 0.99 1.25 1.35
H/D = 2 7913.2 2280.8 2.59 0.93 2.3
Tabla 26. Comparacin de resultados de los depsitos apoyados sobre el terreno.
5
Estudio del depsito
elevado
Para los casos que se van a estudiar se omitirn las vigas horizontales
que aparecen en la fotografa anterior por dos motivos, para simplificar el
problema y a la vista de los resultados obtenidos para el depsito apoyado se
entiende que no sern necesarias para la estabilidad del depsito. Adems
para una elevacin de un metro de altura no es necesario y su instalacin sera
una tarea compleja.
Figura 60. Representacin del depsito de relacin H/D = 0.4 y elevado 1m.
Parmetro Valor
Dimetro del depsito 9m
Altura del agua 3m
Altura del depsito 3.5 m
Espesor de las paredes 5 mm
Espesor del techo 4 mm
Espesor del suelo 5 mm
Altura del soporte 1m
Tabla 27. Geometra del depsito relacin H/D = 0.4 elevado 1m.
Figura 61. Modo 1 del fluido(1 del anlisis) para el depsito relacin H/D = 0.4 elevado 1m.
Figura 62. Modo 2 del fluido(5 del anlisis) para el depsito relacin H/D = 0.4 elevado 1m.
Figura 63.Modo 3 del fluido(13 del anlisis) para el depsito relacin H/D = 0.4 elevado 1m.
Figura 64.Modo 1 del sistema fluido-estructura para el depsito relacin H/D = 0.4 elevado 1m.
Figura 65.Modo 2 del sistema fluido-estructura para el depsito relacin H/D = 0.4 elevado 1m.
Figura 66. tensin de traccin anular para el depsito de relacin H/D = 0.4 y elevado 1m.
Figura 67. tensin de compresin vertical para el depsito de relacin H/D = 0.4 y elevado 1m.
Figura 68.tensin de Von Mises para el depsito de relacin H/D = 0.4 y elevado 1m.
Figura 69. Representacin del depsito de relacin H/D = 0.4 y elevado 5m.
Parmetro Valor
Dimetro del depsito 9m
Altura del agua 3m
Altura del depsito 3.5 m
Espesor de las paredes 5 mm
Espesor del techo 4 mm
Espesor del suelo 5 mm
Altura del soporte 5m
Tabla 29. Geometra del depsito relacin H/D = 0.4 elevado 5m.
Figura 70. Modo 1 del fluido(1 del anlisis) para el depsito relacin H/D = 0.4 elevado 5m.
Figura 71. Modo 2 del fluido(5 del anlisis) para el depsito relacin H/D = 0.4 elevado 5m.
Figura 72. Modo 3 del fluido(12 del anlisis) para el depsito relacin H/D = 0.4 elevado 5m.
Figura 73. Modo 1 del sistema fluido-estructura para el depsito relacin H/D = 0.4 elevado 5m.
Figura 74.Modo 2 del sistema fluido-estructura para el depsito relacin H/D = 0.4 elevado 5m.
Figura 75. Tensin de traccin anular para el depsito de relacin H/D = 0.4 y elevado 5m.
Figura 76. Tensin de compresin vertical para el depsito de relacin H/D = 0.4 y elevado 5m.
Figura 77. Tensin de Von Mises para el depsito de relacin H/D = 0.4 y elevado 5m.
Si ahora se los dividimos al lmite elstico del acero inoxidable AISI 316,
obtenemos que el mnimo coeficiente de seguridad es de 5.6.
Parmetro Valor
Dimetro 6.5 m
Altura 6.5 m
Altura del agua 5.75 m
Espesor del techo 4 mm
Espesor de las paredes 5 mm
Espesor del suelo 5 mm
Altura del soporte 1m
Tabla 31. Geometra del depsito relacin H/D = 1 elevado 1m.
Figura 79. Modo 1 del fluido(1 del anlisis) para el depsito relacin H/D = 1 elevado 1m.
Figura 80. Modo 2 del fluido(5 del anlisis) para el depsito relacin H/D = 1 elevado 1m.
Figura 81. Modo 3 del fluido(12 del anlisis) para el depsito relacin H/D = 1 elevado 1m.
Figura 82. Modo 1 del sistema fluido-estructura para el depsito relacin H/D = 1 elevado 1m.
Figura 83. Modo 2 del sistema fluido-estructura para el depsito relacin H/D = 1 elevado 1m.
Figura 84. Tensin de traccin anular para el depsito de relacin H/D = 1 y elevado 1m.
Figura 85. Tensin de compresin vertical para el depsito de relacin H/D = 1 y elevado 1m.
Figura 86. Tensin de Von Mises para el depsito de relacin H/D = 1 y elevado 1m.
Si ahora se los dividimos al lmite elstico del acero inoxidable AISI 316,
obtenemos que el mnimo coeficiente de seguridad es de 3.2.
Parmetro Valor
Dimetro 6.5 m
Altura 6.5 m
Nivel del agua 5.75 m
Espesor del techo 4 mm
Espesor de las paredes 5 mm
Espesor del suelo 5 mm
Altura del soporte 5m
Tabla 33. Geometra del depsito relacin H/D = 1 elevado 5m.
Figura 88. Modo 1 del fluido(1 del anlisis) para el depsito relacin H/D = 1 elevado 5m.
Figura 89. Modo 1 del fluido(5 del anlisis) para el depsito relacin H/D = 1 elevado 5m.
Figura 90. Modo 3 del fluido(12 del anlisis) para el depsito relacin H/D = 1 elevado 5m.
Figura 91. Modo 1 del sistema fluido-estructura para el depsito relacin H/D = 1 elevado 5m.
Figura 92. Modo 2 del sistema fluido-estructura para el depsito relacin H/D = 1 elevado 5m.
Figura 93. Tensin de traccin anular para el depsito de relacin H/D = 1 y elevado 5m.
Figura 94. Tensin de compresin vertical para el depsito de relacin H/D = 1 y elevado 5m.
Figura 95. Tensin de Von Mises para el depsito de relacin H/D = 1 y elevado 5m.
Si ahora se los dividimos al lmite elstico del acero inoxidable AISI 316,
obtenemos que el mnimo coeficiente de seguridad es de 3.
6
Conclusiones
6. CONCLUSIONES.
7
Trabajos futuros
7. TRABAJOS FUTUROS.
Bibliografa
BIBLIOGRAFA.
[1]. ANSYS v12 Documentation (User guide, UDF guide, Tutorial guide).
[4]. Barbat, A. (1982): Clculo ssmico de las estructuras. Ed: Editores tcnicos
asociados.
[6]. OPS; CEPIS: Reduccin del dao ssmico : gua para las empresas de agua.
[10]. http://www.thyssenkrupp.cl/Documentos/fichasT/Aceros%20Inoxidables/
AISI%20316-316L.pdf.
[11]. http://www.goodfellow.com/S/AceroInoxidable-AISI316.html.
[12]. http://www.icespedes.com/catalog/cubas-est%C3%83%C2%A1ndar-inox-aisi-
3000-adelante-p-33.html.
[13]. http://fluidos.eia.edu.co/fluidos/propiedades/elasticidad_liquidos.html.
Anexos
/FILNAME,DEPOSITODANI,1
/TITLE,DEPOSITODANI
!
! DEPOSITO CON PLANCHA BASAL, CON TECHO Y SIN GAP USANDO UN SOLIDO COMO
EL FLUIDO
/PREP7
CSYS,1 !COORDENADAS CILINDRICAS
*AFUN,DEG !UNIDADES ANGULARES EN GRADOS
/VUP,,Z !EJE Z EN LA VERTICAL
!----------------------------------------------------------------------------------------
!UNIDADES
!----------------------------------------------------------------------------------------
!Longitud m
!Tiempo sg
!Aceleracin m/sg
!Masa Kgf*sg/m
!Fuerza Kg
!Esfuerzo Kg/m
!Densidad Kgf*/m
!----------------------------------------------------------------------------------------
!GEOMETRIA DEL ESTANQUE
!----------------------------------------------------------------------------------------
!Altura del Lquido : 3 m
!Radio del deposito : 4.5 m
!Altura del deposito : 3.5 m
!Espesor de Pared : 5 mm
!Espesor de Fondo : 5 mm
!Espesor de Techo : 4 mm
!Tipo de Fluido : AGUA
!----------------------------------------------------------------------------------------
!CONSTANTES O VARIABLES
!----------------------------------------------------------------------------------------
G = 9.81 !Aceleracin de gravedad
H1 = 3 !Altura del Lquido
H2 = 3.5 !Altura Total del estanque
H = H2-H1 !Altura del deposito sobre el nivel del lquido
DIAM = 9 !Dimetro del fluido
RAD = DIAM/2 !Radio del fluido
NDIV1 = 12 !Nmero de divisiones en un radio del liquido
NDIV2 = 8 !Nmero de divisiones en la profundidad del lquido
NDIV3 = 1 !Nmero de divisiones en la altura sobre el lquido
NDIV4 = 1 !Numero de divisiones en las vigas de techo y columna
ESP1 = 0.005 !Espesor del manto del estanque
ESP2 = 0.004 !Espesor de plancha de techo
ESP3 = 0.005 !Espesor de plancha de fondo
PW = 1000 !Densidad de masa del agua
PS = 7960 !Densidad de masa del acero
DW = PW/G !Densidad de masa del agua
DS = PS/G !Densidad de masa del acero
EW = 224.26e6 !Modulo de Bulk del agua (metido en Kg/m2)
ES = 1.94E10 !Modulo de elasticidad del acero (metido en Kg/m2)
MODOS = 600 !Nmero de modos de vibrar a considerar en el Anlisis Modal
!
ET,1,FLUID80,,1
MP,EX,1,EW
MP,PRXY,1,0
MP,DENS,1,DW
!
ET,2,SHELL63
MP,EX,2,ES
MP,PRXY,2,0.3
MP,DENS,2,DS
!
R,1,ESP1 !CONSTANTE REAL 1 (MANTO DEL ESTANQUE)
R,2,ESP2 !CONSTANTE REAL 2 (PLANCHA DE TECHO)
R,3,ESP3 !CONSTANTE REAL 3 (PLANCHA BASAL)
!
!----------------------------------------------------------------------------------------
!PARAMENTROS DE VIGAS PRINCIPALES (IPE 220)
!----------------------------------------------------------------------------------------
ET,3,BEAM4
R,4,33.4e-4,2.772e-5,2.049e-6,0.22,0.11
!
!----------------------------------------------------------------------------------------
!GENERACION DEL VOLUMEN DE FLUIDO
!----------------------------------------------------------------------------------------
CYL4,0,0,RAD,90,,,-H1
VGEN,2,1,,,,90
VGLUE,ALL
NUMMRG,KP
NUMCMP,ALL
LSEL,S,,,7,10
LESIZE,ALL,,,NDIV2
LSEL,INVE
LESIZE,ALL,,,NDIV1
LSEL,ALL
ALLSEL
!----------------------------------------------------------------------------------------
!GENERACION DE PLANCHA BASAL
!----------------------------------------------------------------------------------------
TYPE,2
MAT,2
REAL,3
WPOFFS,,,-H1
CYL4,0,0,RAD,90
AGEN,2,10,,,,90
ASEL,S,,,10,11,,1
AGLUE,ALL
NUMMRG,KP
NUMCMP,ALL
ALLSEL
WPOFFS,,,H1
WPSTYLE,,,,,,,,0
!----------------------------------------------------------------------------------------
!GENERACION DEL MANTO DEL ESTANQUE
!----------------------------------------------------------------------------------------
TYPE,2
MAT,2
REAL,1
K,,RAD,0,0 !GENERA EL KEYPOINT 13
K,,RAD,90,0 !GENERA EL KEYPOINT 14
K,,RAD,180,0 !GENERA EL KEYPOINT 15
REAL,1
LSEL,S,,,20,22,2
LSEL,A,,,23
LESIZE,ALL,,,NDIV2
LSEL,S,,,15,24,3
LESIZE,ALL,,,NDIV1
LSEL,S,,,25,27,2
LSEL,A,,,28
LESIZE,ALL,,,NDIV3
LSEL,S,,,26,29,3
LESIZE,ALL,,,NDIV1
ASEL,S,,,12,15,1
MSHKEY,1
AMESH,ALL
ASEL,S,,,12,15,1
ESLA
NSLE
CM,MANTO,NODE
ALLSEL
!----------------------------------------------------------------------------------------
!MALLADO DE LA PLANCHA BASAL
!----------------------------------------------------------------------------------------
TYPE,2
MAT,2
REAL,3
LSEL,S,,,15,19,1
LESIZE,ALL,,,NDIV1
ASEL,S,,,10,11,1
AMESH,ALL
ALLSEL
!----------------------------------------------------------------------------------------
!MALLADO DE PLANCHA DE TECHO
!----------------------------------------------------------------------------------------
TYPE,2
MAT,2
REAL,2
ASEL,S,,,16,19,1
AMESH,ALL
ASEL,S,,,10,11,1
ASEL,A,,,12,15,1
ASEL,A,,,16,27,1
ESLA
NSLE
CM,ESTANQUE,NODE
NUMMRG,NODE
!----------------------------------------------------------------------------------------
!MALLADO DEL FLUIDO
!----------------------------------------------------------------------------------------
TYPE,1
MAT,1
MSHAPE,0,3D
MSHKEY,1
VMESH,ALL
VSEL,S,,,ALL
ESLV
NSLE
CM,FLUIDO,ELEM
/VIEW,,0.8,-1,0.4
EPLOT
!----------------------------------------------------------------------------------------
/SOLU
ANTYPE,MODAL
ACEL,,,G
MODOPT,REDUC,MODOS,,,0,OFF
MXPAND,MODOS,,,YES
LUMPM,ON
OUTPR,ALL,NONE
/OUTPUT,DEPOSITODANI,TXT
SOLVE
SAVE
FINISH
!----------------------------------------------------------------------------------------
!EXPANSION DE LOS RESULTADOS
!----------------------------------------------------------------------------------------
/SOLU
EXPASS,ON
MXPAND,MODOS,,,YES
OUTPR,ALL,NONE
OUTRES,ALL,ALL
SOLVE
SAVE
FINISH
!----------------------------------------------------------------------------------------
!ANALISIS ESPECTRAL
!----------------------------------------------------------------------------------------
/SOLU
ANTYPE,SPECTRUM
PSTRES,ON
SPOPT,SPRS
DMPRAT,0.02
SVTYPE,2,1
FREQ,0.2,0.211,0.222,0.235,0.25,0.267,0.400,0.444,0.500
FREQ,0.571,0.667,0.800,1.0,1.25,1.667,1.887,1.923,3.333,7.143
FREQ,8.333,10.00,20.00,100.00,1000.00
SV,0,0.26,0.274,0.289,0.306,0.325,0.347,0.52,0.578,0.65,
SV,0,0.743,0.867,1.04,1.3,1.625,2.167,2.45,2.5,2.5,2.5,
SV,0,2.5,2.385,1.577,1.115,1.012
SED,1.0,0.0,0.0
CQC
SOLVE
FINISH
/OUTPUT
/POST1
RSYS,0
/INPUT,DEPOSITODANI,MCOM
SAVE