第五例杆系结构的静力学分析实例——平面桁架
/clear
L=0.1
A=1e-4
/PREP7
ET,1,BEAM3
R,1,A
MP,EX,1,2E11
MP,PRXY,1,0.3
N,1
N,2,L
N,3,2*L
N,4,L,L
E,1,2
E,2,3
E,1,4
E,2,4
E,3,4
FINISH
/SOLU
D,1,UX
D,1,UY
D,3,UY
F,4,FY,-2000
SOLVE
FINISH
/POST1
ETABLE,FA,SMISC,1
ETABLE,SA,LS,1
PRETAB,FA,SA
FINISH
第六例干系结构的静力学分析——悬臂梁
一工字悬臂梁,分析其在集中力p作用下自由端的变形,一直梁的材料为10号热轧工字钢,其横截面面积为A=14.345cm2,截面高度H=100mm,惯性矩I=245cm4,梁的长度L=1m,集中力p=10000N,钢的弹性模量E=2e11N/m2,泊松比为0.3
/clear
/prep7
Et,1,beam3
R,1,14.345e-4,245e-8,0.1
Mp,ex,1,2e11
Mp,prxy,1,0.3
K,1,0,0,0
Lstr,1,2
Lesize,1,,,50
Lmesh,1
Finish
/solu
Dk,1,ux
Dk,1,uy
Dk,1,rotz
Fk,2,fy,-10000
Solve
Finish
/post1
Pldisp
Finish
第七例平面问题的求解实例——厚壁圆筒问题
一厚壁圆筒,其内半径r1=50mm,外半径r2=100mm,作用在内孔上的压力p=10MPa,无轴向压力,轴向长度很大可视为无穷,要求计算厚壁圆筒的径向应力和切向应力沿半径方向的分布。
/clear
/filname,example7
/prep7
Et,1,plane183,,,2
Mp,ex,1,2e11
Mp,prxy,1,0.3
Pcirc,0.1,0.05,0,90
Lesize,4,,,6
Lesize,3,,,8
Mshape,0
Mshkey,1
Amesh,1
Finish
/solu
Dl,4,,uy
Dl,2,,ux
Sfl,3,pres,10e6
Solve
/post1
Path,p1,2
Ppath,1,30
Ppath,2,1
pdef,sr,s,x
Pdef,st,s,y
Plpath,sr,st
Finish
第八例扳手的受力分析/clear
/filname,example8
/prep7
Et,1,plane42
Et,2,solid45
Mp,ex,1,2e11
Mp,prxy,1,0.3
Rpr4,6,0,0,0.01
K,7,0,0,0
K,8,0,0,0.05
K,9,0,0.1,0.05
Lstr,7,8
Lstr,8,9
Lfillt,7,8,0.015
Lstr,1,4
AsbL,1,10
Lesize,2,,,3
Lesize,3,,,3
Lesize,4,,,3
Lesize,7,0.01
Lesize,8,0.01
Lesize,9,0.01 Mshape,0
Mshkey,1
Amesh,all
Vdrag,all,,,,,,7,9,8 Aclear,2,3,1
Finish
/solu
Da,2,all
Da,3,all
Ksel,s,,,24,29,1
Fk,all,fx,100
Ksel,all
Solve
Finish
/post1
/view,1,1,1,1
Pldisp,2
Plnsol,s,eqv,0,1
Nwpave,159
/type,1,sect
/cplane,1
/replot
finish
第九章圆轴扭转分析
设等直圆轴的圆截面直径D=50mm,长度L=120mm,作用在圆轴两端上的转矩Mn=1.5xe3N.m
/clear
/filname,example9
/prep7
Et,1,plane183
Et,2,solid186
Mp,ex,1,2.08e11
Mp,prxy,1,0.3
Rectng,0,0.025,0,0.12
Lesize,1,,,5
Lesize,2,,,8
Mshape,0
Mshkey,1
Amesh,1
Extopt,esize,5
Extopt,aclear,1
Vrotat,1,,,,,,1,4,360
/view,1,1,1,1
Wprot,0,-90
Cswpla,11,1,1,1
Nsel,s,loc,x,0.025
Rpotat,all
Finish
/solu
D,all,ux
Nsel,r,loc,z,0.12
F,all,fy,1500
Allsel,all
Da,6,all
Da,10,all
Da,14,all
Solve
Finish
/post1
Pldisp,1
Rsys,11
Nsel,s,loc,z,0.045
Nsel,r,loc,y,0
Prnsol,u,y
nsel,s,loc,z,0,0.045
esln,r,1
plesol,s,yz
Finish
第十例模态分析实例——均匀直杆的固有频率分析
一根长度为L的等截面直杆,一端固定,一段自由。已知杆材料的弹性模量为E=2e11,密度为7800kg/m3,杆长L=0.1m,要求计算直杆纵向振动的固有频率。
/CLEAR
/FILNAME,EDAMPLE10
/PREP7
ET,1,SOLID186
MP,EX,1,2E11
MP,PRXY,1,0.3
MP,DENS,1,7800
BLOCK,0,0.01,0,0.01,0,0.01
LESIZE,1,,,3
LESIZE,2,,,3
LESIZE,9,,,15
MSHAPE,0
MSHKEY,1
VMESH,1
FINISH
/SOLU
ANTYPE,MODAL
MODOPT,LANB,5
MXPAND,5
DA,1,UZ
DA,3,UY
DA,5,UX
SOLVE
FINISH
SET,LIST
SET,FIRST
/REPLOT
ANMODE,10,0.5,,0
SET,NEXT
PLDI
ANMODE,10,0.5
FINISH
/CLEAR
/FILNAME,EXAMPLE11
/PREP7
ET,1,SOLID45
MP,EX,1,2E11
MP,PRXY,1,0.3
MP,DENS,1,7800
K,1,21.87E-3
K,2,22.82E-3,1.13E-3
K,3,24.02E-3,1.47E-3
K,4,24.62E-3,1.73E-3
K,5,25.22E-3,2.08E-3
K,6,25.82E-3,2.4E-3
K,7,26.92E-3,3.23E-3
K,8,27.11E-3
BSPLIN,2,3,4,5,6,7
LSYMM,Y,1
LARC,2,9,1
LARC,7,10,8
AL,ALL
CSYS,1
VEXT,1,,,0,8.412,20E-3
VGEN,24,1,,,0,360/24
CSYS,0
CYL4,0,0,10E-3,0,26.37E-3,360,20E-3 BLOCK,-3E-3,3E-3,0,12.8E-3,0,20E-3 VSBV,25,ALL
SMRTSIZE,9
ESIZE,0.002
MSHAPE,1
MSHKEY,0
VMESH,ALL
CSYS,1
NSEL,S,LOC,X,0.01
NROTA T,ALL
D,ALL,UX
ALLSEL,ALL
FINISH
/SOLU
ANTYPE,MODAL
MODOPT,LANB,5
DA,208,UX
NSEL,S,LOC,Z,0
NSEL,A,LOC,Z,20E-3
NSEL,R,LOC,X,0,15E-3
D,ALL,UZ
ALLSEL,ALL
SOLVE
FINISH
/POST1
SET,LIST
FINISH
第十一例模态分析实例——斜齿圆柱齿轮的固有频率分析
/clear
/filname,example11
/prep7
Et,1,SOLID45
MP,EX,1,2E11
MP,PRXY,1,0.3
MP,DENS,1,7800
K,1,0.02187
K,2,0.02282,0.00113
K,3,0.02402,0.00147
K,4,0.02462,0.00173
K,5,0.02522,0.00208
K,6,0.02582,0.0024
K,7,0.02692,0.00323
K,8,0.02711
BSPLIN,2,3,4,5,6,7
LSYMM,Y,1
LARC,2,9,1
LARC,7,10,8
AL,ALL
CSYS,1
VEXT,1,,,0,8.412,0.02
VGEN,24,1,,,0,360/24
CYL4,0,0,0.01,0,0.02637,360,0.02
BLOCK,-3E-3,3E-3,0,12.8E-3,0,20E-3
VSBV,25,ALL
SMRTSIZE,9
ESIZE,0.002
MSHAPE,1
MSHKEY,0
VMESH,ALL
CSYS,1
NSEL,S,LOC,X,0.01
NROTA T,ALL
D,ALL,UX
ALLSEL,ALL
FINISH
/SOLU
ANTYPE,MODAL
MODOPT,LANB,5
MXPAND,5
DA,208,UX
NSEL,S,LOC,Z,0
NSEL,A,LOC,Z,20E-3
NSEL,R,LOC,X,0,15E-3
D,ALL,UZ
ALLSEL,ALL
SOLVE
FINISH
/POST1
SET,LIST
FINISH
第十二例有预应力的模态分析实例——弦的横向振动/CLEAR
/FILNAME,EXAMPLE12
/PREP7
ET,1,LINK1
MP,EX,1,2E11
MP,PRXY,1,0.3
MP,DENS,1,7800
R,1,1E-6
K,1,0,0,0
K,2,1
LSTR,1,2
LESIZE,1,,,50
LMESH,1
/SOLU
DK,1,UX
DK,1,UY
DK,2,UY
FK,2,FX,2000
PSTRES,ON
SOLVE
FINISH
/SOLU
ANTYPE,MODAL
MODOPT,LANB,10
MXPAND,10
DK,2,UX
PSTRES,ON
SOLVE
FINISH
/POST1
SET,LIST
SET,FIRST
PLDI
ANMODE,10,0.5
SET,NEXT
PLDI
ANMODE,10,0.5
FINISH
第十三例谐响应分析实例——单自由度系统的受迫振动/clear
/FILNAME,EXAMPLE13
/PREP7
ET,1,MASS21
ET,2,COMBIN14
R,1,1
R,2,10000,63
N,1
N,2,1
TYPE,2
REAL,2
E,1,2
TYPE,1
REAL,1
E,2
FINISH
/SOLU
ANTYPE,HARMIC
D,1,ALL
D,2,UY
D,2,UZ
D,2,ROTX
D,2,ROTY
D,2,ROTZ
HARFRQ,0,50
NSUBST,25
KBC,1
F,2,FX,2000
SOLVE
FINISH
/POST26
NSOL,2,2,U,X,DISPX
PLV AR,2
PLCPLX,1
PLV AR,2
FINISH
第十四例瞬态动力学分析实例——凸轮从动件运动分析/CLEAR
/FILNAME,EXAMPLE14
/PREP7
ET,1,PLANE42
ET,2,SOLID45
MP,EX,1,2E11
MP,PRXY,1,0.3
MP,DENS,1,7800
K,1
K,2,0.015,0.015
K,3,0.015,0.1
K,4,0,0.1
LSTR,1,2
LSTR,2,3
LSTR,3,4
LSTR,4,1
AL,ALL
LESIZE,1,,,2
LESIZE,3,,,2
LESIZE,2,,,10
AA TT,1,1,1
MSHAPE,0
MSHKEY,1
AMESH,ALL
TYPE,2 EXTOPT,ESIZE,4 EXTOPT,ACLEAR,1 VROTA T,1,,,,,,1,4,360 WPROT,0,-90 CSWPLA,11,CYLIN ASEL,S,,,3,15,4 NSLA,S,1
NROTA T,ALL FINISH
/SOLU
D,ALL,UX ALLSEL,ALL ANTYPE,TRANS OUTRES,ALL,ALL FK,4,FY,-1000 TIME,10 AUTOTS,ON DELTIM,0.5,0.2,1 KBC,0
DK,1,UY,0.02 LSWRITE,1 TIME,20 LSWRITE,2 TIME,35
DK,1,UY,0 LSWRITE,3 TIME,45 LSWRITE,4 LSSOLVE,1,4,1 FINISH
/POST26
NSOL,2,1,U,Y,UY DERIV,3,2,1 DERIV,4,3,1
PLV AR,2,3
PLV AR,2,4
FINISH