- 积分
- 64
- 注册时间
- 2003-3-10
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2010-9-8 12:46:01
|
显示全部楼层
来自 美国
本帖最后由 iomega 于 2010-9-7 20:54 编辑
I raised this question long time ago in below post but nobody answered it so I have to figured it out by myself.
http://forum.simwe.com/viewthread.php?tid=257772&highlight=%2Biomega
For method (1), you can look at my previous post on Mathcad board, where I already gave the source code of analytical solution in mathcad file:
http://forum.simwe.com/viewthread.php?tid=819280&highlight=%2Biomega
For method (2), just refer to Patankar's famous book of numerical heat transfer simulation, assuming the T(x,y,z,t) is in the form of T(x,y,z,omega)*exp(i*omega*t+ phase), then the transient problem becomes steady-state problem.
For method (3), below is an example of a suspended metal bridge with two-sides kept at T0, subjected to a heat generation of Q*sin(omega*t), where omega is the heating frequency and Q is the amplitude. The heating frequency is 100 Hz in this example.
FINISH ! Make sure we are at BEGIN level
/Clear
/PREP7
/PNUM,VOLU,1
/NUMBER,1
/PNUM,MAT,1
/REPLOT
k=90
DEN=3900
Cp=420
/PREP7
um=1e-6
L=20*um
W=2*um
T=2*um
T0=0
Size=0.25*um
ET,1,SOLID70
UIMP,1,DENS, , ,DEN,
UIMP,1,KXX,KYY,KZZ,k,k,k,
UIMP,1,C, , ,Cp,
BLC4, 0, 0, W, L, T
aesize,all,Size,
TYPE,1
MAT,1
mshkey,1
mshape,0,3d
vmesh,1
/SOLU
ANTYPE,HARMIC
HROPT,FULL
OUTRES,ALL,ALL
tunif,T0
TREF,T0
asel,s,LOC,y,0,0
asel,a,LOC,y,L,L
DA,all,TEMP,T0
allsel,all
nsel,all
BF,ALL,HGEN,1e14
allsel,all
HARFRQ,100,
SOLVE
After the solution is done, you have to read the results: set 1 is the real part and set 2 is imaginary part. They together form the complex temperature rise from which you can easily obtain the amplitude and phase delay.
BTW, bulabula, you are coming from engineering mechanics department, right? Is your boss Liang? Just curious... |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|