p53.mw

> restart:#"m11_p53"

Considérese una tarjeta electrónica de 100•150•1 mm3, hecha de fibra de vidrio laminada con epoxy (k0,25 W/(m•K)), con unos componentes que disipan 5 W en total y que se supondrán uniformemente distribuidos. Sabiendo que la tarjeta se va a empotrar por los dos lados mayores en zócalos buenos conductores del calor, que se supone que van a estar a 30 ºC, se pide:
a) Hacer un esquema del perfil de temperatura esperado.
b) Determinar la temperatura máxima que se alcanzaría en la tarjeta despreciando las pérdidas al aire ambiente dentro de la caja electrónica donde va situada, y la radiación térmica (i.e. disipación sólo por conducción hacia los zócalos).
c) Determinar el espesor que debería tener un recubrimiento de cobre por una cara de las de 10•15 cm2 para que la temperatura máxima no superase en más de 40 ºC la de los zócalos, en el caso anterior.

Datos:

> read"../therm_eq.m":read"../therm_proc.m":with(therm_proc):

> su:="Cobre":dat:=[Lw=0.1*m_,Lb=0.15*m_,Lt=1e-3*m_,ke=0.25*W_/(m_*K_),Q=5*W_,Tb=(30+273.15)*K_,DTm=40*K_];

[Lw = `+`(`*`(.1, `*`(m_))), Lb = `+`(`*`(.15, `*`(m_))), Lt = `+`(`*`(0.1e-2, `*`(m_))), ke = `+`(`/`(`*`(.25, `*`(W_)), `*`(m_, `*`(K_)))), Q = `+`(`*`(5, `*`(W_))), Tb = `+`(`*`(303.15, `*`(K_))), ...

Image

> dat:=op(dat),Const,SI2,SI1:Cdat:=get_sol_data(su):eqC:=rho*c*diff(T(x,t),t)=k*diff(T(x,t),x,x)+phi;

`*`(rho, `*`(c, `*`(diff(T(x, t), t)))) = `+`(`*`(k, `*`(diff(diff(T(x, t), x), x))), phi)

a) Hacer un esquema del perfil de temperatura esperado.

Estará más caliente lejos de los anclajes. (Ver arriba.)
b) Determinar la temperatura máxima que se alcanzaría en la tarjeta despreciando las pérdidas al aire ambiente dentro de la caja electrónica donde va situada, y la radiación térmica (i.e. disipación sólo por conducción hacia los zócalos).

Se puede aproximar por un problema unidimensional (a lo largo de la longitud corta, Lw=0,1 m) y en régimen estacionario, integrándose fácilmente la ecuación del calor:

> eqC:=0=k*diff(T(x),x,x)+phi;dsolve(eqC,T(x));eq11_8_1;eq11_8_2;eqA:=A=Lb*Lw;eqA_:=subs(dat,%);eqL:=L=Lw;eqL_:=subs(dat,%);eqphi:=phi=Q/(Lw*Lb*Lt);eqphi_:=subs(dat,%);eqT0:=T[0]=subs(x=0,T[1]=Tb,k=ke,rhs(eq11_8_1));eqT0_:=subs(eqL,eqphi,dat,%);#T[0]=TKC(rhs(%));

0 = `+`(`*`(k, `*`(diff(diff(T(x), x), x))), phi)
T(x) = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(phi, `*`(`^`(x, 2)))), `*`(k))), `*`(_C1, `*`(x)), _C2)
T(x) = `+`(T[1], `/`(`*`(`/`(1, 8), `*`(phi, `*`(`^`(L, 2), `*`(`+`(1, `-`(`/`(`*`(4, `*`(`^`(x, 2))), `*`(`^`(L, 2))))))))), `*`(k)))
Q = `*`(A, `*`(phi, `*`(L)))
A = `*`(Lb, `*`(Lw))
A = `+`(`*`(0.15e-1, `*`(`^`(m_, 2))))
L = Lw
L = `+`(`*`(.1, `*`(m_)))
phi = `/`(`*`(Q), `*`(Lw, `*`(Lb, `*`(Lt))))
phi = `+`(`/`(`*`(333333.3334, `*`(kg_)), `*`(m_, `*`(`^`(s_, 3)))))
T[0] = `+`(Tb, `/`(`*`(`/`(1, 8), `*`(phi, `*`(`^`(L, 2)))), `*`(ke)))
T[0] = `+`(`*`(1969.816667, `*`(K_)))

¡IMPOSIBLE!; sabemos que los equipos electrónicos no pueden funcionar a esas temperaturas tan altas, ni tampoco lo soportarían los materiales inertes poliméricos. Se concluye que, o bien los datos de partida, o las hipótesis de despreciar la convección y la radiación, no son adecuadas (la tarjeta empezaría a echar humo, ¡tóxico!, a unos 150 ºC).

Los datos parecen razonables, aunque el valor de la conductividad térmica de la FR-4 en el plano (i.e. en la dirección de las fibras), puede ser casi el doble (lo que bajaría a la mitad el DT), pero lo que parece más errado es considerar que la conducción de calor por los estrechos bordes de la tarjeta va a ser más importante que la convección y la radiación por toda el área, como puede comprobarse con una sencilla estimación. Supongamos que el ambiente en el que está la tarjeta está a unos 50 ºC (e.g. las paredes de la caja electrónica donde está metida, otras tarjetas adyacentes que no se calienten tanto, el aire circundante si lo hay...). Tomando valores representativos para la emisividad (e=1) y para la convección (h=2 W/(m^2·K) para el aire encerrado), y despreciando ahora la conducción en los anclajes, se obtendría:

