15.07.2014 Views

Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...

Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...

Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>de</strong> ecuaciones), <strong>de</strong> tal forma <strong>de</strong> hacer más sencilla la comprensión <strong>de</strong> la mayor parte <strong>de</strong>l código<br />

para facilitar la ampliación por parte <strong>de</strong> usuarios no-expertos.<br />

El código que se presenta aquí está <strong>de</strong> alguna forma incompleto con el objetivo <strong>de</strong> que el estudiante<br />

genere las partes que pudieran faltar a manera <strong>de</strong> ejercicio. Esas partes faltantes correspon<strong>de</strong>n<br />

a:<br />

Generación automática <strong>de</strong> datos y verificaciones previas.<br />

Generación <strong>de</strong>l vector <strong>de</strong> términos in<strong>de</strong>pendiente <strong>de</strong>bido a condiciones <strong>de</strong> contorno naturales<br />

o <strong>de</strong>bido al término in<strong>de</strong>pendiente <strong>de</strong> la ecuación diferencial.<br />

Evaluación <strong>de</strong>l flujo o estado tensional en los puntos <strong>de</strong> integración una vez conocidas las<br />

variables nodales y evaluación <strong>de</strong> los valores nodales mediante suavizado.<br />

Evaluación <strong>de</strong>l error aproximado y con él los valores <strong>de</strong>l tamaño <strong>de</strong> la malla para generar<br />

una malla adaptable.<br />

Básicamente se han incluido la posibilidad <strong>de</strong> analizar los siguientes problemas físicos:<br />

Ecuación <strong>de</strong> Laplace, que representa una buena cantidad <strong>de</strong> problemas <strong>de</strong>pendiendo <strong>de</strong>l<br />

significado que se le <strong>de</strong>n a los coeficientes y a sus condiciones <strong>de</strong> bor<strong>de</strong> (flujo irrotacional,<br />

flujo en un medio poroso, torsión sin restricción <strong>de</strong> alabeo ya sea usando la función <strong>de</strong> alabeo<br />

o la <strong>de</strong> tensión, transmisión <strong>de</strong>l calor, etc.)<br />

Elasticidad bidimensional, incluyendo típicamente los casos <strong>de</strong>:<br />

1. Tensión Plana<br />

2. Deformación Plana<br />

3. Sólidos Axilsimétricos<br />

Flexión <strong>de</strong> placas planas incluyendo <strong>de</strong>formaciones <strong>de</strong> corte transversales<br />

Base <strong>de</strong> datos<br />

De hecho una parte imprescindible para enten<strong>de</strong>r un programa es saber como están almacenados<br />

los datos. La presente versión <strong>de</strong>l programa incluye sólo dos elementos finitos, un triángulo <strong>de</strong><br />

seis nodos y un cuadrilátero <strong>de</strong> nueve nodos, ambos cuadráticos (en este caso el programa no está<br />

orientado a tener una variada biblioteca <strong>de</strong> elementos sino que se supone que sólo tendrá estos dos).<br />

En base a ello se <strong>de</strong>finen algunas variables escalares, las que pue<strong>de</strong>n ser parámetros in<strong>de</strong>pendientes<br />

<strong>de</strong>l tipo <strong>de</strong> problema, parámetros <strong>de</strong>pendientes <strong>de</strong>l tipo <strong>de</strong> problema, o variables <strong>de</strong>pendientes <strong>de</strong><br />

la características geométricas <strong>de</strong> la malla o <strong>de</strong> las condiciones <strong>de</strong> bor<strong>de</strong> y <strong>de</strong>l material constitutivo.<br />

Se <strong>de</strong>finen las siguientes:<br />

parámetros<br />

DIMEN indica la dimensión espacial <strong>de</strong>l problema, que en este caso está restringido a 2 (bidimensional),<br />

esta variable y otras, aunque no son imprescindibles, se <strong>de</strong>finen por una cuestión <strong>de</strong><br />

or<strong>de</strong>n y <strong>de</strong> claridad <strong>de</strong> lectura <strong>de</strong>l código.<br />

NODET es el número <strong>de</strong> nudos que tiene el triángulo que es 6.<br />

NODEQ es el número <strong>de</strong> nudos que tiene el cuadrilátero que es 9.<br />

TYPEP indica el tipo <strong>de</strong> problema físico que se intenta resolver<br />

164<br />

0 Ecuación <strong>de</strong> Laplace

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!