本帖最后由 小小实习生 于 2011-9-30 15:53 编辑
- Y+ ?+ r8 x0 H. a* m8 D! h0 ? x7 m
这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。. R1 u& g5 x% F! U$ ]7 Z0 h
FINI( a2 I4 D% G/ t( L
/CLE
' |9 O. _- Q( K. v# G/ t!设置参数
% _+ ~7 k2 D& |3 T/CONFIG,NRES,2000
2 }# v* r8 _# p7 P*SET,LONG,40 !梁长度
4 v) I% c v( U# S3 t; Q8 D: g' E*SET,NE,20 !单元数
8 z, @( F5 f ~+ p9 O) R*SET,DL,LONG/NE !单元长度7 ?; f" B; P6 [+ l6 n* Z. y% C
*SET,NN,NE+1 !节点数
& P. V1 G" w& m# h; R5 C*SET,FF,5000 !移动的集中力8 ~- \0 t9 M: E; m2 j; K% V5 V& {
*SET,VV,80 !移动速度# [" ~. I! ~( D' T, v
*SET,DT,DL/VV*3.6 !移动经过一个单元所需时间
7 |+ w5 L2 ~4 x6 i& S/ r, U% ~!定义单元材料实常数: a- X# a' U0 \ c5 D7 r3 b6 S: c
/PREP7) b; x0 |: D* w% E# M
ET,1,BEAM3 !梁单元
- a# B0 l5 t+ A8 C* `8 m5 U' `MP,EX,1,50E9 !弹性模量9 X: [$ ^+ }7 H
MP,DENS,1,2500 !密度
9 n. z6 F3 J0 D! T3 PMP,NUXY,1,0.29 r6 K6 p/ F: L2 P. t' M
R,1,0.5,0.041667,1 !几何常数 Z% }' Y& O; R- }% X
!建立模型
$ r, v. O' L6 a( P*DO,I,1,NN,1! I8 J! W9 ?4 p: G2 e' q( H
N,I,(I-1)*DL !建立
# K, d$ i. D% B; M5 d*ENDDO; q- H6 q, X+ u
*DO,I,1,NE
4 |- u* R# ?6 i, p6 X% x$ k* VE,I,I+1. ~' ~9 U6 U* C# S
*ENDDO8 S0 Z& `. \! S/ \ b: B% M7 P. y
FINI
. O& D; J, f* x, X!加载与求解
+ g5 Y* V0 j& D {8 S/SOL# p" o+ _0 Q& }8 o$ ~
D,1,ALL4 K! v* N+ a1 H: m8 N" l3 u4 e
D,NN,ALL
9 g+ j( M7 T: @* }* z+ nANTYPE,4 !瞬态分析2 [* L3 Q: R$ J: z
TRNOPT,FULL b2 h) r% N& `6 c
TIMINT,OFF
* Z( w. q! z% K3 `9 s1 v7 tOUTRES,BASIC,ALL !输出基本项，每一步都输出
- \7 u, @0 W, `8 t- ^3 S8 P& E, nKBC,1 !指定载荷类型
" _4 ]. L/ i, TTIME,1E-10 !指定载荷步结束时间
& }0 S# k" z) sNSUB,4,8 !指定载荷子步数为4，最大子步数为86 T4 T$ v1 v8 z/ Z3 O6 Z
SSTIF,ON !打开应力刚化效应$ a# f7 R/ r, t
ACEL,0,9.8,0 !施加重力加速度
. q: ^3 e" `6 w. QALLS !选中所有元素
, |+ T/ D/ o3 W/ ] t- DSOLVE !求解第一个载荷步
' c- T6 p# _. V) JTIMINT,ON
2 [4 y) S5 d& l% L) F1 k# R7 W*DO,II,1,NN,1
+ p1 p1 P0 j2 @5 `2 H& i9 }- |5 qOUTRES,ALL,ALL !输出所有加载子步的所有结果
0 ~4 i6 E( c L7 [TIME,II*DT !每个荷载步的求解时间
* H: w1 t' g9 {( ^% zNSUBST,5 !每个荷载步包含五个子步
G& C* T/ Q7 p2 B0 _4 pFDELE,ALL,ALL !删除所有的集中力荷载# T+ X6 Q m- X
NSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点
. @/ u8 R( y! \# ]3 i" OF,ALL,FY,-FF !施加集中荷载
6 i$ z3 F- c2 D4 p. ^ALLSEL !选择所有, r9 t# [/ E: i( m
SOLVE2 Z# C( S- b2 |
*ENDDO3 E8 F' J5 ^5 i/ E0 ~+ b# m9 V
SAVE" h" \3 a3 m3 t1 t
FINI
% F$ o g/ I. |5 }/POST1
. W; r. ]* o# d! `7 b! ESET,NE/2+1,LAST !设定读取的载荷步* @) ]' Q# q7 X
/ESHAPE,1 !显示梁形状7 W; \0 k7 ^' j1 n o
PLNSOL,U,Y,0,1.0 !Y向位移4 f6 G3 K8 G V5 F% l. r0 o
PLNSOL,S,EQV,0,1.0 !等效应力 H2 ^% B& w& h0 s+ ?
FINI
" @7 j/ |" q5 O2 p/POST26
) O& ~% H7 R( u% d+ F7 R9 mNSOL,2,NE/2+1,U,Y,UY_M !取出梁中部节点Y方向的位移2 r- K8 {1 ^) Z2 Z1 @
XVAR,1 !定义时间变量为坐标横轴
1 }% T; Z3 h- j* hPLVAR,2 !绘制位移随时间变化的图形
`9 L* c* g$ `3 q4 [; `4 Q2 J4 {DERIV,3,2,1,,,,,1 !对梁中部节点的Y向位移进行微分运算，得到Y向速度
& S+ U8 y M6 ]* V* d3 WPLVAR,3 !绘制速度随时间变化的图形; `- Z2 G2 j: R7 {+ [+ N
DERIV,4,3,1,,,,,1 !对梁中部节点的Y向速度进行微分运算，得到Y向加速度4 x% [ c% j2 |5 u* h/ A: H
PLVAR,4 !绘制加速度随时间变化的图形, x3 l& `) V; ^# U y# m
: ]1 ?$ l$ W/ p# a, h# W8 y |
发表于 2011-9-30 15:18:32
