找回密码
 注册
Simdroid-非首页
查看: 100|回复: 1

温度场分析求助

[复制链接]
发表于 2012-6-14 19:21:09 | 显示全部楼层 |阅读模式 来自 陕西西安
悬赏20仿真币未解决
哪位大侠能看看这个为什么不收敛啊?急求啊,谢谢啦!命令流如下,附件为加载的温度载荷,
/prep7
et,1,plane77
KEYOPT,1,3,1  
mptemp,1,100,200,300,400,500,600,700
mpdata,kxx,1,,46.91,43.67,40.54,37.11,34.99,32.5,31.1
mpdata,dens,1,,7200,7200,7200,7200,7200,7200,7200
mpdata,c,1,,528,569,610,634,675,740,900
mpdata,ex,1,,1.25e11
mpdata,prxy,1,,0.25
mpdata,alpx,1,,12.1e-6
MPTEMP,,,,,,,,
mptemp,1,33
mptemp,2,113
mptemp,3,793
mptemp,4,1033
mptemp,5,1145
mptemp,6,1185
mptemp,7,1237
mptemp,8,1281
mptemp,9,1333
mptemp,10,1449  
mpdata,kxx,2,,19.6
mpdata,kxx,2,,19.98
mpdata,kxx,2,,22.6
mpdata,kxx,2,,25.6
mpdata,kxx,2,,26.63
mpdata,kxx,2,,28
mpdata,kxx,2,,28.14
mpdata,kxx,2,,27.29
mpdata,kxx,2,,27.1
mpdata,kxx,2,,26.8
mpdata,enth,2,,126.33*8001.41e3
mpdata,enth,2,,162.61*7963.66e3
mpdata,enth,2,,523*7641.63e3
mpdata,enth,2,,668.28*7527.78e3
mpdata,enth,2,,750.6*7347e3
mpdata,enth,2,,852.42*7200.1e3
mpdata,enth,2,,920.12*7115e3
mpdata,enth,2,,982.12*7044e3
mpdata,enth,2,,1081.086*6930e3
mpdata,enth,2,,1199.5*6842e3
mpdata,dens,2,,8001.41
mpdata,dens,2,,7963.66
mpdata,dens,2,,7641.63
mpdata,dens,2,,7527.78
mpdata,dens,2,,7347
mpdata,dens,2,,7200.1
mpdata,dens,2,,7115
mpdata,dens,2,,7044
mpdata,dens,2,,6930
mpdata,dens,2,,6842
mpplot,kxx,2
mpplot,enth,2
MPTEMP,,,,,,,,
mptemp,1,290
mptemp,2,610
mptemp,3,810
mptemp,4,910
mptemp,5,1010
mptemp,6,1110
mptemp,7,1150
mptemp,8,1185
mptemp,9,1191
mptemp,10,1205
mptemp,11,1251
mptemp,12,1271
mpdata,ex,2,,1.725e11
mpdata,ex,2,,1.283e11
mpdata,ex,2,,0.9161e11
mpdata,ex,2,,0.7224e11
mpdata,ex,2,,0.5322e11
mpdata,ex,2,,0.3537e11
mpdata,ex,2,,0.2843e11
mpdata,ex,2,,0.1769e11
mpdata,ex,2,,0.1237e11
mpdata,ex,2,,0.05638e11
mpdata,ex,2,,1.15e9
mpdata,ex,2,,4.031e8
mpdata,prxy,2,,0.2745
mpdata,prxy,2,,0.2745
mpdata,prxy,2,,0.2745
mpdata,prxy,2,,0.2745
mpdata,prxy,2,,0.2745
mpdata,prxy,2,,0.2745
mpdata,prxy,2,,0.2745
mpdata,prxy,2,,0.2877
mpdata,prxy,2,,0.3094
mpdata,prxy,2,,0.3534
mpdata,prxy,2,,0.4029
mpdata,prxy,2,,0.4277
MPTEMP,,,,,,,,
mptemp,1,0
mptemp,2,100
mptemp,3,200
mptemp,4,300
mptemp,5,400
mptemp,6,500
mptemp,7,600
mptemp,8,700
mptemp,9,800
mpdata,alpx,2,,10.4e-6
mpdata,alpx,2,,11.15e-6
mpdata,alpx,2,,11.94e-6
mpdata,alpx,2,,12.6e-6
mpdata,alpx,2,,13.4e-6
mpdata,alpx,2,,14.1e-6
mpdata,alpx,2,,15.3e-6
mpdata,alpx,2,,13.4e-6
mpdata,alpx,2,,10.8e-6
tb,bkin,2,6     !指定材料类型
tbtemp,290       !设置温度点
tbdata,1,185e6,8.6e9   !输入屈服强度和切变模量
tbtemp,610
tbdata,1,102.9e6,6.41e9
tbtemp,810
tbdata,1,144.8e6,0.458e9
tbtemp,910
tbdata,1,87e6,0.361e9
tbtemp,1010                                 
tbdata,1,75.16e6,2.66e8
tbtemp,1150
tbdata,1,50e6,0
MPTEMP,,,,,,,,
mptemp,1,23
mptemp,2,50
mptemp,3,200
mptemp,4,400
mptemp,5,500
mptemp,6,700
mptemp,7,800
mptemp,8,900
mptemp,9,950
mptemp,10,1000
mptemp,11,1100
mptemp,12,1300
mptemp,13,1400  
mpdata,kxx,3,,43.73
mpdata,kxx,3,,43.89
mpdata,kxx,3,,42.66
mpdata,kxx,3,,39.99
mpdata,kxx,3,,39.65
mpdata,kxx,3,,35.5
mpdata,kxx,3,,34.3
mpdata,kxx,3,,35.3
mpdata,kxx,3,,36.66
mpdata,kxx,3,,36.47
mpdata,kxx,3,,37.11
mpdata,kxx,3,,33.42
mpdata,kxx,3,,33.76
MPTEMP,,,,,,,,
mptemp,1,23
mptemp,2,100
mptemp,3,300
mptemp,4,500
mptemp,5,700
mptemp,6,800
mptemp,7,1000
mptemp,8,1100
mptemp,9,1120
mptemp,10,1150
mptemp,11,1200
mptemp,12,1400
mpdata,dens,3,,7125
mpdata,dens,3,,7108
mpdata,dens,3,,7057
mpdata,dens,3,,7001
mpdata,dens,3,,6942
mpdata,dens,3,,6910
mpdata,dens,3,,6801
mpdata,dens,3,,6758
mpdata,dens,3,,6725
mpdata,dens,3,,6600
mpdata,dens,3,,6484
mpdata,dens,3,,6426
MPTEMP,,,,,,,,
mptemp,1,797
mptemp,2,829
mptemp,3,965
mptemp,4,997
mptemp,5,1101
mptemp,6,1149
mptemp,7,1157
mptemp,8,1173
mptemp,9,1181
mptemp,10,1189
mptemp,11,1205
mptemp,12,1221
mptemp,13,1301
mptemp,14,1405
mpdata,enth,3,,378.567*6919.6e3
mpdata,enth,3,,413.892*6894.34e3
mpdata,enth,3,,594.66*6820.25e3
mpdata,enth,3,,619.678*6802.25e3
mpdata,enth,3,,702.889*6756.35e3
mpdata,enth,3,,742.604*6604e3
mpdata,enth,3,,800.072*6585e3
mpdata,enth,3,,1011.86*6547.1e3
mpdata,enth,3,,1026.47*6528.7e3
mpdata,enth,3,,1041.44*6510.3e3
mpdata,enth,3,,1059.6*6481.6e3
mpdata,enth,3,,1065.57*6473.92e3
mpdata,enth,3,,1132.34*6446e3
mpdata,enth,3,,1214.51*6423e3
mpplot,kxx,3
mpplot,enth,3
MPTEMP,,,,,,,,
mptemp,1,56
mptemp,2,176
mptemp,3,236
mptemp,4,356
mptemp,5,536
mptemp,6,656
mptemp,7,836
mptemp,8,1076
mptemp,9,1136
mptemp,10,1142
mptemp,11,1160
mptemp,12,1163
mptemp,13,1196
mpdata,ex,3,,1.3465e11
mpdata,ex,3,,1.3446e11
mpdata,ex,3,,1.335e11
mpdata,ex,3,,1.2695e11
mpdata,ex,3,,1.0460e11
mpdata,ex,3,,0.8522e11
mpdata,ex,3,,5.8353e10
mpdata,ex,3,,2.0229e10
mpdata,ex,3,,1.017e10
mpdata,ex,3,,8.8463e9
mpdata,ex,3,,3.07e8
mpdata,ex,3,,2.43e8
mpdata,ex,3,,0.007e8
mpdata,prxy,3,,0.369
mpdata,prxy,3,,0.369
mpdata,prxy,3,,0.369
mpdata,prxy,3,,0.369
mpdata,prxy,3,,0.369
mpdata,prxy,3,,0.369
mpdata,prxy,3,,0.369
mpdata,prxy,3,,0.36975
mpdata,prxy,3,,0.369889
mpdata,prxy,3,,0.376826
mpdata,prxy,3,,0.39158
mpdata,prxy,3,,0.4659
mpdata,prxy,3,,0.4973
MPTEMP,,,,,,,,
mptemp,1,170
mptemp,2,330
mptemp,3,470
mptemp,4,610
mptemp,5,750
mptemp,6,810
mptemp,7,950
mptemp,8,1130
mptemp,9,1169
mptemp,10,1194
mpdata,alpx,3,,12.677e-6
mpdata,alpx,3,,13.369e-6
mpdata,alpx,3,,13.744e-6
mpdata,alpx,3,,13.983e-6
mpdata,alpx,3,,14.329e-6
mpdata,alpx,3,,9.5466e-6
mpdata,alpx,3,,10.395e-6
mpdata,alpx,3,,10.968e-6
mpdata,alpx,3,,11.457e-6
mpdata,alpx,3,,9.0705e-6
save