> eqBEconv:=Q=2*A*h*(T-Tinf);eqBErad:=Q=2*A*epsilon*sigma*(T^4-Tinf^4);dat_:=h=2*W_/(m_^2*K_),epsilon=1,Tinf=(50+273)*K_;Tconv_:=subs(eqA_,dat_,dat,Tinf+Q/(2*A*h));'Tconv_'=TKC(%);Trad_:=subs(eqA_,dat_,dat,SI0,(Tinf^4+Q/(2*A*epsilon*sigma))^(1/4))*K_;'Trad_'=TKC(%);

Q = `+`(`*`(2, `*`(A, `*`(h, `*`(`+`(T, `-`(Tinf)))))))
Q = `+`(`*`(2, `*`(A, `*`(epsilon, `*`(sigma, `*`(`+`(`*`(`^`(T, 4)), `-`(`*`(`^`(Tinf, 4))))))))))
h = `+`(`/`(`*`(2, `*`(W_)), `*`(`^`(m_, 2), `*`(K_)))), epsilon = 1, Tinf = `+`(`*`(323, `*`(K_)))
`+`(`*`(406.3333334, `*`(K_)))
Tconv_ = `+`(`*`(133.1833334, `*`(?C)))
`+`(`*`(342.8927748, `*`(K_)))
Trad_ = `+`(`*`(69.7427748, `*`(?C))) (1)

i.e. la conducción no sería el modo dominante; simplemente por radiación ya se ve que la temperatura máxima sería unos 20 ºC superior a la del ambiente (o 37 ºC si eps=0.5). Y la suma de conveccción y radiación todavía la disminuiría un poco más. De todas formas, cuando la tarjeta está junto a otras tarjetas que también disipan, la conducción por los extremos cobra importancia, y por eso se añade (durante la fabricación) una película de cobre por el lado opuesto al de montaje de los componentes.

c) Determinar el espesor que debería tener un recubrimiento de cobre por una cara de las de 10•15 cm2 para que la temperatura máxima no superase en más de 40 ºC la de los zócalos, en el caso anterior.

Volviendo al estudio de la conducción exclusiva, podemos ahora empezar despreciando la conducción longitudinal por la resina, frente a la del cobre.

> eqphi:=phi=Q/(Lw*Lb*Lc);eqphi_:=subs(dat,%);eqT0:=T[0]=subs(x=0,eqL,T[1]=Tb,k=kc,rhs(eq11_8_1));eqT0:=DTm=subs(eqphi,rhs(%))-Tb;eqLc:=Lc=solve(%,Lc);subs(kc=k,Cdat,dat,%)*1e6;

phi = `/`(`*`(Q), `*`(Lw, `*`(Lb, `*`(Lc))))
phi = `+`(`/`(`*`(333.3333334, `*`(kg_)), `*`(`^`(s_, 3), `*`(Lc))))
T[0] = `+`(Tb, `/`(`*`(`/`(1, 8), `*`(phi, `*`(`^`(Lw, 2)))), `*`(kc)))
DTm = `+`(`/`(`*`(`/`(1, 8), `*`(Q, `*`(Lw))), `*`(Lb, `*`(Lc, `*`(kc)))))
Lc = `+`(`/`(`*`(`/`(1, 8), `*`(Q, `*`(Lw))), `*`(DTm, `*`(Lb, `*`(kc)))))
`+`(`*`(0.1e7, `*`(Lc))) = `+`(`*`(26.50551315, `*`(m_)))

i.e. basta una capa de 27 micrómetros (<0,03 mm) de cobre para que la tarjeta disipe adecuadamente.

Si se considera la conducción combinada del cobre y la FR-4, vemos que el resultado apenas cambia.

> eqT0_combi:=subs(kc=keff,Lc=Lc+Lt,eqT0);keff:=(kc*Lc+ke*Lt)/(Lc+Lt);eqLc_:=Lc=solve(eqT0_combi,Lc);eq_Lc__:=subs(dat,kc=k,Cdat,dat,%);%*1e6;

DTm = `+`(`/`(`*`(`/`(1, 8), `*`(Q, `*`(Lw))), `*`(Lb, `*`(`+`(Lc, Lt), `*`(keff)))))
`/`(`*`(`+`(`*`(Lc, `*`(kc)), `*`(Lt, `*`(ke)))), `*`(`+`(Lc, Lt)))
Lc = `+`(`-`(`/`(`*`(`/`(1, 8), `*`(`+`(`*`(8, `*`(DTm, `*`(Lb, `*`(Lt, `*`(ke))))), `-`(`*`(Lw, `*`(Q)))))), `*`(DTm, `*`(Lb, `*`(kc))))))
Lc = `+`(`*`(0.2586938084e-4, `*`(m_)))
`+`(`*`(0.1e7, `*`(Lc))) = `+`(`*`(25.86938084, `*`(m_))) (2)

CONCLUSIÓN. La conducción térmica a través de la resina es siempre despreciable; la radiación y la convección asegurarán que no suba la temperatura muchísimo, pero en la mayoría de los casos es necesario usar PCB con un baño de cobre (o varias láminas de cobre intermedias, o disponer un disipador térmico especial sobre los chips más potentes, como se hace en la CPU de un PC) para que no sea tan grande la temperatura máxima.

Para una simulación numérica del proceso transitorio puede usarse un método de diferencias finitas como el que se muestra en "HT_1D.m" para Matlab..

>