|
2#

樓主 |
發(fā)表于 2011-9-30 15:37:04
|
只看該作者
本帖最后由 小小實習生 于 2011-9-30 15:53 編輯 . k# e: ^2 W- W* _4 V
k1 p. A+ P; T' g5 i; i這個是經典的簡支梁的分析,,移動載荷,,但我做出來的結果不是怎么的理想。; E" J* G4 n- N8 ?: f
FINI& i8 P! Y) K1 W! n' b
/CLE) a, V$ v* n* k& v7 A5 h2 @
!設置參數(shù)
" Q3 f$ t8 z# y: W0 y1 {/CONFIG,NRES,20003 j! G& n0 W) Z5 H6 W3 k' v% |
*SET,LONG,40 !梁長度. \/ L9 _8 ?/ E# [
*SET,NE,20 !單元數(shù)3 ^5 x% Q8 ]1 Z) ` R' w4 i
*SET,DL,LONG/NE !單元長度
5 J4 m {, R. [/ J, A# j*SET,NN,NE+1 !節(jié)點數(shù)
# ?7 x' c _1 t*SET,FF,5000 !移動的集中力! T. Y+ Z$ `, z9 d5 w6 ^# F+ t% A5 U1 ]' ?
*SET,VV,80 !移動速度 }1 |% I% p. k% z( D. K
*SET,DT,DL/VV*3.6 !移動經過一個單元所需時間+ B( b* `4 h# z6 R3 ^
!定義單元材料實常數(shù) h# s: e# V) ]4 y3 T
/PREP7
& Q6 a% R6 e/ d+ ~4 f. }ET,1,BEAM3 !梁單元1 w+ Q! b& u! j5 r
MP,EX,1,50E9 !彈性模量
9 n) I, Q6 y) O! {+ ]4 G( N! NMP,DENS,1,2500 !密度, T5 G' j1 H, G- ]! h$ Z$ V) t
MP,NUXY,1,0.2
2 V G7 @ U b7 B* _R,1,0.5,0.041667,1 !幾何常數(shù)7 w% ^3 R2 ^+ r! |5 o7 q# F" q' u' U
!建立模型
2 n3 i8 a$ z( t) h- `) ~*DO,I,1,NN,1! ^5 w9 T/ e' K( T
N,I,(I-1)*DL !建立
$ r# L1 D( a N# l& l/ z9 ]*ENDDO+ P/ x2 M$ a, |
*DO,I,1,NE/ ^% ^1 A" V" M2 n: |
E,I,I+1
. w3 {3 B1 O3 v& H) e& c*ENDDO3 ^2 ^7 _/ e5 O$ W: Z' I' V
FINI, I5 d% A' K: C; W' {
!加載與求解
% k- p5 K9 H- Q5 N+ e. n- A/SOL
, ^( ]" s5 ^7 Y2 m' c$ Z7 u9 a% f, V" {D,1,ALL2 q" U: ^. q, h, b- v
D,NN,ALL$ H; u( h- b y" }0 y$ w: a
ANTYPE,4 !瞬態(tài)分析& ]0 L7 ^$ a* }+ ]' Z
TRNOPT,FULL+ r0 |" k' A3 A0 v" U$ ~4 }
TIMINT,OFF
& m2 _1 w& G4 h3 I" [- V+ TOUTRES,BASIC,ALL !輸出基本項,,每一步都輸出( q5 f2 k `7 [1 M
KBC,1 !指定載荷類型4 Y y2 {- S R/ L
TIME,1E-10 !指定載荷步結束時間: Z$ w L5 d) ~& s; p( D
NSUB,4,8 !指定載荷子步數(shù)為4,,最大子步數(shù)為8
( P0 H! S( |2 y. y: z1 FSSTIF,ON !打開應力剛化效應
$ K% g9 G0 i2 [# KACEL,0,9.8,0 !施加重力加速度
3 ?: z- l3 D% P1 G h* sALLS !選中所有元素
9 y& f: s" G& p- a1 @ H8 g- Z# ^SOLVE !求解第一個載荷步9 ]- G& m( m: O+ _1 w3 R5 O' Q
TIMINT,ON- F. r* c, Z( U+ ?- v9 Y; c9 Z
*DO,II,1,NN,1
, G# N+ d" \$ q& u. M S1 m% ]OUTRES,ALL,ALL !輸出所有加載子步的所有結果$ b* {- b1 a0 e! m/ K
TIME,II*DT !每個荷載步的求解時間3 p* ^# I5 }* T" Q: o
NSUBST,5 !每個荷載步包含五個子步6 s2 @+ t" J3 ?2 g5 G
FDELE,ALL,ALL !刪除所有的集中力荷載
* v1 O& i) ^; D9 X" a; |) eNSEL,S,LOC,X,(II-1)*DL !選擇集中力作用節(jié)點* h' j: C( f7 T; n; i' T0 Y
F,ALL,FY,-FF !施加集中荷載
3 U( j3 n. b6 s+ ?' Y/ P% KALLSEL !選擇所有$ Y. U! j! G& P: |" e
SOLVE' I1 ^8 W8 f, k. L) W) H0 ]
*ENDDO
# i0 u9 h& X0 A0 O5 m' U j. ZSAVE
5 y7 e4 ^( c, f5 q6 lFINI
) u9 O% d9 O* A( J/POST14 `: _5 L" {: Q. C3 @5 G
SET,NE/2+1,LAST !設定讀取的載荷步
" E% p# k5 }/ g/ESHAPE,1 !顯示梁形狀4 r2 U0 b6 z. M9 s) c
PLNSOL,U,Y,0,1.0 !Y向位移
3 V% {8 ?5 j3 |PLNSOL,S,EQV,0,1.0 !等效應力3 b) c$ f% X | j: m0 f
FINI
. o# @* O2 _# [/POST26 C+ i! j9 D, U
NSOL,2,NE/2+1,U,Y,UY_M !取出梁中部節(jié)點Y方向的位移
7 r* d0 ^8 j+ ^' L" J5 vXVAR,1 !定義時間變量為坐標橫軸
8 Y$ d+ f7 b4 ~PLVAR,2 !繪制位移隨時間變化的圖形
8 C0 F4 X4 y7 K# IDERIV,3,2,1,,,,,1 !對梁中部節(jié)點的Y向位移進行微分運算,得到Y向速度: D3 `3 F4 Q# o& W. t- h6 c3 j$ ~1 Z' s
PLVAR,3 !繪制速度隨時間變化的圖形7 o* K* X/ }: `1 J$ P; y' {3 v
DERIV,4,3,1,,,,,1 !對梁中部節(jié)點的Y向速度進行微分運算,,得到Y向加速度
' v2 K( C @4 p# v6 u3 KPLVAR,4 !繪制加速度隨時間變化的圖形
; H/ e) x! |( \; M% K
2 A6 g! a4 f2 c9 E" d/ R |
-
4.JPG
(28.91 KB, 下載次數(shù): 116)
下載附件
2011-9-30 15:39 上傳
這個是我做的一個懸臂梁中點的位移時間曲線,。
-
526.jpg
(30.62 KB, 下載次數(shù): 171)
下載附件
2011-9-30 15:52 上傳
這個是簡支梁的中點位移時間曲線,,哪位大俠指正下。
|