Elemento Finito Analisis Termico 2d
Elemento Finito Analisis Termico 2d
Elemento Finito Analisis Termico 2d
FACULTAD DE INGENIERÍA
PROYECTO PARCIAL 3
Donde:
∇ ∗ (−𝛼∇𝑇) = 𝑄𝑖
Las condiciones de frontera están dadas por las siguientes ecuaciones:
𝑇| = 500 𝐾
𝑥=0
𝑇| = 500 𝐾
𝑦=1
𝑊
𝑞| = 10 (𝑇𝑖𝑛𝑓 − 𝑇)
𝑥 = 1, 𝑦 = 1 𝑚 𝐾
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 𝑇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
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 𝑇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 𝑇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
11. Gráfica 2D de la función solución para la variable primaria (Para 16elementos) (ANSYS)
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
[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
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
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.