【原创】+【讨教高手】循环加载热流密度（经典命令流）

asandylau

我做的是制动盘和摩擦片的温度场分析，我做的是让制动盘移动，每次让制动盘移动一个距离，然后加载热流密度和对流进行计算，但发现循环的时候，第3步以后的结果却都是一样的，也就是说相当于加载热流密度和对流没有起到作用，我调试了一个多星期，各个方面都考虑了，我的命令流应该说是没有问题的，但实在找不除原因，希望高手帮忙解决。QQ:1157806538（可以加QQ聊）
              TEL：13109806861

FINI
/CLEAR
/FILN,Brake_thermal_1   
/TITLE,THERMAL ANALYSIS OF DISC BRAKE   
WPSTYLE,,,,,,,,1
/TRIAD,OFF  
/PREP7  
CYL4,,,70E-3,,440E-3,,17.5E-3  !摩擦盘   
BLOCK,200E-3,400E-3,-150E-3,150E-3,0,-13E-3
asel,,loc,z,-13E-3  
VEXT,all,,,,,-20E-3
ALLSEL  
VPLOT   
VGLUE,2,3   
ET,1,70
ET,2,151
!Q345D制动盘材料特性
MP,PRXY,1,0.31      !泊松比  
MP,DENS,1,7870     !制动盘密度
MP,MU,1,0.38       !摩擦系数
MPTEMP  
MPTEMP,1,260,373,473,573,673,1773
MPDATA,EX,1,,2.12e11,2.08e11,2.03e11,1.97e11,1.85e11,1.82e11  !弹性模量
MPDATA,KXX,1,,43.6,42.3,40.9,39.4,37.8,36.1  !热导率
MPDATA,C,1,,410 , 460,507,548,582,605  !比热容   
MPDATA,ALPX,1,,12.2E-6, 12.8E-6, 13.4E-6, 13.9E-6, 14.5E-6, 14.7E-6 !线膨胀系数
!烧结金属材料特性   
MP,PRXY,2,0.32      !泊松比  
MP,DENS,2,8410     !制动盘密度
MP,MU,2,0.38       !摩擦系数
MP,C, 2,384  !比热容
MPTEMP  
MPTEMP,1,260,373,473,573,673,1773
MPDATA,EX,2,,0.853e11,0.834e11,0.795e11,0.726e11,0.697e11,0.668e11  !弹性模量   
MPDATA,KXX,2,,0.39e3, 0.38e3,0.368e3,0.352e3,0.344e3,0.339e3  !热导率
MPDATA,ALPX,2,,17.7E-6, 17.8E-6, 18.0E-6, 18.5E-6, 19.2E-6, 19.9E-6 !线膨胀系数
!ZG340-640制动衬片基体材料特性
MP,PRXY,3,0.3      !泊松比   
MP,DENS,3,7820    !制动盘密度  
MP,ALPX,3,11.3E-6 !线膨胀系数  
MPTEMP  
MPTEMP,1,260,373,473,573,673,1773
MPDATA,EX,3,,1.95e11,1.92e11,1.87e11,1.85e11,1.82e11,1.80e11  !弹性模量
MPDATA,KXX,3,,0.036e3,0.038e3,0.0424e3,0.0413e3,0.0396e3,0.0377e3  !热导率   
MPDATA,C,3,,470,486,494,511,517,557 !比热容  
VSEL,S,,,1  
VATT,1,1,1  
VSEL,S,,,2  
VATT,2,1,1  
VSEL,S,,,3  
VATT,3,1,1  
ALLSEL  
ESIZE,20e-3
WPROTA,,90,
VSBW,1  
WPROTA,,,-90
VSBW,ALL
MSHAPE,0,3D
MSHKEY,1
VMESH,ALL   
WPROTA,,,90
WPROTA,,-90,
ALLSEL  
!NUMMRG,ALL
!NUMCMP,ALL
!DEFINE FULX NODE   
VSEL,S,LOC,Z,-13E-3,0   
ASLV
ASEL,R,LOC,Z,0  
CM,PADFLU,AREA       !定义摩擦片热流面  
NSLA
CM,NPADFLU,NODE          !定义热流节点
VSEL,S,LOC,Z,0,17.5E-3  
VPLOT   
NSLV,S,1
NPLOT   
NSEL,R,LOC,X,200E-3,400E-3  
NSEL,R,LOC,Y,-150E-3,150E-3
NSEL,R,LOC,Z,0  
CM,DISFLU,NODE        !定义制动盘热流节点  
ALLSEL  
VSEL,S,LOC,Z,0,17.5E-3  
ASLV,s,1
CSYS,1  
ASEL,S,LOC,X,70E-3  
ASEL,A,LOC,X,440E-3
CSYS,4  
ASEL,A,LOC,Z,17.5E-3
!ASEL,A,,,1
!ASEL,A,,,27
!ASEL,A,,,32
!ASEL,A,,,36
CM,ADIS,AREA !定义圆盘对流面
ALLSEL  
NSEL,S,LOC,Z,0  
CMSEL,U,DISFLU  
CMSEL,U,NPADFLU
CM,NCONT,NODE   !定义接触面对流点组件
ALLSEL  
VSEL,S,LOC,Z,0,-33E-3   
ASLV
ASEL,R,LOC,Y,-150E-3,150E-3
ASEL,U,LOC,Z,0  
CM,APAD,AREA !定义衬片及烧结金属圆柱体  
FINISH  
/SOL
csys,4  
ANTYPE,TRANS
!NLGEOM,1   
LUMPM,0                              !关闭集中质量矩阵计算格式  
NROPT,FULL  
!TIMINT,0,STRUCT                  !关闭结构时间积分选项   
!TOFFST,460
!TINTP,,,,1                   ！定义瞬态积分参数   
DELTIM,0.001,0.001,0.003
AUTOTS,1
LNSRCH,1
PRED,ON
KBC,0   
TIME,0.01  
!SSTIF,0                 !关闭应力刚化影响  
OUTRES,ERASE
OUTRES,ALL,ALL  
EQSLV, , ,  
PRECISION,0              !定义双精度求解
!MSAVE,0
!IC,ALL,TEMP,298          !定义初始温度   
!SFA, ACYL,,CONV,5.3,298
TUNIF,298   
SFA, APAD,,CONV,5.3,298
SFA, ADIS,,CONV,5.3,298
SF,NCONT,CONV,5.3,298   
csys,4  
!CSWPLA,11,0
!CSYS,11
allsel  
CPINTF,TEMP,1E-8
*SET,FW,0.2505  
!定义制动盘的热流密度函数   
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'tdisflu'
*SET,_FNCCSYS,0
! /INPUT,disflu.func,,,1
*DIM,%_FNCNAME%,TABLE,6,16,1,,,,%_FNCCSYS%  
!   
! Begin of equation: 0.2505*623800*(219.8-49.15*{TIME})*sqrt({X}^2+{Y}^2)   
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0.2505, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 623800, 0, 0, -1
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 49.15, 0, 0, 1  
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 219.8, 0, 0, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -4, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 1, -3, 3, -4
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 2, 0, 0, 2  
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -2   
*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 1, 3, 17, -2   
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 1, -3, 1, -4   
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 16, 1, -2, 0, 0   
*SET,%_FNCNAME%(0,15,1), 0.0, -3, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,16,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 0.2505*623800*(219.8-49.15*{TIME})*sqrt({X}^2+{Y}^2)
!-->
!定义摩擦片的热流密度函数   
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'tpadflu'
*SET,_FNCCSYS,0
! /INPUT,padflu.func,,,1
*DIM,%_FNCNAME%,TABLE,6,16,1,,,,%_FNCCSYS%  
!   
! Begin of equation: 0.7495*623800*(219.8-49.15*{TIME})*sqrt({X}^2+{Y}^2)   
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0.7495, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 623800, 0, 0, -1
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 49.15, 0, 0, 1  
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 219.8, 0, 0, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -4, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 1, -3, 3, -4
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 2, 0, 0, 2  
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -2   
*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 1, 3, 17, -2   
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 1, -3, 1, -4   
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 16, 1, -2, 0, 0   
*SET,%_FNCNAME%(0,15,1), 0.0, -3, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,16,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 0.7495*623800*(219.8-49.15*{TIME})*sqrt({X}^2+{Y}^2)
!-->
SF,DISFLU,HFLUX,%tdisflu%   
SFA,PADFLU,,HFLUX,%tpadflu%
ALLSEL  
SOLVE   

