UNIVERSIDAD AUTÓNOMA DE SAN LUIS POTOSÍ

FACULTAD DE INGENIERÍA

ÁREA MECÁNICA ELÉCTRICA

INGENIERIA MECANICA ELÉCTRICA

MATERIA: ELEMENTO FINITO

PROFESOR :ALBERTO TORRES CRUZ

PROYECTO PARCIAL 3

NOMBRE : ALAN ACOSTA MARTINEZ

FECHA DE ENTREGA : 16 DE ABRIL DEL 2023

1. Diagrama donde muestre el sistema de referencia a usar e indique las condiciones de frontera

Donde:

Los elementos se enumeran con números romanos, siendo 16.

Los nodos se enumeran en un círculo, son 15.

Es una placa cuadrada de 1mx1m.

Los elementos triangulares tienen las siguientes características, a=.25m y b=.5m.

La línea de la izquierda y de abajo tienen una temperatura de 500 K.

La línea de la derecha y la de arriba tienen una convención de 10 W/m2K.

Se da una generación de calor interna de 𝑄(𝑥,𝑦,𝑑)=10000∙(4𝑥−4𝑥2)−400𝑑[𝑊/𝑚3].

 2. Forma fuerte

La forma fuerte es la siguiente:

∇ ∗ (−𝛼∇𝑇) = 𝑄𝑖
Las condiciones de frontera están dadas por las siguientes ecuaciones:

𝑇| = 500 𝐾
𝑥=0

𝑇| = 500 𝐾
𝑦=1
𝑊
𝑞| = 10 (𝑇𝑖𝑛𝑓 − 𝑇)
𝑥 = 1, 𝑦 = 1 𝑚 𝐾

3. Función de aproximación para la variable primaria

𝑥 𝑦 𝑥 𝑦
𝑇(𝑥) = 1 − − 𝑇𝑖 + 𝑇𝑗 + 𝑇𝑘
. 25 . 5 . 25 .5
4. Gráfico de la función de generación de calor

5. Matrices elementales
𝐾 𝑎 +𝑏 −𝑏 −𝑎 𝑇1 𝑞1 1 𝑄𝑖 ∗ 𝑎 ∗ 𝑏
−𝑏 𝑏 0 𝑇2 = 𝑞2 + 1
2∗𝑎∗𝑏 𝑞3 6
−𝑎 0 𝑎 𝑇3 1

Elemento 1
 𝐾 . 25 +. 5 −. 5 −. 25 𝑇1 1 (𝑄(0) + 𝑄(0) + 𝑄(1/4)) ∗ 𝑎 ∗ 𝑏
−. 5 .5 0 𝑇2 = + 1
2 ∗ .5 ∗ .25 6∗3
−. 25 0 . 25 𝑇3 1
Elemento 2

𝐾 . 25 +. 5 −. 5 −. 25 𝑇2 1 (𝑄(0) + 𝑄(1/4) + 𝑄(1/4)) ∗ 𝑎 ∗ 𝑏

−. 5 .5 0 𝑇3 = + 1
. 25 6∗3
−. 25 0 . 25 𝑇4 1
Elemento 3

𝐾 . 25 +. 5 −. 5 −. 25 𝑇2 1 (𝑄(1/2) + 𝑄(1/4) + 𝑄(1/4)) ∗ 𝑎 ∗ 𝑏

−. 5 .5 0 𝑇4 = + 1
. 25 6∗3
−. 25 0 . 25 𝑇5 1
Elemento 4

𝐾 . 25 +. 5 −. 5 −. 25 𝑇4 1 (𝑄(1/2) + 𝑄(1/2) + 𝑄(1/4)) ∗ 𝑎 ∗ 𝑏

−. 5 .5 0 𝑇5 = + 1
. 25 6∗3
−. 25 0 . 25 𝑇6 1
Elemento 5

𝐾 . 25 +. 5 −. 5 −. 25 𝑇5 1 (𝑄(1/2) + 𝑄(1/2) + 𝑄(3/4)) ∗ 𝑎 ∗ 𝑏

−. 5 .5 0 𝑇6 = + 1
. 25 6∗3
−. 25 0 . 25 𝑇7 1
Elemento 6

𝐾 . 25 +. 5 −. 5 −. 25 𝑇6 1 (𝑄(1/2) + 𝑄(3/4) + 𝑄(3/4)) ∗ 𝑎 ∗ 𝑏

−. 5 .5 0 𝑇7 = + 1
. 25 6∗3
−. 25 0 . 25 𝑇8 1
Elemento 7