tb,bkin,3,6     !指定材料类型
tbtemp,56       !设置温度点
tbdata,1,186.4e6,6.72e9   !输入屈服强度和切变模量
tbtemp,236
tbdata,1,169.9e6,6.6e9
tbtemp,449
tbdata,1,144.8e6,5.472e9
tbtemp,499
tbdata,1,131e6,5.419e9
tbtemp,1153                                 
tbdata,1,13.16e6,2.02e8
tbtemp,1160  
tbdata,1,5e6,0
MPTEMP,,,,,,,,
mptemp,1,300
mptemp,2,500
mptemp,3,800
mptemp,4,900
mptemp,5,950
mptemp,6,1000
mptemp,7,1100
mptemp,8,1200
mpdata,kxx,4,,0.65
mpdata,kxx,4,,0.55
mpdata,kxx,4,,0.38
mpdata,kxx,4,,0.45
mpdata,kxx,4,,0.55
mpdata,kxx,4,,0.6
mpdata,kxx,4,,0.7
mpdata,kxx,4,,0.8
mpdata,c,4,,430
mpdata,c,4,,590
mpdata,c,4,,710
mpdata,c,4,,755
mpdata,c,4,,795
mpdata,c,4,,805
mpdata,c,4,,810
mpdata,c,4,,820
mp,dens,4,2800
mp,dens,5,5000
mp,c,5,1000
mp,kxx,5,0.5
mp,dens,6,1700
mp,c,6,960
MPTEMP,,,,,,,,
mptemp,1,25
mptemp,2,100
mptemp,3,200
mptemp,4,300
mptemp,5,400
mptemp,6,500
mptemp,7,600
mptemp,8,700
mpdata,kxx,6,,1.1
mpdata,kxx,6,,1.2
mpdata,kxx,6,,1.3
mpdata,kxx,6,,1.4
mpdata,kxx,6,,1.5
mpdata,kxx,6,,1.6
mpdata,kxx,6,,1.7
mpdata,kxx,6,,1.7
mpdata,ex,4,,1.25e11
mpdata,prxy,4,,0.25
mpdata,alpx,4,,12.1e-6
mpdata,ex,5,,1.25e11
mpdata,prxy,5,,0.25
mpdata,alpx,5,,12.1e-6
mpdata,ex,6,,1.25e11
mpdata,prxy,6,,0.25
mpdata,alpx,6,,12.1e-6
!建立模型
rectng,0.125,0.185,-0.306,0.306
aplot
rectng,0.185,0.187,-0.306,0.306
aplot
rectng,0.187,0.267,-0.306,0.306
aplot
rectng,0.267,0.335,-0.31,0.31
aplot
rectng,0.125,0.267,-0.31,-0.306
alpot
rectng,0.125,0.267,0.306,0.31
aplot
rectng,0.1,0.185,0.31,0.96
rectng,0.1,0.185,-0.86,-0.31
rectng,0,0.185,-0.94,-0.86
rectng,0,0.125,-0.31,0.31
rectng,0,0.1,-0.86,0.96
aplot
/pnum,area,1
aplot
aadd,3,4
aplot
aadd,8,9
aadd,10,11
aplot
numcmp,all
aglue,all
numcmp,all

