本帖最后由 小小實習生 于 2011-9-30 15:53 編輯
# M2 G# R+ H( v$ m2 V( R! B" {3 ]2 F' _1 s' p/ o( J# C7 W M8 I
這個是經典的簡支梁的分析，移動載荷，但我做出來的結果不是怎么的理想。
, ^7 f2 [: ^# f0 D7 G. R1 GFINI
' v: l# Q$ u9 J4 g/CLE
: M+ B. r. h# d# A* u% P/ `!設置參數
! W5 \" F( T Z$ V, }& W; {( I/CONFIG,NRES,2000
+ ~6 }: A+ Z; V( }. S, w# J*SET,LONG,40 !梁長度6 `5 F" J4 m% G a8 _8 D4 p
*SET,NE,20 !單元數8 F* g8 H* U6 J9 A: T$ ^+ E
*SET,DL,LONG/NE !單元長度6 Q, e3 F" G# H& e& s/ L8 K
*SET,NN,NE+1 !節點數
$ j: c/ ]# m5 Z: ?9 x+ [5 j3 t3 B* t*SET,FF,5000 !移動的集中力; S0 x! w! h" ]4 \6 {/ t
*SET,VV,80 !移動速度
; _: p; Q! z$ U t+ S- Q- v*SET,DT,DL/VV*3.6 !移動經過一個單元所需時間
8 c& ?1 r: Q) s0 ~, f Q( |/ R' A!定義單元材料實常數; l$ @, P1 J; \9 E- r2 y* }
/PREP7
7 u# G$ c) M; g% ~/ qET,1,BEAM3 !梁單元
4 l* v9 c7 n* S2 [$ P- J' eMP,EX,1,50E9 !彈性模量3 f# E/ x {' [7 s, Y
MP,DENS,1,2500 !密度
: d. }8 v; |* L' Y2 a# iMP,NUXY,1,0.2
% }0 x$ o: w- M; M, F: h" MR,1,0.5,0.041667,1 !幾何常數
+ F w( N* w8 U7 u# F!建立模型 \$ C4 m# U9 b. z# u
*DO,I,1,NN,1
d$ H" Q5 N( f4 {4 dN,I,(I-1)*DL !建立: t* t! N" n+ [/ Y- z- B
*ENDDO
0 S& D7 k1 Y6 H" f3 c" M/ }*DO,I,1,NE
/ R* X; n+ I: K# VE,I,I+1
, w) z z+ `5 k) C. l*ENDDO
+ }7 z" D" \! q8 A. V$ e4 z# m: ZFINI
" B! L* s* L: Y! ^!加載與求解6 f9 v9 u( J7 ^, n+ h5 b
/SOL1 [" Q6 V; a/ D8 Z
D,1,ALL m# k& j7 r, ?. c8 W
D,NN,ALL
. ` U g# J( U) S1 P% nANTYPE,4 !瞬態分析
$ ?1 }$ c9 o6 m# w9 V, k% z$ JTRNOPT,FULL
6 \1 L$ M$ L5 @5 c2 vTIMINT,OFF7 k8 e; h8 t: W( r5 I/ U
OUTRES,BASIC,ALL !輸出基本項，每一步都輸出" P) @: |3 s& H) ], |
KBC,1 !指定載荷類型
8 ^0 a5 ]& v0 ^$ nTIME,1E-10 !指定載荷步結束時間
& H( @, g; ^+ KNSUB,4,8 !指定載荷子步數為4，最大子步數為8* P$ X) S- A% r* [" ]
SSTIF,ON !打開應力剛化效應/ |. J4 v( p; B. w9 H. d- |# j
ACEL,0,9.8,0 !施加重力加速度. h2 N) V$ {0 l
ALLS !選中所有元素
n2 I: \' T0 Z7 F1 W: USOLVE !求解第一個載荷步
6 L2 F# E L- g8 D0 r- {TIMINT,ON6 R+ q, H) }: X, y8 B% g
*DO,II,1,NN,1
: _# N* b2 a+ Q+ T/ s4 oOUTRES,ALL,ALL !輸出所有加載子步的所有結果
; d+ X6 k5 I& S% \TIME,II*DT !每個荷載步的求解時間8 Q1 H1 R K! c3 x6 K1 w2 J- v$ h
NSUBST,5 !每個荷載步包含五個子步1 M: W1 s' V* L8 T$ G
FDELE,ALL,ALL !刪除所有的集中力荷載; ]! b. D3 M1 P1 o' g1 i6 o
NSEL,S,LOC,X,(II-1)*DL !選擇集中力作用節點
# y- {3 X( l; k! [; }2 [F,ALL,FY,-FF !施加集中荷載
% U9 _; s% B9 {& C" T, wALLSEL !選擇所有3 e- J* w, b2 }0 |$ G( b' d3 s
SOLVE
4 _6 e( B, B/ q# {! _2 ]( b*ENDDO
]; d$ V$ v0 }, Y! h# s2 @' ~SAVE
! y' I/ w* q( T, N6 @+ r8 MFINI4 D2 Y! S9 {7 g8 y1 j3 X
/POST15 \! ^% Q0 ^3 m. q% _
SET,NE/2+1,LAST !設定讀取的載荷步
& y- Z4 ^" R8 c, o" x s7 E6 P/ESHAPE,1 !顯示梁形狀/ p7 f) x# x4 g6 L9 F* Y
PLNSOL,U,Y,0,1.0 !Y向位移# u5 ~; h) G; i
PLNSOL,S,EQV,0,1.0 !等效應力; E, B7 \; f' ^. q
FINI5 O) g+ t) F8 X" U
/POST266 K; J a- o: ~1 h7 Z$ ~
NSOL,2,NE/2+1,U,Y,UY_M !取出梁中部節點Y方向的位移2 l! T* Q. A5 m5 T" ^/ N* b5 `4 h5 I% D
XVAR,1 !定義時間變量為坐標橫軸
/ Z' G1 c7 Z! Y# d6 ]; ^PLVAR,2 !繪制位移隨時間變化的圖形
( v6 \: p/ ~( ?1 s2 [0 g3 g4 @9 g# bDERIV,3,2,1,,,,,1 !對梁中部節點的Y向位移進行微分運算，得到Y向速度
" W% L6 @ T6 c1 r% p- C IPLVAR,3 !繪制速度隨時間變化的圖形
0 O9 H0 l# P% w- U4 Y& y U8 ^ z* bDERIV,4,3,1,,,,,1 !對梁中部節點的Y向速度進行微分運算，得到Y向加速度+ q/ C- E- V2 a5 G/ z# w
PLVAR,4 !繪制加速度隨時間變化的圖形
* \# B3 ^! F6 ~9 j1 H: c9 \* D1 D* W$ W# G, A
|
發表于 2011-9-30 15:18:32
QQ好友和群
QQ空間
