Archive for the ‘código’ category

Distribución de tensiones alrededor de un agujero en una placa a tracción.

febrero 14, 2010

Si es cierto que una imagen vale más que mil palabras, una animación ‘.gif’ es una formidable verborrea que sacia hasta los ojos del más curioso. Las animaciones realizadas en este blog para la simulación de fluidos parecen haber gustado a los lectores y, en concreto, han inspirado a los administradores de un blog amigo que han decidido hacernos un encargo. Por eso interrumpimos brevemente las publicaciones sobre maximización de potencia fotovoltaica para presentar los resultados de dicha petición.

Se trata del estudio del estado tensional alrededor de un agujero en una placa sometida a una carga lateral uniforme. La realización de agujeros en una laja sometida a tracción es una necesidad habitual en todo tipo de estructuras (depósitos, elementos de estructuras civiles, fuselajes de aviación, etc.) y es importante conocer ciertos datos como el hecho de que el orificio actúa como concentrador de tensiones y el valor de tensión máxima que se produce (σθ = 3p, siendo p la tensión lateral aplicada)

(más…)

Integración en el tiempo de campos 2D (II)

diciembre 9, 2009

En este artículo vamos a presentar un ejemplo de simulaciones hechas con el código del artículo anterior. (Con animaciones incluidas)

Se trata de un problema conocido y sencillo de dinámica de fluidos: movimiento de un fluido alrededor de una placa plana finita. Para empezar, se ha implementado una función con las derivadas temporales de las variables: densidad, velocidad “Vx”, velocidad “Vy”, presión y entropía. Es decir, las ecuaciones de Navier Stokes (en su variante que utiliza la ecuación de la entropía para gases en lugar de la energía interna) más la ecuación de estado del fluido derivada con respecto al tiempo (considerando que éste tiene comportamiento de gas ideal).  (más…)

Una pequeña herramienta

diciembre 4, 2009

Al subir código al blog “codigospoloestable” para el artículo anterior, me he dado cuenta de una serie de problemas que da el editor de html: Elimina las tabulaciones, interpreta como definición de una etiqueta todo lo que esté entre un < y un > (por lo que desaparecen grandes trozos del código), no representa los “&” si no van seguidos de espacio y la etiqueta <code> sólo funciona por párrafos, por lo que hay que ponerla en cada párrafo para conservar la fuente en todo el código. Además, hay que adaptar el ancho de línea para que se ajuste al ancho de la entrada, para eso hay que escribir tres puntos al final de las líneas de instrucciones y, si se trata de Matlab, un “%” al principio del nuevo renglón sie estamos dividiendo un comentario… con todo, la apariencia de los comentarios queda destrozada (lleno de líneas cortas y largas intercaladas).

Para solucionarlo (más…)

Integración en el tiempo de campos 2D (I)

diciembre 2, 2009

En este artículo se presenta una solución sencilla (aunque limitada) para integrar sistemas de ecuaciones diferenciales con variables definidas en un campo 1D o 2D y un código en Matlab.

No hay cosa más satisfactoria que ver una animación resultado de la simulación de un sistema y comprobar que “tiene buena pinta”, es decir, los resultados son lógicos y verosímiles. Los que nos dedicamos a la ingeniería de sistemas tenemos asociada la idea de simulación a la de integración a lo largo del tiempo y, aunque no necesariamente es así, la gran mayoría de las simulaciones que hacemos incluyen integraciones de sistemas cuyo comportamiento viene descrito por un sistema de ecuaciones diferenciales. Para nosotros no hay mucha diferencia entre simular e integrar usando la famosa “ode45” de Matlab o entre el sistema a simular y la función “odefun” que se le pasa a la ode.

Sin embargo, en algunos de los estudios más apasionantes que existen (problemas de elasticidad, ecuaciones de ondas, ecuación del calor, dinámica de fluídos, …)  se manejan variables definidas en el espacio mediante campos escalares¹. (más…)