把出现以上问题的命令流贴上来,希望有人看下. 想做轴径的优化. 不知道这样做对不对.
!file:mode.txt
!初始化设计变量参数
*set,r1,0.375
*set,r2,0.3
*set,r3,0.05
*set,r4,0.16
*set,r5,0.18
*set,d1,0.125
*set,d2,0.35
*set,d3,0.415
*set,d4,0.35
*set,d5,0.08
*set,d6,0.07
*set,L1,0.3
*set,L2,0.366
*set,L3,0.365
*set,L4,0.186
*set,L5,0.18
*set,L6,0.12
*set,L7,0.047
*set,L8,0.1
*set,L9,0.38
*set,L10,0.3344
*set,L11,1.13 !-------增加
*set,L12,0.085 !锤头尺寸
*set,L13,0.08 !锤头尺寸
*set,L14,0.125 !锤头
*set,L15,0.0675 !锤头
*set,h_r1,0.05
*set,h_d1,0.05
*set,h_d2,0.045
*set,shaft_r1,0.157
*set,shaft_r2,0.16
*set,shaft_r3,0.135
*set,shaft_r4,0.13
*set,shaft_r5,0.125
*set,shaft_r6,0.119
*set,shaft_r7,0.116
*set,shaft_l1,0.144
*set,shaft_l2,0.42
*set,shaft_l3,0.073
*set,shaft_l4,0.216
*set,shaft_l5,0.045
*set,shaft_l6,0.12
/prep7
et,1,solid92
mp,ex,1,2.06e11
mp,prxy,1,0.3
mp,dens,1,7800
!建立锤体模型
k,1
k,2,,-d1
k,3,l6,d5-d3,
k,4,l7,d5-d3,
k,5,l7,-r1-d1
k,6,l6,-d1
circle,1,r2
circle,2,r1
l,3,4
l,4,5
l,3,6
lovlap,10,8
lovlap,11,4
lovlap,2,6
ldele,13,14,1,1
ldele,10,11,1,1
ldele,17,18,1,1
ldele,3,5,2,1
kdele,2
al,all !由线生成面
voffst,1,L2/2 !由面偏移生成体
agen,2,1,,,,,L3/2-L5/2+L2/2
wpoffs,,,L3/2-L5/2+L2/2
csys,1
agen,2,12,,,,180
csys,0
agen,2,12,13,,,,L4
adele,12,,,1
aadd,14,15
voffst,13,L5/2-L3/2
voffst,13,L4/2
voffst,12,-L4/2
voffst,12,L3/2-L4/2
vadd,all
kwpave,1
cyl4,0,0,,,r5,,L1/2
vsbv,6,1
cyl4,0,0,,,shaft_r2,,1
vsbv,2,1
!锥体部分,在锤体上开键槽
k,1000,-d6/2,L10-shaft_r2
k,1001,d6/2,L10-shaft_r2
k,1002,d6/2,0.13
k,1003,-d6/2,0.13
k,1004,L10-shaft_r2,d6/2
k,1005,L10-shaft_r2,-d6/2
k,1006,0.13,-d6/2
k,1007,0.13,d6/2
a,1000,1001,1002,1003
a,1004,1005,1006,1007
voffst,1,-L11/2
voffst,2,-L11/2
vsbv,3,1
vsbv,4,2
!在锤体上开与锤头连接用的槽
k,200,-L8,-d4
k,201,-L8,-0.2-d4
k,202,-L8-0.3,-d4
l,200,201
l,200,202
lovlap,7,103,104
ldele,164,171,7,1 !删除分割线后得到的多余线段
ldele,105,106,1,1
al,172,163,160
agen,2,2,,,,,L2/4
agen,2,2,,,,,L3/2-L5/2+L2/2+L4
csys,1
agen,2,2,,,,180
csys,0
agen,2,36,,,,,L2/2+L3/2-L5/2
agen,2,36,,,,,L2/2+L3/2-L5/2+L4
adele,2,,,1
adele,36,,,1
voffst,24,-0.06 !由面偏移生成体
voffst,24,0.06
voffst,34,-0.06
voffst,34,0.06
voffst,37,-0.06
voffst,37,0.06
voffst,38,-0.06
voffst,38,0.06
vsbv,1,2
vsbv,10,3
vsbv,1,4
vsbv,2,5
vsbv,1,6
vsbv,2,7
vsbv,1,8
vsbv,2,9
!开螺栓孔
k,300,-L8,-d3,L2/4
kwpave,300
wprota,,,90
cyl4,0,0,,,0.105/2,,0.052
wpoff,,,0.052
cyl4,0,0,,,0.038/2,,0.043 !0.038为锤头小端内孔直径
wpoff,,,0.043
cyl4,0,0,,,0.105/2,,0.052
vadd,2,3,4
kwpave,1
WPROTA,,,-90
vgen,2,5,,,,,L2/4+L3/2-L5/2+L4 !复制生成全部体积
csys,1
vgen,2,5,,,,180
csys,0
vgen,2,3,,,,,L2/4+L3/2-L5/2
vgen,2,3,,,,,L2/4+L3/2-L5/2+L4
vdele,3,,,1
vsbv,1,5
vsbv,3,2
vsbv,1,4
vsbv,2,6
!建立锤头的模型
kwpave,300
WPROTA,,,90 !工作平面绕y轴顺时针方向旋转90度
wpoff,,,L7+L8
cyl4,0,0,,,h_r1,,-h_d1
cyl4,0,0,,,0.038/2,,-h_d1 !0.038为锤头小端内孔直径
vsbv,2,3
rectng,L5/2,-L5/2,-L12,L13
voffst,55,L14
vadd,2,4
cyl4,0,0,,,0.038/2,,h_d2
vsbv,3,2
lfillt,389,390,0.025
lfillt,391,387,0.025
l,265,266
adrag,397,,,,,,406
vsba,4,55
vdele,2,,,1
kwpave,210
csys,4
k,600,,,L15
l,600,267
adrag,386,,,,,,410
vsba,3,55
vdele,2,,,1
kwpave,300 !锤头上的通孔
wpoff,,,L7+L8+h_d2
csys,4
cyl4,0,0,,,0.06/2,,0.08
vsbv,4,2
!建立螺栓组建
cyl4,0,0,,,0.054/2,,0.036
cyl4,0,0,,,0.036/2,,-0.17
wpoff,,,-0.14
rpr4,6,0,0,0.055/2,,-0.03
vadd,2,4,5
wpoff,,,0.156
rpr4,6,0,0,0.03085/2,,0.02
vsbv,6,2
vgen,2,3,4,1,-(L2/4+L3/2-L5/2+L4) !复制生成所有体
kwpave,1
WPROTA,,,-90
csys,1
vgen,2,3,4,1,,180
csys,0
vgen,2,6,7,1,,,L2/4+L3/2-L5/2
vgen,2,6,7,1,,,L2/4+L3/2-L5/2+L4
vdele,6,7,1,1
!建立键模型
lgen,2,166,,,,,0.06
lsbl,161,577,,,keep
lsbl,162,577,,,keep
al,579,577,580,165
kl,577,0.5
kwpave,398
WPROTA,,,-90
csys,4
cyl4,0,0,,,d6/2
aadd,207,208
voffst,209,0.0364 !键的厚度0.0346为一关键参数
kwpave,1
WPROTA,,,90
csys,1
vgen,2,6,,,,90
!建立主轴模型,需注意参数化建模
cyl4,0,0,,,shaft_r1,,shaft_l1
wpoff,,,shaft_l1
cyl4,0,0,,,shaft_r2,,shaft_l2
wpoff,,,shaft_l2
cyl4,0,0,,,shaft_r3,,shaft_l3
wpoff,,,shaft_l3
cyl4,0,0,,,shaft_r4,,shaft_l4
wpoff,,,shaft_l4
cyl4,0,0,,,shaft_r5,,shaft_l5
wpoff,,,shaft_l5
cone,shaft_r6,shaft_r7,0,shaft_l6
wpoff,,,shaft_l6
vsbv,13,7,,,keep !生成轴上键槽1
vsbv,18,6,,,keep !生成轴上键槽2
vadd,12,13,14,15,16,17 !体相加生成主轴
save
!完成建模,glue体
vglue,1,11,10
vglue,14,8,9
vglue,11,2,5
vglue,14,3,4
vglue,11,6
vglue,4,7
vglue,3,18
vglue,7,6
!选择锤头
vsel,s,volu,,1,2,1
vsel,a,volu,,8,12,4
esize,0.04
vmesh,all
vsel,all
!螺栓
vsel,s,volu,,10,13,3
vsel,a,volu,,5,9,4
esize,0.04
vmesh,all
vsel,all
!选择键及轴进行网格划分
vsel,s,volu,,3,4,1
vsel,a,volu,,14
esize,0.04
vmesh,all
vsel,all
!选择锤体进行网格划分
vsel,s,volu,,11
esize,0.04
vmesh,all
allsel,all
nummrg,all,,,,low
numcmp,all
finish
!求解
/solu
da,13,symm !施加对称约束
da,221,symm !施加对称约束
da,242,ux,0 !轴端面施加约束
da,242,uy,0
da,242,uz,0
da,242,rotx,0
da,242,roty,0
/input,loads,lod !读入loads.lod载荷步文件,循环求解
*do,i,1,3 !3代表loads文件中有3个载荷步,需根据实际情况修改
lsread,i
irlf,1
solve
*enddo
save
finish
!进入后处理提取参数
/post1
vsel,s,,,10 !选择轴
eslv,s !选择与轴相连的单元
etable,EVOLUME,volu, !evolume=每个单元的体积,将轴上每个单元的体积放入表格
ssum !将单元表中每列的数据相加
*get,VTOT,ssum,,item,evolume !参数VTOT=总体积
*set,wt,7800*vtot
nsle,s,,,
nsort,s,eqv,1,1
*get,MAXSTRESS,sort,0,max !maxsrress=最大应力
nsort,U,Y,,1 !以UY为基准对节点排序
*get,DEFL,sort,,max !DEFL=节点的竖向最大位移
DEFL=ABS(DEFL) !定义竖向变形的绝对值
esel,s,,,all
nsel,s,,,all
allsel,all |