这是一段ansys apdl代码,请为我逐行进行中文注释,注意不要更改代码,且每一行都有注释,不要忽略任意一行
FINISH
/CLEAR
/UNITS,MPA
/TITLE,THREE_GEARS
! ---------- 优化参数 ----------
rm1=1 ! 齿轮轴孔半径(轮毂内径)
rm2=1.5
rm3=1
rh1=1.75 ! 轮毂外径
rh2=2
rh3=1.75
b=2.5 ! 齿轮齿宽
bh13=2.2 ! 轮毂宽
bh2=3
d1=0.5-(b-bh13)/2 ! 轮毂z正向偏移距离(随齿宽/轮毂宽变动)
d2=0
d3=d1
z1=20 ! 齿轮齿数
z2=21
z3=21
m=0.25 ! 模数
ha=1 ! 齿顶高系数
c=0.25 ! 顶隙系数
cf=0.38 ! 齿根圆角系数
!ca=0.2 ! 齿顶圆角系数
alphak=20 ! 分度圆压力角
xx1=0.5426 ! 1号轮变位系数
a2=5.25 ! 齿轮1/2中心距(由xx1/2无法显式计算a2,故a2替代xx2作为变量)
pi=acos(-1)
alphab=acos((z1+z2)*m*cos(20*pi/180)/(2*a2))! 实际啮合角
invafb=tan(alphab)-alphab
inva=tan(20*pi/180)-20*pi/180 ! 渐开线函数值(分度圆上展角)
xxx=(z1+z2)*(invafb-inva)/(2*tan(20*pi/180))! 齿轮1/2的总变位系数
xx2=xxx-xx1
xx3=0.542561273-xxx-xx2 ! 3号齿轮变位系数随动
! ---------- 计算齿1齿廓线参数 ----------
ang1=360/z1 ! 单齿圆心角
halfang1=ang1/2
r1=m*z1/2 ! 分度圆半径
rf1=(z1-2*ha-2*c+2*xx1)*m/2 ! 齿根圆半径
ra1=(z1+2*ha+2*xx1)*m/2 ! 齿顶圆半径
rff=m*cf ! 齿根倒角半径
!raf=m*ca ! 齿顶倒角半径
*AFUN,DEG ! 角度单位为度
s1=pi*m/2+2*xx1*m*tan(alphak) ! 分度圆齿厚
rb1=r1*COS(alphak) ! 基圆半径
*AFUN,RAD ! 角度单位为弧度
N=12
*DIM,x1,,N ! 定义数组存放齿廓线离散点
*DIM,y1,,N
*IF,rf1,GT,rb1,THEN ! 如果齿根圆半径大于基圆,计算齿根圆上的压力角
af1=ACOS(rb1/rf1) ! 齿根圆上的压力角
invaf1=TAN(af1)-af1 ! 齿根圆的展角
phi1=s1/r1-2*(invaf1-inva) ! 齿根圆对应的圆心角
x1(1)=rf1*COS(phi1/2) ! 计算齿根圆的坐标值
y1(1)=rf1*SIN(phi1/2)
delta1=(ra1-rf1)/(N-1) ! 计算任意圆上的半径增量
*DO,i,1,N-1 ! 计算任意圆上的坐标
ri=rf1+i*delta1
j=i+1
ai=ACOS(rb1/ri) ! 计算任意圆上的压力角
invai=TAN(ai)-ai
phi=s1/r1-2*(invai-inva) ! 任意圆对应齿厚的圆心角
x1(j)=ri*COS(phi/2) ! 任意圆对应渐开线的X坐标值
y1(j)=ri*SIN(phi/2) ! 任意圆对应渐开线的Y坐标值
*ENDDO
*ELSE ! 齿根圆小于基圆,取af=0
phi1=s1/r1+2*inva ! 计算齿根圆齿厚对应的圆心角
x1(1)=rf1*COS(phi1/2) ! 计算齿根圆的坐标值
y1(1)=rf1*SIN(phi1/2)
x1(2)=rb1*COS(phi1/2) ! 计算基圆的坐标值
y1(2)=rb1*SIN(phi1/2)
delta1=(ra1-rb1)/(N-2)
*DO,i,1,N-2 ! 计算任意圆上的坐标
ri=rb1+i*delta1
j=i+2
ai=ACOS(rb1/ri) ! 计算任意圆上的压力角
invai=TAN(ai)-ai
phi=s1/r1-2*(invai-inva) ! 任意圆对应齿厚的圆心角
x1(j)=ri*COS(phi/2) ! 任意圆对应曲线的X坐标值
y1(j)=ri*SIN(phi/2) ! 任意圆对应曲线的Y坐标值
*ENDDO
*ENDIF
! ---------- 计算齿2齿廓线参数 ----------
ang2=360/z2 ! 单齿圆心角
halfang2=ang2/2
r2=m*z2/2 ! 分度圆半径
rf2=(z2-2*ha-2*c+2*xx2)*m/2 ! 齿根圆半径
ra2=(z2+2*ha+2*xx2)*m/2 ! 齿顶圆半径
*AFUN,DEG ! 角度单位为度
s2=pi*m/2+2*xx2*m*tan(alphak) ! 分度圆齿厚
rb2=r2*COS(alphak) ! 基圆半径
*AFUN,RAD ! 角度单位为弧度
N=12
*DIM,x2,,N ! 定义数组存放齿廓线离散点
*DIM,y2,,N
*IF,rf2,GT,rb2,THEN ! 如果齿根圆半径大于基圆,计算齿根圆上的压力角
af2=ACOS(rb2/rf2) ! 齿根圆上的压力角
invaf2=TAN(af2)-af2 ! 齿根圆的展角
phi2=s2/r2-2*(invaf2-inva) ! 齿根圆对应的圆心角
x2(1)=rf2*COS(phi2/2) ! 计算齿根圆的坐标值
y2(1)=rf2*SIN(phi2/2)
delta2=(ra2-rf2)/(N-1) ! 计算任意圆上的半径增量
*DO,i,1,N-1 ! 计算任意圆上的坐标
ri=rf2+i*delta2
j=i+1
ai=ACOS(rb2/ri) ! 计算任意圆上的压力角
invai=TAN(ai)-ai
phi=s2/r2-2*(invai-inva) ! 任意圆对应齿厚的圆心角
x2(j)=ri*COS(phi/2) ! 任意圆对应渐开线的X坐标值
y2(j)=ri*SIN(phi/2) ! 任意圆对应渐开线的Y坐标值
*ENDDO
*ELSE ! 齿根圆小于基圆,取af=0
phi2=s2/r2+2*inva ! 计算齿根圆齿厚对应的圆心角
x2(1)=rf2*COS(phi2/2) ! 计算齿根圆的坐标值
y2(1)=rf2*SIN(phi2/2)
x2(2)=rb2*COS(phi2/2) ! 计算基圆的坐标值
y2(2)=rb2*SIN(phi2/2)
delta2=(ra2-rb2)/(N-2)
*DO,i,1,N-2 ! 计算任意圆上的坐标
ri=rb2+i*delta2
j=i+2
ai=ACOS(rb2/ri) ! 计算任意圆上的压力角
invai=TAN(ai)-ai
phi=s2/r2-2*(invai-inva) ! 任意圆对应齿厚的圆心角
x2(j)=ri*COS(phi/2) ! 任意圆对应曲线的X坐标值
y2(j)=ri*SIN(phi/2) ! 任意圆对应曲线的Y坐标值
*ENDDO
*ENDIF
! ---------- 计算齿3齿廓线参数 ----------
ang3=360/z3 ! 单齿圆心角
halfang3=ang3/2
r3=m*z3/2 ! 分度圆半径
rf3=(z3-2*ha-2*c+2*xx3)*m/2 ! 齿根圆半径
ra3=(z3+2*ha+2*xx3)*m/2 ! 齿顶圆半径
*AFUN,DEG ! 角度单位为度
s3=pi*m/2+2*xx3*m*tan(alphak) ! 分度圆齿厚
rb3=r3*COS(alphak) ! 基圆半径
*AFUN,RAD ! 角度单位为弧度
N=12
*DIM,x3,,N ! 定义数组存放齿廓线离散点
*DIM,y3,,N
*IF,rf3,GT,rb3,THEN ! 如果齿根圆半径大于基圆,计算齿根圆上的压力角
af3=ACOS(rb3/rf3) ! 齿根圆上的压力角
invaf3=TAN(af3)-af3 ! 齿根圆的展角
phi3=s3/r3-2*(invaf3-inva) ! 齿根圆对应的圆心角
x3(1)=rf3*COS(phi3/2) ! 计算齿根圆的坐标值
y3(1)=rf3*SIN(phi3/2)
delta3=(ra3-rf3)/(N-1) ! 计算任意圆上的半径增量
*DO,i,1,N-1 ! 计算任意圆上的坐标
ri=rf3+i*delta3
j=i+1
ai=ACOS(rb3/ri) ! 计算任意圆上的压力角
invai=TAN(ai)-ai
phi=s3/r3-2*(invai-inva) ! 任意圆对应齿厚的圆心角
x3(j)=ri*COS(phi/2) ! 任意圆对应渐开线的X坐标值
y3(j)=ri*SIN(phi/2) ! 任意圆对应渐开线的Y坐标值
*ENDDO
*ELSE ! 齿根圆小于基圆,取af=0
phi3=s3/r3+2*inva ! 计算齿根圆齿厚对应的圆心角
x3(1)=rf3*COS(phi3/2) ! 计算齿根圆的坐标值
y3(1)=rf3*SIN(phi3/2)
x3(2)=rb3*COS(phi3/2) ! 计算基圆的坐标值
y3(2)=rb3*SIN(phi3/2)
delta3=(ra3-rb3)/(N-2)
*DO,i,1,N-2 ! 计算任意圆上的坐标
ri=rb3+i*delta3
j=i+2
ai=ACOS(rb3/ri) ! 计算任意圆上的压力角
invai=TAN(ai)-ai
phi=s3/r3-2*(invai-inva) ! 任意圆对应齿厚的圆心角
x3(j)=ri*COS(phi/2) ! 任意圆对应曲线的X坐标值
y3(j)=ri*SIN(phi/2) ! 任意圆对应曲线的Y坐标值
*ENDDO
*ENDIF
!***************************前处理
/PREP7 ! 前处理
ET,1,PLANE182 ! 指定单元
ET,2,PLANE183
ET,3,SOLID185
ET,4,SOLID186
KEYOPT,3,2,3 ! Simplified enhanced strain formulation
MP,EX,1,2.0e5 ! 钢弹性模量MPa
MP,PRXY,1,0.285 ! 指定材料的泊松比
MP,DENS,1,7.85e-9 ! 钢密度t/mm3
MP,MU,1,0.1 ! 摩擦系数
! ---------- 中心点设置 ----------
N,1,0,0,d1 ! 创建主动轮中心点,编号1
N,2,a2,0,d2 ! 中间轮中心点变动
N,3,2*m*z2,0,d3
! ---------- 齿廓1建模 ----------
*AFUN,DEG
*DO,i,1,N !生成N=12个关键点
K,i,x1(i),y1(i)
*ENDDO
K,N+1,ra1,0 !生成齿顶圆上的关键点
LSTR,N,N+1 !直线连接近似替代齿顶圆圆弧
KSEL,S,KP,,1,N
BSPLIN,ALL !样条拟合生成齿廓
CYL4,,,rm1,0,rf1,halfang1 !生成一个圆环
ADELE,1 !删除生成的圆环面
LDELE,3 !删除一条与齿廓相连的线
KSEL,ALL
LSTR,14,13 !生成两条线
LSTR,15,1
KL,4,(rf1-rh1)/(rf1-rm1) !在线上按比率生成一个关键点
KL,6,(rf1-rh1)/(rf1-rm1)
LCOMB,3,6,0 !线相加
LFILLT,7,2,rff,, !对齿根进行倒角
!LFILLT,2,1,raf,, !对齿顶进行倒角
LANG,3,20,90, !由倒角半径上的一点生成一条垂线
LARC,18,19,17,rh1
LSBL,4,10,,DELETE,KEEP !通过生成的垂线对另一条线进行分割
LSBL,8,10,,DELETE,KEEP
AL,3,1,2,9 !形成一个齿形面
AL,6,7,11,10,13,9
AL,4,5,12,10
ALLSEL,ALL
NUMCMP,ALL
! ---------- 齿廓2建模 ----------
WPOFFS,1.5*M*Z2
CSYS,4
WPSTYL,,,,,,0
*GET,K_2,KP,0,NUM,MAX
*GET,L_2,LINE,0,NUM,MAX
*GET,A_2,AREA,0,NUM,MAX
*DO,i,1,N !生成N=12个关键点
K,i+K_2,x2(i),y2(i)
*ENDDO
K,K_2+N+1,ra2,0 !生成齿顶圆上的关键点
LSTR,K_2+N,K_2+N+1 !直线连接近似替代齿顶圆圆弧
KSEL,S,KP,,K_2+1,K_2+N
BSPLIN,ALL !样条拟合生成齿廓
CYL4,,,rm2,0,rf2,halfang2 !生成一个圆环
ADELE,A_2+1 !删除生成的圆环面
LDELE,L_2+3 !删除一条与齿廓相连的线
KSEL,ALL
LSTR,K_2+14,K_2+13 !生成两条线
LSTR,K_2+15,K_2+1
KL,L_2+4,(rf2-rh2)/(rf2-rm2) !在线上按比率生成一个关键点
KL,L_2+6,(rf2-rh2)/(rf2-rm2)
LCOMB,L_2+3,L_2+6,0 !线相加
LFILLT,L_2+7,L_2+2,rff,, !对齿根进行倒角
!LFILLT,L_2+2,L_2+1,raf,, !对齿顶进行倒角
LANG,L_2+3,K_2+20,90, !由倒角半径上的一点生成一条垂线
LARC,K_2+18,K_2+19,K_2+17,rh2
LSBL,L_2+4,L_2+10,,DELETE,KEEP !通过生成的垂线对另一条线进行分割
LSBL,L_2+8,L_2+10,,DELETE,KEEP
AL,L_2+3,L_2+1,L_2+2,L_2+9 !形成一个齿形面
AL,L_2+6,L_2+7,L_2+11,L_2+10,L_2+13,L_2+9
AL,L_2+4,L_2+5,L_2+12,L_2+10
ALLSEL,ALL
NUMCMP,ALL
! ---------- 齿廓3建模 ----------
WPOFFS,1.5*M*Z2
CSYS,4
WPSTYL,,,,,,0
*GET,K_3,KP,0,NUM,MAX
*GET,L_3,LINE,0,NUM,MAX
*GET,A_3,AREA,0,NUM,MAX
*DO,i,1,N !生成N=12个关键点
K,i+K_3,x3(i),y3(i)
*ENDDO
K,K_3+N+1,ra3,0 !生成齿顶圆上的关键点
LSTR,K_3+N,K_3+N+1 !直线连接近似替代齿顶圆圆弧
KSEL,S,KP,,K_3+1,K_3+N
BSPLIN,ALL !样条拟合生成齿廓
CYL4,,,rm3,0,rf3,halfang3 !生成一个圆环
ADELE,A_3+1 !删除生成的圆环面
LDELE,L_3+3 !删除一条与齿廓相连的线
KSEL,ALL
LSTR,K_3+14,K_3+13 !生成两条线
LSTR,K_3+15,K_3+1
KL,L_3+4,(rf3-rh3)/(rf3-rm3) !在线上按比率生成一个关键点
KL,L_3+6,(rf3-rh3)/(rf3-rm3)
LCOMB,L_3+3,L_3+6,0 !线相加
LFILLT,L_3+7,L_3+2,rff,, !对齿根进行倒角
!LFILLT,L_3+2,L_3+1,raf,, !对齿顶进行倒角
LANG,L_3+3,K_3+20,90, !由倒角半径上的一点生成一条垂线
LARC,K_3+18,K_3+19,K_3+17,rh3
LSBL,L_3+4,L_3+10,,DELETE,KEEP !通过生成的垂线对另一条线进行分割
LSBL,L_3+8,L_3+10,,DELETE,KEEP
AL,L_3+3,L_3+1,L_3+2,L_3+9 !形成一个齿形面
AL,L_3+6,L_3+7,L_3+11,L_3+10,L_3+13,L_3+9
AL,L_3+4,L_3+5,L_3+12,L_3+10
ALLSEL,ALL
NUMCMP,ALL
! ---------- 齿1稀疏面网格划分 ----------
LESIZE,5,,,3,,,,,1 ! 对线指定划分单元1的等分数
LESIZE,9,,,3,,,,,1
LESIZE,11,,,4,,,,,1
LESIZE,4,,,4,,,,,1
LESIZE,10,,,4,,,,,1
LESIZE,12,,,5,,,,,1
LESIZE,8,,,3,,,,,1
LESIZE,7,,,1,,,,,1
LESIZE,6,,,2,,,,,1
!LESIZE,2,,,6,,,,,1
!LESIZE,3,,,7,,,,,1
!LESIZE,8,,,1,,,,,1
!LESIZE,1,,,1,,,,,1
CSYS,1
ASEL,S,LOC,X,0,rh1
AGEN,2,ALL,,,,2*ang1,,,0
ALLSEL,ALL
ASEL,S,LOC,X,rh1,ra1
AGEN,2,ALL,,,,2*ang1,,,0
ALLSEL,ALL
ADELE,3,,,1
WPOFFS,-3*M*Z2
WPROTA,2.5*ang1
CSYS,4
WPSTYL,,,,,,1
ASEL,S,LOC,Y,-halfang1,0
ASEL,R,LOC,X,0,rf1+rff
TYPE,1
MSHAPE,0,2D ! 指定为2D结构四边形网格
MSHKEY,2 ! 指定为映射+自由网格
AMESH,ALL
ASEL,S,LOC,Y,-halfang1,0
ASEL,R,LOC,X,rf1+rff,ra1
TYPE,2
ESIZE,0.1
AMESH,ALL
ASEL,S,LOC,Y,-halfang1,0
WPSTYL,,,,,,0
ARSYM,Y,ALL,,,,0,0
ALLSEL,ALL
WPSTYL,,,,,,1
ASEL,S,LOC,X,0,rh1
AGEN,Z1,ALL,,,,ang1,,,0
ALLSEL,ALL
ASEL,S,LOC,X,rh1,ra1
ASEL,R,LOC,Y,-halfang1,halfang1
AGEN,Z1-3,ALL,,,,ang1,,,0
ALLSEL,ALL
! ---------- 齿1稠密面 ----------
WPROTA,-2*ang1
ASEL,S,LOC,Y,-halfang1,0
ASEL,R,LOC,X,rh1,ra1
!LESIZE,3,,,10,,,,,1
LESIZE,12,,,7,,,,,1
LESIZE,7,,,1,,,,,1
LESIZE,6,,,6,,,,,1
!LESIZE,2,,,23,,,,,1
!LESIZE,8,,,4,,,,,1
!LESIZE,1,,,3,,,,,1
LESIZE,8,,,6,0.6,,,,1
LESIZE,9,,,4,,,,,1
LESIZE,10,,,6,,,,,1
TYPE,2
ESIZE,0.03
AMESH,ALL
WPSTYL,,,,,,0
ARSYM,Y,ALL,,,,0,0
WPSTYL,,,,,,1
AGEN,2,ALL,,,,-ang1,,,0
ASEL,R,LOC,Y,-halfang1,halfang1
AGEN,2,ALL,,,,ang1,,,0
ALLSEL,ALL
! ---------- 齿2稀疏面网格划分 ----------
WPROTA,-halfang1
WPOFFS,1.5*M*Z2
LESIZE,17,,,3,,,,,1 ! 对线指定划分单元1的等分数
LESIZE,21,,,3,,,,,1
LESIZE,23,,,4,,,,,1
LESIZE,16,,,4,,,,,1
LESIZE,22,,,4,,,,,1
LESIZE,24,,,5,,,,,1
LESIZE,20,,,3,,,,,1
LESIZE,19,,,1,,,,,1
LESIZE,18,,,2,,,,,1
!LESIZE,14,,,6,0.8,,,,1
!LESIZE,15,,,7,,,,,1
!LESIZE,21,,,1,,,,,1
!LESIZE,13,,,1,,,,,1
WPSTYL,,,,,,1
ASEL,S,LOC,X,0,rh2
AGEN,2,ALL,,,,ang2,,,0
ALLSEL,ALL
ASEL,S,LOC,X,rh2,ra2
AGEN,2,ALL,,,,ang2,,,0
ALLSEL,ALL
ADELE,6,,,1
WPROTA,1.5*ang2
CSYS,4
WPSTYL,,,,,,1
ASEL,S,LOC,Y,-halfang2,0
ASEL,R,LOC,X,0,rf2+rff
TYPE,1
MSHAPE,0,2D ! 指定为2D结构四边形网格
MSHKEY,2 ! 指定为映射+自由网格
AMESH,ALL
ASEL,S,LOC,Y,-halfang2,0
ASEL,R,LOC,X,rf2+rff,ra2
TYPE,2
ESIZE,0.2
AMESH,ALL
ASEL,S,LOC,Y,-halfang2,0
WPSTYL,,,,,,0
ARSYM,Y,ALL,,,,0,0
ALLSEL,ALL
WPSTYL,,,,,,1
ASEL,S,LOC,X,0,rh2
AGEN,Z2,ALL,,,,ang2,,,0
ALLSEL,ALL
ASEL,S,LOC,X,rh2,ra2
ASEL,R,LOC,Y,-halfang2,halfang2
AGEN,Z2-3,ALL,,,,ang2,,,0
ALLSEL,ALL
! ---------- 齿2稠密面 ----------
WPROTA,-ang2
ASEL,S,LOC,Y,-halfang2,0
ASEL,R,LOC,X,rh2,ra2
!LESIZE,15,,,11,,,,,1
LESIZE,24,,,4,,,,,1
LESIZE,19,,,3,,,,,1
LESIZE,18,,,5,,,,,1
!LESIZE,14,,,22,,,,,1
!LESIZE,21,,,4,,,,,1
!LESIZE,13,,,3,,,,,1
LESIZE,20,,,6,0.6,,,,1
LESIZE,21,,,4,,,,,1
LESIZE,22,,,4,,,,,1
TYPE,2
ESIZE,0.03
MSHAPE,0,2D
MSHKEY,0
AMESH,ALL
WPSTYL,,,,,,0
ARSYM,Y,ALL,,,,0,0
WPSTYL,,,,,,1
AGEN,3,ALL,,,,-ang2,,,0
ALLSEL,ALL
WPROTA,-halfang2 ! 调整面
ASEL,S,LOC,X,rh2,ra2
ASEL,R,LOC,Y,180-2*ang2,180+ang2
ACLEAR,ALL
ADELE,ALL,,,1
ALLSEL,ALL
ASEL,S,LOC,X,rh2,ra2
ASEL,R,LOC,Y,-halfang2,halfang2
AGEN,2,ALL,,,,180-1.5*ang2,,,0
ALLSEL,ALL
ASEL,S,LOC,X,rh2,ra2
ASEL,R,LOC,Y,180-2*ang2,180-ang2
AGEN,3,ALL,,,,ang2,,,0
ALLSEL,ALL
! ---------- 齿3稀疏面网格划分 ----------
WPOFFS,1.5*M*Z2
LESIZE,29,,,3,,,,,1 ! 对线指定划分单元1的等分数
LESIZE,33,,,3,,,,,1
LESIZE,35,,,4,,,,,1
LESIZE,28,,,4,,,,,1
LESIZE,34,,,4,,,,,1
LESIZE,36,,,5,,,,,1
LESIZE,32,,,3,,,,,1
LESIZE,31,,,1,,,,,1
LESIZE,30,,,2,,,,,1
!LESIZE,26,,,6,0.8,,,,1
!LESIZE,27,,,7,,,,,1
!LESIZE,34,,,1,,,,,1
!LESIZE,25,,,1,,,,,1
ASEL,S,LOC,X,0,rh3
AGEN,2,ALL,,,,180-1.5*ang3,,,0
ALLSEL,ALL
ASEL,S,LOC,X,rh3,ra3
AGEN,2,ALL,,,,180-1.5*ang3,,,0
ALLSEL,ALL
ADELE,9,,,1
WPROTA,180-1.5*ang3
ASEL,S,LOC,Y,0,halfang3
ASEL,R,LOC,X,0,rf3+rff
TYPE,1
ESIZE,0.2
MSHAPE,0,2D ! 指定为2D结构四边形网格
MSHKEY,2 ! 指定为映射+自由网格
AMESH,ALL
ASEL,S,LOC,Y,0,halfang3
ASEL,R,LOC,X,rf3+rff,ra3
TYPE,2
ESIZE,0.2
AMESH,ALL
ASEL,S,LOC,Y,0,halfang3
ASEL,R,LOC,X,0,ra3
WPSTYL,,,,,,0
ARSYM,Y,ALL,,,,0,0
ALLSEL,ALL
WPSTYL,,,,,,1
ASEL,S,LOC,X,0,rh3
AGEN,Z3,ALL,,,,ang3,,,0
ALLSEL,ALL
WPROTA,1.5*ang3-180
ASEL,S,LOC,X,rh3,ra3
ASEL,R,LOC,Y,0,halfang3
AGEN,,ALL,,,,180-halfang3,,,,1
ALLSEL,ALL
! ---------- 齿3稠密面 ----------
ASEL,S,LOC,Y,180-halfang3,180
ASEL,R,LOC,X,rh3,ra3
!LESIZE,27,,,11,,,,,1
LESIZE,36,,,6,,,,,1
LESIZE,31,0.03,,,,,,,1
LESIZE,30,,,5,,,,,1
!LESIZE,26,,,22,,,,,1
!LESIZE,34,,,4,,,,,1
!LESIZE,25,,,3,,,,,1
LESIZE,32,,,6,0.6,,,,1
LESIZE,33,,,4,,,,,1
LESIZE,34,,,6,,,,,1
TYPE,2
ESIZE,0.03
AMESH,ALL
WPROTA,180-halfang3
WPSTYL,,,,,,0
ARSYM,Y,ALL,,,,0,0
WPSTYL,,,,,,1
AGEN,3,ALL,,,,ang3,,,0
ALLSEL,ALL
ASEL,S,LOC,X,rh3,ra3
ASEL,R,LOC,Y,-1.5*ang3,-halfang3
AGEN,Z3-3,ALL,,,,-ang3,,,0
ALLSEL,ALL
WPROTA,halfang3-180
! ---------- 分区建3D齿轮 ----------
! -MERGE-
CSYS,1 ! 各面积区域分别建立集合
ASEL,S,LOC,X,rh1,ra1
ASEL,R,LOC,Y,-ang1,2*ang1
CM,A_CON1,AREA
ALLSEL,ALL
ASEL,S,LOC,X,0,rh1
CM,A_HUB1,AREA
CSYS,4
WPSTYL,,,,,,1
ASEL,S,LOC,X,rh3,ra3
ASEL,R,LOC,Y,180-ang3,180+2*ang3
CM,A_CON3,AREA
ALLSEL,ALL
ASEL,S,LOC,X,0,rh3
CM,A_HUB3,AREA
ALLSEL,ALL
WPOFFS,-1.5*m*z2
ASEL,S,LOC,Y,-2*ang2,ang2
ASEL,A,LOC,Y,180-2*ang2,180+ang2
ASEL,R,LOC,X,rh2,ra2
CM,A_CON2,AREA
ALLSEL,ALL
ASEL,S,LOC,X,0,rh2
CM,A_HUB2,AREA
ALLSEL,ALL
ASEL,ALL ! 齿轮外圈非接触区/齿轮2内轮毂各自合并节点/线/关键点/单元
CMSEL,U,A_CON1
CMSEL,U,A_CON2
CMSEL,U,A_CON3
CMSEL,U,A_HUB1
CMSEL,U,A_HUB2
CMSEL,U,A_HUB3
NSLA,S,1
LSLA,S,1
KSLL,S,1
NUMMRG,ALL,,,,LOW
NUMCMP,ALL
ALLSEL,ALL
CMSEL,S,A_CON1 ! 齿轮外圈接触区/齿轮2内轮毂各自合并节点/线/关键点/单元
CMSEL,A,A_CON2
CMSEL,A,A_CON3
NSLA,S,1
LSLA,S,1
KSLL,S,1
NUMMRG,ALL,,,,LOW
NUMCMP,ALL
ALLSEL,ALL
CMSEL,S,A_HUB1 ! 轮毂区/齿轮2内轮毂各自合并节点/线/关键点/单元
CMSEL,A,A_HUB2
CMSEL,A,A_HUB3
NSLA,S,1
LSLA,S,1
KSLL,S,1
NUMMRG,ALL,,,,LOW
NUMCMP,ALL
ALLSEL,ALL
NUMCMP,ALL
! -EXTENSION-
CMSEL,U,A_CON1 ! 非接触区拉伸
CMSEL,U,A_CON2
CMSEL,U,A_CON3
CMSEL,U,A_HUB1
CMSEL,U,A_HUB2
CMSEL,U,A_HUB3
CSYS,0
TYPE,3
ESIZE,,6
VEXT,ALL,,,0,0,b,0,0,0
VSEL,ALL
CM,V_NOC,VOLU
ALLSEL,ALL
CMSEL,S,A_HUB1
CMSEL,A,A_HUB3
VEXT,ALL,,,0,0,bh13,0,0,0
VSEL,ALL
CMSEL,U,V_NOC
CM,V_HUB13,VOLU
ALLSEL,ALL
CMSEL,S,A_HUB2
VEXT,ALL,,,0,0,bh2,0,0,0
VSEL,ALL
CMSEL,U,V_NOC
CMSEL,U,V_HUB13
CM,V_HUB2,VOLU
ALLSEL,ALL
CMSEL,S,A_CON1 ! 接触区拉伸
CMSEL,A,A_CON2
CMSEL,A,A_CON3
TYPE,4
ESIZE,,18
VEXT,ALL,,,0,0,b,0,0,0
VSEL,ALL
CMSEL,U,V_NOC
CMSEL,U,V_HUB13
CMSEL,U,V_HUB2
CM,V_CON,VOLU
ALLSEL,ALL
ACLEAR,ALL ! 删除2D单元
ALLSEL,ALL
NUMCMP,ALL ! 压缩实体编号
CMSEL,S,V_CON ! 将齿轮移至中心
CMSEL,A,V_NOC
VGEN,,ALL,,,,,-b/2,,,1
CMSEL,S,V_HUB13
VGEN,,ALL,,,,,-bh13/2+d1,,,1
CMSEL,S,V_HUB2
VGEN,,ALL,,,,,-bh2/2+d2,,,1
ALLSEL,ALL
! ---------- 选择接触面1 ----------
CSYS,1 ! 切换到柱坐标系
ASEL,s,loc,x,rf1,ra1
*DO,i,1,z1
ASEL,u,loc,y,i*ang1
*ENDDO
ASEL,u,loc,x,cos(halfang1)*(rf1+rff),rf1+rff
ASEL,u,loc,z,b/2
ASEL,u,loc,z,-b/2
ASEL,a,loc,x,ra1
ASEL,u,loc,x,rf1
ASEL,R,LOC,Y,-ang1,2*ang1
CM,CONTA_SURF1,AREA
ALLSEL,ALL
! ---------- 选择目标面1 ----------
CSYS,4
WPSTYL,,,,,,1
ASEL,s,loc,x,rf2,ra2
*DO,i,1,z2
ASEL,u,loc,y,i*ang2
*ENDDO
ASEL,u,loc,x,cos(halfang2)*(rf2+rff),rf2+rff
ASEL,U,LOC,Z,-B/2
ASEL,U,LOC,Z,B/2
ASEL,a,loc,x,ra2
ASEL,u,loc,x,rf2
ASEL,R,LOC,Y,180-2*ang2,180+ang2
CM,TARGE_SURF1,AREA
ALLSEL,ALL
! ---------- 选择接触面2 ----------
ASEL,s,loc,x,rf2,ra2
*DO,i,1,z2
ASEL,u,loc,y,i*ang2
*ENDDO
ASEL,u,loc,x,cos(halfang2)*(rf2+rff),rf2+rff
ASEL,U,LOC,Z,-B/2
ASEL,U,LOC,Z,B/2
ASEL,a,loc,x,ra2
ASEL,u,loc,x,rf2
ASEL,R,LOC,Y,-2*ang2,ang2
CM,CONTA_SURF2,AREA
ALLSEL,ALL
! ---------- 选择目标面2 ----------
WPOFFS,1.5*m*z2
WPSTYL,,,,,,1
ASEL,s,loc,x,rf3,ra3
*DO,i,1,z3
ASEL,u,loc,y,i*ang3
*ENDDO
ASEL,u,loc,x,cos(halfang3)*(rf3+rff),rf3+rff
ASEL,u,loc,z,b/2
ASEL,u,loc,z,-b/2
ASEL,a,loc,x,ra3
ASEL,u,loc,x,rf3
ASEL,R,LOC,Y,180-ang3,180+2*ang3
CM,TARGE_SURF2,AREA
ALLSEL,ALL
! ---------- 从动轮复位 ----------
WPOFFS,-1.5*m*z2
VSEL,s,loc,x,0,ra2
WPSTYL,,,,,,0
VGEN,,ALL,,,-1.5*m*z2+a2,,,,,1 ! 齿轮2移至啮合位
WPOFFS,1.5*m*z2
WPSTYL,,,,,,1
VSEL,s,loc,x,0,ra3
WPSTYL,,,,,,0
VGEN,,ALL,,,-m*z2,,,,,1 ! 齿轮3移至啮合位
WPOFFS,-3*m*z2
ALLSEL,ALL
! ---------- 齿轮1接触区与其他区域BOND1 ----------
ET,5,174 ! Define CONTACT175 element
ET,6,170 ! Define TARGET170 element
KEYOPT,5,12,5 ! Bonded Contact
KEYOPT,5,2,2 ! MPC style contact
CMSEL,S,V_CON
NSLV,S,1
CM,N_CON,NODE
ALLSEL,ALL
CSYS,1
NSEL,S,LOC,Y,-ang1
NSEL,R,LOC,X,rh1,ra1
CMSEL,R,N_CON
TYPE,5
REAL,2
ESURF
ALLSEL,ALL
CMSEL,S,V_NOC ! 目标面
NSLV,S,1
CM,N_NOC,NODE
ALLSEL,ALL
NSEL,S,LOC,Y,-ang1
NSEL,R,LOC,X,rh1,ra1
CMSEL,R,N_NOC
ESLN
TYPE,6
REAL,2
ESURF ! 创建轴与凸轮之间的接触对
ALLSEL,ALL
! ---------- 齿轮1接触区与其他区域BOND2 ----------
CMSEL,S,V_CON
NSLV,S,1
CM,N_CON,NODE
ALLSEL,ALL
CSYS,1
NSEL,S,LOC,Y,2*ang1
NSEL,R,LOC,X,rh1,ra1
CMSEL,R,N_CON
TYPE,5
REAL,3
ESURF
ALLSEL,ALL
CMSEL,S,V_NOC ! 目标面
NSLV,S,1
CM,N_NOC,NODE
ALLSEL,ALL
NSEL,S,LOC,Y,2*ang1
NSEL,R,LOC,X,rh1,ra1
CMSEL,R,N_NOC
ESLN
TYPE,6
REAL,3
ESURF ! 创建轴与凸轮之间的接触对
ALLSEL,ALL
! ---------- 齿轮1内外圈BOND ----------
CMSEL,S,V_HUB13
NSLV,S,1
CM,N_HUB13,NODE
NSEL,S,LOC,X,rh1
CMSEL,U,N_HUB13
TYPE,5
REAL,4
ESURF
ALLSEL,ALL
NSEL,S,LOC,X,rh1 ! 选择齿轮外圈轮毂侧节点
CMSEL,R,N_HUB13
ESLN
TYPE,6
REAL,4
ESURF ! 创建轮毂与齿轮之间的接触对
ALLSEL,ALL
! ---------- 齿轮2接触区与非接触区BOND1 ----------
WPOFFS,a2
CSYS,4
WPSTYL,,,,,,1
ASEL,S,LOC,Y,180+ang2
ASEL,R,LOC,X,rh2,rf2
ASEL,U,LOC,X,rf2
NSLA,S,1
CMSEL,R,N_CON
TYPE,5
REAL,5
ESURF
ALLSEL,ALL
ASEL,S,LOC,Y,180+ang2 ! 目标面
ASEL,R,LOC,X,rh2,rf2
ASEL,U,LOC,X,rf2
NSLA,S,1
CMSEL,R,N_NOC
ESLN
TYPE,6
REAL,5
ESURF
ALLSEL,ALL
! ---------- 齿轮2接触区与非接触区BOND2 ----------
CSYS,4
WPSTYL,,,,,,1
ASEL,S,LOC,Y,180-2*ang2
ASEL,R,LOC,X,rh2,rf2
ASEL,U,LOC,X,rf2
NSLA,S,1
CMSEL,R,N_CON
TYPE,5
REAL,6
ESURF
ALLSEL,ALL
ASEL,S,LOC,Y,180-2*ang2 ! 目标面
ASEL,R,LOC,X,rh2,rf2
ASEL,U,LOC,X,rf2
NSLA,S,1
CMSEL,R,N_NOC
ESLN
TYPE,6
REAL,6
ESURF
ALLSEL,ALL
! ---------- 齿轮2接触区与非接触区BOND3 ----------
ASEL,S,LOC,Y,ang2
ASEL,R,LOC,X,rh2,ra2
ASEL,U,LOC,X,ra2
NSLA,S,1
CMSEL,R,N_CON
TYPE,5
REAL,7
ESURF
ALLSEL,ALL
ASEL,S,LOC,Y,ang2 ! 目标面
ASEL,R,LOC,X,rh2,ra2
ASEL,U,LOC,X,ra2
NSLA,S,1
CMSEL,R,N_NOC
ESLN
TYPE,6
REAL,7
ESURF
ALLSEL,ALL
! ---------- 齿轮2接触区与非接触区BOND4 ----------
ASEL,S,LOC,Y,-2*ang2
ASEL,R,LOC,X,rh2,ra2
ASEL,U,LOC,X,ra2
NSLA,S,1
CMSEL,R,N_CON
TYPE,5
REAL,8
ESURF
ALLSEL,ALL
ASEL,S,LOC,Y,-2*ang2 ! 目标面
ASEL,R,LOC,X,rh2,ra2
ASEL,U,LOC,X,ra2
NSLA,S,1
CMSEL,R,N_NOC
ESLN
TYPE,6
REAL,8
ESURF
ALLSEL,ALL
! ---------- 齿轮2内外圈BOND ----------
CMSEL,S,V_HUB2
NSLV,S,1
CM,N_HUB2,NODE
NSEL,S,LOC,X,rh2
CMSEL,U,N_HUB2
TYPE,5
REAL,9
ESURF
ALLSEL,ALL
NSEL,S,LOC,X,rh2 ! 选择齿轮外圈轮毂侧节点
CMSEL,R,N_HUB2
ESLN
TYPE,6
REAL,9
ESURF ! 创建轮毂与齿轮之间的接触对
ALLSEL,ALL
! ---------- 齿轮3接触区与其他区域BOND1 ----------
WPOFFS,-a2+2*m*z2
ASEL,S,LOC,Y,180+2*ang3
ASEL,R,LOC,X,rh3,rf3
ASEL,U,LOC,X,rf3
NSLA,S,1
CMSEL,R,N_CON
TYPE,5
REAL,10
ESURF
ALLSEL,ALL
ASEL,S,LOC,Y,180+2*ang3 ! 目标面
ASEL,R,LOC,X,rh3,rf3
ASEL,U,LOC,X,rf3
NSLA,S,1
CMSEL,R,N_NOC
ESLN
TYPE,6
REAL,10
ESURF
ALLSEL,ALL
! ---------- 齿轮3接触区与其他区域BOND2 ----------
ASEL,S,LOC,Y,180-ang3
ASEL,R,LOC,X,rh3,rf3
ASEL,U,LOC,X,rf3
NSLA,S,1
CMSEL,R,N_CON
TYPE,5
REAL,11
ESURF
ALLSEL,ALL
ASEL,S,LOC,Y,180-ang3 ! 目标面
ASEL,R,LOC,X,rh3,rf3
ASEL,U,LOC,X,rf3
NSLA,S,1
CMSEL,R,N_NOC
ESLN
TYPE,6
REAL,11
ESURF
ALLSEL,ALL
! ---------- 齿轮3内外圈BOND ----------
CMSEL,S,N_HUB13
NSEL,S,LOC,X,rh3
CMSEL,U,N_HUB13
TYPE,5
REAL,12
ESURF
ALLSEL,ALL
NSEL,S,LOC,X,rh3 ! 选择齿轮外圈轮毂侧节点
CMSEL,R,N_HUB13
ESLN
TYPE,6
REAL,12
ESURF ! 创建轮毂与齿轮之间的接触对
ALLSEL,ALL
! ---------- 摩擦接触对1 ----------
ET,7,CONTA174 ! 接触单元
ET,8,TARGE170 ! 目标单元
KEYOPT,7,9,0 ! default,包括初始穿透/间隙
KEYOPT,7,12,0 ! 标准单向接触模型
KEYOPT,7,5,3 ! 自动关闭初始间隙或减小初始穿透
KEYOPT,7,2,0 ! 增强拉格朗日default
KEYOPT,7,10,2 ! 子步自动更新FKN和FKT
CMSEL,S,TARGE_SURF1 ! -创建接触面单元- (主动轮和从动轮互换,这里反选targe面)
NSLA,s,1 ! 选择依附于面的节点
TYPE,7
REAL,13
ESURF
ALLSEL,ALL
CMSEL,S,CONTA_SURF1 ! -创建目标面单元-
NSLA,s,1 ! 选择依附于面的节点
ESLN,s,0
TYPE,8
REAL,13
ESURF
ALLSEL,ALL
! ---------- 摩擦接触对2 ----------
ET,9,CONTA174 ! 接触单元
ET,10,TARGE170 ! 目标单元
KEYOPT,9,9,0 ! default,包括初始穿透/间隙
KEYOPT,9,12,0 ! 标准单向接触模型
KEYOPT,9,5,3 ! 自动关闭初始间隙或减小初始穿透
KEYOPT,9,2,0 ! 增强拉格朗日default
KEYOPT,9,10,2 ! 子步自动更新FKN和FKT
CMSEL,S,TARGE_SURF2 ! -创建接触面单元-
NSLA,s,1 ! 选择依附于面的节点
TYPE,9
REAL,14
ESURF
ALLSEL,ALL
CMSEL,S,CONTA_SURF2 ! -创建目标面单元-
NSLA,s,1 ! 选择依附于面的节点
ESLN,s,0
TYPE,10
REAL,14
ESURF
ALLSEL,ALL
! -------- 创建齿轮1旋转副 --------
ET,11,184 ! Defines MPC184 element
KEYOPT,11,1,6 ! Selects revolute joint type
KEYOPT,11,4,1 ! Selects Z axis as revolute axis
LOCAL,11,0,0,0,0 ! Creates local coordinate system such that axis of shaft becomes the local Z axis
ET,12,175 ! Defines CONTACT175 element
ET,13,170 ! Defines TARGET170 element
KEYOPT,12,12,5 ! Bonded Contact
KEYOPT,12,4,2 ! Rigid CERIG style load
KEYOPT,12,2,2 ! MPC style contact
KEYOPT,13,2,1 ! Don't fix the pilot node
KEYOPT,13,4,0 ! Activates all DOF's due to CERIG style
CSYS,1
NSEL,S,loc,x,rm1
TYPE,12 ! Creating bonded contact pairs for defining joints
REAL,15
NN=0
*GET,NUMN,NODE,,COUNT
*DO,I,1,NUMN
NN=NDNEXT(NN)
E,NN
*ENDDO
ALLSEL,ALL
TYPE,13
REAL,15
TSHAP,PILO ! 定义目标单元形状为导向节点
*GET,N_PILO,ELEM,0,NUM,MAX
N_PILO1=N_PILO+1
EN,N_PILO1,1 ! 创建目标单元,使用导向节点1
TSHAP ! 重置目标单元形状为默认
CSYS,11
NROTAT,1 ! 对齐点坐标系与当前激活坐标系(冗余)
CSYS,0
ALLSEL,ALL
SECTYPE,1,JOINT,REVO, ! Creates section type for joint elements
SECJOINT,,11,11 ! Defines local coordinate systems at joint element nodes
REAL,16
TYPE,11
SECNUM,1
*GET,N_REVO,ELEM,0,NUM,MAX
N_REVO1=N_REVO+1
EN,N_REVO1,,1 ! Defines body to ground joint (revolute joint)
ALLSEL,ALL
! -------- 创建齿轮2旋转副 --------
ET,14,184 ! Defines MPC184 element
KEYOPT,14,1,6 ! Selects revolute joint type
KEYOPT,14,4,1 ! Selects Z axis as revolute axis
LOCAL,12,0,a2,0,0 ! Creates local coordinate system at one end of the shaft
ET,15,175 ! Defines CONTACT175 element
ET,16,170 ! Defines TARGET170 element
KEYOPT,15,12,5 ! Bonded Contact
KEYOPT,15,4,2 ! Rigid CERIG style load
KEYOPT,15,2,2 ! MPC style contact
KEYOPT,16,2,1 ! Don't fix the pilot node
KEYOPT,16,4,0 ! Activates all DOF's due to CERIG style
WPOFFS,-2*m*z2+a2
CSYS,4
WPSTYL,,,,,,1
NSEL,S,loc,x,rm2
TYPE,15 ! Creating bonded contact pairs for defining joints
REAL,17
NN=0
*GET,NUMN,NODE,,COUNT
*DO,I,1,NUMN
NN=NDNEXT(NN)
E,NN
*ENDDO
ALLSEL,ALL
TYPE,16
REAL,17
TSHAP,PILO ! 定义目标单元形状为导向节点
*GET,N_PILO,ELEM,0,NUM,MAX
N_PILO2=N_PILO+1
EN,N_PILO2,2 ! 创建目标单元999,使用导向节点99
TSHAP ! 重置目标单元形状为默认
CSYS,12
NROTAT,2 ! 对齐点坐标系与当前激活坐标系
CSYS,4
ALLSEL,ALL
SECTYPE,2,JOINT,REVO, ! Creates section type for joint elements
SECJOINT,,12,12 ! Defines local coordinate systems at joint element nodes
REAL,18
TYPE,14
SECNUM,2
*GET,N_REVO,ELEM,0,NUM,MAX
N_REVO2=N_REVO+1
EN,N_REVO2,,2 ! Defines body to ground joint (revolute joint)
ALLSEL,ALL
! -------- 创建齿轮3旋转副 --------
ET,17,184 ! Defines MPC184 element
KEYOPT,17,1,6 ! Selects revolute joint type
KEYOPT,17,4,1 ! Selects Z axis as revolute axis
LOCAL,13,0,2*m*z2,0,0 ! Creates local coordinate system at one end of the shaft
ET,18,175 ! Defines CONTACT175 element
ET,19,170 ! Defines TARGET170 element
KEYOPT,18,12,5 ! Bonded Contact
KEYOPT,18,4,2 ! Rigid CERIG style load
KEYOPT,18,2,2 ! MPC style contact
KEYOPT,19,2,1 ! Don't fix the pilot node
KEYOPT,19,4,0 ! Activates all DOF's due to CERIG style
WPOFFS,-a2+2*m*z2
CSYS,4
WPSTYL,,,,,,1
NSEL,S,loc,x,rm3
TYPE,18 ! Creating bonded contact pairs for defining joints
REAL,19
NN=0
*GET,NUMN,NODE,,COUNT
*DO,I,1,NUMN
NN=NDNEXT(NN)
E,NN
*ENDDO
ALLSEL,ALL
TYPE,19
REAL,19
TSHAP,PILO ! 定义目标单元形状为导向节点
*GET,N_PILO,ELEM,0,NUM,MAX
N_PILO3=N_PILO+1
EN,N_PILO3,3 ! 创建目标单元N_PILO3,使用导向节点3
TSHAP ! 重置目标单元形状为默认
CSYS,13
NROTAT,3 ! 对齐点坐标系与当前激活坐标系
CSYS,4
ALLSEL,ALL
SECTYPE,3,JOINT,REVO, ! Creates section type for joint elements
SECJOINT,,13,13 ! Defines local coordinate systems at joint element nodes
REAL,20
TYPE,17
SECNUM,3
*GET,N_REVO,ELEM,0,NUM,MAX
N_REVO3=N_REVO+1
EN,N_REVO3,,3 ! Defines body to ground joint (revolute joint)
ALLSEL,ALL
! -------- 载荷与约束 ----------
*AFUN,RAD
ESEL,S,TYPE,,17 ! Selects joint elements
DJ,ALL,OMGZ,-50*PI ! Applies load on joint elements
ALLSEL,ALL
ESEL,S,TYPE,,11 ! Selects joint elements
FJ,ALL,MZ,20 ! Applies load on joint elements
ALLSEL,ALL
CNCHECK,ADJUST
CNCHECK,AUTO
!CNCHECK,DETAIL
FINISH
!***************************进入求解器设置
/SOLU
ALLSEL,all ! 恢复选择全部对象
ANTYPE,0 ! 瞬态分析类型(0静态4瞬态)
TIMINT,OFF !No transient (mass or inertia) effects
!TRNOPT,FULL ! 使用完全法瞬态算法
KBC,1
NLGEOM,ON ! 激活大变形
NROPT,FULL,,OFF ! Newton-Raphson option to FULL
ERESX,no ! 积分点结果不外推(适用非线性)
OUTRES,ALL,all ! 存储所有结果
NLDIAG,CONT,SUBS ! 接触分析/子步/.cnd
AUTOTS,ON ! 启用自动时间步
DELTIM,0.00001,0.00001,0.0001! 时间步:初始子步s,最小子步s,最大子步s
TIME,1/(25*z3) ! 第一载荷步时间(齿轮1走1齿)
SOLVE ! 第一次求解
SAVE ! 保存数据库
FINISH ! 退出求解器
!***************************进入后处理设置
/POST1
ESEL,S,TYPE,,7 ! 选择接触单元
ESEL,A,TYPE,,9
! --- 计算第一载荷步的能量 ---
SET,1,last
*GET,N_SUB_1,ACTIVE,,SOLU,NCMSS
*DIM,E_LST_1,ARRAY,N_SUB_1
*DO,i,1,N_SUB_1 ! 循环计算每个时间步的能量
SET,1,i
ETABLE,CON_ARE1,VOLU ! 提取接触面积
ETABLE,E_DISS_I,SMISC,18 ! 提取DISSIPATION方法的能量
ETABLE,E_DISS_J,SMISC,19
ETABLE,E_DISS_K,SMISC,20
ETABLE,E_DISS_L,SMISC,21
*GET,T_STEP,ACTIVE,,SOLU,DTIME
SADD,E_DISS1,E_DISS_I,E_DISS_J,T_STEP,T_STEP
SADD,E_DISS2,E_DISS_K,E_DISS_L,T_STEP,T_STEP
SADD,E_DISS3,E_DISS1,E_DISS2,0.25,0.25
SMULT,E_DISS4,E_DISS3,CON_ARE1
SSUM
*GET,E_DISS_1,SSUM,,Item,E_DISS4
E_LST_1(i)=E_DISS_1
*ENDDO
*STATUS,E_LST_1
*VSCFUN,SUM_E_1,sum,E_LST_1(1)
fricEnergy=SUM_E_1
最新发布