Re:求助:如何使用谐响应分析求解这样的热源
ANSYS可以解热的谐响应分析吗？？？ 请教请教。。如果能求解正弦变化的热源的话，对于1+sint的形式其实很简单：
得到热源为1的稳态解，然后和sint的谐响应解叠加。
Re:求助:如何使用谐响应分析求解这样的热源
谢谢，我去试试。因该没有将稳态解和谐响应解叠加的命令，一般只有自己手动或编程做了。。。
Re:求助:如何使用谐响应分析求解这样的热源
看了下帮助文件，好像热源是不能作为谐响应载荷的，只有温度可以。。。试了下，也不对。。。
Re:求助:如何使用谐响应分析求解这样的热源
谢谢，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可以用复数解来谐响应分析，不过没试过。
Re:求助:如何使用谐响应分析求解这样的热源
好，那我就提一个具体问题：一根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
Re:求助:如何使用谐响应分析求解这样的热源
谢谢，自己搞出来了，和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
页:
[1]