FINISH

FINISH
*do,I,2,10
!I=2
/FILNAME,Brake_thermal_%I%,1
/SOLU   
ALLSEL  
SFDELE,DISFLU,HFLUX
SFDELE,NCONT,CONV   
/REPLO  
FINISH  
/PREP7  
vsel,s,loc,z,0,17.5e-3  
vplot   
csys,1  
alph1=2.198*(I-1)-0.5*49.15*1E-4*(I-1)**2
alph2=2.198*(I-2)-0.5*49.15*1E-4*(I-2)**2
PI=4*ATAN(1)
dis=(alph1-alph2)*180/PI
vgen,,all,,,,dis,,,,1   
vplot   
allsel  
vplot   
csys,4  
CMDELE,DISFLU   
CMDELE,NCONT
VSEL,S,LOC,Z,0,17.5E-3  
VPLOT   
NSLV,S,1
NPLOT   
NSEL,R,LOC,X,200E-3,400E-3  
NSEL,R,LOC,Y,-150E-3,150E-3
NSEL,R,LOC,Z,0  
CM,DISFLU,NODE        !定义制动盘热流节点  
ALLSEL  
NSEL,S,LOC,Z,0  
CMSEL,U,DISFLU  
CMSEL,U,NPADFLU
CM,NCONT,NODE   !定义接触面对流点组件
FINISH  
/SOL
!ALLSEL
!ANTYPE,,REST,1,5,1
LDREAD,TEMP,,,0.01,1,'Brake_thermal_%I-1%','rth',' '
DELTIM,0.001,0.001,0.003
TIME,0.01   
SFA, APAD,,CONV,5.3,298
SFA, ADIS,,CONV,5.3,298
SF,NCONT,CONV,5.3,298   
csys,4  
allsel  
CPINTF,TEMP,1E-8
!定义制动盘的热流密度函数   
*DEL,tdisflu
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'tdisflu'
*SET,_FNCCSYS,0
! /INPUT,disflu.func,,,1
*DIM,%_FNCNAME%,TABLE,6,18,1,,,,%_FNCCSYS%  
!   
! Begin of equation:  0.2505*623800*(219.8-49.15*({TIME}+0.03))*sqrt({X}^2+{Y}^
! 2)
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0.2505, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 623800, 0, 0, -1
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 0.0%I%, 0, 0, 1   
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, 1, 1, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 49.15, 0, 0, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -4, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 219.8, 0, 0, -4
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 1, -1, 2, -4
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -3, 3, -2   
*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, 2, 17, -2   
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, 3, 17, -2   
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, -3, 1, -4   
*SET,%_FNCNAME%(0,16,1), 0.0, -2, 16, 1, -2, 0, 0   
*SET,%_FNCNAME%(0,17,1), 0.0, -3, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,18,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation:  0.2505*623800*(219.8-49.15*({TIME}+0.03))*sqrt({X}^2+{Y}^2)
!-->
!定义摩擦片的热流密度函数   
*DEL,tpadflu
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'tpadflu'
*SET,_FNCCSYS,0
! /INPUT,padflu.func,,,1
*DIM,%_FNCNAME%,TABLE,6,18,1,,,,%_FNCCSYS%  
!   
! Begin of equation: 0.7495*623800*(219.8-49.15*({TIME}+0.01))*sqrt({X}^2+  
! {Y}^2)
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0.7495, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 623800, 0, 0, -1
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 0.0%I%, 0, 0, 1   
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, 1, 1, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 49.15, 0, 0, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -4, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 219.8, 0, 0, -4
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 1, -1, 2, -4
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -3, 3, -2   
*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, 2, 17, -2   
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, 3, 17, -2   
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, -3, 1, -4   
*SET,%_FNCNAME%(0,16,1), 0.0, -2, 16, 1, -2, 0, 0   
*SET,%_FNCNAME%(0,17,1), 0.0, -3, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,18,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 0.7495*623800*(219.8-49.15*({TIME}+0.01))*sqrt({X}^2+{Y}^2)  
!-->

SF,DISFLU,HFLUX,%tdisflu%   
SFA,PADFLU,,HFLUX,%tpadflu%
ALLSEL  
SOLVE
finish
*enddo

辰星小c

全世界最爱美的女性到底是谁？

爱惜自己的名声胜过外表的那些人.......





减肥必看:快速减肥不反弹,吃什么减肥最快,怎么瘦身,怎样减肥最有效和减肥方法食谱