𝐾 . 25 +. 5 −. 5 −. 25 𝑇7 1 (𝑄(1) + 𝑄(3/4) + 𝑄(3/4)) ∗ 𝑎 ∗ 𝑏

−. 5 .5 0 𝑇8 = + 1
. 25 6∗3
−. 25 0 . 25 𝑇9 1
Elemento 8

𝐾 . 25 +. 5 −. 5 −. 25 𝑇8
−. 5 .5 0 𝑇9
. 25
−. 25 0 . 25 𝑇10
1 (𝑄(1) + 𝑄(1) + 𝑄(3/4)) ∗ 𝑎 ∗ 𝑏 0
=+ 1 + 0 ∗ ℎ ∗ (𝑇 − 𝑇10) ∗ 𝑎
6∗3
1 1
Elemento 9

𝐾 . 25 +. 5 −. 5 −. 25 𝑇3 1 (𝑄(0) + 𝑄(0) + 𝑄(1/4)) ∗ 𝑎 ∗ 𝑏

−. 5 .5 0 𝑇4 = + 1
2 ∗ .5 ∗ .25 6∗3
−. 25 0 . 25 𝑇11 1
Elemento 10
 𝐾 . 25 +. 5 −. 5 −. 25 𝑇4
−. 5 .5 0 𝑇11
. 25
−. 25 0 . 25 𝑇12
1 (𝑄(0) + 𝑄(1/4) + 𝑄(1/4)) ∗ 𝑎 ∗ 𝑏 0
=+ 1 + 0 ∗ ℎ ∗ (𝑇 − 𝑇14) ∗ 𝑏
6∗3
1 1
Elemento 11

𝐾 . 25 +. 5 −. 5 −. 25 𝑇4 1 (𝑄(1/2) + 𝑄(1/4) + 𝑄(1/4)) ∗ 𝑎 ∗ 𝑏

−. 5 .5 0 𝑇6 = + 1
. 25 6∗3
−. 25 0 . 25 𝑇12 1

Elemento 12

𝐾 . 25 +. 5 −. 5 −. 25 𝑇6
−. 5 .5 0 𝑇12
. 25
−. 25 0 . 25 𝑇13
1 (𝑄(1/2) + 𝑄(1/2) + 𝑄(1/4)) ∗ 𝑎 ∗ 𝑏 0
=+ 1 + 0 ∗ ℎ ∗ (𝑇 − 𝑇13) ∗ 𝑏
6∗3
1 1
Elemento 13

𝐾 . 25 +. 5 −. 5 −. 25 𝑇6 1 (𝑄(1/2) + 𝑄(1/2) + 𝑄(3/4)) ∗ 𝑎 ∗ 𝑏

−. 5 .5 0 𝑇8 = + 1
. 25 6∗3
−. 25 0 . 25 𝑇13 1
Elemento 14

𝐾 . 25 +. 5 −. 5 −. 25 𝑇8
−. 5 .5 0 𝑇13
. 25
−. 25 0 . 25 𝑇14
1 (𝑄(1/2) + 𝑄(3/4) + 𝑄(3/4)) ∗ 𝑎 ∗ 𝑏 0
=+ 1 + 0 ∗ ℎ ∗ (𝑇 − 𝑇114) ∗ 𝑏
6∗3
1 1
Elemento 15

𝐾 . 25 +. 5 −. 5 −. 25 𝑇8 1 (𝑄(1) + 𝑄(3/4) + 𝑄(3/4)) ∗ 𝑎 ∗ 𝑏

−. 5 .5 0 𝑇10 = + 1
. 25 6∗3
−. 25 0 . 25 𝑇14 1
Elemento 16

𝐾 . 25 +. 5 −. 5 −. 25 𝑇10
−. 5 .5 0 𝑇14
. 25
−. 25 0 . 25 𝑇15
1 (𝑄(1) + 𝑄(1) + 𝑄(3/4)) ∗ 𝑎 ∗ 𝑏 0
=+ 1 + 0 ∗ ℎ ∗ (𝑇 − 𝑇15) ∗ 𝑎
6∗3
1 1
6. Matrices globales

La matriz global es la siguiente

Término fuente, condición de frontera Dirichlet y convección

Suma total del término fuente, condición de frontera Dirichlet y convección

 7. Matriz que se usa para resolver el sistema(Para 16 elementos)
