温度场做完了，应力场老是不成功，求助
/PREP7/VIEW,1,1
ET,1,SOLID70
mptemp,,20,200,500,800,1000,1500,1800,2500 !温度
mpdata,dens,1,,7820,7700,7610,7540,7490,7350,7290,7090
mpdata,kxx,1,,50,48.5,40,25.4,30,35.8,140,142
mpdata,c,1,,460,485,530,676,670,660,790,820
mpdata,ex,1,,2.09,1.9,1.72,1.0,0.05,0.01,0.001,0.001 !弹性模量
mpdata,alpx,1,,1.1,1.22,1.39,1.48,1.34,1.33,1.32,1.31 ! 热膨胀系数
mpdata,prxy,1,,0.27,0.27,0.28,0.28,0.29,0.3,0.31,0.32 !泊松比
!屈服强度?
MPDATA,HF,1,,100,350,520,1050,1500,3000,,3100,3500 !定义对应温度下的对流换热系数数值
mpcopy,,1,2
!******************************定义应力分析参数***************************************
mpdata,ex,1,1,2.05e11,2.05e11,2.05e11,2.05e11,2e11,2e11 !定义弹性模量
mpdata,ex,1,,1.95e11,1.9e11,1.8e11,1.65e11,1.5e11,1.45e11
mpdata,ex,1,,1e11,0.8e11,0.5e11,0.5e11,0.5e11,0.5e11
tb,bkin,1,6 !定义屈服强度和切变模量
tbtemp,25,1
tbdata,1,345e6,0.8e11
tbtemp,200,2
tbdata,1,345e6,0.05e11
tbtemp,500,3
tbdata,1,300e6,0.045e11
tbtemp,1000,4
tbdata,1,250e6,0.01e11
tbtemp,1500,5
tbdata,1,150e6,0.005e11
tbtemp,2000,6
tbdata,1,100e5,0.004e10
mpdata,alpx,1,1,12e-6,12e-6,12e-6,12e-6,12e-6,12e-6 !定义线膨胀系数(考虑凝固收缩中的固液区）
mpdata,alpx,1,,12e-6,12e-6,12e-6,12e-6,12e-6,36e-6
mpdata,alpx,1,,36e-6,36e-6,12e-6,12e-6,12e-6,12e-6
mpdata,prxy,1,1,0.3,0.3,0.3,0.3,0.3,0.3 !定义泊松比
mpdata,prxy,1,,0.3,0.3,0.3,0.3,0.3,0.3
mpdata,prxy,1,,0.3,0.3,0.3,0.3,0.3,0.3
mpcopy,,1,2
!*********************************************
k,1,0,0.000,0.000
k,2,0,0.000,0.003
k,3,0,0.017,0.003
k,4,0,0.017,0.005
k,5,0,0.019,0.007
k,6,0,0.021,0.007
k,7,0,0.050,0.030
k,8,0,0.050,0.025
k,9,0,0.025,0.000
k,10,0,0.020,0.000
k,11,0,0.020,0.002
k,12,0,0.021,0.003
k,13,0,0.022,0.005
k,14,0,0.050,0.000
k,15,0,0.018,0.002
k,16,0,0.017,0
k,17,0,0.026,0.006
LSTR, 1, 2
LSTR, 2, 3
BSPLIN,3,4,5,6
BSPLIN,11,12,13,6
BSPLIN,3,15,11
LARC,7,6,14,0.03,
LARC,10,11,14,0.03,
LARC,8,17,14,0.025,
LARC,17,9,14,0.025,
LSTR, 7, 8
LSTR, 1, 16
LSTR, 16, 10
LSTR, 9, 10
LSTR, 16, 3
LSTR, 6, 17
al,3,4,5
al,4,15,9,13,7
al,6,10,8,15
al,5,14,12,7
al,1,2,14,11
vext,all,,,0.11,
vglue,all
rb=0.0011
LESIZE,3, , ,6, , , , ,1
LESIZE,4, , ,5, , , , ,1
LESIZE,5, , ,4, , , , ,1
LESIZE,14, , ,4, , , , ,1
LESIZE,12, , ,4, , , , ,1
LESIZE,7, , ,4, , , , ,1
LESIZE,13, , ,4, , , , ,1
LESIZE,9, , ,8, , , , ,1
LESIZE,15, , ,4, , , , ,1
LESIZE,8, , ,20, , , , ,1
LESIZE,6, , ,20, , , , ,1
LESIZE,10, , ,4, , , , ,1
LESIZE,1, , ,4, , , , ,1
LESIZE,11, , ,10, , , , ,1
LESIZE,2, , ,10, , , , ,1
Lesize,40,4*rb
Lesize,36,4*rb
Lesize,28,4*rb
Lesize,26,4*rb
Lesize,32,4*rb
Lesize,41,4*rb
Lesize,19,4*rb
Lesize,32,4*rb
TYPE, 1
MAT, 1
VSWEEP,1
TYPE, 1
MAT, 2
VSWEEP,2
VSWEEP,3
VSWEEP,4
VSWEEP,5
！*****************************初试条件**********************
finish
/SOLU
AnTYPE,TRANS
TUNIF,100！初始温度
kbc,1
timint,on！瞬时效应
CONV1=110
allsel,all !选择所有外面的节点
nsel,s,ext
sf,all,conv,CONV1,25
allsel
！******************稳态分析************************************************
vsel,s,volume,,1
eslv,s,1
ekill,all
allsel
delt=0.1 ！步长随便
*do,i,1,100
time,i*delt
vsel,s,volume,,1
nslv,s,1
nsel,r,loc,x,(i-1)*0.0011-0.000001,(i-1)*0.0011+0.000001
esln,s,0
ealive,all
allsel
nsel,s,loc,x,(i-1)*0.0011-0.000001,(i-1)*0.0011+0.000001
esln,s,0
esel,r,mat,,1
bfe,all,hgen,1,8e+10 ！热源温度自己随便
allsel
*if,i,gt,1,then
nsel,s,loc,x,-0.000001,(i-2)*0.0011+0.000001
esln,s,0
esel,r,mat,,1
BFEDELE,ALL,HGEN！删除单元体上的载荷
allsel
*endif
nsubst，1
autots,on
outres,all,all
solve
*enddo
!************************************冷却***************************************
DDELE,ALL,TEMP ! 删除温度载荷
TIME,180 ! 定义求解时间
DELTIME,2,1,10 ! 定义时间子步
SOLVE ! 求解
!************************************应力场***************************************
/com,structural
/prep7
alls
eplot
etchg,tts
lsclear,solid
!************************************施加位移约束***************
asel,s,area,,15
nsla,s,1
d,all,all,0
alls
eplot
!***********非线性求解选项**************
/solu
antype,4
trnopt,full !完全瞬态积分法
lumpm,0 !使用与单元相关的质量矩阵公式
nlgeom,on !打开大变形选项
pred,on !打开预测校正
sstif,1 !打开应力刚化功能
solcontrol,on !自动控制
nropt,full,,on !完全牛顿-拉普森方法
cnvtol,u,,0.5 !设置收敛值，位移
!cnvtol,f,,0.05 !设置收敛值，力
cutcontrol,plslimit,10,!设置在一个非线性求解中时间步的减少
eqslv,sparse,,2 !求解器类型
lnsrch,on !线性搜索
pstres,1 !打开预应力选项
autots,on !自动步长
tintp,,,,1,, !向后积分
timint,on !瞬态效果
tref,25 !参考温度
kbc,0 !连续性载荷
neqit,50 !规定每个子步中最大迭代次数,默认为25
!*******************************死亡(为使对应载荷步内未焊焊缝金属不参与应力计算)********************
/solu
/auto,1
/replot
alls
eplot
esel,s,mat,,1
eplot
ekill,all
alls
esel,s,live
eplot
alls
eplot
!****************************应力场********************
delt=0.1 ！步长随便
*do,i,1,100
time,i*delt
vsel,s,volume,,1
nslv,s,1
nsel,r,loc,x,(i-1)*0.0011-0.000001,(i-1)*0.0011+0.000001
esln,s,0
ealive,all
allsel
nsel,s,loc,x,(i-1)*0.0011-0.000001,(i-1)*0.0011+0.000001
esln,s,0
esel,r,mat,,1
antype,,rest
ldread,temp,i,,,,,rth!
nsubst，1
autots,on
outres,all,all
solve
*enddo 大侠，能帮我看看上面的程序吗？
温度场是用生热率做的，应力场老是出不来
能否帮我改改 楼主，可以交流么，你的问题解决没啊 参数好多，是不是可以尝试简化模型算，再逐步细化条件 我也有同样的问题要问 你的命令流
写的太杂乱无章了
例如，热应力分析时
传热膜系数的材料定义 应直接定义为边界条件比较合适
结构分析中
应该是模型中单元划分有一些问题
警告有100多个
页:
[1]