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 ...
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