Both solid model and finite element model boundary conditions have been applied to this model.As solid loads are transferred to the nodes or elements,they can overwritte diretly applied loads.
我是搞一个花键轴的优化,当加把轴向力,径向力,弯矩(斜齿)转化为节点力加在花键上后,把转矩加在花键齿上(面力),就弹出了上面的警告!
这是命令流:(请大家帮我分析一下)
/filename,APDLaxis1
/title,model
/prep7
et,1,solid92
mp,ex,1,2.06e11
mp,prxy,1,0.3
/RGB,INDEX,100,100,100,0
/RGB,INDEX,80,80,80,13
/RGB,INDEX,0,0,0,15
/REPLOT
r1=0.0175 !最左端轴径范围(0.01-0.0175M)长度0.017M
r2=0.021 !第二段和花键底径,轴径范围(0.0125-0.021M)长度0.017M
r3=0.024 !花键大径,轴径范围(0.01-0.024M)长度0.017M
r4=0.02 !最右端轴径范围(0.01-0.0175M)长度0.017M
CYLIND,r1,0,0,0.017,0,360,
CYLIND,r2,0,0.017,0.085,0,360,
/VIEW,1,1,1,1
WPSTYLE,,,,,,,,1
wpoff,0,0,0.085
PCIRC,r2,0,0,360,
RECTNG,-0.004,0.004,-0.025,0.025,
/pnum,area,1
/replot
ASBA, 10, 9
numcmp,all
CSYS,1
asel,s,,,9,10,1
AGEN,4,all, , , ,45, , ,0
asel,s,,,9,10,1
asel,a,,,11,12,1
asel,a,,,13,14,1
asel,a,,,15,16,1
asel,a,,,6,,1
aadd,all
numcmp,all
PCIRC,0.027,r3,0,360,
asba,9,10
voffst,11,0.017
numcmp,all
ALLSEL,ALL
APLOT
VSEL,S, , ,3
VGEN,2,all, , , , ,0.110, ,0
VPLOT
wpoff,0,0,0.110
CYLIND,r3,0,0,-0.093,0,360,
CYLIND,r4,0,0.017,0.039,0,360,
allsel,all
/replot
numcmp,all
vsel,s,,,1,2,1
vsel,a,,,3,4,1
vsel,a,,,5,6,1
vadd,all
allsel,all
smrt,5
mshaoe,1,3d
mshkey,0
vmesh,7
finish
/solu
antype,static
aplot
!斜齿
!加载力Pr2=525
allsel,all
ASEL,S,LOC,Z,0.195,0.212
CSYS,1
ASEL,R,LOC,Y,90,270
ASEL,R,LOC,X,0.021
NSLA,R,1
CSYS,0
*SET,number1,ndinqr(0,13)
F,all,FX,525/number1
!加载力矩M=30n.m
allsel,all
CSYS,0
NSEL,S,LOC,Z,0.195,0.212
NSEL,R,LOC,Y,-0.004,0.004
NSEL,R,LOC,X,-0.024,-0.021
*SET,number2,ndinqr(0,13)
F,ALL,FZ,666.6/number2
NSEL,S,LOC,Z,0.195,0.212
NSEL,R,LOC,Y,-0.004,0.004
NSEL,R,LOC,X,0.021,0.024
*SET,number3,ndinqr(0,13)
F,ALL,FZ,-666.6/number3
!加载力P2=1375
allsel,all
ASEL,S,LOC,Z,0.195,0.212
CSYS,1
ASEL,R,LOC,Y,180,360
ASEL,R,LOC,X,0.021
NSLA,R,1
CSYS,0
*SET,number4,ndinqr(0,13)
F,all,Fy,1375/number4
!加载力Pa2=429
allsel,all
NSEL,S,LOC,Z,0.2035
*SET,number5,ndinqr(0,13)
F,all,Fz,429/number5
!加载扭矩M=97,每个面上产生的扭矩为97/8=12.125,每个面的面积为A=0.017*0.003,压力为BF=12.125/0.0225
ALLSEL,ALL
A=0.017*0.003
BF=12.125/0.0225
asel,s,,,67,71,4
asel,a,,,58,62,4
asel,a,,,49,53,4
asel,a,,,76,80,4
SFA,all,1,PRES,BF/A
!直齿
!加载力P1y=706
allsel,all
ASEL,S,LOC,Z,0.085,0.102
CSYS,1
ASEL,R,LOC,Y,0,180
ASEL,R,LOC,X,0.021
NSLA,R,1
CSYS,0
*SET,number6,ndinqr(0,13)
F,all,Fy,-706/number6
!加载力P1z=1776
allsel,all
ASEL,S,LOC,Z,0.085,0.102
CSYS,1
ASEL,R,LOC,Y,-90,90
ASEL,R,LOC,X,0.021
NSLA,R,1
CSYS,0
*SET,number7,ndinqr(0,13)
F,all,FX,-1776/number7
!加载扭矩M=97,每个面上产生的扭矩为97/8=12.125,每个面的面积为A=0.017*0.003,压力为allsel,all
BF=12.125/0.0225
A=0.017*0.003
BF=12.125/0.0225
allsel,all
asel,s,,,31,35,4
asel,a,,,40,44,4
asel,a,,,13,17,4
asel,a,,,22,26,4
SFA,all,1,PRES,BF/A
asel,s,,,3,4,1
DA,all,ALL,0
ASEL,s,LOC,Z,0.217,0.234
NSLA,R,1
nSEL,R,LOC,Z,0.217,0.234
D,all,,,,,,UX,UY
allsel,all
/replot |