本帖最后由 小小实习生 于 2011-9-30 15:53 编辑 % A; e# v( }/ s& H7 Q
2 V# E3 \+ M$ Q! G5 R$ d8 x
这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。
: C' {7 f" b e3 f6 L9 m6 N' GFINI# c. c+ b/ v g4 r% T+ L
/CLE3 `7 A2 Q# {; H- ]4 j- X+ I; c! v
!设置参数% n9 a9 T6 Z$ K
/CONFIG,NRES,2000
3 C2 X5 Z6 f/ ]( f ]. h*SET,LONG,40 !梁长度; S6 R5 X' F9 J R0 M" r% E
*SET,NE,20 !单元数; ?6 M; t2 \. P3 D/ H! C
*SET,DL,LONG/NE !单元长度5 b" {4 V, I3 f5 i* w5 L
*SET,NN,NE+1 !节点数
. j* U( h1 j* |! W*SET,FF,5000 !移动的集中力: ^2 o; n4 R" D5 P1 E
*SET,VV,80 !移动速度
2 ^% o5 s* k; q' c! q4 a @6 z; h*SET,DT,DL/VV*3.6 !移动经过一个单元所需时间
0 U" b4 H4 o& R8 n6 Y' B' V" B! P6 x!定义单元材料实常数; @. E: q( n1 w5 `8 z, p: l! _
/PREP7
# C. ~2 F! W& [8 r; a4 t' yET,1,BEAM3 !梁单元; Q$ B% ]. s; E- o
MP,EX,1,50E9 !弹性模量
+ T# n9 i4 I. a) e- {MP,DENS,1,2500 !密度
4 [, S) h2 \8 x" P2 ~MP,NUXY,1,0.2
1 d2 i$ r- ]: p7 R L0 b. z$ T- UR,1,0.5,0.041667,1 !几何常数1 `4 I4 c9 H$ f+ }
!建立模型& X8 B$ e# z) U
*DO,I,1,NN,1
$ T! \& Q& o" [N,I,(I-1)*DL !建立- g; j2 G8 t9 w+ A$ u' ^. D V5 y
*ENDDO( x' w* `# B6 C& @! |( }
*DO,I,1,NE
) Z1 F0 B# J3 `# V* W% Z( Z/ wE,I,I+1
3 Y. X( I! C6 u*ENDDO [1 ?# o# `) I
FINI
1 c2 {" Z/ g2 g E!加载与求解( ^# @/ f8 h! K3 A; V9 v
/SOL
. z0 K& ]3 z5 U5 w# e5 ~" oD,1,ALL
2 S, U4 z- U* v- wD,NN,ALL
( @1 S! w& D' {" _ANTYPE,4 !瞬态分析
6 x) s% \/ F1 x. o, A; P6 P/ g% WTRNOPT,FULL
" h- ~" v3 u( m$ g+ a9 Q# R9 ?* s9 GTIMINT,OFF$ ]1 z8 X' R' r& k! W! r4 K! v' ?5 p
OUTRES,BASIC,ALL !输出基本项，每一步都输出
" U$ [' H6 Z5 B1 G- R j5 l' ^2 `KBC,1 !指定载荷类型- B; ?8 {3 j' `5 D
TIME,1E-10 !指定载荷步结束时间/ C6 A/ {% L3 a4 z t P6 g' U/ d
NSUB,4,8 !指定载荷子步数为4，最大子步数为8
* ~4 @1 M- K5 @$ k2 ZSSTIF,ON !打开应力刚化效应
1 f2 G8 z; M/ j/ l: `1 AACEL,0,9.8,0 !施加重力加速度0 }- O" v3 p; D. X) l! N* I$ A
ALLS !选中所有元素
' q4 p$ F& Y8 }3 P% m5 {SOLVE !求解第一个载荷步
: P, f* w7 w! E6 D; vTIMINT,ON
5 t; n3 D2 r: o* `- f* ^3 E6 J*DO,II,1,NN,1
% N$ {0 D! g' K/ a5 q* V4 R7 jOUTRES,ALL,ALL !输出所有加载子步的所有结果
; B: e, l+ |0 L4 N. u$ ~# ^7 _# oTIME,II*DT !每个荷载步的求解时间
' @& t( ]% J/ XNSUBST,5 !每个荷载步包含五个子步
! c6 T0 t; z) {) n4 k8 U0 HFDELE,ALL,ALL !删除所有的集中力荷载) U1 v4 H2 a* R! R
NSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点
: U6 G; V! U) i7 ^' i! GF,ALL,FY,-FF !施加集中荷载6 d) I5 k: r* ?6 k' l
ALLSEL !选择所有! d0 L$ B$ U/ @
SOLVE
! ^1 T8 g7 K& w*ENDDO( |2 n2 F# h9 I1 h q5 ~" q
SAVE: u1 y% g1 D! E. \; c3 l. q; L
FINI
1 r+ ~! J' H+ o' ]/POST1
/ h5 ~+ |/ {' Z: d# cSET,NE/2+1,LAST !设定读取的载荷步+ c [( ?- H% y; {3 r
/ESHAPE,1 !显示梁形状
5 ~3 z, @6 y) i( O: bPLNSOL,U,Y,0,1.0 !Y向位移$ P' g( n4 f6 U
PLNSOL,S,EQV,0,1.0 !等效应力# e5 ^+ X2 g8 K4 Q; e
FINI; ^$ D& p' X9 p* e" w+ G% m- u6 U8 |
/POST26
1 G2 R, d# L# q' p! N% H) FNSOL,2,NE/2+1,U,Y,UY_M !取出梁中部节点Y方向的位移- r' W4 C; e, z4 B4 J
XVAR,1 !定义时间变量为坐标横轴7 W' H: F5 W; A. r, |% x
PLVAR,2 !绘制位移随时间变化的图形* {3 {7 T/ b7 R1 L) U
DERIV,3,2,1,,,,,1 !对梁中部节点的Y向位移进行微分运算，得到Y向速度
2 U. z- Y/ t/ h& lPLVAR,3 !绘制速度随时间变化的图形
0 Y" D7 i( X6 T# M/ T/ c) x& g6 nDERIV,4,3,1,,,,,1 !对梁中部节点的Y向速度进行微分运算，得到Y向加速度4 K+ ~( ?4 n) g" H ?6 c2 A
PLVAR,4 !绘制加速度随时间变化的图形
v7 d1 \. t; R1 w. |
/ ?0 b. Y( B( B; ^' v5 I9 y! M |
发表于 2011-9-30 15:18:32
