焊接过程怎么总是出现与Newton-Raphson有关的警告
fini/clear
!进入前处理器
/PREP7
!选择单元
ET,1,PLANE55
ET,2,SOLID70
!定义材料模型
MPTEMP,1,20,500,1000,1500,2000
MPDATA,EX,1,1,1.93E11,1.5E11,0.7E11,0.1E11,0.01E11
MPDATA,EX,2,1,1.02E11,0.5E11,0.08E11,0.001E11,0.0001E11
MPDATA,EX,3,1,1.17E11,0.9E11,0.3E11,0.05E11,0.005E1
!定义材料密度
MP,DENS,1,8030
MP,DENS,2,4850
MP,DENS,3,8900
!定义热传导系数
MP,KXX,1,16.3
MP,KXX,2,7.44
MP,KXX,3,393
!定义比热
MP,C,1,502
MP,C,2,544
MP,C,3,385
!建立关键点
K,,,,,
K,,,0.002,,
K,,,0.006,,
K,,0.0029,0.006,,
K,,0.1008,0.006,,
K,,0.1008,0.002,,
K,,0.1008,0.00,,
K,,0.1008,-0.01,,
K,,0,-0.01,,
K,,0.0008,-0.01,,
K,,0.0008,0,,
K,,0.0008,0.002,,
!由关键点建立面
A,1,2,12,11
A,2,3,4,12
A,11,7,6,12
A,12,4,5,6
A,1,9,10,11
A,10,8,7,11
!拷贝关键点
KGEN,2,9,,,,,-0.1,,0
!生成直线
LSTR,9,13
!进行布尔操作，将所有面粘合
AGLUE,1,2,3,4,5,6
!定义面的属性
AATT,3, , 1, 0,1
AATT,3, , 1, 0,2
AATT,1, , 1, 0,3
AATT,1, , 1, 0,4
AATT,2, , 1, 0,5
AATT,2, , 1, 0,6
!对线定义单元数量
LESIZE,1,,,5,,,,,1
LESIZE,3,,,5,,,,,1
LESIZE,5,,,5,,,,,1
LESIZE,7,,,5,,,,,1
LESIZE,9,,,5,,,,,1
LESIZE,12,,,5,,,,,1
LESIZE,13,,,3,,,,,1
LESIZE,15,,,3,,,,,1
LESIZE,17,,,3,,,,,1
LESIZE,2,,,5,,,,,1
LESIZE,4,,,5,,,,,1
LESIZE,6,,,5,,,,,1
LESIZE,14,,,5,,,,,1
LESIZE,8,,,40,10,,,,1
LESIZE,11,,,40,10,,,,1
LESIZE,16,,,40,10,,,,1
LESIZE,10,,,40,0.1,,,,1
!划分网格
AMESH,1,6,1
!选择单元类型
TYPE,2
!定义单元分数
EXTOPT,ESIZE,60,0,
!选择建立模型后删除面
EXTOPT,ACLEAR,1
!体的单元属性和面相同
EXTOPT,ATTR,1,0,0
!延伸面
VDRAG,1,2,3,4,5,6,18
/SOLU
!选择所有
ALLSEL,ALL
!选择外表面
ASEL,R,EXT
!选择面
asel,s,area,,1,4,1
asel,a,area,,11
asel,a,area,,13
asel,a,area,,15
asel,a,area,,17
asel,a,area,,19,22,1
!选择面上所有节点
NSLA,S,1
!施加对流系数
SF,ALL,CONV,200,25
!选择所有
ALLSEL,ALL
!选择外表面
ASEL,S,EXT
!选择面
asel,s,area,,26,29,1
asel,a,area,,5
asel,a,area,,6
asel,a,area,,24
!选择面上所有节点
NSLA,S,1
!施加载荷
SF,ALL,CONV,100,25
!选择所有
ALLSEL,ALL
!获取参数
*GET,EMAX,ELEM,,COUNT, , , ,
*GET,NMAX,NODE,,COUNT, , , ,
!选择体1、2
VSEL,S,VOLU,,1,2
!在选择的体下编辑
ALLSEL,BELOW,VOLU
!获取参数
*GET,AEMAX,ELEM,,NUM,MAX, , , ,
!定义列
*DIM,ANEO,ARRAY,AEMAX,1,1, , ,
*DIM,ANEZ,ARRAY,AEMAX,1,1, , ,
*DIM,ANEL,ARRAY,AEMAX,1,1, , ,
!获取单元号码
*SET,II,0
*DO,I,1,EMAX
*IF,ESEL(I),EQ,1,THEN
*SET,II,II+1
*SET,ANEO(II),I
*ENDIF
*ENDDO
!获取单元的Z坐标
*DO,I,1,AEMAX
*GET,ANEZ(I),ELEM,ANEO(I),CENT,Z
*ENDDO
!设置开始温度
TUNIF,25,
!杀死单元
*DO,I,1,AEMAX
EKILL,ANEO(I)
*ENDDO
!选择求解器
ANTYPE,4
TRNOPT,FULL
NROPT,FULL
!调用宏文件PAR_FUNCT，输入参数
Z0=-0.1
V=0.00277
R0=0.1/60
MU=0.6
U=22
EI=160
T=0
DT=0.1/60/V
GENQ=U*EI*MU/7.2E-8
!输出结果文件
LUMPM,0
OUTRES,ALL,ALL
!调用宏文件，激活单元、施加载荷
*DO,I,1,60
*SET,T,T+DT
!激活单元宏文件EAELMENT
ZMAX=Z0+V*T
ZMIN=Z0+V*T-R0
J=0
*DO,II,1,AEMAX
*IF,ANEZ(II),LE,ZMAX,THEN
*IF,ANEZ(II),GT,ZMIN,THEN
EALIVE,ANEO(II)
J=J+1
ANEL(J)=ANEO(II)
*ENDIF
*ENDIF
*ENDDO
ANSEL=J
ESEL,S,LIVE
EPLOT
ESEL,ALL
!施加载荷宏文件DQGEN
*DO,K,1,ANSEL
BFE,ANEL(K),HGEN,,GENQ
*ENDDO
TIME,T
AUTOTS,1
NSUBST,4,,,1
KBC,1
TSRES,ERASE
!求解
SOLVE
!删除载荷DEQGEN
*DO,JJ,1,ANSEL
BFEDELE,ANEL(JJ),HGEN
*ENDDO
*ENDDO
FINISH
页:
[1]