求助:如何使用谐响应分析求解这样的热源

iomega

ANSYS可以解热的谐响应分析吗？？？ 请教请教。。
如果能求解正弦变化的热源的话，对于1+sint的形式其实很简单：
得到热源为1的稳态解，然后和sint的谐响应解叠加。

iomega

谢谢，我去试试。

因该没有将稳态解和谐响应解叠加的命令，一般只有自己手动或编程做了。。。

iomega

看了下帮助文件，好像热源是不能作为谐响应载荷的，只有温度可以。。。

试了下，也不对。。。

iomega

谢谢，BTW,能不能帖一个例子，让我学习学习。
我一直以来自己写程序来进行热谐响应分析，用finite difference method 加上复数形式解 T(x,y,z,t) = TT(x,y,z,omega)*exp(i*omega*t), 这里TT(x,y,z,omega)本身是复数。用复数的好处是把transient传热方程变成了steady-state方程， TT(x,y,z,omega)的模就是在频率omega处温度的振幅。  问题是自己写的程序叠代的太慢。。。

好像FEMLAB可以用复数解来谐响应分析，不过没试过。

iomega

好，那我就提一个具体问题：
一根beam,长200um,宽2um,厚2um, 两端温度为0,在梁上施加uniform Q(t):
Q(t)=1e12*sin(omega*t), 求在omega=100,200,300Hz时 beam上各点温度变化的amplitude和相对与Q(t)的相位delay. 我已经解出了analytical solution,可以和ANSYS结果做比较。

我贴出了模型部分， 但不知道如何加载Q(t)和求解。 谢谢帮助。。。

FINISH  ! Make sure we are at BEGIN level   
/Clear
/PREP7
um=1e-6
Length=20*um               !Length of the block in z direction
Width=2*um  
Thickness=2*um
Size=0.5*um                !Mesh size
T0=0

Den=3900
Cp=450
k=100                     
UIMP,1,DENS, , ,Den,          !Define density of material
UIMP,1,KXX,KYY,KZZ,k,k,k,  !Define thermal conductivity of material
UIMP,1,C, , ,Cp,                  !Define heat capacity of material  

ET,1,SOLID70

BLC4,0,0,Width,Length,Thickness
aesize,all,Size,
TYPE,1
MAT,1
mshkey,1
mshape,0,3d
vmesh,ALL   

/solu
tunif,T0                 !Ambient temperature
TREF,T0
asel,s,LOC,y,0,0
asel,a,LOC,y,Length,Length
DA, ALL, TEMP, T0
asel,all

!volume num =1 , 如何加载Q(t)=1e12*sin(omega*t)

!加载热源之后.
/SOLU
ANTYPE,HARMIC ! HARMONIC ANALYSIS
hropt,full
harfrq,100,300
kbc,1
nsubst,3
OUTPR,,NONE
OUTRES,,1
allsel
solve
FINISH

iomega

谢谢，自己搞出来了，和analytical solution的解完全一样。

另外，对你的问题进行一下讨论：

一般情况下热源不可能为负，即都为q = q0+  q1*sin(w*t).  (w为频率）
对于这种热源，物体内任意一点的温度
T(x,y,z,t)= T_steady(x,y,z)+  T(x,y,w)*sin(w*t+D)  (D: phase delay)

T(x,y,w）的大小随w的增大而减小，一般遵循一阶信号系统的规律decay.
其thermal time constant t=RC , R:thermal resistance, C: volume heat capacity