本帖最后由 小小實習生 于 2011-9-30 15:53 編輯
% d V+ \8 s( b. Z" X6 K+ o9 h3 [- H* J: ^! p9 D& e$ l- x
這個是經(jīng)典的簡支梁的分析，移動載荷，但我做出來的結(jié)果不是怎么的理想。 y3 q6 M& K: f1 F* \0 C
FINI4 `6 x% ]1 D. z- | X6 R( e9 ]
/CLE3 Z8 x8 V, I! ?1 ~) Y' S/ n, q
!設(shè)置參數(shù); q \' e" W2 h
/CONFIG,NRES,2000
: R1 I4 X. h, O7 D# D*SET,LONG,40 !梁長度/ [8 D9 `' ?! _; R* {8 f
*SET,NE,20 !單元數(shù)
& @: [+ e3 E, D1 o* b0 {*SET,DL,LONG/NE !單元長度 P) r3 O' A1 \' o
*SET,NN,NE+1 !節(jié)點數(shù)
1 m5 P2 R. B0 p* Y8 S/ g*SET,FF,5000 !移動的集中力8 a4 M( ^7 l( u! o) s* G
*SET,VV,80 !移動速度# U* X1 q2 F1 T5 F4 S/ k0 F
*SET,DT,DL/VV*3.6 !移動經(jīng)過一個單元所需時間
# V" H; ]+ f. w. [$ T5 z8 i4 N2 v!定義單元材料實常數(shù)4 `3 J- o3 a; `& O4 b
/PREP77 `+ q/ X% U- U- |% M5 {! C% j& w3 ?
ET,1,BEAM3 !梁單元) o9 _; G2 n( {7 `. Z6 w
MP,EX,1,50E9 !彈性模量9 F9 |' N+ W# i8 E9 z" S6 O
MP,DENS,1,2500 !密度
- T( |. J, b/ a6 @1 eMP,NUXY,1,0.2
4 p c2 B7 n, v4 M; tR,1,0.5,0.041667,1 !幾何常數(shù)
X; a+ o7 p' F# P- |!建立模型
# i# `8 C. W! D3 r/ `*DO,I,1,NN,1# I7 U4 ^2 w* A3 ^# N+ d. d
N,I,(I-1)*DL !建立4 Q7 |* c* J! P! ]7 a6 g
*ENDDO
, g+ f/ W8 G( G" I/ c$ u; ?*DO,I,1,NE* A C% }2 ]' |/ c4 l& q; }
E,I,I+1
7 y4 {/ n l7 A- m4 X0 l9 m*ENDDO
. v9 O5 i" b* nFINI9 ]. E1 Q, `) N
!加載與求解* r3 M1 g6 l; B3 y% a) p
/SOL, B: ~/ S; k6 E( O0 O
D,1,ALL
- @% l9 _" \! ], h- R2 m8 C4 gD,NN,ALL
2 q. k! k: [4 y, l( L$ JANTYPE,4 !瞬態(tài)分析5 n( A2 r6 B$ V8 J$ p6 w V; o
TRNOPT,FULL8 Y/ l: ^( P$ c' ^, Y) Z/ l0 u
TIMINT,OFF
; O. n+ t5 Q8 }& D3 ?OUTRES,BASIC,ALL !輸出基本項，每一步都輸出) I5 S6 y! B/ n! G& G4 }; ~
KBC,1 !指定載荷類型
* i: }. t; P/ `- s2 ^3 b3 `TIME,1E-10 !指定載荷步結(jié)束時間" X9 i* c+ \: H4 j" H6 D0 o
NSUB,4,8 !指定載荷子步數(shù)為4，最大子步數(shù)為8
# A2 `9 W' Y1 U7 Y9 ^/ bSSTIF,ON !打開應(yīng)力剛化效應(yīng)1 _* B3 g) ~; t7 k, W
ACEL,0,9.8,0 !施加重力加速度
/ f% X+ }6 O' T# kALLS !選中所有元素
# q8 c w; G; F2 l5 R( V9 jSOLVE !求解第一個載荷步2 B' n- u1 V" p R* T
TIMINT,ON
5 k* y |$ v; L, P+ d& Y( B*DO,II,1,NN,1: T" }" W3 F4 g
OUTRES,ALL,ALL !輸出所有加載子步的所有結(jié)果
4 X/ v# u, a! `" \! B" P% NTIME,II*DT !每個荷載步的求解時間' B V6 ]4 z" e
NSUBST,5 !每個荷載步包含五個子步0 W% Q; P+ _' T7 m0 X) K7 \% L W* C
FDELE,ALL,ALL !刪除所有的集中力荷載( b$ C$ C8 L6 a: j$ p k
NSEL,S,LOC,X,(II-1)*DL !選擇集中力作用節(jié)點
- ?# Q$ } K# l* A/ I: I- BF,ALL,FY,-FF !施加集中荷載+ J1 G. y9 r) E- B1 E
ALLSEL !選擇所有8 N* K" O- Q* U6 o- T, y' b
SOLVE
) e( t7 H/ P' R* a*ENDDO0 _( N$ d9 e, |+ b
SAVE
6 ^, F" R( F5 ]4 G+ eFINI7 e9 M+ ~4 p) h6 }! s
/POST1% S# \/ E: }/ N; x! C8 v5 x
SET,NE/2+1,LAST !設(shè)定讀取的載荷步
5 T5 b6 w, t- d; ? q- S- O/ESHAPE,1 !顯示梁形狀
4 [/ Z R! C: K9 ?# q' R8 d: f6 F% GPLNSOL,U,Y,0,1.0 !Y向位移4 Y; G* U/ W, c3 h" M0 x4 ]
PLNSOL,S,EQV,0,1.0 !等效應(yīng)力
6 z3 A$ Q# W( @$ }7 H9 ~9 DFINI
: N6 W1 X: y% T! U/POST267 ^& Q$ m3 a/ {
NSOL,2,NE/2+1,U,Y,UY_M !取出梁中部節(jié)點Y方向的位移) I5 T9 B3 a3 O' U
XVAR,1 !定義時間變量為坐標橫軸1 m/ g% A; K2 S f. i% `- t7 ?0 O; T4 E
PLVAR,2 !繪制位移隨時間變化的圖形
# O$ C0 e0 C ?' ?" F! R1 v! GDERIV,3,2,1,,,,,1 !對梁中部節(jié)點的Y向位移進行微分運算，得到Y(jié)向速度
$ }0 E" a* [! W, E4 OPLVAR,3 !繪制速度隨時間變化的圖形
: o* u! E8 F% g/ T5 R. R7 X- p7 _DERIV,4,3,1,,,,,1 !對梁中部節(jié)點的Y向速度進行微分運算，得到Y(jié)向加速度
, k' i& w m* }5 A/ @6 O4 l( cPLVAR,4 !繪制加速度隨時間變化的圖形
- D* d% ^; I- D
5 F a1 a; r9 o9 q4 w, }. z2 s1 m! }7 A |
發(fā)表于 2011-9-30 15:18:32
QQ好友和群
QQ空間