!划分网格
/pnum,area,1
replot
mat,2
lesize,1,0.001
lesize,3,0.001
lesize,2,0.003
lesize,4,0.003
amesh,1
mat,4
lesize,2,0.003
lesize,5,0.003
lesize,21,0.001
lesize,22,0.001
amesh,2,,3
mat,1
esize,0.01
amesh,6
mat,5
lesize,31,0.005
lesize,23,0.005
lesize,3,0.005
lesize,30,0.005
lesize,29,0.01
lesize,10,0.01
amesh,5,,3
mat,5
lesize,3,0.005
lesize,26,0.005
lesize,27,0.005
lesize,28,0.005
lesize,25,0.01
lesize,9,0.01
amesh,4,,3
mat,3
esize,0.005
amesh,8
mat,6
esize,0.005
amesh,7
mat,6
esize,0.005
amesh,3
aplot
/pnum,area,0
*dim,t,table,48275,1
*tread,t,nodetemp,txt
/SOLU
ANTYPE,TRANS             ! 设置为瞬态求解
NROPT,FULL
timint,off
/pnum,line,1
replot
esel,s,mat,,3
nsle,s
d,all,temp,1400
esel,s,mat,,6
nsle,s
d,all,temp,140
nsel,s,node,,1,48275
*GET,N_TOTAL,NODE,,COUNT   
*DO,I,1,N_TOTAL                     !循环开始
D,I,temp,T(I,)  
*ENDDO  
allsel
lplot
sfl,7,CONV,15,,35
sfl,33,CONV,15,,35
sfl,32,CONV,15,,35  
outres,all,all
outres,all,all
NLGEOM,ON !打开大位移效果
NROPT,FULL !设定牛顿-拉夫森选项
time,0.001
deltim,0.001   
outres,all,all
/psf,conv,2
allsel
solve
time,0.005
autots,on
timint,on
deltim,0.002,0.001,0.2
ddele,all,temp
outres,all,all
solve

time,1
autots,on
timint,on
deltim,0.002,0.001,0.2
outres,all,all
solve

time,10
autots,on
timint,on
ddele,all,temp
deltim,1,0.01,3
outres,all,all
solve
time,50
autots,on
timint,on
deltim,5,0.1,10
outres,all,all
solve

time,3500
autots,on
timint,on
deltim,50,30,80
outres,all,all
solve

附件: 您需要 登录 才可以下载或查看,没有账号?注册
 楼主| 发表于 2012-6-17 11:39:14 | 显示全部楼层 来自 陕西西安
Simdroid开发平台
材料3热焓有问题
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-9-18 01:22 , Processed in 0.030403 second(s), 9 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表