𝑇1
⎡ 𝑇2 ⎤ ⎡ ⎤
⎢ 𝑇3 ⎥ ⎢ ⎥
⎡ ⎤ ⎢ 𝑇4 ⎥ ⎢ ⎥
⎢ ⎥ ⎢ 𝑇5 ⎥ ⎢ ⎥
⎢ ⎥ ⎢ 𝑇6 ⎥ ⎢ ⎥
⎢ ⎥ ⎢ 𝑇7 ⎥ ⎢ ⎥
⎢ ⎥ ⎢ 𝑇8 ⎥ = ⎢ ⎥
⎢ ⎥ ⎢ 𝑇9 ⎥ ⎢ ⎥
⎢ ⎥ ⎢𝑇10⎥ ⎢ ⎥
⎢ ⎥ ⎢𝑇11⎥ ⎢ ⎥
⎢ ⎥ ⎢𝑇12⎥ ⎢ ⎥
⎣ ⎦ ⎢𝑇13⎥ ⎢ ⎥
⎢𝑇14⎥ ⎢ ⎥
⎣𝑇15⎦ ⎣ ⎦
8. Temperatura en el punto (0.5, 0.5)(Para 16elementos) (Matlab)

Como se realizaron 16 elementos, simplemente tomamos la temperatura correspondiente a ese punto,

la cual es

9. Temperatura en el punto (0.5, 0.5) (Para 16elementos) (Ansys)

10. Gráfica2Dde la función solución para la variable primaria (Para 16elementos) (Matlab)

11. Gráfica 2D de la función solución para la variable primaria (Para 16elementos) (ANSYS)

Matlab Ansys Error

507.14 507.58 0.08%
 12. Código
clc
close all
clear all
%MEF 2-D
%Datos
ancho=1;
largo=1;
x=0:1/11:1;
NL=1;
Q=10000*(4*x-4*x.^2)-400*NL;
alfa=50;
a=1/4;
b=0.5;
h=10;
T_inf=300;

%%Funcion de calor ??????

Qf=@(x) 10000*(4*x-4*x.^2)-400*NL;
Qav1= (Qf(0)+Qf(0)+Qf(1/4))/3
Qav9= (Qf(0)+Qf(0)+Qf(1/4))/3

Qav2= (Qf(1/4)+Qf(1/4)+Qf(0))/3
Qav10=(Qf(1/4)+Qf(1/4)+Qf(0))/3

Qav3= (Qf(1/2)+Qf(1/4)+Qf(1/4))/3
Qav11= (Qf(1/2)+Qf(1/4)+Qf(1/4))/3

Qav4= (Qf(1/2)+Qf(1/2)+Qf(1/4))/3
Qav12= (Qf(1/2)+Qf(1/2)+Qf(1/4))/3

Qav5=(Qf(1/2)+Qf(1/2)+Qf(3/4))/3
Qav13=(Qf(1/2)+Qf(1/2)+Qf(3/4))/3

Qav6= (Qf(1/2)+Qf(1/2)+Qf(3/4))/3
Qav14= (Qf(1/2)+Qf(1/2)+Qf(3/4))/3

Qav7= (Qf(3/4)+Qf(3/4)+Qf(1))/3
Qav15= (Qf(3/4)+Qf(3/4)+Qf(1))/3

Qav8= (Qf(1)+Qf(1)+Qf(3/4))/3
Qav16= (Qf(1)+Qf(1)+Qf(3/4))/3

%Grafica del calor interno

