侵权不好
/TITLE,Buckling
/config,nres,2500
/PREP7
ET,1,PLANE42,,,3,,2
ET,2,COMBIN14,,,2
R,1,4
R,2,1E6
MP,EX,1,25E6
MP,PRXY,1,0.17
MP,DENS,1,2500
K,1,-2
K,2,-2,100
K,3,2
K,4,2,100
A,1,2,4,3
LSEL,S,LOC,X,-1,1
LESIZE,ALL,,,4
ALLSEL
LSEL,S,LOC,Y,1,99
LESIZE,ALL,,,50
ALLSEL
AMESH,ALL
N,1000,-5,100
TYPE,2
REAL,2
E,1000,2
save,bucklingmodel
/SOLU
NSEL,S,LOC,Y,100
NSEL,R,LOC,X,-1.5,1.5
F,ALL,FY,-25
ALLSEL
NSEL,S,LOC,Y,100
NSEL,R,LOC,X,-2
F,ALL,FY,-12.5
ALLSEL
NSEL,S,LOC,Y,100
NSEL,R,LOC,X,2
F,ALL,FY,-12.5
ALLSEL
NSEL,S,LOC,Y
D,ALL,ALL
ALLSEL
D,1000,ALL
PSTRES,ON
SOLVE
finish
save,bucklingrst1
/solu
ANTYPE,BUCKLE
BUCOPT,SUBSP,1
MXPAND,1
SOLVE
/POST1
SET,,,1,,1
finish
save,bucklingrst2
/PREP7
TB,MISO,1
TBPT,,0.0002,5000
TBPT,,0.0004,9250
TBPT,,0.0006,13000
TBPT,,0.0008,16250
TBPT,,0.001,19000
TBPT,,0.0012,21250
TBPT,,0.0014,23000
TBPT,,0.0016,24250
TBPT,,0.0018,25000
TBPT,,0.0020,25250
TBPT,,0.0022,24312.5
TBPT,,0.0024,23375
TBPT,,0.0026,22437.5
TBPT,,0.0028,21500
TBPT,,0.003,20562.5
/SOLU
fac=10664
NSEL,S,LOC,Y,100
NSEL,R,LOC,X,-1.5,1.5
F,ALL,FY,-25*fac
ALLSEL
NSEL,S,LOC,Y,100
NSEL,R,LOC,X,-2
F,ALL,FY,-12.5*fac
ALLSEL
NSEL,S,LOC,Y,100
NSEL,R,LOC,X,2
F,ALL,FY,-12.5*fac
ALLSEL
/solu
UPCOORD,0.1
SAVE,bucklingrst3
ANTYPE,STATIC
NLGEOM,ON
SSTIF,ON
NEQIT,10
NSUBST,80
OUTRES,ALL,ALL
ARCLEN,ON,2
SOLVE
FINISH
SAVE
/POST1
SET,,,1,,,,875
PLDISP,1
PLNSOL,S,1,0,1
PLNSOL,U,X,0,1
PLNSOL,NL,SEPL,0,1
PLNSOL,EPEL,EQV,0,1
SAVE,bucklingrst4
FINISH |