cfrfosh 发表于 2010-3-8 21:03:13

温度与结构分析中出现了下面问题?

*** WARNING ***                         CP =      13.047   TIME= 21:11:52
The ETCHG command has changed some element types.Please examine not   
only the KEYOPTs, but also the real constants, material properties,   
boundary conditions, and loadings, and then modify as needed.         
*** WARNING ***                         CP =      14.156   TIME= 21:11:54
Constraint on area 2 with inactive load label TEMP not transferred to   
nodes.

请教大侠,是什么问题造成的呢?(还有我下面瞬态分析有无问题呢?有些点出现了离散状态??)
/filename,wendu
/title,the anaslsys of xiangliang thermal
/config,nres,100000
/prep7                        !进入处理器
         
ET,1,solid70                  !定义单元类型
mp,ex,1,3.50e10               !指定混凝土材料与温度有关的属性(弹性模量)
mp,dens,1,2500                !混凝土材料属性( 质量密度)
mp,C,1,0.880                  !混凝土材料属性( 比热容)
mp,KXX,1,3.3                  !混凝土材料属性( 热导率)
                     
k,1,0,0                   !建立半截面模型(关键点)
k,2,1,0
k,3,1,1
k,4,0,1
k,11,0,0,2                  !建立半截面模型(关键点)
k,12,1,0,2
k,13,1,1,2
k,14,0,1,2
v,1,2,3,4,11,12,13,14
allsel
mat,1
vsweep,1
asel,s,loc,y,0                !选择顶板顶面0
cm,mian1,area               !生成顶板顶面元件               

asel,s,loc,y,1                !选择顶板顶面0
cm,mian2,area
asel,s,loc,x,1                !选择顶板顶面0
cm,mian3,area
asel,s,loc,z,0                !选择顶板顶面0
cm,genbu,area               
allsel
*dim,a1,arrab,13      !!!!!!!!!定义综合温度下顶板数组为一维整数列
a1(1)=28
a1(2)=29
a1(3)=31
a1(4)=33
a1(5)=34
a1(6)=35
a1(7)=37
a1(8)=37
a1(9)=37
a1(10)=37
a1(11)=36
a1(12)=36
a1(13)=34
*dim,a2,arrab,13      
a2(1)=28
a2(2)=29
a2(3)=31
a2(4)=33
a2(5)=34
a2(6)=35
a2(7)=37
a2(8)=37
a2(9)=37
a2(10)=37
a2(11)=36
a2(12)=36
a2(13)=34
*dim,a3,arrab,13      
a3(1)=28
a3(2)=29
a3(3)=31
a3(4)=33
a3(5)=34
a3(6)=35
a3(7)=37
a3(8)=37
a3(9)=37
a3(10)=37
a3(11)=36
a3(12)=36
a3(13)=34

*dim,xwdq,arrab,13       !!!!!!!!!定义箱外大气平均温度变化
xwdq(1)=32
xwdq(2)=32
xwdq(3)=33
xwdq(4)=34
xwdq(5)=34
xwdq(6)=34
xwdq(7)=34
xwdq(8)=35
xwdq(9)=35
xwdq(10)=35
xwdq(11)=35
xwdq(12)=35
xwdq(13)=35
/solu                     !加载
antype,transient,new       !做一次新的瞬态分析
timint,off               !先作稳态(静态)分析,确定初始条件
time,0.01                  !设定只有一个子步的时间很小的载荷步
deltim,0.01                !设定时间步长
allsel                      !选择全部
sfa,genbu,,hflux,0
                  
               !选择下顶板元件
sfa,mian1,,CONV,9.4,30
da,mian1,temp,28               
sfa,mian2,,CONV,19.6,30
da,mian2,temp,28               
sfa,mian3,,CONV,21.1,30
da,mian3,temp,28
            
                  
SOLVE                     !得到初始温度分布


sfa,genbu,,hflux,0
*do,i,1,13,1                !!!用do循环来进行瞬态求解
time,i                     !设置每一步荷载结束时的时间
timint,on                  !打开时间积分效果
                     
allsel                        !删除前一阶段边界条件
               !选择下顶板元件
sfa,mian1,,CONV,9.4,xwdq(i)   !在线上施加对流面荷载(19.6为膜系数,a1(i)为平均温度)
da,mian1,temp,a1(i)         
sfa,mian2,,CONV,19.6,xwdq(i)
da,mian2,temp,a2(i)         
sfa,mian3,,CONV,21.1,xwdq(i)
da,mian3,temp,a3(i)         
allsel
TRNOPT,FULL                  !指定完全法为瞬态分析方法
deltim,0.05,0.005,0.1       !使用自动时间跟踪设定时间步长(0.05为子步开始时间,0.005为最小时间步,0.1为最大)
autots,on                   !瞬态分析自动时间步长跟踪
KBC,0
outres,all,1                !将荷载步每隔1个子步写入到数据库中
sftran                      !转换面荷载
SOLVE                     !求解
sfadele,all,all
ddele,all,temp
*enddo                      !do命令运行结束

save

/prep7
/replot
sfadele,all,all
ddele,all,temp
etchg,tts
save
et,1,solid45,1,1
mp,ex,1,3.50e10               !指定混凝土材料与温度有关的属性(弹性模量)
mp,prxy,1,0.167               !混凝土材料属性(泊松比)
mp,dens,1,2500                !混凝土材料属性( 质量密度)
mp,alpx,1,1.0e-5            !混凝土材料属性( 热膨胀系数)
allsel
mat,1
da,genbu,all                  !根部自由度全约束
finish
/solu
Allsel,all
Ldread,temp,0.01,,,2,‘wendu.rth’
*do,i,1,13,1
   
Ldread,temp,i,,,0,‘wendu.rth’!读取temp.rth温度场
time,i                     !设置每一步荷载结束时的时间
timint,on                  !打开时间积分效果
TRNOPT,FULL               !指定完全法为瞬态分析方法
deltim,0.05,0.04,0.1       !使用自动时间跟踪设定时间步长(0.05为子步开始时间,0.005为最小时间步,0.1为最大)
autots,on                   !瞬态分析自动时间步长跟踪
KBC,0
outres,all,5                !将荷载步每隔1个子步写入到数据库中
sftran
SOLVE                     !求解
sfadele,all,all
ddele,all,temp
*enddo                      !do命令运行结束
finish
/post1
plnsol,s,eqv

cfrfosh 发表于 2010-3-9 21:13:27

版主啊,求救啊
页: [1]
查看完整版本: 温度与结构分析中出现了下面问题?