CAPITULO 8

8.0 Introducción.

En general los problemas en Ingeniería, son modelados por ecuaciones diferenciales

parciales de segundo orden. La solución analíticas asociada a este tipo de problemas es
complicados y en casos especiales se necesita mucho desarrollo matemático para obtener
una representación analítica, de esta manera hace que este tipo de procedimientos a usar
para la solución de estos problemas este restringido para ciertas personas, ahora con el
desarrollo tecnológico y la creación de las computadoras, cada vez de mayor capacidad y
velocidad, los métodos numéricos se toman cada día mas en cuenta para dar solución de
problemas que son muy complicados. Los métodos numéricos que se trabajan en este
capitulo son el método de las diferencias finitas y el método de el elemento finito. Para el
régimen estacionario y transitorio lineal.

Aplicación 1.

Considere una barra de sección transversal cuadrada como se muestra en la figura (8.1)
determinar la temperatura en el punto p que tiene coordenadas x=0.2, y=0.2 usando el
método de las diferencias finitas, y el método de elemento finito compáralos con la
solución exacta:

La solución de este problema se plantea de la siguiente forma:

Sea
T − T1
θ= (8.1.1)
T2 − T1

El modelo matemático queda definido por la ecuación diferencial:

304
 ∂ 2θ ∂ 2θ
+ =0 (8.1.2)
∂x 2 ∂y 2
y sus condiciones de frontera:
θ (0, y ) = 0 (8.1.3)
θ (w, y ) = 0 (8.1.4)
θ (x, y ) = 0 (8.1.5)
θ (x, L) = 1 (8.1.6)

Usando el método de variables separables se llega al siguiente resultado:

⎡ ⎛ nπ y ⎞ ⎤
⎢ senh ⎜ ⎟
w ⎠ ⎥⎥
n +1
θ=
T − T1 2
= ∑
(− 1) + 1
α
⎢ ⎛ nπ x ⎞
sen⎜ ⎝ (8.1.7)
⎟
T2 − T1 π n =1 n ⎢ ⎝ w ⎠ ⎛ nπL ⎞ ⎥
⎢ senh ⎜ ⎟
⎣ ⎝ w ⎠ ⎥⎦