[X,Y]=meshgrid(0:1/4:1,0:0.5:1);
qi_gra=[Qf(0) Qf(1/4) Qf(1/2) Qf(3/4) Qf(1);
Qf(0) Qf(1/4) Qf(1/2) Qf(3/4) Qf(1) ;
Qf(0) Qf(1/4) Qf(1/2) Qf(3/4) Qf(1) ]
%surf(X,Y,qi_gra)
n=15;
%Numero de nodos
%Conectividad, da el numero de nodo y sus esquinas con sus coordenadas
% N i j k xi yi xj yj xk yk alfa
conec=[1 1 2 3 0 0 1/4 0 0 0.5 alfa
2 4 3 2 1/4 .5 0 0.5 1/4 0 alfa
3 2 5 4 1/4 0 2/4 0 1/4 0.5 alfa
4 6 4 5 2/4 .5 1/4 0.5 2/4 0 alfa
5 5 7 6 2/4 0 3/4 0 2/4 0.5 alfa
6 8 6 7 3/4 .5 2/4 0.5 3/4 0 alfa
7 7 9 8 3/4 0 4/4 0 3/4 0.5 alfa
8 10 8 9 4/4 .5 3/4 0.5 4/4 0 alfa
9 3 4 11 0 .5 1/4 0.5 0 1 alfa
10 12 11 4 1/4 1 0 1 1/4 0.5 alfa
11 4 6 12 1/4 .5 2/4 .5 1/4 1 alfa
12 13 12 6 2/4 1 1/4 1 2/4 .5 alfa
13 6 8 13 2/4 .5 3/4 .5 2/4 1 alfa
14 14 13 8 3/4 1 2/4 1 3/4 .5 alfa
15 8 10 14 3/4 .5 4/4 .5 3/4 1 alfa
16 15 14 10 1 1 3/4 1 4/4 0.5 alfa]
%Matrices elementales
%Rigidez
k1=[a^2+b^2 -b^2 -a^2
-b^2 b^2 0
-a^2 0 a^2]*alfa/(2*a*b);
k2=k1;
k3=k1;
k4=k1;
k5=k1;
k6=k1;
k7=k1;
k8=k1;
k9=k1;
k10=k1;
k11=k1;
k12=k1;
k13=k1;
k14=k1;
k15=k1;
k16=k1;
%Termino fuente (generacion de calor)
f1=Qav1*a*b*[1;1;1]/6
f2=Qav2*a*b*[1;1;1]/6
f3=Qav3*a*b*[1;1;1]/6
f4=Qav4*a*b*[1;1;1]/6
f5=Qav5*a*b*[1;1;1]/6
f6=Qav6*a*b*[1;1;1]/6
f7=Qav7*a*b*[1;1;1]/6
f8=Qav8*a*b*[1;1;1]/6
f9=Qav9*a*b*[1;1;1]/6
f10=Qav10*a*b*[1;1;1]/6
f11=Qav11*a*b*[1;1;1]/6
f12=Qav12*a*b*[1;1;1]/6
f13=Qav13*a*b*[1;1;1]/6
f14=Qav14*a*b*[1;1;1]/6
f15=Qav15*a*b*[1;1;1]/6
f16=Qav16*a*b*[1;1;1]/6

%Matrices elementales con dimension de la matriz global

kG1=zeros(n,n);
conec_1=conec(1,[2 3 4]);
kG1(conec_1,conec_1)=k1;

kG2=zeros(n,n);
conec_2=conec(2,[2 3 4]);
kG2(conec_2,conec_2)=k2;

kG3=zeros(n,n);
conec_3=conec(3,[2 3 4]);
kG3(conec_3,conec_3)=k3;

kG4=zeros(n,n);
conec_4=conec(4,[2 3 4]);
kG4(conec_4,conec_4)=k4;

kG5=zeros(n,n);
conec_5=conec(5,[2 3 4]);
kG5(conec_5,conec_5)=k5;

kG6=zeros(n,n);
conec_6=conec(6,[2 3 4]);
kG6(conec_6,conec_6)=k6;

kG7=zeros(n,n);
conec_7=conec(7,[2 3 4]);
kG7(conec_7,conec_7)=k7;

kG8=zeros(n,n);
conec_8=conec(8,[2 3 4]);
kG8(conec_8,conec_8)=k8;

kG9=zeros(n,n);
conec_9=conec(9,[2 3 4]);
kG9(conec_9,conec_9)=k9;

kG10=zeros(n,n);
conec_10=conec(10,[2 3 4]);
kG10(conec_10,conec_10)=k10;

kG11=zeros(n,n);
conec_11=conec(11,[2 3 4]);
kG11(conec_11,conec_11)=k11;

kG12=zeros(n,n);
conec_12=conec(12,[2 3 4]);
kG12(conec_12,conec_12)=k12;

kG13=zeros(n,n);
conec_13=conec(13,[2 3 4]);
kG13(conec_13,conec_13)=k13;

kG14=zeros(n,n);
conec_14=conec(14,[2 3 4]);
kG14(conec_14,conec_14)=k14;

kG15=zeros(n,n);
conec_15=conec(15,[2 3 4]);
kG15(conec_15,conec_15)=k15;

kG16=zeros(n,n);
conec_16=conec(16,[2 3 4]);
kG16(conec_16,conec_16)=k16;

%Fuente
fG1=zeros(n,1);
fG1(conec_1)=f1;

fG2=zeros(n,1);
fG2(conec_2)=f2;

fG3=zeros(n,1);
fG3(conec_3)=f3;

