本帖最后由 小小實習生 于 2011-9-30 15:53 編輯
3 ~, q6 ]( X' {8 g
4 h: `# _6 X y) t N7 s3 J這個是經典的簡支梁的分析，移動載荷，但我做出來的結果不是怎么的理想。4 r# p* H; p C3 i; x& @/ Z# R8 G( T
FINI- l# e6 p4 C$ t0 }( b
/CLE! S% M0 k6 b( r% J* @+ I
!設置參數
+ \+ f* P5 }; J! P ]/CONFIG,NRES,2000% M) R6 ?/ b0 |. A- E2 @
*SET,LONG,40 !梁長度
3 o: q! N; H# j/ z*SET,NE,20 !單元數
. V% a' ]8 @% x A( T* z*SET,DL,LONG/NE !單元長度
* _- Q: w b3 d" L2 a4 B( _0 `*SET,NN,NE+1 !節點數
2 ~2 C7 g% p# _9 x" o2 n/ n: k*SET,FF,5000 !移動的集中力) y# W. `2 T1 U4 A. O. d8 t
*SET,VV,80 !移動速度7 Y; K0 M, o( y
*SET,DT,DL/VV*3.6 !移動經過一個單元所需時間
8 ~7 u. h |; N5 j' E' r3 j!定義單元材料實常數
* q/ W! Z& R9 F6 \/PREP7
3 @5 J' A' ?2 _ c3 Q6 c- cET,1,BEAM3 !梁單元' Q# U& h9 i2 @3 ?/ v- Q
MP,EX,1,50E9 !彈性模量
% l: L; ?" e8 LMP,DENS,1,2500 !密度
( @. T0 ~8 }3 s3 p2 |, nMP,NUXY,1,0.2
7 g/ Q) E0 S5 E. QR,1,0.5,0.041667,1 !幾何常數
1 B1 m; p ^% P+ q!建立模型
O' o) Z& B4 }*DO,I,1,NN,1
* a, W+ ~+ A7 [$ g/ E0 VN,I,(I-1)*DL !建立
9 H7 R+ `( x0 e8 D% s E! c*ENDDO* Z+ R; }" ~1 }7 I+ Y5 Y
*DO,I,1,NE2 A. \: D- j9 B& I
E,I,I+1
$ \- X9 z4 c" ~2 j u. C*ENDDO
8 x/ n) M4 P+ y$ K# w9 ?FINI' ~2 @. m7 b; ~; ^! [& X
!加載與求解
, P( o" z* J) U V0 O5 m/SOL
0 T- e4 N' _: b0 m1 s( zD,1,ALL7 p! }; e: M" e+ D
D,NN,ALL+ |: h" ^! u/ u7 v1 G9 B
ANTYPE,4 !瞬態分析
* i. x, r1 N* k, t( ITRNOPT,FULL* c$ `( c$ m- J4 t0 M6 y2 _! r
TIMINT,OFF6 \" D* f1 L( C9 o
OUTRES,BASIC,ALL !輸出基本項，每一步都輸出9 Y3 [$ X3 [: s$ ]
KBC,1 !指定載荷類型
% [, u2 D" `$ ?5 z) n! M& TTIME,1E-10 !指定載荷步結束時間2 Z6 r+ L4 Q* b
NSUB,4,8 !指定載荷子步數為4，最大子步數為8
5 \9 K) b- z8 }/ Q5 U1 y9 ^/ ZSSTIF,ON !打開應力剛化效應
$ W- h" V* F! Z: uACEL,0,9.8,0 !施加重力加速度" x) D. ]& @) R% b A5 Z+ \# h+ i
ALLS !選中所有元素
; ?& P- l! t# |" Y- V$ }SOLVE !求解第一個載荷步
* U0 g: B. B5 U" }0 h0 L' MTIMINT,ON9 g9 S; E. H: f
*DO,II,1,NN,1" Y5 q3 n' v2 B% G% t: Z4 u, z3 R- J
OUTRES,ALL,ALL !輸出所有加載子步的所有結果' V% `1 r+ B4 o9 M9 L: W2 o' `
TIME,II*DT !每個荷載步的求解時間, m; m7 @( V, k+ ~
NSUBST,5 !每個荷載步包含五個子步8 N$ o# ~$ I1 e! \, F3 l, @
FDELE,ALL,ALL !刪除所有的集中力荷載& z: h. s# S8 e* @& E
NSEL,S,LOC,X,(II-1)*DL !選擇集中力作用節點2 Z7 h+ C* [+ Z$ h' W
F,ALL,FY,-FF !施加集中荷載
! |6 n4 \8 Y9 Q. WALLSEL !選擇所有
# C( V' ~5 a/ b" Y$ tSOLVE
8 ^) r, u3 C8 B. Q% _ p*ENDDO
1 v8 V" N/ c! v4 zSAVE3 {* N- x0 D: {
FINI. y9 }/ E/ u( I+ V* ]
/POST1
0 O' @8 z/ p" ySET,NE/2+1,LAST !設定讀取的載荷步
8 [: [9 F6 c4 C( H0 v1 B+ d5 d/ESHAPE,1 !顯示梁形狀
% a/ ]+ Z" p4 n* u& C0 C1 J# z) aPLNSOL,U,Y,0,1.0 !Y向位移$ k8 m% j8 }% s! F5 x$ v
PLNSOL,S,EQV,0,1.0 !等效應力2 Z+ \" Y+ l6 H1 b$ e
FINI0 ]" M5 q. W) c5 a g
/POST26
4 V$ A9 R- O3 n' u# g3 Q. j" mNSOL,2,NE/2+1,U,Y,UY_M !取出梁中部節點Y方向的位移
; H4 d) {/ g8 ^: ^2 bXVAR,1 !定義時間變量為坐標橫軸! }0 g4 K. _! s8 ]
PLVAR,2 !繪制位移隨時間變化的圖形 C- ]5 y% s; _ s2 M# K0 E
DERIV,3,2,1,,,,,1 !對梁中部節點的Y向位移進行微分運算，得到Y向速度; ^) x; H3 T, P3 e8 y# @% H
PLVAR,3 !繪制速度隨時間變化的圖形
* L5 r% k( w0 y* j! V1 x4 }, Q- kDERIV,4,3,1,,,,,1 !對梁中部節點的Y向速度進行微分運算，得到Y向加速度
. s C! t% c5 {( ?/ o3 E, APLVAR,4 !繪制加速度隨時間變化的圖形. B' p7 V* R- o( P3 @' Y
' f; G0 p1 E- c. B% W- o6 w7 ?% p: n
|
發表于 2011-9-30 15:18:32
QQ好友和群
QQ空間
