本帖最后由 小小实习生 于 2011-9-30 15:53 编辑 : g/ \% Q9 S7 m O$ F
/ \2 o1 a) l$ ]* z& Z8 E; Y这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。
* _" l: Y9 [- I4 H+ z! W( H' b+ JFINI
# r9 ?9 C% x! l' g2 C: b- }0 L/CLE$ A* d, H" [. T! m/ u0 y7 s
!设置参数7 }% `# {6 Q) w [* m$ }% j
/CONFIG,NRES,20008 R/ F' i, i0 ]' y
*SET,LONG,40 !梁长度! p# X/ p! @. S. o; c
*SET,NE,20 !单元数6 e% P6 j! I& H- J- l+ D3 x3 E
*SET,DL,LONG/NE !单元长度7 `; ~- @ x. G; C$ L1 q+ g
*SET,NN,NE+1 !节点数
8 E% i& p: W; \5 Q& H; ~*SET,FF,5000 !移动的集中力7 G$ j8 @* g3 x
*SET,VV,80 !移动速度
8 x! y: H% L* V*SET,DT,DL/VV*3.6 !移动经过一个单元所需时间/ b% B4 n# q9 r
!定义单元材料实常数& k3 l. ?* C( [% b3 K2 R; A+ Y
/PREP7
8 T* i4 f2 d% F! W' q QET,1,BEAM3 !梁单元
/ x* M+ G( A) A, w% F! }MP,EX,1,50E9 !弹性模量3 {4 K2 t: B4 D% E9 v. Z& o
MP,DENS,1,2500 !密度
J' p( e3 ^' O* PMP,NUXY,1,0.2
9 o7 _/ b- q, y) s# G" ^8 l2 s4 pR,1,0.5,0.041667,1 !几何常数/ c3 _$ {( P' W/ o
!建立模型
; K5 V$ m9 I# K" Z8 N*DO,I,1,NN,1
4 R2 x% w# M/ \ X% h! W1 M# ZN,I,(I-1)*DL !建立
& `% V- H( p# M*ENDDO
7 j4 {$ V0 ^4 g7 D*DO,I,1,NE
4 r: n j% U+ _7 |7 nE,I,I+1* o u1 J t. p: N. L
*ENDDO+ d, O+ S5 ^: h: L3 b+ f5 N
FINI& W' S2 y, X, c1 J8 j. |1 L0 r
!加载与求解' }7 j8 C" |1 ]5 C9 C
/SOL
9 v+ h2 ~3 M TD,1,ALL
/ z$ R# B N4 K' ]( o2 ]D,NN,ALL
& R. u- r+ I4 t5 k xANTYPE,4 !瞬态分析, w' b% P; q6 w+ h3 t. r' C
TRNOPT,FULL
- h6 W6 y3 x! \6 }* _/ uTIMINT,OFF/ w3 o4 O& g7 H3 |9 x
OUTRES,BASIC,ALL !输出基本项，每一步都输出
9 c2 ?) Q4 e: C7 B7 W- U$ MKBC,1 !指定载荷类型
4 ^6 Q u- G" g8 c5 o8 O. j1 bTIME,1E-10 !指定载荷步结束时间
% e- p0 g& Q' WNSUB,4,8 !指定载荷子步数为4，最大子步数为8
# @( j P/ p8 t4 h. Z6 BSSTIF,ON !打开应力刚化效应4 M4 {' e" q. R% R
ACEL,0,9.8,0 !施加重力加速度1 z4 `4 w3 o, Z1 |2 i
ALLS !选中所有元素) l+ [( R' H* g$ T+ ~/ C
SOLVE !求解第一个载荷步
9 y, Q, Y S7 i* d+ s& tTIMINT,ON
2 q- Q2 X) d7 P& d% O6 I+ ]*DO,II,1,NN,1: n, w7 W0 a& Z
OUTRES,ALL,ALL !输出所有加载子步的所有结果3 B! D2 b/ n, Z* {
TIME,II*DT !每个荷载步的求解时间
/ D* N6 [2 F' ]5 K1 j) BNSUBST,5 !每个荷载步包含五个子步' o$ j; E& y% W4 @9 _0 Y3 t
FDELE,ALL,ALL !删除所有的集中力荷载3 `, O) s& J* N+ K& w, v
NSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点
" F) w( M: s# u0 wF,ALL,FY,-FF !施加集中荷载* v* J1 _6 j7 ^7 e* P" l5 Y0 a
ALLSEL !选择所有
4 @. ^$ j2 M9 l" Q9 cSOLVE* N4 c* q1 Y. c$ g9 N
*ENDDO$ ^4 [+ Z( j0 [$ q
SAVE
4 e. d2 H/ Y. p5 y! N* ?+ B' m* V4 }FINI
% T' T8 H/ T4 |# P1 {/POST1. g$ B$ ]2 G# P; K& D0 q9 r9 |" F
SET,NE/2+1,LAST !设定读取的载荷步
5 b8 D$ c1 n1 e' I- z! C/ESHAPE,1 !显示梁形状9 d0 p8 f, q0 D4 z" ~
PLNSOL,U,Y,0,1.0 !Y向位移# [8 M8 `2 p; ^# f7 f5 a
PLNSOL,S,EQV,0,1.0 !等效应力
3 `( F1 A3 z# M, y' y3 b2 n. kFINI# c/ n x; o5 b" o- Q( _# C7 H
/POST26
" v. b0 a; k3 GNSOL,2,NE/2+1,U,Y,UY_M !取出梁中部节点Y方向的位移% W6 k# X& W1 n) d5 d; z" y
XVAR,1 !定义时间变量为坐标横轴
) m! r. k; M: q0 p3 o4 c# o' n# n# zPLVAR,2 !绘制位移随时间变化的图形5 d4 q" u3 S2 y/ P8 G$ d# v
DERIV,3,2,1,,,,,1 !对梁中部节点的Y向位移进行微分运算，得到Y向速度4 s( U& }' J5 y* D* c
PLVAR,3 !绘制速度随时间变化的图形
1 Z6 X* r5 m3 X2 s& VDERIV,4,3,1,,,,,1 !对梁中部节点的Y向速度进行微分运算，得到Y向加速度3 ]- y0 E" H/ [+ `
PLVAR,4 !绘制加速度随时间变化的图形( I5 s2 L( G7 J# i+ q
# n8 f/ @8 ` G4 ]0 ^, l n j |
发表于 2011-9-30 15:18:32
