Placas Delgadas Diferencias Finitas
Placas Delgadas Diferencias Finitas
Placas Delgadas Diferencias Finitas
5 (2015)
Resumen: La determinacin de cargas crticas de inestabilidad de placas planas delgadas, constituye un problema que
encuentra gran aplicacin en diversas ramas de la ingeniera tales como la ingeniera civil, mecnica, aeronutica y
naval. Con el objeto de investigar su comportamiento y aplicabilidad, en el presente trabajo, se emplea el Mtodo de
Diferencias Finitas para la determinacin de las cargas de inestabilidad por corte puro en una placa de acero de
pequeo espesor constante y relacin de aspecto igual a uno. Se analiza la influencia sobre los resultados numricos
obtenidos de la discretizacin espacial utilizada. Se tienen en cuenta para esto, dos condiciones de borde de la placa
como son la correspondiente a las cuatro aristas simplemente apoyadas y las mismas perfectamente empotradas. Se
presentan conclusiones generales relativas al trabajo desarrollado, a partir de la comparacin de los resultados
obtenidos con valores de cargas crticas publicados en la bibliografa especializada.
Palabras claves: inestabilidad, placas, corte, diferencias finitas
1. INTRODUCCIN
En la actualidad, se encuentra ampliamente difundida la utilizacin de placas planas de pequeo
espesor en la construccin de estructuras y dispositivos diversos. Este tipo estructural, que comunmente se
utiliza como parte integrante de grandes maquinarias, vigas de acero de gran altura, estructuras
componentes de naves y transportes aeronuticos, martimos y fluviales, entre otros, en los ltimos aos se
ha comenzado a utilizar como sistema de disipacin de energa en estructuras sismorresistentes y en
sistemas estructurales diseados para soportar acciones dinmicas debidas al viento atmosfrico. Los
mismos se encuentran sometidos a una determinada distribucin de acciones externas que los solicitan
generando de esta manera un campo de tensiones internas sobre las placas, especfico para cada caso.
De acuerdo a lo mencionado, un aspecto clave a considerar en el estudio y anlisis de problemas relativos a
placas de pequeo espesor lo constituye la posibilidad cierta de aparicin de fenmenos de inestabilidad
elstica derivados precisamente de los grandes valores de esbeltez que las mismas poseen. Es as que la
determinacin de cargas crticas de inestabilidad, es decir las cargas mnimas a partir de las cuales se
produce el pandeo de la placa bajo una determinada configuracin de acciones exteriores, es de
importancia fundamental para asegurar el adecuado comportamiento de la estructura analizada [1].
La evaluacin de las cargas de bifurcacin mencionadas requiere la resolucin de un modelo matemtico
constituido por ecuaciones diferenciales de cuarto orden. Diversos mtodos tanto numricos como semianalticos permiten, a partir de su adecuada utilizacin al caso en estudio, hallar la solucin buscada, [2].
En este trabajo se investigan las posibilidades de aplicacin del Mtodo de Diferencias Finitas [3] al
problema planteado, considerando un estado de corte puro sobre la placa. La discretizacin de la ecuacin
diferencial conduce a un problema de valores y vectores propios, a partir del cual se obtiene la carga crtica
de pandeo por corte. La solucin obtenida se compara con resultados publicados en la bibliografa
especializada para dos tipos de condiciones de borde, estableciendo de esta manera criterios de utilizacin
del mtodo propuesto para diversas discretizaciones espaciales de la placa.
2. GENERALIDADES
Se considera una placa de acero plana, rectangular de ancho a y altura b, cuyo espesor t es pequeo en
relacin al resto de sus dimensiones, de acuerdo a lo indicado en la Figura 1. Para el caso en estudio la
placa se encuentra sometida a un estado de solicitacin de corte puro en el plano de la misma, distribuido
en forma uniforme sobre la longitud de sus cuatro aristas, tal como puede observarse en dicha figura, donde
Nzy y Nyz, representan las fuerzas de corte por unidad de longitud, actuantes en las aristas mencionadas.
163
, E es el mdulo de elasticidad
(3)
(4)
Por otro lado, a partir de la derivada primera:
(5)
Se deduce la siguiente aproximacin para la derivada segunda cruzada:
(6)
164
Reemplazando las ecuaciones (2), (3), (4) y (6) en (1), y considerando y= z, se obtiene:
(7)
La expresin (7), aplicada sucesivamente en cada uno de los nodos de la discretizacin bidimensional de la
placa, conduce a la siguiente ecuacin matricial:
(8)
Donde la matriz
est
asociada a los cambios en la rigidez debido a la no linealidad cinemtica. La resolucin del problema de
valores y vectores propios dado en (8) permite hallar la menor carga de corte que produce el fenmeno de
inestabilidad por pandeo de la placa en estudio.
4. RESULTADOS OBTENIDOS
Se analiza una placa de acero con relacin de aspecto unitaria,
, t=1 cm,
a=b=300cm. La tensin de corte crtica de inestabilidad de referencia [4], [5], se encuentra dada por:
(9)
siendo ks un coeficiente que vara para las distintas condiciones de borde y en este est caso, adopta valores
de ks=9.34 para la placa con cuatro bordes simplemente apoyados (SA-SA-SA-SA) y ks=14.58 para la
placa con cuatro bordes empotrados (E-E-E-E). De esta manera se obtiene un valor de tensin crtica de
inestabilidad
en el primer caso y de
en el segundo. La Tabla 1
indica los valores obtenidos a partir de la aplicacin del mtodo propuesto en el presente trabajo, para
distintas discretizaciones empleadas:
[kg/cm2]
[%]
rel =
[kg/cm2]
SA-SA-SA-SA
5
10
15
20
25
30
35
40
45
50
181.71
187.39
191.89
193.82
194.79
195.34
195.68
195.90
196.06
196.17
Tabla 1: Valores de
[%]
rel =
E-E-E-E
7.75
4.86
2.58
1.60
1.10
0.83
0.66
0.54
0.46
0.41
251.56
274.88
290.61
297.80
301.50
303.63
304.95
305.83
306.43
306.87
18.18
10.60
5.49
3.15
1.94
1.25
0.82
0.54
0.34
0.20
En los grficos de la Figura 2 es posible observar la evolucin de los errores relativos y la convergencia de
los valores hallados de tensin crtica hacia los valores de referencia mencionados.
165
20
300
275
325
E-E-E-E
250
SA-SA-SA-SA
225
200
E-E-E-E
15
SA-SA-SA-SA
10
5
0
175
0
20
ny
40
60
10
20
30
40
50
ny
SA-SA-SA-SA
E-E-E-E
1.374
1.560
1.540
1.493
1.439
1.383
1.964
2.411
2.686
2.949
3.237
3.573
5. CONCLUSIONES GENERALES
A partir de una adecuada seleccin de las aproximaciones numricas de las derivadas parciales
presentes en la ecuacin diferencial que gobierna el problema tratado, es posible emplear el mtodo de
Diferencias Finitas para la determinacin de cargas crticas de inestabilidad por corte en placas delgadas.
El procedimiento utilizado, brinda excelentes resultados ya que permite obtener adecuados valores de dicha
carga con valores relativamente pequeos de intervalos de discretizacin.
Las diversas condiciones de borde analizadas, muestran diferencias apreciables en los exponentes que
definen la funcin potencial utilizada en una aproximacin por mnimos cuadrados, lo cual indica la
tendencia y evolucin de dichos errores a medida que crece el nmero de intervalos de discretizacin
utilizados.
AGRADECIMIENTOS
Los autores agradecen a la Secretara de Ciencia, Tcnica y Posgrado de la Universidad Nacional de
Cuyo por el financiamiento de los Proyectos de Investigacin Bianuales 2013-2015 nro. 06/B324 y B027.
REFERENCIAS
[1] D.O. BRUSH, B.O. ALMROTH, Buckling of bars, plates and shells, Mc Graw Hill, Kogakusha (1975).
[2] J. K. PAIK AND A. K. THAYAMBALLI, Ultimate limit state design of steel-plated structures, John Wiley
and Sons, 2006.
[3] J. D. HOFFMAN, Numerical methods for engineers and scientists, CRC Press, 2001.
[4] M.M. ALINIA, AND M. DASTFAN, Behaviour of thin steel plate shear walls regarding frame members,
Journal of Construccional Steel Research, 62 (2006), pp.730-738.
[5] R. D. ZIEMIAN, Guide to Stability design criteria for metal structures, 6th Edition, John Wiley & Sons.
INC., (2010).
166