- 积分
- 0
- 注册时间
- 2010-6-19
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2010-10-11 20:17:06
|
显示全部楼层
来自 福建泉州
补充一下 ,我是先编写的材料特性参数文件,以后用命令调用:
*dim,t_s,,1440 !定义两个小时内,时间间隔为5s的1440个点的钢管温度数组
*dim,t_deta,,1440!定义两个小时内,1440个时间点的钢管温度增量数组
*dim,t_pipe,,360 !定义两个小时内,时间间隔为20s的360个点的钢管温度数组
desi=7850 !定义钢管材料的密度
f_v=30 !定义钢管的截面形状系数
deta_t=5 !定义计算钢管1440个时间点的温度值时的时间间隔为5s
c_f=45 !定义钢管材料的热传到系数
*set,t_s(1),20 !设置钢管温度的初值(即常温下钢管的温度)为20摄氏度
*do,i,2,1440,1 !利用循环计算其余1439个时间点的钢管温度值
!计算当前时刻空气温度值(基于ISO834标准升温曲线)
Tg=20+345*LOG10(8*(i-1)*5/60+1)
!计算以辐射方式由空气向构件表面传热的传热系数
a_r=0.5*5.67e-8*((tg+273)**4-(t_s(i-1)+273)**4)/(tg-t_s(i-1))
!计算综合传热系数
k_f=a_r+25
!计算当前钢管温度的温度增量
*set,t_deta(i-1),k_f*f_v*(tg-t_s(i-1))*deta_t/(desi*c_f)
!计算当前的钢管温度
*set,t_s(i),t_s(i-1)+t_deta(i-1)
*enddo
!提取间隔20s的钢管温度值并赋值给t_pipe数组
*do,i,1,360,1
*set,t_pipe(i),t_s((i-1)*4+1)
*enddo
*cfopen,E:\BOOK_APDL\chapter_7\temp_pipe,txt,,append
*vwrite,t_pipe(1)
(f16.8)
*cfclos
:) 没加载也没求解,这样怎么办呢,再请教高手:handshake |
|