本帖最后由 小小实习生 于 2011-9-30 15:53 编辑 , x) w, c+ F9 r- k7 S! u
) ]% y$ Y" E$ Z H# f9 p& p* S
这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。- E6 |1 y! u2 E, g( a, e+ p2 s7 B
FINI
7 o4 _5 E0 J) D+ R: B" A" X/ O6 G/CLE
, q0 P) S! m- |2 j!设置参数
6 G0 f0 C0 q2 ?, e" d/CONFIG,NRES,2000: L) c* ]; I( N6 R
*SET,LONG,40 !梁长度
, C+ G6 w+ p+ F9 z0 V4 i6 I6 ^*SET,NE,20 !单元数
8 f# u2 I( C, |# H*SET,DL,LONG/NE !单元长度; S a5 A! z9 U5 M. B- A
*SET,NN,NE+1 !节点数0 x. }8 x5 u9 @% T( w' |
*SET,FF,5000 !移动的集中力; k9 p/ R R3 G, S% a2 @# m3 V8 L( x
*SET,VV,80 !移动速度
, P1 L$ R! i0 b6 e*SET,DT,DL/VV*3.6 !移动经过一个单元所需时间; ^. R6 \7 i+ z2 S: A
!定义单元材料实常数$ j9 p- x; M" {0 x2 O
/PREP73 |" y. w8 w0 q$ A: m5 h
ET,1,BEAM3 !梁单元6 f3 D4 G1 w6 p4 y4 R1 O5 {( N
MP,EX,1,50E9 !弹性模量
6 S# r- t. z0 n5 G# s9 fMP,DENS,1,2500 !密度8 i# [! g, M! D7 j6 B
MP,NUXY,1,0.2
. i3 x2 D7 _6 X: eR,1,0.5,0.041667,1 !几何常数
' u( T2 i5 C n1 {!建立模型( f4 a9 K0 d) ^- ?" t7 y6 q
*DO,I,1,NN,1
# o/ s1 v9 O% B: S: p! PN,I,(I-1)*DL !建立
- d* P. I/ X d9 P*ENDDO
2 x8 D9 `* @5 d8 Q7 y5 n$ N0 I*DO,I,1,NE
4 {$ B \# D# e, UE,I,I+18 J4 _6 F, ^9 R
*ENDDO
4 ~. Q# g G7 G. }. Y4 m/ t$ S0 `FINI3 R% _1 ]1 P& v
!加载与求解
6 a# P+ x! ~$ N: U/SOL' U5 n1 a! @& g$ S# N8 W1 b; @
D,1,ALL
% R2 Z& r* y9 B4 v8 w# w# n& TD,NN,ALL
4 s7 E6 Y& U7 b( vANTYPE,4 !瞬态分析
* q& D" ~, T, a, RTRNOPT,FULL8 B$ x% f( H i
TIMINT,OFF
9 o5 y! @$ H! N IOUTRES,BASIC,ALL !输出基本项，每一步都输出
6 x. [, Y0 h8 @KBC,1 !指定载荷类型
, l- j5 H6 Z, g( g. `+ q T! e5 eTIME,1E-10 !指定载荷步结束时间9 b \* f8 F5 f& B
NSUB,4,8 !指定载荷子步数为4，最大子步数为8& z, o& F6 _$ B5 A/ V, D
SSTIF,ON !打开应力刚化效应5 R, s- T; Y8 a7 f
ACEL,0,9.8,0 !施加重力加速度
+ ?, m: k2 ]$ \0 PALLS !选中所有元素* k/ F6 K. r% W
SOLVE !求解第一个载荷步
4 S! i9 s) G& V" v, M: i; m# DTIMINT,ON
" m. R9 X l2 z: v! g*DO,II,1,NN,1
2 E& O$ Y: J0 L4 f+ @5 G+ V" GOUTRES,ALL,ALL !输出所有加载子步的所有结果
6 R& z5 j. Q$ x m* ?TIME,II*DT !每个荷载步的求解时间
9 X. \0 q: h( U L2 v/ ]NSUBST,5 !每个荷载步包含五个子步
# |# s j& S9 LFDELE,ALL,ALL !删除所有的集中力荷载& E) x6 s8 L/ G7 l( L
NSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点1 i! u! H2 B( |2 f- |
F,ALL,FY,-FF !施加集中荷载$ g" p. A3 @! F6 e3 I9 {( F( ^, P
ALLSEL !选择所有
0 n! V' ?# x6 pSOLVE
2 _# F f. i! W2 S*ENDDO
. l" n `+ B* Z) MSAVE
8 Z8 R I |# O& j/ t' FFINI" S% O0 I7 w* Z1 l1 Q1 V1 n
/POST1
1 i3 m3 \- ^% T0 u9 [, s2 z2 ESET,NE/2+1,LAST !设定读取的载荷步
) `* n s8 Q2 J- r/ESHAPE,1 !显示梁形状
( t/ x9 a! n* l' GPLNSOL,U,Y,0,1.0 !Y向位移0 n4 D' R/ b3 H
PLNSOL,S,EQV,0,1.0 !等效应力
H& D7 J( I7 V5 { G7 kFINI5 v9 w" M; K5 F4 P9 a6 i# g: P# J
/POST26! c, ]# D [! @
NSOL,2,NE/2+1,U,Y,UY_M !取出梁中部节点Y方向的位移
1 G/ C4 B; h+ {/ B+ e; E6 X* V {XVAR,1 !定义时间变量为坐标横轴
; E# Z! z5 D& S; U# ePLVAR,2 !绘制位移随时间变化的图形
! m/ \3 R4 w+ r2 tDERIV,3,2,1,,,,,1 !对梁中部节点的Y向位移进行微分运算，得到Y向速度
6 |' l* j, v' tPLVAR,3 !绘制速度随时间变化的图形0 r- b9 t/ T- Y V
DERIV,4,3,1,,,,,1 !对梁中部节点的Y向速度进行微分运算，得到Y向加速度
! ~" v G9 ?. I: CPLVAR,4 !绘制加速度随时间变化的图形
* G' ?7 P9 l f, I+ Z5 L! O' o# s$ L/ A) f* x8 E
|
发表于 2011-9-30 15:18:32
楼主
