请版主帮我的模型指点一下
在设置螺栓预紧力计算的时候出现问题
谢谢
/PREP7
ET,1,SOLID45
MP,EX,1,2e5
MP,PRXY,1,0.3
BLC4,0,0,1200,12
!BLC4,1195,0,10,3528
BLC4,0,3512,1200,16
BLC4,170,3232,8,296
BLC4,170,3232,260,8
BLC4,422,3232,8,296
aadd,3,4,5
agen,2,6,,,600
VOFFST,1,-1000, ,
VOFFST,2,-1000, ,
VDRAG,3, 6, , , , , 46
vovlap,all
a,5,67,16,9,14,15,75,72,31,32,23,24,65,6,3,4
VOFFST,2,5, ,
allsel
vsel,s,,,1,2,1
vsel,a,,,5,15,1
VSYMM,z,all,,, ,0,0
vglue,all
allsel
!vadd,all
allsel
!建立加劲板
wpoff,740,3528 !移动工作平面
wprot,0,0,-90 !旋转工作平面
BLC4,65,-16,312,-20 !建立矩形
BLC4,5,-76,20,-480 !建立矩形
CYL4,65,-76,40,90,60,180 !建立半圆
aadd,7,2,3
VOFFST,8,280, ,
wpoff,0,0,20 !移动工作平面
FLST,2,2,5,ORDE,2
FITEM,2,15
FITEM,2,13
ASBW,P51X
wpoff,0,0,240 !移动工作平面
FLST,2,2,5,ORDE,2
FITEM,2,28
FITEM,2,26
ASBW,P51X
wpoff,0,-511,0 !移动工作平面
wprot,0,-90,0 !旋转工作平面
FLST,2,3,5,ORDE,3
FITEM,2,13
FITEM,2,30
FITEM,2,25
ASBW,P51X
wpoff,0,260,511 !移动工作平面
wprot,0,90,0 !旋转工作平面
k,158+1000,kx(158),ky(158),kz(158)+65
k,118+1000,kx(118),ky(118)-76,kz(118)
a,158,138,117,118,118+1000,158+1000
VOFFST,13,-20, ,
k,154+1000,kx(154),ky(154),kz(154)+65
k,148+1000,kx(148),ky(148)-76,kz(148)
a,154,145,147,148,148+1000,154+1000
VOFFST,59,20, ,
vglue,1,2,3
!VGEN,2,5,7,1,1200,,,,
allsel
VSYMM,z,4,6,1, ,0,0
allsel
aovlap,9,299
aovlap,115,296
allsel
!建立上部吊耳
wprot,0,-90,0 !旋转工作平面
BLC4,380,0,-760,-280,30
wpoff,0,-127.5,0 !移动工作平面
wprot,0,90,0 !旋转工作平面
ASBW, 133
CYL4,0,400,250
wpoff,0,476,0 !移动工作平面
wprot,0,-90,0 !旋转工作平面
ASBW, 133
a,254,253,260,259
aadd,143,133,142
wpoff,0,0,-76 !移动工作平面
wprot,0,90,0 !旋转工作平面
VOFFST,144,-25, ,
CYL4,0,0,230, , , ,-16
wpoff,0,0,25 !移动工作平面
CYL4,0,0,230, , , ,16
wpoff,0,0,-41 !移动工作平面
CYL4,0,0,100, , , ,57
vglue,8,9,10
vadd,9,10,14
vsbv,8,12
allsel
vglue,7,9
allsel
wpoff,230,-370,28.5 !移动工作平面
wprot,0,0,90 !旋转工作平面
asbw,141
asbw,149
k,290+1000,kx(290),ky(290)+190,kz(290)
k,290+1000+1000,kx(290)-60,ky(290)+190,kz(290)
k,289+1000,kx(289),ky(289)+40,kz(289)
a,289,290,290+1000,290+1000+1000,289+1000
k,288+1000,kx(288),ky(288)+190,kz(288)
k,288+1000+1000,kx(288)+60,ky(288)+190,kz(288)
k,287+1000,kx(287),ky(287)+40,kz(287)
a,287,288,288+1000,288+1000+1000,287+1000
VOFFST,149,20, ,
VOFFST,156,-20, ,
VGEN,2,7,10,3,,,-480,,
vglue,7,8,9,10,12,14
wpoff,55,0,70 !移动工作平面
wprot,0,-90,0 !旋转工作平面
CYL4,0,0,13, , , ,-66
CYL4,0,120,13, , , ,-66
CYL4,0,200,13, , , ,-66
CYL4,0,400,13, , , ,-66
CYL4,0,480,13, , , ,-66
CYL4,0,600,13, , , ,-66
CYL4,-110,0,13, , , ,-66
CYL4,-110,120,13, , , ,-66
CYL4,-110,200,13, , , ,-66
CYL4,-110,400,13, , , ,-66
CYL4,-110,480,13, , , ,-66
CYL4,-110,600,13, , , ,-66
vsel,s,,,8,10,1
vsel,a,,,12,14,2
vsel,a,,,20,26,1
cm,kuai5,volu
vsel,s,,,3,6,3
vsel,a,,,18
vsel,a,,,40,44,4
cm,kuai6,volu
allsel
vsbv,kuai6,kuai5
allsel
!建立螺栓
CYL4,0,0,12, , , ,-66!螺杆
CYL4,0,0,25, , , ,25!螺母
VGEN,2,6,,,,-91,,,
vglue,3,6,8
CYL4,-110,0,12, , , ,-66!螺杆
CYL4,-110,0,25, , , ,25!螺母
VGEN,2,8,,,,-91,,,
vglue,6,8,12
allsel
vsel,s,,,3,9,3
vsel,a,,,10,18,4
cm,kuai8,volu
allsel
VGEN,2,kuai8,,,,,-120
VGEN,2,kuai8,,,,,-200
VGEN,2,kuai8,,,,,-400
VGEN,2,kuai8,,,,,-480
VGEN,2,kuai8,,,,,-600
allsel
vsel,s,,,3,6,3
vsel,a,,,8,10,1
vsel,a,,,12,14,2
vsel,a,,,18,20,2
vsel,a,,,21,26,1
vsel,a,,,40,44,4
vsel,a,,,56,74,1
cm,kuai9,volu
allsel
KWPAVE, 1
wprot,0,90,0
WPSTYLE,,,,,,,,0
CSYS,0
allsel
!划分网格
VATT,1,1,1
SMRT,8
MSHAPE,1,3D
MSHKEY,0
VMESH,all
cmsel,s,kuai8
SMRT,5
MSHAPE,1,3D
MSHKEY,0
VMESH,all
allsel
aplot
!* 设置接触对和目标对
ET,2,TARGE170 !建立目标单元
ET,3,CONTA174 !建立接触单元
KEYOPT,3,9,0
mp,mu,2,0.5 !定义摩擦系数
!1
!* 建立第1对接触单元
!* 建立接触对常数
R,1
REAL,1
!* 建立目标面1
ASEL,S,,,381
NSLA,S,1
NSEL,r,LOC,X,460,740
NSEL,r,LOC,Y,3512,3512
NSEL,r,LOC,Z,65,377
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面1
ASEL,S,,,394
NSLA,S,1
CM,_CONTACT,NODE
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!* 建立第2对接触单元
!* 建立接触对常数
R,2
REAL,2
!* 建立目标面2
ASEL,S,,,366
NSLA,S,1
NSEL,r,LOC,X,460,740
NSEL,r,LOC,Y,3512,3512
NSEL,r,LOC,Z,-65,-377
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面2
ASEL,S,,,408
NSLA,S,1
CM,_CONTACT,NODE
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!2
!* 建立第3对接触单元
!* 建立接触对常数
R,3
REAL,3
!* 建立目标面3
ASEL,S,,,367,380,13
NSLA,S,1
NSEL,r,LOC,X,460,740
NSEL,r,LOC,Y,3528,3528
NSEL,r,LOC,Z,-380,380
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面3
ASEL,S,,,440
NSLA,S,1
CM,_CONTACT,NODE
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!3
!上部螺栓接触
!* 建立第4对接触单元
!* 建立接触对常数
R,4
REAL,11
!* 建立目标面4
ASEL,S,,,434
NSLA,S,1
NSEL,r,LOC,X,630,680
NSEL,r,LOC,Y,3558,3558
NSEL,r,LOC,Z,275,325
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面4
ASEL,S,,,165
NSLA,S,1
CM,_CONTACT,NODE
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!* 建立第5对接触单元
!* 建立接触对常数
R,5
REAL,5
!* 建立目标面5
ASEL,S,,,438
NSLA,S,1
NSEL,r,LOC,X,630,680
NSEL,r,LOC,Y,3558,3558
NSEL,r,LOC,Z,-180-25,180+25
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面5
ASEL,S,,,327,351,24
ASEL,a,,,462,486,24
NSLA,S,1
CM,_CONTACT,NODE
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!* 建立第6对接触单元
!* 建立接触对常数
R,6
REAL,6
!* 建立目标面6
ASEL,S,,,435
NSLA,S,1
NSEL,r,LOC,X,630,680
NSEL,r,LOC,Y,3558,3558
NSEL,r,LOC,Z,-275,-325
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面6
ASEL,S,,,510
NSLA,S,1
CM,_CONTACT,NODE
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!* 建立第7对接触单元
!* 建立接触对常数
R,7
REAL,7
!* 建立目标面7
ASEL,S,,,437
NSLA,S,1
NSEL,r,LOC,X,520,570
NSEL,r,LOC,Y,3558,3558
NSEL,r,LOC,Z,275,325
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面7
ASEL,S,,,311
NSLA,S,1
CM,_CONTACT,NODE
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!* 建立第8对接触单元
!* 建立接触对常数
R,8
REAL,8
!* 建立目标面8
ASEL,S,,,439
NSLA,S,1
NSEL,r,LOC,X,520,570
NSEL,r,LOC,Y,3558,3558
NSEL,r,LOC,Z,-180-25,180+25
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面8
ASEL,S,,,335,446,111
ASEL,a,,,470,494,24
NSLA,S,1
CM,_CONTACT,node
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!* 建立第9对接触单元
!* 建立接触对常数
R,9
REAL,9
!* 建立目标面9
ASEL,S,,,436
NSLA,S,1
NSEL,r,LOC,X,520,570
NSEL,r,LOC,Y,3558,3558
NSEL,r,LOC,Z,-275,-325
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面9
ASEL,S,,,518
NSLA,S,1
CM,_CONTACT,node
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!4螺栓下部
!* 建立第10对接触单元
!* 建立接触对常数
R,10
REAL,10
!* 建立目标面10
ASEL,S,,,395
NSLA,S,1
NSEL,r,LOC,X,520,570
NSEL,r,LOC,Y,3492,3492
NSEL,r,LOC,Z,75,325
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面10
ASEL,S,,,310,332,22
ASEL,a,,,443
NSLA,S,1
CM,_CONTACT,node
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!* 建立第11对接触单元
!* 建立接触对常数
R,11
REAL,11
!* 建立目标面11
ASEL,S,,,395
NSLA,S,1
NSEL,r,LOC,X,630,680
NSEL,r,LOC,Y,3492,3492
NSEL,r,LOC,Z,75,325
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面11
ASEL,S,,,145
ASEL,a,,,324,348,24
NSLA,S,1
CM,_CONTACT,node
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!* 建立第12对接触单元
!* 建立接触对常数
R,12
REAL,12
!* 建立目标面12
ASEL,S,,,409
NSLA,S,1
NSEL,r,LOC,X,520,570
NSEL,r,LOC,Y,3492,3492
NSEL,r,LOC,Z,-75,-325
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面12
ASEL,S,,,515
ASEL,a,,,467,491,24
NSLA,S,1
CM,_CONTACT,node
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
!* 建立第13对接触单元
!* 建立接触对常数
R,13
REAL,13
!* 建立目标面13
ASEL,S,,,409
NSLA,S,1
NSEL,r,LOC,X,630,680
NSEL,r,LOC,Y,3492,3492
NSEL,r,LOC,Z,-75,-325
CM,_TARGET,NODE
type,2
ESLN,S,0
esurf,all,
ALLSEL,ALL
!*
!* 建立接触面13
ASEL,S,,,507
ASEL,a,,,459,483,24
NSLA,S,1
CM,_CONTACT,node
type,3
ESLN,S,0
esurf,all,
ALLSEL,ALL
allsel
aplot
!!!!!!!!!!!!高强螺栓psmesh
!建立螺栓预紧力单元
!shuxiang
psmesh,1,fbolt1,,volu,3,0,y,3525,,,,fbolt1
psmesh,2,fbolt2,,volu,8,0,y,3525,,,,fbolt2
psmesh,3,fbolt3,,volu,24,0,y,3525,,,,fbolt3
psmesh,4,fbolt4,,volu,59,0,y,3525,,,,fbolt4
psmesh,5,fbolt5,,volu,65,0,y,3525,,,,fbolt5
psmesh,6,fbolt6,,volu,69,0,y,3525,,,,fbolt6
psmesh,7,fbolt7,,volu,6,0,y,3525,,,,fbolt7
psmesh,8,fbolt8,,volu,12,0,y,3525,,,,fbolt8
psmesh,9,fbolt9,,volu,25,0,y,3525,,,,fbolt9
psmesh,10,fbolt10,,volu,58,0,y,3525,,,,fbolt10
psmesh,11,fbolt11,,volu,64,0,y,3525,,,,fbolt11
psmesh,12,fbolt12,,volu,70,0,y,3525,,,,fbolt12
allsel
eplot
/solu
lsclear,all
/SOLU
ANTYPE,0 !分析类型Specifies the analysis type and restart status,静力
NLGEOM,1 !打开大变形分析Includes large-deflection effects in a static or full transient analysis
NSUBST,5,0,0 !Specifies the number of substeps to be taken this load step
OUTRES,all,1 !Controls the solution data written to the database
NEQIT,5
LNSRCH,0
RESCONTRL,DEFINE,ALL,ALL,3
TIME,10
allsel
eplot
/solu
ASEL,S,LOC,X,0,0
DA,all,symm
ASEL,s,LOC,X,1200,1200
DA,all,symm
ASEL,S,LOC,Z,-1000,-1000
DA,all,symm
ASEL,s,LOC,Z,1000,1000
DA,all,symm
ASEL,s,LOC,y,0,0
DA,all,all
allsel
!!!在螺栓中施加螺栓预紧力
sload,1,9,SLID,FORC,240000,1,2
sload,2,9,SLID,FORC,240000,1,2
sload,3,9,SLID,FORC,240000,1,2
sload,4,9,SLID,FORC,240000,1,2
sload,5,9,SLID,FORC,240000,1,2
sload,6,9,SLID,FORC,240000,1,2
sload,7,9,SLID,FORC,240000,1,2
sload,8,9,SLID,FORC,240000,1,2
sload,9,9,SLID,FORC,240000,1,2
sload,10,9,SLID,FORC,240000,1,2
sload,11,9,SLID,FORC,240000,1,2
sload,12,9,SLID,FORC,240000,1,2
allsel
SFA,162,1,PRES,220
allsel
solve |