/FILNAME,Whee-rail contact,0 /PREP7 !定义单元 ET,1,SOLID185 KEYOPT,1,2,1 !定义材料参数 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2e11 MPDATA,PRXY,1,,0.3 MP,DENS,1,7800 TB,BISO,1,1,2, TBTEMP,0 TBDATA,,500e6,0,,,, !建立模型 BLOCK,0,1,0,0.06,0,0.15, CYL4,0.1,0.16,0.03,0,0.1,360,0.05 VGEN, ,2, , , , ,0.05, , ,1 !划分网格 ESIZE,0.012 VSWEEP,ALL !定义接触对 !* CM,_NODECM,NODE CM,_ELEMCM,ELEM CM,_KPCM,KP CM,_LINECM,LINE CM,_AREACM,AREA CM,_VOLUCM,VOLU /GSAV,cwz,gsav,,temp ! Define surface-based constraint type of pair MAT,1 R,3 REAL,3 ET,2,170 ET,3,174 KEYOPT,3,12,5 KEYOPT,3,4,2 KEYOPT,3,2,2 KEYOPT,2,2,1 KEYOPT,2,4,111111 TYPE,2 ! Create a pilot node ! At center of mass of contact geometric entities N,8881, 0.1,0.16,0.075 TSHAP,PILO E,8881 ! Generate the contact surface ASEL,S,,,11 ASEL,A,,,12 CM,_CONTACT,AREA TYPE,3 NSLA,S,1 ESLN,S,0 NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) ESURF ALLSEL ESEL,ALL ESEL,S,TYPE,,2 ESEL,A,TYPE,,3 ESEL,R,REAL,,3 /PSYMB,ESYS,1 /PNUM,TYPE,1 /NUM,1 EPLOT ESEL,ALL ESEL,S,TYPE,,2 ESEL,A,TYPE,,3 ESEL,R,REAL,,3 CMSEL,A,_NODECM CMDEL,_NODECM CMSEL,A,_ELEMCM CMDEL,_ELEMCM CMSEL,S,_KPCM CMDEL,_KPCM CMSEL,S,_LINECM CMDEL,_LINECM CMSEL,S,_AREACM CMDEL,_AREACM CMSEL,S,_VOLUCM CMDEL,_VOLUCM /GRES,cwz,gsav CMDEL,_TARGET CMDEL,_CONTACT !* !* CM,_NODECM,NODE CM,_ELEMCM,ELEM CM,_KPCM,KP CM,_LINECM,LINE CM,_AREACM,AREA CM,_VOLUCM,VOLU /GSAV,cwz,gsav,,temp MP,MU,1,0.25 MAT,1 MP,EMIS,1,7.88860905221e-031 R,4 REAL,4 ET,4,170 ET,5,174 R,4,,,1.0,0.1,0, RMORE,,,1.0E20,0.0,1.0, RMORE,0.0,0,1.0,,1.0,0.5 RMORE,0,1.0,1.0,0.0,,1.0 KEYOPT,5,4,0 KEYOPT,5,5,0 KEYOPT,5,7,0 KEYOPT,5,8,0 KEYOPT,5,9,0 KEYOPT,5,10,2 KEYOPT,5,11,0 KEYOPT,5,12,0 KEYOPT,5,2,0 KEYOPT,4,5,0 ! Generate the target surface ASEL,S,,,9 ASEL,A,,,10 CM,_TARGET,AREA TYPE,4 NSLA,S,1 ESLN,S,0 ESLL,U ESEL,U,ENAME,,188,189 NSLE,A,CT2 ESURF CMSEL,S,_ELEMCM ! Generate the contact surface ASEL,S,,,4 CM,_CONTACT,AREA TYPE,5 NSLA,S,1 ESLN,S,0 NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) ESURF ALLSEL ESEL,ALL ESEL,S,TYPE,,4 ESEL,A,TYPE,,5 ESEL,R,REAL,,4 /PSYMB,ESYS,1 /PNUM,TYPE,1 /NUM,1 EPLOT ESEL,ALL ESEL,S,TYPE,,4 ESEL,A,TYPE,,5 ESEL,R,REAL,,4 CMSEL,A,_NODECM CMDEL,_NODECM CMSEL,A,_ELEMCM CMDEL,_ELEMCM CMSEL,S,_KPCM CMDEL,_KPCM CMSEL,S,_LINECM CMDEL,_LINECM CMSEL,S,_AREACM CMDEL,_AREACM CMSEL,S,_VOLUCM CMDEL,_VOLUCM /GRES,cwz,gsav CMDEL,_TARGET CMDEL,_CONTACT /MREP,EPLOT )/GOP ! Resume printing after UNDO process )! We suggest a save at this point !进入求解 /SOL ANTYPE,4 !定义分析类型 TRNOPT,FULL !使用完全法进行瞬态动力学分析 NLGEOM,1 !激活大变形分析 EQSLV,SPAR !设置求解器类型 !定义边界条件 DA,3,ALL, D,8881, , , , , ,ALL, , , , , D,8881, ,-0.00039, , , ,UY, , , , , !设置计算结果的输出控制 OUTRES,ALL,ALL, !打开接触优化控制 SOLCONTROL,ON,1,NOPL !关闭结构瞬态效应 TIMINT,0 TINTP,0.005, , , , , , TIME,1e-3 AUTOTS,-1 NSUBST,1,1,1,1 KBC,0 NEQIT,40, !设置每一个子步的最大迭代次数 SOLVE !删除平动和转动位移约束 DDELE,8881,UX DDELE,8881,ROTZ !施加一个新的转动位移约束 D,8881, ,0.6, , , ,UX, , , , TIMINT,1 !激活结构瞬态效应 TIME,2.5 NSUBST,200,225 ,100 ,1 AUTOTS,1 SOLVE |