本帖最后由 小小实习生 于 2011-9-30 15:53 编辑 ) w+ ^) c: @: Z' F% b
8 G& h3 }* F5 F2 K! @' a
这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。$ Z6 J$ _( E1 b, e5 T
FINI6 T- K/ Y: f5 V \
/CLE
1 |# x3 }; w Y w" `5 M1 C3 a!设置参数$ x6 f6 ^, |( M' j, v; `
/CONFIG,NRES,2000: w+ }1 Q t! Y% P
*SET,LONG,40 !梁长度
% d# ?3 O9 J) v*SET,NE,20 !单元数
6 g/ F: S8 L2 X*SET,DL,LONG/NE !单元长度6 \& g4 o+ T, r/ Z* v! l
*SET,NN,NE+1 !节点数9 Z% p9 I& X, v$ D
*SET,FF,5000 !移动的集中力
\( g$ Q6 I) G5 u*SET,VV,80 !移动速度; E+ f$ u5 r$ |3 r6 S
*SET,DT,DL/VV*3.6 !移动经过一个单元所需时间4 m, k/ A/ g& ?! v" w& v
!定义单元材料实常数
+ A D/ G8 b7 h/PREP7
1 f& g& Z; \6 S# AET,1,BEAM3 !梁单元
" l6 d' m! K( e8 M4 |" pMP,EX,1,50E9 !弹性模量; X- R9 [4 p. ]/ d; Y
MP,DENS,1,2500 !密度
- h5 u; f$ Q% o, _MP,NUXY,1,0.2
! r, P9 [8 z; K. H! RR,1,0.5,0.041667,1 !几何常数
+ e3 c: R# O& @9 S. O! a!建立模型
9 E' v( j, ~" f7 D4 w2 M*DO,I,1,NN,1( T6 s7 T C6 a9 h, a- P8 a V
N,I,(I-1)*DL !建立9 A6 {3 Y0 t* t
*ENDDO$ u+ I. N" U% L4 A) _: |
*DO,I,1,NE$ L: d4 r' u) n, L* \' u
E,I,I+1
' \ H. B5 b# W8 T5 Y) p: S& w*ENDDO2 g' n( s8 d! @3 ^
FINI
U! U* ]$ n* K* S" G!加载与求解+ _2 v2 g" l0 S; G' J; o& n. ]
/SOL: [7 u( T; z0 ^3 o# M4 B/ ?
D,1,ALL9 {4 j5 S# J. K% n+ Q
D,NN,ALL: Y4 k6 y4 z) O
ANTYPE,4 !瞬态分析% L! k* a3 Q0 `7 X; C4 s+ n; V
TRNOPT,FULL7 e. B. Z- I! p+ R6 Y( t
TIMINT,OFF; W) O' P- g/ f+ p" t
OUTRES,BASIC,ALL !输出基本项，每一步都输出
# o# P7 u& d1 T& K: SKBC,1 !指定载荷类型; W' S( x0 F0 @) d, q
TIME,1E-10 !指定载荷步结束时间
1 E' |4 T; P: \5 M# iNSUB,4,8 !指定载荷子步数为4，最大子步数为84 C6 J/ h' B1 {* I6 L. a- m' S1 O
SSTIF,ON !打开应力刚化效应
; }9 d1 r4 f8 S, Q7 XACEL,0,9.8,0 !施加重力加速度
; V2 \' y o- w) i0 WALLS !选中所有元素/ [: _2 A( [4 b# J# M2 R( s
SOLVE !求解第一个载荷步 O. G# r( d% K- w( x$ W7 F
TIMINT,ON
% }$ [$ u4 e% S6 I% D |; Q' p$ g o*DO,II,1,NN,1- a O$ V |1 U* Z; {. @3 L
OUTRES,ALL,ALL !输出所有加载子步的所有结果
6 Y* Z: v6 d7 h. F( x* ]; cTIME,II*DT !每个荷载步的求解时间: T2 p* M) R* a* B4 I7 b
NSUBST,5 !每个荷载步包含五个子步& a* V( |+ ~6 N3 ~: Q
FDELE,ALL,ALL !删除所有的集中力荷载% O& |5 ]" f% H. o' f
NSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点; Y4 F- h6 X) K
F,ALL,FY,-FF !施加集中荷载/ T; w' W* R# \, r# M
ALLSEL !选择所有
' | m4 S0 H" r+ J3 ?3 fSOLVE% C! z7 B) h% L
*ENDDO
7 V8 V8 f' l' }1 P7 jSAVE
2 i3 y' I V! f; _: lFINI# Z/ v/ C( j# ^6 g7 i
/POST1
% ?3 F% F+ C" |" F1 jSET,NE/2+1,LAST !设定读取的载荷步/ g6 J) f0 b0 O, L I' u
/ESHAPE,1 !显示梁形状5 P& ^( H/ e1 |3 u8 b
PLNSOL,U,Y,0,1.0 !Y向位移
, P9 s. A9 E: f0 l: m3 I4 d UPLNSOL,S,EQV,0,1.0 !等效应力( D9 X( d3 r4 Y. I; z2 v( _
FINI5 ` T) }1 H6 |
/POST269 q% w+ J# S- f* l- x. Q
NSOL,2,NE/2+1,U,Y,UY_M !取出梁中部节点Y方向的位移/ Z x& v( j. R
XVAR,1 !定义时间变量为坐标横轴' P6 H* a: r0 p
PLVAR,2 !绘制位移随时间变化的图形: x5 @# u3 b. q7 o& h
DERIV,3,2,1,,,,,1 !对梁中部节点的Y向位移进行微分运算，得到Y向速度6 X! q5 b$ _/ z& E5 l, o
PLVAR,3 !绘制速度随时间变化的图形 z* f8 [0 U( v2 }% a
DERIV,4,3,1,,,,,1 !对梁中部节点的Y向速度进行微分运算，得到Y向加速度1 N7 m+ c4 |4 V* j
PLVAR,4 !绘制加速度随时间变化的图形
% E3 d7 D( C2 d: w; t5 l
0 o% u( v% N! V3 ~! a( o+ f3 L |
发表于 2011-9-30 15:18:32
