本帖最后由 小小实习生 于 2011-9-30 15:53 编辑
* S+ X' {3 b4 U+ x! E% w9 i
' p- R1 I y$ Y9 g这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。
4 R" I# I1 i, k m$ s7 FFINI
' m9 V+ _- B- O0 b! l7 a/CLE' J7 M0 X5 a, h7 r7 {( z
!设置参数
/ f/ x% K3 P+ Q7 f# [% P& Z: X. I. F/CONFIG,NRES,2000
5 P% D7 C6 [- [! r) o, G: l*SET,LONG,40 !梁长度9 ^. ~( N( F9 J! r
*SET,NE,20 !单元数
8 J7 U T& [# f( D*SET,DL,LONG/NE !单元长度) ^4 }, \2 X, w; h) O3 w
*SET,NN,NE+1 !节点数
2 e0 D p+ v% U8 ~- o*SET,FF,5000 !移动的集中力
+ k2 c# E. Z6 K+ P" M*SET,VV,80 !移动速度$ c: P. l( y7 f
*SET,DT,DL/VV*3.6 !移动经过一个单元所需时间
6 a) g& t2 V3 g- k w |!定义单元材料实常数2 X0 U1 h: G4 t! K
/PREP7( _# a4 h- `" Z
ET,1,BEAM3 !梁单元
+ N! v5 }9 l1 {$ V# K# d G3 ^MP,EX,1,50E9 !弹性模量6 q6 m" h6 f6 j7 ^' [9 @
MP,DENS,1,2500 !密度/ n7 V" c# t9 `, U0 M
MP,NUXY,1,0.2
5 \- _" a' d( t7 pR,1,0.5,0.041667,1 !几何常数
4 u, L7 d" Z3 b!建立模型
5 \1 \+ P+ p* v0 O* p. k7 F*DO,I,1,NN,1
* B! c6 ]4 H6 q5 z# _6 fN,I,(I-1)*DL !建立
# f5 n' }$ N/ X*ENDDO
! _3 i6 v3 n$ v*DO,I,1,NE- s3 c* ~* l& J+ E- n
E,I,I+1
% H! `0 ~ V% ?' c* E0 `*ENDDO
. `0 i7 b( A4 d1 r' M- e7 hFINI
, h# G6 D; T3 U1 U!加载与求解
: \8 @/ t: N# N+ H7 t/SOL' q# U' o$ y$ D. C: T) H1 L" l
D,1,ALL/ m1 U* d& m" Q
D,NN,ALL4 f4 u! A0 X1 U% c; w! G
ANTYPE,4 !瞬态分析
' y k4 V' _8 z" b. P7 m. VTRNOPT,FULL
/ k8 @4 c# M, L( xTIMINT,OFF
7 F2 D2 N+ \* c3 G' \3 \( HOUTRES,BASIC,ALL !输出基本项，每一步都输出
1 {# Y7 w( i$ U4 l( u+ w: e' ?KBC,1 !指定载荷类型
# g1 G# A+ R+ ^TIME,1E-10 !指定载荷步结束时间
' U% s/ z4 r2 Z% a1 RNSUB,4,8 !指定载荷子步数为4，最大子步数为8! p6 X4 i* C5 V2 R* X
SSTIF,ON !打开应力刚化效应
$ U. l2 R. b5 e D2 MACEL,0,9.8,0 !施加重力加速度3 _, w- F1 |5 ], J# r( r
ALLS !选中所有元素
# v% b3 f3 L! W; ^' p# T4 D6 Y- BSOLVE !求解第一个载荷步
2 l, e2 v3 X% P% e; ~- u/ B8 |TIMINT,ON" z; w4 {) h- e. w" ]: o% {
*DO,II,1,NN,1
! L4 t Z7 N0 P( V+ k* s% P; C( b* I* g6 TOUTRES,ALL,ALL !输出所有加载子步的所有结果3 ^7 Z/ W- h$ Y8 b; N
TIME,II*DT !每个荷载步的求解时间+ Z; `3 \1 S1 _6 x
NSUBST,5 !每个荷载步包含五个子步
1 C) f& h" ^ W2 ]- _& [* \$ zFDELE,ALL,ALL !删除所有的集中力荷载& y; Y1 n& y8 H: f6 o$ r$ r
NSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点
9 d7 X" `0 K! a- x3 H; q. }1 E ]F,ALL,FY,-FF !施加集中荷载6 t4 @9 r/ U# _/ j1 q9 N
ALLSEL !选择所有
, D5 V- n) V8 @+ ]% fSOLVE/ d# M& c! O2 ?
*ENDDO
4 k" a. _! f/ x- v' Q; z; c7 _SAVE
+ J/ L+ X/ x9 Y3 a) r2 w* x. o: XFINI
9 ?& d3 u- X5 |1 R7 E f4 y/POST1. V; X+ q2 E8 C5 ~! C' K9 v
SET,NE/2+1,LAST !设定读取的载荷步9 V* c& I0 T( j/ \. ]: O# p6 c" J
/ESHAPE,1 !显示梁形状# p7 a8 ~, `% w. _: p/ y
PLNSOL,U,Y,0,1.0 !Y向位移
" v2 q5 K% X' c& w$ s0 f; ~PLNSOL,S,EQV,0,1.0 !等效应力" s3 O1 U) Q& ^6 Z5 C! ]( u6 e1 N
FINI
! w: u& d! k3 Z: Z4 v; }/POST26
( K3 ?3 O) p ?$ w& v1 ~5 uNSOL,2,NE/2+1,U,Y,UY_M !取出梁中部节点Y方向的位移' P$ _7 B+ O8 j
XVAR,1 !定义时间变量为坐标横轴
+ ^& K" J2 J# ]PLVAR,2 !绘制位移随时间变化的图形
' _; K( ~" P$ [3 U' R aDERIV,3,2,1,,,,,1 !对梁中部节点的Y向位移进行微分运算，得到Y向速度" T( [! V" e' M4 k% k$ n! o
PLVAR,3 !绘制速度随时间变化的图形
3 e7 d; W# M1 b) s0 U2 ZDERIV,4,3,1,,,,,1 !对梁中部节点的Y向速度进行微分运算，得到Y向加速度
1 h+ I6 m6 C8 i. k1 r7 aPLVAR,4 !绘制加速度随时间变化的图形& n' H3 O$ D, ^4 y& Z
L6 U) O" V3 W/ x( |
|
发表于 2011-9-30 15:18:32
楼主