⎡ 2
⎢ sen⎛⎜ 0.2π ⎞⎟
(
senh 0.2π )
0.3 + 2 sen⎛ 0.6π
(
⎞ senh
0.6π )
0.3
⎤
⎜ ⎟ + ⎥
⎢ 1
2 ⎢ ⎝ 0.3 ⎠ senh(0.π 3 ⎝ 0.3 ⎠ senh(0 .9π ⎥
0.3 0.3
θ (0.2,0.2 ) = ⎥
π ⎢
⎢ 2 sen⎛⎜ 1.0π ⎞⎟
(
senh 1.0π )
0.3 + ...
⎥
⎥
⎢ 5 ⎝ 0.3 ⎠ senh( 1 .5π ⎥
⎣ 0.3 ⎦

2
=
π
(0.60 + 0 − 1.84 x10 −3
)
+ ...... =0.381 (8.1.8)
por lo tanto:
T −0
= 0.381 (8.1.9)
100 − 0

T=38.1 oC (8.1.10)

Ahora usando el método de diferencias finitas:

Como se muestra en la figura(8.1) se han marcado los nodos con 1,2,3,4 que corresponden
a las temperaturas T1, T2, T3 y T4 las ecuaciones para cada uno de los nodos están dadas en
la siguiente forma:
Nodo 1 : -4T1+T2+T3+100=0
Nodo 2 : T1-4T2+T4+100=0
Nodo 3 : T1-4T3+T4=0
Nodo 4 : T2+T3-4T4=0

Estas ecuaciones se pueden resolver usando cualquier método numérico así llegamos a la
solución:

305
 T1=37.5 oC
T2=37.5 oC
T3=12.5 oC
T4=12.5 oC

Por el método de elemento finito se llega al siguiente resultado que se muestra en la tabla
8.1, la solución gráfica es mostrada en la figura (8.1.3):

Tabla 8.1 solución por el método de Elemento finito para la ecuación (8.1.1)

Nodo Temperatura 6 0 12 0
1 0 7 0 13 11.11
2 0 8 50 14 38.889
3 0 9 100 15 11.111
4 0 10 100 16 38.889
5 50 11 0

306
307
Aplicación 2. Se tiene un cilindro Hueco como se muestra en la figura (8.1.4) con radios
r1=1 y r2=2 con temperaturas T1=100 y T2=40, obtener la distribución de temperatura.

La solución analítica a este problema esta dado por la ecuación (8.2.1):

⎛ r ⎞
ln⎜⎜ ⎟⎟
r ⎛ ln(r ) ⎞
T (r ) = T1 − (T1 − T2 ) ⎝ 1 ⎠ = 100 − 60⎜⎜ ⎟⎟ (8.2.1)
⎛ r2 ⎞ ⎝ ln (r )
2 ⎠
ln⎜⎜ ⎟⎟
r
⎝ 1 ⎠

La solución usando el Método de Elemento Finito se da en la tabla (8.2), en la figura

(8.1.5) se muestra la distribución de temperatura, la tabla 8.3 muestra una comparación de
los resultados obtenidos por el método de elemento finito con la solución exacta, en la tabla
8.4 se muestran las coordenadas de cada uno de los nodos y el la figura (8.1.5') se observan
la distribución de los nodos dentro de la placa.

308
309
310
Tabla 8.2 temperatura en nodos usando el Método de elemento finito

Nodo Temperatura 48 100.000 96 46.664 144 91.968

1 40.000 49 100.000 97 45.029 145 84.682
2 40.000 50 100.000 98 52.510 146 92.142
3 40.000 51 100.000 99 48.910 147 91.632
4 40.000 52 100.000 100 47.911 148 83.706
5 40.000 53 100.000 101 56.138 149 91.395
6 40.000 54 100.000 102 47.508 150 83.330
7 40.000 55 100.000 103 48.819 151 89.629
8 40.000 56 100.000 104 46.284 152 82.493
9 40.000 57 100.000 105 52.578 153 89.987
10 40.000 58 100.000 106 49.947 154 91.882
11 40.000 59 100.000 107 50.375 155 91.526
12 40.000 60 100.000 108 50.063 156 91.628
13 40.000 61 100.000 109 51.601 157 92.299
14 40.000 62 100.000 110 47.056 158 83.615
15 40.000 63 100.000 111 44.889 159 94.694
16 40.000 64 100.000 112 87.921 160 83.666
17 40.000 65 100.000 113 80.602 161 91.503
18 40.000 66 100.000 114 90.262 162 89.900
19 40.000 67 100.000 115 80.744 163 72.646
20 40.000 68 100.000 116 90.729 164 89.332
21 40.000 69 100.000 117 92.628 165 81.013
22 40.000 70 100.000 118 80.548 166 89.740
23 40.000 71 100.000 119 92.172 167 80.726
24 40.000 72 100.000 120 85.107 168 60.869
25 40.000 73 100.000 121 91.875 169 63.741
26 40.000 74 100.000 122 83.454 170 54.905
27 40.000 75 100.000 123 91.626 171 54.669
28 40.000 76 100.000 124 92.330 172 85.316
29 40.000 77 100.000 125 91.821 173 78.902
30 40.000 78 100.000 126 90.528 174 80.326
31 40.000 79 100.000 127 87.238 175 90.038
32 40.000 80 100.000 128 81.658 176 53.548
33 40.000 81 48.844 129 88.726 177 70.180
34 40.000 82 50.364 130 79.797 178 61.110
35 40.000 83 94.577 131 89..734 179 72.254
36 40.000 84 91.679 132 90.980 180 71.409
37 40.000 85 91.909 133 81.300 181 72.902
38 40.000 86 46.802 134 90.868 182 63.283
39 40.000 87 58.105 135 80.892 183 60.531
40 40.000 88 46.838 136 90.429 184 70.293
41 100.000 89 47.309 137 74.338 185 59.440
42 100.000 90 48.310 138 91.037 186 57.739
43 100.000 91 50.248 139 93.319 187 55.125
44 100.000 92 49.847 140 92.000 188 70.057
45 100.000 93 47.538 141 92.905 189 52.716
46 100.000 94 49.413 142 93.533 190 61.813
47 100.000 95 48.468 143 84.591 191 53.247

311
192 54.141 213 83.749 234 66.957 255 49.496
193 69.085 214 75.476 235 61.562 256 49.706
194 57.191 215 71.472 236 72.260 257 87.873
195 50.522 216 70.309 237 66.966 258 86.146
196 55.961 217 74.055 238 65.284 259 79.315
197 61.003 218 68.456 239 59.761 260 86.456
198 57.142 219 67.893 240 64.419 261 85.334
199 59.337 220 73.339 241 69.405 262 84.226
200 59.741 221 81.239 242 61.104 263 81.815
201 64.090 222 81.085 243 69.969 264 65.195
202 48.717 223 71.008 244 62.633 265 63.548
203 56.637 224 76.369 245 65.324 266 58.164
204 59.161 225 75.816 246 65.937 267 55.744
205 57.570 226 76.207 247 67.044 268 53.496
206 57.266 227 76.739 248 71.853 269 52.197
207 60.689 228 78.915 249 45.921 270 46.589
208 68.770 229 82.211 250 47.040 271 47.281
209 76.821 230 83.978 251 45.734 272 58.718
210 74.339 231 75.754 252 47.970 273 48.414
211 81.672 232 73.721 253 47.170 274 45.002
212 85.623 233 62.940 254 44.365

Tabla 8.3 comparación de la solución analítica con el método de elemento finito

Numero Solución Usando Solución Error Absoluto Error

De Nodo M.E.F. Analítica Porcentual
81 48.844 48.8307472 0.0132528 0.027140276
82 50.364 50.3873312 0.0233312 0.046303702
83 94.577 94.47297728 0.10402274 0.110108459
89 47.309 47.34259561 0.03359561 0.070962754
94 49.413 49.44943169 0.03643169 0.073674638
253 47.170 47.22904356 0.05904356 0.124972253
207 60.689 60.56969768 0.11930232 0.196967005
265 63.548 63.56141737 0.01341737 0.021106549
86 46.802 47.5528452 0.750845201 1.578970084
234 66.957 66.89144193 0.06555807 0.098006662
128 81.658 81.74673125 0.08873125 0.108544095
226 76.207 76.15962278 0.04737722 0.062207792
220 73.339 73.13695677 0.20204323 0.276253263
212 85.623 85.50176504 0.12123494 0.141792324
228 78.915 78.87622796 0.03877204 0.049155545
238 65.284 65.330866522 0.04686522 0.071735188
219 67.893 67.93765493 0.044654993 0.06572937
244 62.633 62.59156068 0.04143932 0.066205922
261 85.334 85.19427781 0.13972219 0.164004195
142 93.533 93.43117934 0.10182066 0.108979315

312
Tabla 8.4 Coordenadas de los nodos
Nodo X Y 49 0.8910 -0.4540
1 4.0000 0.0000 50 0.9511 -0.3090
2 0.0000 4.0000 51 0.9877 -0.1564
3 3.9508 0.6257 52 -1.0000 0.0000
4 3.8042 1.2361 53 -0.9877 -0.1564
5 3.5640 1.8160 54 -0.9511 -0.3090
6 3.2361 2.3511 55 -0.8910 -0.4540
7 2.8284 2.8284 56 -0.8090 -0.5878
8 2.3511 3.2361 57 -0.7071 -0.7071
9 1.8160 3.5640 58 -0.5878 -0.8090
10 1.2361 3.8042 59 -0.4540 -0.8910
11 0.6257 3.9508 60 -0.3090 -0.9511
12 -4.0000 0.0000 61 -0.1564 -0.9877
13 -0.6257 3.9508 62 0.0000 1.0000
14 -1.2361 3.8042 63 -0.1564 0.9877
15 -1.8160 3.5640 64 -0.3090 0.9511
16 -2.3511 3.2361 65 -0.4540 0.8910
17 -2.8284 2.8284 66 -0.5878 0.8090
18 -3.2361 2.3511 67 -0.7071 0.7071
19 -3.5640 1.8160 68 -0.8090 0.5878
20 -3.8042 1.2361 69 -0.8910 0.4540
21 -3.9508 0.6257 70 -0.9511 0.3090
22 0.0000 -4.0000 71 -0.9877 0.1564
23 -3.9508 -0.6257 72 0.9877 0.1564
24 -3.8042 -1.2361 73 0.9511 0.3090
25 -3.5640 -1.8160 74 0.8910 0.4540
26 -3.2361 -2.3511 75 0.8090 0.5878
27 -2.8284 -2.8284 76 0.7071 0.7071
28 -2.3511 -3.2361 77 0.5878 0.8090
29 -1.8160 -3.5640 78 0.4540 0.8910
30 -1.2361 -3.8042 79 0.3090 0.9511
31 -0.6257 -3.9508 80 0.1564 0.9877
32 0.6257 -3.9508 81 3.2614 -0.4705*10-01
33 1.2361 -3.8042 82 3.0837 0.6324
34 1.8160 -3.5640 83 1.1221 -0.1785
35 2.3511 -3.2361 84 1.1747 -0.3115
36 2.8284 -2.8284 85 1.2040 -0.6765*10-01
37 3.2361 -2.3511 86 3.3588 -0.6439
38 3.5640 -1.8160 87 2.6218 -0.1962
39 3.8042 -1.2361 88 3.0326 1.5540
40 3.9508 -0.62574 89 2.7384 1.9742
41 0.0000 -1.0000 90 1.6770 2.8479
42 1.0000 0.0000 91 1.1539 2.9403
43 0.1564 -0.9877 92 -0.1328 3.1873
44 0.3090 -0.9511 93 -0.6760 3.2927
45 0.4540 -0.8910 94 -1.4190 2.8854
46 0.5878 -0.8090 95 -2.0056 2.6044
47 0.7071 -0.7071 96 -2.5052 2.3396
48 0.8090 -0.5878 97 -2.8684 2.1144

313
98 -2.8288 0.9704 151 -0.2680 1.2420
99 -3.2564 -0.1585 152 -0.5229 1.4098
100 -3.1333 -1.1293 153 -0.5218*10-01 1.2609
101 -2.5476 -1.0643 154 0.9331*10-01 1.2059
102 -2.4826 -2.2654 155 0.2562 1.1602
103 -1.9700 -2.5927 156 0.4479 1.1269
104 -0.9676 -3.3254 157 0.6399 1.0065
105 -0.5552 -2.9427 158 0.6093 1.3308
106 0.3568 -3.1591 159 0.6919 0.8954
107 0.9610 -3.0005 160 0.9082 1.1380
108 1.9286 -2.5195 161 0.8481 0.8708
109 2.2654 -2.0507 162 1.0210 0.7411
110 3.1188 -1.3590 163 1.5209 1.1270
111 3.4111 -1.0781 164 1.1589 0.5422
112 1.1855 -0.5838 165 1.4201 0.6345
113 1.3110 -0.8541 166 1.2265 0.3173
114 0.9761 -0.7872 167 1.5332 0.2935
115 1.1634 -1.0352 168 2.4354 0.4149
116 0.8081 -0.9376 169 2.1115 0.9581
117 0.6444 -0.9976 170 2.5309 1.2864
118 0.9387 -1.2534 171 2.3134 1.6607
119 0.4880 -1.0955 172 1.3886 -0.2244
120 0.5632 -1.2936 173 1.5417 -0.5270
121 0.3049 -1.1666 174 1.5714 -0.3392*10-01
122 0.3308 -1.4308 175 1.2539 0.8259*10-01
123 0.1079 -1.2073 176 2.8121 -0.8245
124 -0.7636*10-01 -1.1942 177 1.9714 -0.3193
125 -0.2725 -1.1786 178 2.1976 -1.1043
126 -0.4972 -1.1424 179 1.6457 -0.9389
127 -0.7687 -1.0963 180 1.9246 0.2652
128 -0.6697 -1.3697 181 1.7280 0.7394
129 -0.8780 -0.9516 182 1.8979 1.3673
130 -1.0030 -1.2361 183 1.7034 1.8154
131 -0.9649 -0.8212 184 1.4972 -1.3042
132 -1.0284 -0.6803 185 1.1696 2.2659
133 -1.2287 -0.9265 186 -0.3326 2.6347
134 -1.1118 -0.5355 187 1.6221 2.3115
135 -1.4031 -0.6785 188 1.2049 1.5942
136 -1.1942 -0.3334 189 0.5009 2.9336
137 -1.7497 -0.4233 190 -0.8873 2.2496
138 -1.2226 -0.1369 191 -2.3183 1.8111
139 -1.1675 0.5737*10-01 192 -0.8287 2.7656
140 -1.1754 0.2474 193 -0.8724 1.8489
141 -1.0956 0.4351 194 -1.7784 2.0105
142 -1.0105 0.5775 195 -3.1225 0.3538
143 -1.2941 0.6097 196 -2.6907 -0.6638
144 -0.9445 0.7470 197 -2.1255 1.2418
145 -1.1156 0.8857 198 -2.6725 -0.3339
146 -0.8051 0.8884 199 -1.6147 -1.9983
147 -0.6489 1.0264 200 -0.9669 -2.3482
148 -0.9442 1.1103 201 -2.1211 -0.8749
149 -0.4795 1.1243 202 -1.4414 -2.9288
150 -0.7232 1.2839 203 0.5022 -2.6758

314
204 1.5887 -2.0191 245 0.6352 -2.1383
205 -0.7732*10-01 -2.6655 246 -0.3587 -2.1732
206 0.8507 -2.5452 247 0.9855 -1.9034
207 1.8872 -1.6196 248 0.6730*10-01 -1.9189
208 1.2692 -1.6174 249 3.2972 1.1393
209 0.7334 -1.5468 250 2.4396 2.3670
210 0.4417 -1.7587 251 -1.1103 3.3274
211 0.1127*10-02 -1.5299 252 -2.8839 1.6703
212 -0.1408 -1.3908 253 -2.9573 -1.6464
213 -0.3786 -1.4058 254 -0.6007 -3.5693
214 -0.3903 -1.7237 255 0.4880 -2.8496
215 -0.8350 -1.7418 256 2.6882 -1.7175
216 -1.2795 -1.5208 257 -1.1469 0.6680
217 -1.5381 -0.9725 258 0.1692 1.3725
218 -1.9585 -0.6723 259 -1.1908 -1.0800
219 -2.0855 -0.2251 260 0.7064 -1.1736
220 -1.8467 0.7241*10-01 261 -1.4052 0.8682*10-01
221 -1.5409 -0.1150 262 -1.2128 -0.7828
222 -1.5070 0.3493 263 1.2274 0.9087
223 -1.8135 0.7303 264 -0.4712 2.1882
224 -1.4014 1.0131 265 -2.2640 -0.5102
225 -1.1090 1.3571 266 1.2466 -2.3142
226 -0.8354 1.5203 267 -2.2476 -1.6322
227 -0.6367 1.5912 268 2.0987 2.0400
228 -0.2696 1.6067 269 1.9386 2.3168
229 0.5216*10-01 1.5143 270 -0.2089 -3.4319
230 0.3111 1.4178 271 2.1049 2.6503
231 0.2848 1.7361 272 -1.3283 2.2255
232 0.7098 1.7003 273 -3.2307 -0.6536
233 0.9653 2.2441 274 -2.8920 -2.0877
234 0.3828 2.1145
235 0.1086 2.4278
236 -0.5581 1.8170
237 -1.2354 1.7537
238 -1.6371 1.5110
239 -2.5307 0.4036*10-01
240 -2.2422 0.3666
241 -1.6361 -1.1867
242 -2.3585 -0.6963
243 -0.9430*10-01 2.0036
244 0.2733 -2.3576

315
Aplicación 3. Se tiene el problema mostrado en la figura (8.1.6) con r1=1, r2=2 y
temperaturas T1=100, T2=0, y sus cara opuestas aisladas, por método del elemento finito
calcular la distribución de temperatura.

La solución analítica para este problema esta dado por:

ln(r )
T (r ) = T1 − (T1 − T2 )
ln(r2 )
ln(r )
T (r ) = 100 − 100
ln(r2 )

usando las coordenadas del punto X=.91845E-16 y Y=1.5 se obtiene que R=1.5 sustituida
en la ecuación anterior se tiene que T(1.5)=41.50374988oC. Que comparado con el valor
obtenido por el método de elemento finito da un error porcentual del 0.393338239. Tabla
8.5 se muestran los valores de las temperaturas en los puntos nodales, en la tabla 8.6 son
las coordenadas de los nodos, y la figura 8.1.7 muestra las isotermas.

Tabla 8.5 Temperatura en los Nodos

Numero Temperatu 11 0 23 100

De Nodo ra 12 100 24 41.667
1 0 13 41.667 25 41.667
2 0 14 100 26 41.667
3 0 15 100 27 41.667
4 0 16 100 28 41.667
5 0 17 100 29 41.667
6 0 18 100 30 41.667
7 0 19 100 31 41.667
8 0 20 100 32 41.667
9 0 21 100 33 41.667
10 0 22 100

316
317
Tabla 8.6 Coordenadas de los nodos

Nodo X Y 12 .6123E-16 1 24 1.5 0

1 2 0 13 .91845E-16 1.5 25 1.4815 .23465
2 .224E-15 2 14 1 0 26 1.4266 .46353
3 1.9754 .31287 15 .98769 .15643 27 1.3365 .68099
4 1.9021 .61803 16 .95106 .30902 28 1.2135 .88168
5 1.7820 .90798 17 .89101 .45399 29 1.0607 1.0607
6 1.6180 1.1756 18 .80902 .58779 30 .88168 1.2135
7 1.4142 1.4142 19 .70711 .70711 31 .68099 1.3365
8 1.1756 1.6180 20 .58779 .80902 32 .46353 1.4266
9 .90798 1.7820 21 .45399 .89101 33 .23465 1.4815
10 .61809 1.9021 22 .30902 .95106
11 .31287 1.9754 23 .15643 .98769

Aplicación 4. Se tiene la siguiente placa semicircular con radio interior ri=1 y el radio
exterior re=4, una temperatura T=100oC en toda la parte superior del radio exterior, en el
resto de la frontera una temperatura de T=0oC, usar el método de elemento finito para
calcular la distribución de temperatura. La figura 8.1.8 muestra gráficamente el problema a
resolver.

La solución de este esta dado por la tabla 8.7,donde se indican los nodos y la temperatura
que tienen , la tabla 8.8 muestra las coordenadas de estos nodos, la figura 8.1.9 muestra
gráficamente la distribución de temperatura en la placa.

Tabla 8.7 Distribución de temperaturas en los nodos

Nodo Temperatura 7 0 14 0
1 50 8 0 15 0
2 0 9 0 16 0
3 0 10 0 17 0
4 0 11 0 18 0
5 0 12 0 19 0
6 0 13 0 20 0

318
21 0 72 4.8638 123 38.961
22 0 73 8.0500 124 43.297
23 0 74 10.499 125 28.213
24 0 75 17.862 126 29.791
25 0 76 10.182 127 22.321
26 0 77 11.632 128 21.025
27 0 78 9.9236 129 17.436
28 0 79 10.846 130 14.723
29 0 80 12.707 131 24.596
30 0 81 21.009 132 20.205
31 0 82 25.198 133 14.265
32 50 83 14.301 134 22.115
33 0 84 15.491 135 31.380
34 0 85 12.449 136 48.446
35 0 86 10.382 137 61.484
36 0 87 10.088 138 52.987
37 0 88 8.5142 139 50.336
38 0 89 6.8566 140 52.503
39 0 90 6.3621 141 76.208
40 0 91 5.0101 142 77.182
41 0 92 12.683 143 71.882
42 100 93 9.1492 144 68.152
43 100 94 2.9346 145 64.428
44 100 95 7.7496 146 58.431
45 100 96 11.498 147 55.836
46 100 97 16.257 148 46.310
47 100 98 31.695 148 53.210
48 100 99 40.384 150 58.206
49 100 100 29.306 151 47.333
50 100 101 54.556 152 34.562
51 100 102 25.702 153 33.640
52 100 103 73.118 154 27.894
53 100 104 74.744 155 37.933
54 100 105 76.192 156 36.648
55 100 106 76.019 157 64.384
56 100 107 85.964 158 27.219
57 100 108 91.015 159 21.433
58 100 109 85.964 160 87.837
59 100 110 84.207 161 82.113
60 100 111 78.567 162 19.066
61 52.567 112 74.456 163 22.365
62 31.401 113 53.688 164 5.3401
63 59.160 114 39.869 165 53.920
64 31.330 115 39.391 166 47.116
65 23.265 116 29.391 167 63.510
66 16.769 117 20.452 168 79.617
67 42.427 118 11.810
68 11.070 119 25.203
69 6.5116 120 33.517
70 30.109 121 24.396
71 2.8638 122 32.897

319
Tabla 8.8 Coordenadas de los nodos
Nodo X Y 52 -.62574 3.9508
1 -4 0 53 -1.2361 3.8042
2 -1 0 54 -1.8160 3.5640
3 -3.7 0 55 -2.3511 3.2361
4 -3.4 0 56 -2.8284 2.8284
5 -3.1 0 57 -3.2361 2.3511
6 -2.8 0 58 -3.5640 1.8160
7 -2.5 0 59 -3.8042 1.2361
8 -2.2 0 60 -3.9508 .62574
9 -1.9 0 61 -3.5456 .45715
10 -1.6 0 62 -3.1207 .42383
11 -1.3 0 63 -3.2098 .85694
12 0 1 64 -2.8218 .53808
13 -.15643 .98769 65 -2.4596 .50268
14 -.30902 .95106 66 -2.1236 .44024
15 -.45399 .89101 67 -2.4161 .99812
16 -.58779 .80902 68 -1.8160 .35229
17 -.70711 .70711 69 -1.5093 .27133
18 -.80902 .58779 70 -2.0587 .82975
19 -.89101 .45399 71 -1.2208 .20416
20 -.95106 .30902 72 -1.1280 .40689
21 -.98769 .15643 73 -1.0684 .61789
22 1 0 74 -.95608 .80067
23 .98769 .15643 75 -1.2986 .81002
24 .95106 .30902 76 -.80673 .90025
25 .89101 .45399 77 -.66050 1.0300
26 .80902 .58779 78 -.51261 1.0629
27 .70711 .70711 79 -.34825 1.1379
28 .58779 .80902 80 -.16110 1.2084
29 .45399 .89101 81 -.42943 1.3275
30 .30902 .95106 82 -.19412 1.4621
31 .15643 .98769 83 .54416E-01 1.2458
32 4 0 84 .26291 1.2444
33 1.3 0 85 .40784 1.1510
34 1.6 0 86 .56212 1.0484
35 1.9 0 87 .71309 .96058
36 2.2 0 88 .85924 .81751
37 2.5 0 89 .98338 .65056
38 2.8 0 90 1.0990 .51098
39 3.1 0 91 1.2001 .34502
40 3.4 0 92 1.3111 .61227
41 3.7 0 93 1.4274 .40696
42 0 4 94 1.2572 .18737
43 3.9508 .62574 95 1.7817 .25163
44 3.8042 1.2361 96 2.0957 .30524
45 3.5640 1.8160 97 2.4293 .35606
46 3.2361 2.3511 98 2.3655 .73968
47 2.8284 2.8224 99 2.7945 .75463
48 2.3511 3.2361 100 3.1511 .38231
49 1.8160 3.5640 101 2.7360 1.1202
50 1.2361 3.8042 102 3.3800 .27419
51 .62574 3.9508 103 3.0659 1.3722

320
104 2.7186 1.7729 137 3.2809 .83709
105 2.2952 2.1810 138 1.8882 1.6093
106 1.7695 2.5033 139 1.3545 1.7913
107 .76960 3.2777 140 .90847 2.0280
108 .48269 3.5448 141 1.1384 2.7796
109 -.44044 3.3252 142 .24793 2.9956
110 -.97069 3.1546 143 -.24467 2.8004
111 -2.3486 2.2325 144 -.72850 2.5986
112 -2.6773 1.8000 145 -1.2080 2.3249
113 -2.2281 1.4272 146 -1.9645 1.7455
114 -1.9209 1.1793 147 -1.5248 1.88899
115 -1.5372 1.3786 148 -.82615 1.8705
116 -16169 1.0109 148 -.49839 2.1536
117 -1.7230 .68044 150 -.8051E-01 2.3494
118 -1.4003 .52888 151 .43424 1.9893
119 -1.1819 1.1133 152 .73343 1.5656
120 -1.0104 1.4390 153 1.0216 1.4308
121 -.70441 1.3199 154 1.1811 1.1889
122 -.58774 1.5695 155 1.3667 1.3989
123 -.30181 1.7818 156 1.6091 1.2357
124 .33308E-01 1.9146 157 .42834 2.5251
125 .81324E-01 1.5377 158 -3.3511 .30031
126 .40411 1.5325 159 2.7868 .37463
127 .65014 1.2830 160 .71531E-01 3.4339
128 .84324 1.1674 161 -1.4777 2.8902
129 1.0135 .97343 162 -.55898 1.2461
130 1.1704 .78174 163 .48752 1.3409
131 1.3716 .98474 164 1.5116 .22064
132 1.5475 .74941 165 3.5888 .44025
133 1.6760 .49597 166 -1.8164 1.4570
134 1.9982 .62445 167 -2.9199 1.2058
135 1.8556 .95691 168 -1.9240 2.5562
136 2.2941 1.2293

321
322
Aplicación 5. Se tiene el mismo problema anterior solo que ahora la temperatura T=1000C
lo tiene en otra parte de la frontera como lo muestra la figura 8.1.10. Usando el método de
elemento finito calcular la distribución de temperatura.

La solución para este caso esta dado en la tabla 8.9, las coordenadas en los nodos se
encuentran en la tabla 8.8. Las Isotermas son mostradas en la figura 8.1.11.

Tabla 8.9 temperatura en los nodos

Nodo Temperatura 27 0 54 0
1 50 28 0 55 0
2 50 29 0 56 0
3 100 30 0 57 0
4 100 31 0 58 0
5 100 32 0 59 0
6 100 33 0 60 0
7 100 34 0 61 46.417
8 100 35 0 62 65.896
9 100 36 0 63 37.050
10 100 37 0 64 63.163
11 100 38 0 65 68.299
12 0 39 0 66 71.169
13 0 40 0 67 43.865
14 0 41 0 68 72.753
15 0 42 0 69 70.387
16 0 43 0 70 49.958
17 0 44 0 71 53.875
18 0 45 0 72 25.276
19 0 46 0 73 18.221
20 0 47 0 74 12.784
21 0 48 0 75 28.778
22 0 49 0 76 8.2320
23 0 50 0 77 5.5299
24 0 51 0 78 3.3669
25 0 52 0 79 2.5435
26 0 53 0 80 2.0454

323
81 4.8238 111 10.523 141 .82346
82 3.6602 112 16.168 142 1.6877
83 1.5386 113 28.327 143 2.9781
84 1.1536 114 34.605 144 4.9695
85 .68892 115 23.317 145 8.4823
86 .42814 116 34.41 146 19.850
87 .31833 117 50.827 147 15.121
88 .19866 118 44.928 148 8.9544
89 .12326 119 19.742 149 5.6273
90 .95103E-01 120 14.019 150 3.4318
91 .64485E-01 121 8.2211 151 2.1118
92 .16751 122 7.4060 152 1.2135
93 .10161 123 4.9768 153 .79942
94 .348E-01 124 3.3088 154 .51820
95 .59837 E-01 125 2.5710 155 .57226
96 .62394 E-01 126 1.6871 156 .40981
97 .56569 E-01 127 .88635 157 1.9312
98 .12432 128 .61082 158 71.558
99 .79093 E-01 129 .37026 159 .41762 E-01
100 .27064 E-01 130 .23992 160 1.1043
101 .11605 131 .34699 161 4.7242
102 .12466 E-01 132 .22091 162 5.7152
103 .77351 E-01 133 .12817 163 1.1266
104 .14861 134 .13991 164 .53155 E-01
105 .26124 135 .24952 165 .10275 E-01
106 .47262 136 .21372 166 25.507
107 .75435 137 .40541 E-01 167 29.027
108 .62635 138 .39355 168 7.5620
109 1.8080 138 .74238
110 2.9097 140 1.2837

324
325
Aplicación 6. Se tiene el siguiente problema una placa cuadrada de 4X4 con un semicirculo
con una temperatura T=100 en la parte semicircular y cero en el resto de la frontera como
se muestra en la figura 8.1.12 . Usar el método de elemento finito para calcular la
distribucion de temperatura.

En la tabala 8.10 se muestra la temperatura en los nodos, la tabla 8.10 contiene sus
respectivas coordenadas para cada uno de los nodos. En la figura 8.1.13 se muestran la
solución grafica, en la figura 8.113´ muestra la distribución de los nodos dentro de la placa.

Tabla 8.10 Muestra los nodos con sus temperaturas

Nodo Temperatura 29 60 58 0 87 49.413
1 0 30 60 59 0 88 47.740
2 30 31 60 60 0 89 45.754
3 0 32 0 61 0 90 44.449
4 0 33 0 62 0 91 43.401
5 0 34 0 63 0 92 41.478
6 0 35 0 64 0 93 37.913
7 0 36 0 65 0 94 21.239
8 0 37 0 66 0 95 5.9635
9 0 38 0 67 0 96 6.6686
10 0 39 0 68 0 97 4.5819
11 0 40 0 69 0 98 3.1429
12 60 41 0 70 0 99 7.5054
13 60 42 0 71 4.4957 100 11.183
14 60 43 0 72 8.9409 101 11.099
15 60 44 0 73 9.8991 102 12.326
16 60 45 0 74 12.382 103 14.789
17 60 46 0 75 11.226 104 26.643
18 60 47 0 76 35.684 105 20.850
19 60 48 0 77 42.652 106 22.549
20 60 49 0 78 45.326 107 7.8828
21 60 50 0 79 46.173 108 5.3116
22 30 51 0 80 47.034 109 3.5621
23 60 52 0 81 47.874 110 1.0082
24 60 53 0 82 48.607 111 13.168
25 60 54 0 83 49.412 112 1.7777
26 60 55 0 84 49.933 113 3.5848
27 60 56 0 85 50.181 114 2.4754
28 60 57 0 86 50.076 115 3.0374

326
116 3.4299 170 24.890
117 3.2135 171 22.332
118 2.1430 172 27.707
119 0.91107 173 30.559
120 2.9685 174 32.082
121 3.9447 175 32.547
122 4.6804 176 32.618
123 7.0326 177 33.438
124 11.468 178 24.644
125 7.4893 179 69.893
126 17.956 180 12.871
127 26.178 181 8.5468
128 10.926 182 15.781
129 20.450 183 12.310
130 21.765 184 18.929
131 16.381 185 22.800
132 24.606 186 7.6204
133 34.646 187 2.9953
134 34.971 188 2.5668
135 36.056 189 5.9869
136 37.625 190 2.0825
137 39.414 191 1.1178
138 40.686 192 6.8524
139 41.067 193 6.8641
140 40.848 194 33.642
141 40.101 195 34.843
142 37.458 196 19.417
143 32.857 197 9.8119
144 30.828 198 25.638
145 41.621 199 19.028
146 26.463 200 14.071
147 11.639 201 7.8005
148 14.071
149 5.4131
150 8.6877
151 16.484
152 10.704
153 5.9947
154 12.730
155 4.6387
156 6.0891
157 7.9191
158 7.0944
159 12.296
160 11.882
161 7.9596
162 28.463
163 24.458
164 18.206
165 25.097
166 18.1054
167 10.554
168 17.185
169 14.452

327
Tabla 8.10´ coordenadas de los nodos.

Nodo X Y 50 2 3.2 100 1.6720 .88392

1 -2 0 51 2 3.6 101 1.6393 1.2081
2 -1 0 52 -2 4 102 1.4889 1.6084
3 -1.9 0 53 1.6 4 103 1.5656 .60535
4 -1.8 0 54 1.2 4 104 1.3131 .68963
5 -1.7 0 55 .8 4 105 1.3645 1.1002
6 -1.6 0 56 .4 4 106 1.2145 1.3314
7 -1.5 0 57 0 4 107 1.5311 2.0845
8 -1.4 0 58 -.4 4 108 1.5380 2.5237
9 -1.3 0 59 -.8 4 109 1.4977 3.0018
10 -1.2 0 60 -1.2 4 110 1.5805 3.6288
11 -1.1 0 61 -1.6 4 111 1.1202 2.1833
12 0 1 62 -2 3.6 112 1.1943 3.6452
13 1 .98769 63 -2 3.2 113 1.1900 3.3109
14 .98769 .95106 64 -2 2.8 114 .80676 3.6411
15 .95106 .89101 65 -2 2.4 115 .38772 3.6297
16 .89101 .80902 66 -2 2 116 -.68175E-01 3.5974
17 .80902 .70711 67 -2 1.6 117 -.56370 3.5855
18 .70711 .58779 68 -2 1.2 118 -.99869 3.6642
19 .58779 .45399 69 -2 .8 119 -1.5549 3.6624
20 -.95106 .30902 70 -2 .4 120 -1.3805 3.2846
21 -.98769 .15643 71 -1.8280 .41842 121 -1.5526 2.8003
22 1 0 72 -1.5719 .2593 122 -1.6408 2.3992
23 .98769 .15643 73 -1.4534 .18207 123 -1.6207 1.9753
24 .95106 .30902 74 -13487 .16848 124 -1.5590 1.5182
25 .89101 .45399 75 -1.2878 .12055 125 -1.7592 1.2181
26 .80902 .58779 76 -1.09998 .16588 126 -1.4163 .37142
27 .70711 .70711 77 -1.0894 .32747 127 -1.2613 .32459
28 .58779 .80902 78 -1.0454 .51037 128 -1.6246 .44471
29 .45399 .89101 79 -.96384 .69759 129 -1.2162 .16598
30 .30902 .95106 80 -.84993 .85311 130 -1.3959 .56760
31 .15643 .98769 81 -.68987 .99831 131 -1.5273 .97823
32 2 0 82 -.52308 1.0971 132 -1.3265 .91143
33 1.9 0 83 -.33682 1.1628 133 -1.1306 .81049
34 1.8 0 84 -.14684 1.1969 134 -1.0093 1.0416
35 1.7 0 85 .44678E-01 1.2009 135 -.82804 1.2295
36 1.6 0 86 .22558 1.1800 136 -.60336 1.3569
37 1.5 0 87 .39439 1.1433 137 -.35519 1.4203
38 1.4 0 88 .56708 1.0942 138 -.11900 1.4359
39 1.3 0 89 .72834 1.0211 139 .11746 1.4253
40 1.2 0 90 .86615 .90871 140 .32572 1.3895
41 1.1 0 91 .98120 .76003 141 .50161 1.3411
42 2 4 92 1.1038 .38104 142 .67846 1.3115
43 2 .4 93 1.1039 .19794 143 1.0232 1.1101
44 2 .8 94 1.2455 .20293 144 1.1664 .93978
45 2 1.2 95 1.4242 .10525 145 1.0728 .58979
46 2 1.6 96 1.5112 .15368 146 1.2890 .40510
47 2 2 97 1.6273 .14783 147 1.3613 .16207
48 2 2.4 98 1.7589 .18420 148 1.4832 .33455
49 2 2.8 99 1.6666 .32094 149 1.1927 3.0049

328
150 1.1597 2.5951 168 -1.1268 1.8567 186 -1.5485 .18838
151 .75659 2.2392 169 -1.5244 1.2510 187 1.7071 .13340
152 .82789 2.6526 170 -1.0096 1.5210 188 1.8618 .28246
153 .92361 3.1118 171 -1.2521 1.2894 189 1.7942 .51331
154 .44920 2.6579 172 -.69814 1.6681 190 1.5670 3.2854
155 .83315 3.3360 173 -.37264 1.7157 191 -1.2392 3.7652
156 .41229 3.2802 174 -.70073E-01 1.7170 192 -1.7266 1.5923
157 -.11009 3.1376 175 .23733 1.6767 193 -1.7892 .79626
158 -.77669 3.0666 176 .46594 1.5997 194 -1.2010 .54846
159 -.25368 2.7346 177 .59471 1.5102 195 .85592 1.2439
160 -.64999 2.6503 178 .13995E-01 2.0230 196 1.4477 .90480
161 -1.1791 2.6750 179 -.29721 2.4016 197 .14953 2.9536
162 .78922 1.5668 180 .72235E-01 2.7138 198 .98125 1.5067
163 .68894 1.8285 181 .48644 3.0095 199 1.0696 1.8011
164 .40923 2.2960 182 -.57238 2.3880 200 -1.3896 1.6599
165 .34628 1.9549 183 -.98595 2.3984 201 -1.3867 2.4070
166 .54195E-01 2.3589 184 -.75385 2.0839
167 -1.3402 2.1467 185 -.36714 2.0574

329
330
331
Aplicación 7. Se tiene el mismo problema anterior solo que ahora se tiene aislada la placa y
temperatura en la parte semicircular T=60, y T=0 en los nodos 1,2,3,4,5,6,12,17,18,19,20,21 sus
dimensiones son de 8X4, el radio de la parte semicircular es r=1 . Como lo muestra la figura 8.1.14,
calcular por medio de elemento finito la distribucion de temperatura.

La solución a esta aplicación se muestra en la tabla 8.11 la cual contiene la temperatura en los
nodos, la tabla 8.12 contiene la coordenadas de los nodos. La figura 8.1.15 muestra la distribución
de temperatura en la placa. En la figura 8.1.15’ se da la distribuciòn de los nodos en la placa.

Tabla 8.11 temperatura en los nodos

Nodo Temperatura 35 32.5510 70 20.7990
1 0.0000 36 33.2790 71 23.8900
2 30.0000 37 33.5620 72 26.8350
3 0.0000 38 33.3270 73 28.7900
4 0.0000 39 32.6330 74 30.2930
5 0.0000 40 31.5660 75 33.0130
6 0.0000 41 30.2370 76 33.7450
7 60.0000 42 28.809 77 33.9450
8 60.0000 43 27.515 78 33.6560
9 60.0000 44 26.260 79 32.8040
10 60.0000 45 25.324 80 31.4910
11 60.0000 46 24.799 81 28.9520
12 30.0000 47 23.884 82 25.4960
13 60.0000 48 21.000 83 23.0140
14 60.0000 49 16.168 84 20.8760
15 60.0000 50 9.0358 85 26.4860
16 0.0000 51 9.7602 86 26.4330
17 0.0000 52 17.811 87 32.4820
18 0.0000 53 11.3180 88 40.5860
19 0.0000 54 25.1670 89 44.1420
20 0.0000 55 14.5360 90 46.3950
21 0.0000 56 25.7200 91 47.0570
22 24.6900 57 35.8120 92 46.9120
23 8.8840 58 37.5940 93 41.1620
24 16.2880 59 48.0850 94 27.1310
25 21.1090 60 51.4630 95 21.2150
26 23.5810 61 53.0470 96 28.2310
27 24.6750 62 53.6620 97 30.3770
28 24.9160 63 46.8090 98 29.7240
29 25.3350 64 46.7590 99 34.7900
30 26.1550 65 31.5630 100 34.2500
31 27.3720 66 17.4000 101 34.9880
32 28.7850 67 17.6000 102 34.6720
33 30.1940 68 9.9791 103 34.7930
34 31.4790 69 16.5260 104 36.0560

332
105 36.3770 115 41.3540 125 27.6800
106 36.2860 116 40.2460 126 31.7830
107 33.7920 117 34.6800 127 30.2600
108 35.9880 118 36.9510 128 24.9000
109 35.8140 119 38.4080 129 24.3370
110 31.2560 120 38.2580 130 32.6240
111 39.1260 121 37.6670 131 22.7480
112 41.5220 122 53.1590 132 26.6550
113 41.0280 123 14.1460
114 41.9860 124 25.2510

Tabla 8.12 coordenadas de los nodos.

Nodo X Y 38 -0.4000 4.0000
1 -4.0000 0.0000 39 -0.8000 4.0000
2 -1.0000 0.0000 40 -1.2000 4.0000
3 -3.4000 0.0000 41 -1.6000 4.0000
4 -2.8000 0.0000 42 -2.0000 4.0000
5 -2.2000 0.0000 43 -2.4000 4.0000
6 -1.6000 0.0000 44 -2.8000 4.0000
7 0.0000 1.0000 45 -3.2000 4.0000
8 -0.3090 0.9511 46 -3.6000 4.0000
9 -0.5877 0.8090 47 -4.0000 3.2000
10 -0.8090 0.5878 48 -4.0000 2.4000
11 -0.9510 0.3090 49 -4.0000 1.6000
12 1.0000 0.0000 50 -4.0000 0.8000
13 0.9510 0.3090 51 -3.3155 0.7927
14 0.8090 0.5878 52 -3.2247 1.6318
15 0.5877 0.8090 53 -2.6565 0.7162
16 0.3090 0.9511 54 -1.4764 0.4878
17 4.0000 0.0000 55 -2.0356 0.5971
18 3.4000 0.0000 56 -1.8285 1.1618
19 2.8000 0.0000 57 -1.3333 0.9419
20 2.2000 0.0000 58 -1.1997 1.3285
21 1.6000 0.0000 59 -0.4996 1.3843
22 4.0000 4.0000 60 -0.1828 1.3277
23 4.0000 0.8000 61 0.1498 1.2547
24 4.0000 1.6000 62 0.4230 1.1425
25 4.0000 2.4000 63 0.7830 1.2097
26 4.0000 3.2000 64 1.0051 0.8862
27 -4.0000 4.0000 65 1.4118 0.6992
28 3.6000 4.0000 66 2.0890 0.8358
29 3.2000 4.0000 67 3.0137 1.5291
30 2.8000 4.0000 68 3.3898 0.8566
31 2.4000 4.0000 69 3.4467 1.5512
32 2.0000 4.0000 70 3.3834 2.2191
33 1.6000 4.0000 71 3.2529 2.9458
34 1.2000 4.0000 72 2.5013 3.3922
35 0.8000 4.0000 73 1.9980 3.5125
36 0.4000 4.0000 74 1.6065 3.4443
37 0.2220 4.0000 75 0.7559 3.4489

333
76 0.3620 3.4795 105 -0.3723*10-01 2.7208
77 -0.3486*10-01 3.5143 106 -0.3966 2.6671
78 -0.4458 3.5006 107 -0.8926 2.9607
79 -0.8653 3.4757 108 -0.7273 2.5311
80 -1.2921 3.4658 109 -1.0429 2.1449
81 -1.9208 3.3557 110 -1.4895 2.8343
82 -2.9277 3.2707 111 -0.6541 2.0710
83 -3.1693 2.5931 112 -0.6353 1.8142
84 -2.4749 1.4749 113 -0.3234 2.0226
85 -2.5340 3.0865 114 0.6297*10-01 1.9689
86 -2.2174 2.3379 115 0.4395 1.9068
87 -1.4856 1.6091 116 0.8230 1.7070
88 -0.8579 1.6792 117 1.3773 1.4504
89 -0.3699 1.7215 118 0.5288 2.4557
90 -0.4091*10-01 1.6239 119 0.5336*10-01 2.3805
91 0.2876 1.5307 120 -0.3085 2.3578
92 0.5869 1.3771 121 -0.5767 2.3200
93 1.0063 1.3170 122 0.7173 0.9631
94 1.9236 1.5977 123 2.7613 1.0173
95 2.8820 1.9658 124 2.8199 2.9633
96 2.0868 2.9930 125 2.2986 3.6126
97 1.6441 2.8777 126 1.1836 3.4229
98 1.7641 2.3176 127 -1.6097 3.5910
99 1.1459 2.1646 128 -3.2074 3.4033
100 0.7016 2.9316 129 -3.5031 3.3199
101 -0.5556*10-01 3.0770 130 1.1671 2.8454
102 -0.4847 3.0397 131 2.4639 1.77870
103 0.3331 3.0607 132 2.2670 2.5490
104 0.2749 2.7457

334
335
Aplicación 8. Se tiene el siguiente problema una placa circular con un circulo de radio R=8, dentro
de este otro circulo de radio r=1, con distancia entre centros de d=3, este se encuentra sobre el
diámetro del circulo mayor, con temperatura T1=30oC y T2=100oC como se muestra en la figura
8.1.16, por el método de elemento finito calcular la distribución de temperatura.

La solución a esta aplicación esta dada en la tabla 8.13 muestra la distribución de temperatura en
los nodos, mientras que la tabla 8.14 muestra las coordenadas de los nodos. La figura 8.1.17
muestra la distribución de temperatura en la placa.

Tabla 8.13 temperatura en los nodos.

Nodo Temperatura 22 30 44 77.868

1 100 23 30 45 93.373
2 100 24 30 46 94.028
3 100 25 30 47 84.76
4 100 26 30 48 93.536
5 100 27 30 48 88.573
6 100 28 33.342 50 74.656
7 100 29 30 51 480706
8 100 30 30 52 46.29
9 100 31 30 53 47.638
10 100 32 30 54 51.028
11 100 33 30 55 48.664
12 100 34 34.74 56 43.295
13 100 35 30 57 48.663
14 100 36 32.599 58 66.78
15 100 37 50.596 59 69.112
16 100 38 69.811 60 64.922
17 30 39 49.983 61 80.442
18 30 40 46.678 62 62.116
19 30 41 69.774 63 71.851
20 30 42 86.496 64 67.803
21 30 43 89.712 65 84.133

336
 66 82.603 70 49.545 74 89.722
67 91.644 71 65.283 75 95.583
68 94 72 71.559
69 59.804 73 38.863

Tabla 8.14 coordenadas de los nodos

Nodo X Y 39 -1.8448 -.30016
1 4 0 40 1.3889 -.75409
2 -4 8 41 -.81192 2.8436
3 3.391 3.0615 42 -2.7397 4.5934
4 1.6569 5.6569 43 -5.2580 4.4242
5 -.93853 7.391 44 -3.772 2.4640
6 -12 0 45 -8.5338 1.5156
7 -7.0615 7.391 46 -8.6409 -1.8210
8 -9.6569 5.6569 47 -6.1041 -.33595
9 -11.391 3.0615 48 -7.3591 -3.8089
10 -4 -8 48 -4.1355 -4.7528
11 -11.391 -3.0615 50 -1.3561 -3.1414
12 -9.6569 -5.6569 51 -1.7711 .36268
13 -7.0615 -7.3910 52 -1.4107 .90848
14 -.93853 -7.3910 53 -1.1918 1.2896
15 1.6569 -5.6569 54 -.73025 1.7956
16 3.391 -3.0615 55 -1.0165 -1.5649
17 0 1 56 -.90776 -1.1854
18 -1 0 57 -1.4189 -1.0252
19 -.38268 .92388 58 -3.1146 .76144
20 -.70711 .70711 59 -3.4485 -.39080
21 -.92388 .38268 60 -2.5928 1.4501
22 1 .0 61 -4.9104 1.5238
23 .92388 .38268 62 -2.0869 1.7884
24 ..70711 .70711 63 -2.3632 2.6371
25 38268 .92388 64 -27998 -1.7794
26 0 -1 65 -5.0996 -2.6693
27 .25882 -.96593 66 .22231 3.9856
28 .5 -.86603 67 -7.0588 3.2484
29 .70711 -.70711 68 -6.1517 -5.2502
30 .86603 -.5 69 2.0229 -.69122
31 .96593 -.25882 70 1.2776 -1.0711
32 -.96593 -.25882 71 1.5862 -1.7390
33 -.86603 -.5 72 .55773 -2.8759
34 -.70711 -.70711 73 -.99889 -.85921
35 -.5 -.86603 74 -1.1520 5.3233
36 -.25882 -.96593 75 -9.7081 -.77718E-01
37 1.4133 1.0683
38 1.2567 2.2789

337
336
8.2 APLICACIÓN DE METODOS NUMERICOS A ESTADO TRANSITORIO

APLICACIÓN
Una placa sólida rectangular de Lxl como se muestra en la figura (8.2.1), con una
temperatura uniforme de To=100 oC , tiene sus cuatro cantos o superficies laterales
térmicamente aisladas, la temperatura de una de sus caras expuesta se eleva a T=100 oC
mientras que la otra se mantiene constante a T(0,t)=0oC obténgase la distribución de
temperatura en función del tiempo.

Como las superficies están aisladas, el flujo de calor sólo se efectúa en una dirección, así
que el modelo matemático que gobierna este fenómeno esta dado de la siguiente manera:

∂ 2T ∂T
α2 = (8.2.1)
∂x 2 ∂t

sujeto a las siguientes condiciones de frontera:

T(x=0,t)=To=0, T(x=L,t)=TL=100 (8.2.2)

Donde:
X: esta en cm
α: difusividad térmica cm2/seg=10
T: temperatura en oC
t: tiempo en segundos:

Usando el método de variables separables se llega al siguiente resultado:

Este problema es no Homogéneo por lo que se encuentra una solución para el caso
estacionario:
(T − T0 )
Te ( x, t ) = T0 + L x (8.2.3)
L
se define una nueva variable Z de la siguiente manera:
(T − T0 )
Z = T − T0 − 1 x (8.2.4)
L
sustituyendo en la ecuación diferencial (8.2.1) se obtiene la nueva ecuación.

338
 ∂ 2 Z ∂Z
2
α = (8.2.5)
∂x 2 ∂t
Ahora las condiciones de frontera son:

X=0, Z=0
X=L, Z=0 (8.2.6)

En esta aplicación la condición final es la que temperatura sea la inicial To=0, así que la Z
se convierte en:
T
Z =− 1 x (8.2.7)
L
ahora la solución del problema (8.2.5) esta dado por:
∞ ⎛ nπx ⎞ − n 2π 2α t / L2
Z = ∑ A sen ⎜ ⎟e
n (8.2.8)
n =1 ⎝ ⎠
L
donde:
L
⎛ nπx ⎞ TL
An = ∫ f o (x ) sen⎜ ⎟dx donde f o (x ) = − x (8.2.9)
0 ⎝ L ⎠ L
realizando las operaciones respectivas se obtiene:
T n
An = L (− 1) (8.2.10)
nπ
así sustituyendo este resultado en (8.28) se tiene:
∞ 2
( )n ⎛ nπx ⎞ − n 2π 2α t / L2
Z = ∑ − 1 TL sen ⎜ ⎟e (8.2.11)
n =1 nπ ⎝ L ⎠
regresando a la variable original utilizando (8.2.4)

T x
∞ 2 n ⎛ nπx ⎞ − n 2π 2α t / L2
( )
= + ∑ − 1 TL sen ⎜ ⎟e (8.2.12)
TL L n =1 nπ ⎝ L ⎠
⎛ x ∞ 2 n
( ) ⎛ nπx ⎞ − n 2π 2α t / L2 ⎞
T = TL ⎜ + ∑ − 1 TL sen ⎜ ⎟e ⎟ (8.2.13)
⎝ l n =1 nπ ⎝ ⎠
L ⎠
sustituyendo valores y tomando solo los dos primeros términos:
⎛ 10 2 ⎛ πx ⎞ 2 ⎞
T (x, t ) = 100⎜ − sen⎜ ⎟e −π t / 10 ⎟ (8.2.14)
⎝ x π ⎝ 10 ⎠ ⎠
ahora dividiendo la placa en 10 partes de longitud de un centímetro y evaluando para t=2
se tiene la tabla(8.2.1):

Tabla 8.2.1 Temperatura en las coordenadas T(x,y)

339
T(1,2) 7.26725 T(5,2) 41.15664471 T(9,2) 87.26725293
T(2,2) 14.80200 T(6,2) 51.58946 T(10,2) 100.00
T(3,2) 22.88455 T(7,2) 62.8455752
T(4,2) 31.5844 T(8,2) 74.80200618

Se resolverá el mismo problema usando el método explícito, con

a) Δt = 0.02
b) Δt = 0.055 (8.2.15)

Solución para el inciso (a):

La aproximación por diferencias finitas del método explícito para la ecuación (8.2.1) es:

(
Ti ( n +1) = Ti ( n ) + γ Ti −( n1) − 2Ti ( n ) + Ti +( n1) ) (8.2.16)
donde
α Δt
γ = (8.2.17)
Δx 2
es el criterio de estabilidad que debe satisfacer Δt para que γ ≤ 0.5 , es
2
0.5(Δx )
Δt ≤ = 0.05 seg (8.2.18)
α
los resultados computacionales con Δt = 0.02 ,( γ = 0.2) aparecen en la tabla (8.2.2)

Tabla (8.2.2) Cálculo de la distribución de temperatura usando el Método explícito

n t(s) i=0 1 2 3 4 5 6 7 8 9 10
1 0.02 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 20.0 100.0
2 0.04 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 4.0 32.0 100.0
3 0.06 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.8 8.8 40.0 100.0
4 0.08 0.0 0.0 0.0 0.0 0.0 0.0 0.2 2.2 13.4 45.8 100.0
5 0.10 0.0 0.0 0.0 0.0 0.0 0.0 0.5 4.1 17.7 50.1 100.0
6 0.12 0.0 0.0 0.0 0.0 0.0 0.1 1.1 6.1 21.4 53.6 100.0
7 0.14 0.0 0.0 0.0 0.0 0.0 0.3 1.9 8.2 24.8 56.5 100.0
8 0.16 0.0 0.0 0.0 0.0 0.1 0.6 2.9 10.2 27.8 58.8 100.0
9 0.18 0.0 0.0 0.0 0.0 0.2 0.9 3.9 12.3 30.5 60.9 100.0
10 0.20 0.0 0.0 0.0 0.0 0.3 1.4 5.0 14.2 32.9 62.6 100.0
20 0.40 0.0 0.1 0.5 1.4 3.5 8.0 16.2 29.5 48.5 72.7 100.0
30 0.60 0.0 0.8 2.1 4.4 8.5 15.2 25.2 39.0 56.7 77.5 100.0
40 0.80 0.0 1.9 4.4 8.1 13.5 21.4 32.0 45.6 61.9 80.4 100.0
50 1.00 0.0 3.1 6.7 11.5 18.0 26.5 37.3 50.4 65.6 82.4 100.0
60 1.20 0.0 4.2 8.9 14.7 21.8 30.7 41.5 54.2 68.4 83.9 100.0
70 1.40 0.0 5.2 10.8 17.3 25.0 34.2 44.9 57.1 70.6 85.0 100.0
80 1.60 0.0 6.0 12.4 19.5 27.7 37.0 47.6 59.4 72.3 86.0 100.0
90 1.80 0.0 6.7 13.8 21.4 29.9 39.3 49.9 61.4 73.7 86.7 100.0
100 2.00 0.0 7.3 14.9 22.9 31.7 41.3 51.7 62.9 74.9 87.3 100.0

340
Para el inciso (b). En la tabla (8.2.3) se muestran los resultados con Δt = 0.055 ( γ = 0.55) ,
que es mayor que el criterio de estabilidad. Las oscilaciones de magnitud creciente tanto en
el espacio como en el tiempo son síntomas de inestabilidad:

Tabla (8.2.3) Cálculo de la distribución de temperatura usando el Método Explícito

n t(s) i=0 1 2 3 4 5 6 7 8 9 10
1 0.055 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 55.0 100.0
2 0.11 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 30.3 49.5 100.0
3 0.17 0.0 0.0 0.0 0.0 0.0 0.0 0.0 16.6 24.2 66.7 100.0
4 0.22 0.0 0.0 0.0 0.0 0.0 0.0 9.2 11.6 43.4 61.6 100.0
5 0.28 0.0 0.0 0.0 0.0 0.0 5.0 5.5 27.7 36.0 72.7 100.0
6 0.33 0.0 0.0 0.0 0.0 2.8 2.5 17.5 20.0 51.7 67.5 100.0
7 0.38 0.0 0.0 0.0 1.5 1.1 10.9 10.7 36.0 43.0 76.7 100.0
8 0.44 0.0 0.0 0.8 0.5 6.7 5.4 24.7 25.9 57.7 71.0 100.0
9 0.50 0.0 0.5 0.2 4.1 2.5 16.8 14.7 42.7 47.5 79.6 100.0
10 0.55 0.0 0.0 2.5 1.1 11.2 7.8 31.2 30.0 62.5 73.2 100.0
15 0.82 0.0 4.2 0.1 15.9 4.4 34.7 20.0 59.5 53.0 86.4 100.0
20 1.10 0.0 -2.2 19.7 -3.7 41.8 4.9 64.6 30.7 84.3 74.3 100.0
25 1.38 0.0 18.5 -15.0 53.7 -19.1 82.1 -2.2 98.3 40.0 101.4 100.0
30 1.65 0.0 -21.6 66.4 -54.2 117.4 -56.9 140.2 -18.8 131.0 55.6 100.0
35 1.92 0.0 63.8 -93.5 172.3 - 228.5 - 216.6 -37.8 146.5 100.0
145.8 128.5
40 2.20 0.0 - 232.0 - - 385.2 - 407.7
296.2 -28.1 100.0
105.6 273.8 328.1 237.9
45 2.47 0.0 234.0 - 617.8 - 778.9 - 661.7 - 316.5 100.0
412.6 662.0 664.4 356.6

Se usara el método implícito para resolver el mismo problema con γ = 2 :

Las ecuaciones en diferencias son:
− γ Ti −( n1+1) + (2γ + 1) Ti ( n +1) − γ Ti +( n1+1) = Ti ( n ) i=1,2,..,n-1, (8.2.19)

donde i=10, T0=0 y T10=100, de esta parte se tiene un sistema de ecuaciones simultaneas.
La solución a este se muestra en la tabla (8.2.3), claramente se nota que a pesar que a pesar
que se viola el criterio de estabilidad el método implícito es estable.

341
Tabla 8.2.3 Cálculo de la distribución de temperatura usando el Método Implícito

n t(s) i=0 1 2 3 4 5 6 7 8 9 10
0 .00 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 100.0
1 .2 0.0 0.2 0.4 0.8 1.6 3.1 6.3 12.5 25.0 50.0 100.0
2 .4 0.0 0.6 1.3 2.5 4.7 8.3 14.6 25.0 41.7 66.7 100.0
3 .6 0.0 1.2 2.8 5.0 8.6 14.1 22.4 34.7 51.8 74.1 100.0
4 .8 0.0 2.1 4.5 7.9 12.7 19.5 29.1 41.9 58.4 78.2 100.0
5 1 0.0 3.0 6.4 10.8 16.6 24.3 34.5 47.4 63.0 80.8 100.0
6 1.2 0.0 3.9 8.2 13.5 20.1 28.5 38.9 51.5 66.3 82.7 100.0
7 1.4 0.0 4.7 9.9 15.9 23.2 32.0 42.5 54.8 68.8 84.1 100.0
8 1.6 0.0 5.5 11.4 18.1 25.8 34.9 45.5 57.5 70.8 85.1 100.0
9 1.8 0.0 6.2 12.8 20.0 28.1 37.4 47.9 59.6 72.4 86.0 100.0
10 2.0 0.0 6.8 13.9 21.6 30.0 39.4 49.9 61.4 73.7 86.7 100.0
Se usara el método de elemento finito para resolver el problema usando, Δt = 0.02 ,
Δt = 0.55 , Δt = 0.6 la tabla 8.2.4 muestra los resultados:

Tabla 8.2.4 Transitoria elemento finito

Nodo Tiempo=2se Tiempo=2.4 Tiempo=6 Teórica Teórica Teórica

g 7seg n=300 t=2 t=2.47 t=6
1 0 0 0 0 0 0
2 3.1557 3.9772 7.2562 7.2742 8.631 9.9473
3 6.6742 8.2965 14.679 14.813 17.396 19.9
4 10.931 13.304 22.433 22.857 26.415 29.862
5 16.329 19.355 30.687 31.596 35.786 39.838
6 23.303 26.823 39.604 41.157 45.568 49.829
7 32.331 36.108 49.352 51.583 55.785 59.838
8 43.935 47.638 60.100 62.834 66.414 69.862
9 58.678 61.841 72.017 74.791 77.394 79.9
10 77.158 79.080 85.249 87.26 88.63 89.94
11 100 100 100 100 100 100

342
343
344
345