fG4=zeros(n,1);
fG4(conec_4)=f4;

fG5=zeros(n,1);
fG5(conec_5)=f5;

fG6=zeros(n,1);
fG6(conec_6)=f6;

fG7=zeros(n,1);
fG7(conec_7)=f7;

fG8=zeros(n,1);
fG8(conec_8)=f8;

fG9=zeros(n,1);
fG9(conec_9)=f9;

fG10=zeros(n,1);
fG10(conec_10)=f10;

fG11=zeros(n,1);
fG11(conec_11)=f11;

fG12=zeros(n,1);
fG12(conec_12)=f12;

fG13=zeros(n,1);
fG13(conec_13)=f13;
fG14=zeros(n,1);
fG14(conec_14)=f14;

fG15=zeros(n,1);
fG15(conec_15)=f15;

fG16=zeros(n,1);
fG16(conec_16)=f16;

%Ensamble
%Rigidez
kG=kG1+kG2+kG3+kG4+kG5+kG6+kG7+kG8+kG9+kG10+kG11+kG12+kG13+kG14+kG15+kG16
%fuente
fG=fG1+fG2+fG3+fG4+fG5+fG6+fG7+fG8+fG9+fG10+fG11+fG12+fG13+fG14+fG15+fG16

%Condicin de Frontera Tipo Neuman

%Simetria
%sim=zeros(15,1)
%Aislamiento
%aisl=zeros(15,1)
%Extraccin de calor
%CFN=[0;0;0;0;0;0;-300;-300]*b;
%Condicin de Frontera Tipo Robin en fG
CFR_fG=[0;0;0;0;0;0;0;0;0;b;0;a;a;a;b]*h*T_inf;
%Condicin de Frontera Tipo Robin en kG
kG(10,10) = kG(10,10) + h*b;
kG(12,12) = kG(12,12) + h*a;
kG(13,13) = kG(13,13) + h*a;
kG(14,14) = kG(14,14) + h*a;
kG(15,15) = kG(15,15) + h*b;

%Condiciones de frontera tipo Dirichlet

T1=500;
T2=500;
T5=500;
T7=500;
T9=500;
T3=500;
T11=500;
CFD1=-kG(:,1)*T1;
CFD2=-kG(:,2)*T2;
CFD3=-kG(:,5)*T5;
CFD4=-kG(:,7)*T7;
CFD5=-kG(:,9)*T9;
CFD6=-kG(:,3)*T3;
CFD7=-kG(:,11)*T11;

%Reacondicionamiento de la matriz de rigidiez

kG(1,:)=0;
kG(2,:)=0;
kG(5,:)=0;
kG(7,:)=0;
kG(9,:)=0;
kG(3,:)=0;
kG(11,:)=0;
kG(:,1)=0;
kG(:,2)=0;
kG(:,5)=0;
kG(:,7)=0;
kG(:,9)=0;
kG(:,3)=0;
kG(:,11)=0;

kG(1,1)=1;
kG(2,2)=1;
kG(5,5)=1;
kG(7,7)=1;
kG(9,9)=1;
kG(11,11)=1;
kG(3,3)=1;

kG
fG
%CFN
b=fG+CFD1+CFD2+CFD3+CFD4+CFD5+CFD6+CFD7+CFR_fG
%Imponer la temperatura de 500 en los nodos correspondientes
b(1,1)=500;
b(2,1)=500;
b(5,1)=500;
b(7,1)=500;
b(9,1)=500;
b(3,1)=500;
b(11,1)=500;

%Solucion
T=inv(kG)*b

[X,Y]=meshgrid(0:1/4:1,0:0.5:1);
T_graf=[T(1) T(2) T(5) T(7) T(9);
T(3) T(4) T(6) T(8) T(10);
T(11) T(12) T(13) T(14) T(15)]
surf(X,Y,T_graf)

13. Conclusión

Mediante este trabajo se pudo corroborar como se puede realizar una simulación térmica y a su vez, un
método meramente numérico que nos puede ayudar a resolver sistemas de elemento finito complicados,
en esta ocasión se trató de un esquema 2D y de una pieza bastante sencilla con geometrías simples, pero
nos da una idea de lo que podríamos hacer si llegamos a toparnos con un problema de este tipo más
adelante en nuestra carrera.

También se observó la importancia de la convergencia de resultados, que si bien se realizaron de manera

numérica ambos y hay una cierta discrepancia, se aproximan bastante y esto nos dice que ambos métodos
son confiables.