求助：懸臂梁移動載荷分析。

小小實習生

屬性：長度為20mm，直徑為1mm，材料為45號鋼，楊氏模量為2.1×1011pa，泊松比為0.3，密度為7.8×103kg/m³
5 m  j" I  Z/ o0 i' q6 G力：大小1.37N，移動速度為：3mm/s.從一端移動到另一端
求梁中點處的位移相應
. J6 g6 A7 Y8 S/ }3 c

小小實習生

5 T3 }; a5 V" ]0 l1 o: n* b: E0 J
這個是經典的簡支梁的分析，移動載荷，但我做出來的結果不是怎么的理想。
. Y. i  X2 J$ k# QFINI
/CLE
3 m0 j4 ?6 {5 `! v) j2 K" w5 N!設置參數
/CONFIG,NRES,2000
  n* ^# y1 [3 }! c  R9 ~" r*SET,LONG,40    !梁長度
+ E5 D/ J, J/ [8 }/ _; S*SET,NE,20    !單元數
3 R' U  ?1 K) S$ T! ^( u: k* b*SET,DL,LONG/NE    !單元長度
- @7 \2 R: v: X0 c* }& G" v*SET,NN,NE+1    !節(jié)點數
3 K( e9 B- p1 g( U" v+ U' a*SET,FF,5000    !移動的集中力
*SET,VV,80    !移動速度
*SET,DT,DL/VV*3.6   !移動經過一個單元所需時間
!定義單元材料實常數
/PREP7
8 _" L0 K, ?0 r+ e6 q& bET,1,BEAM3    !梁單元
: G7 W5 D' \: W: aMP,EX,1,50E9    !彈性模量
% Z6 u4 R, I. n. ZMP,DENS,1,2500    !密度
MP,NUXY,1,0.2
: e, O% O6 b4 h/ QR,1,0.5,0.041667,1   !幾何常數
% b4 ]: r# y, O; y3 t: U!建立模型
5 e: D" `2 y: [2 ?1 v' n, k*DO,I,1,NN,1
3 T) _4 C; a- X4 T. EN,I,(I-1)*DL    !建立
, `4 b" R9 Z; m9 G) n* Y( l*ENDDO
*DO,I,1,NE
* V- f7 ?9 b3 s; Y& V4 L% R, qE,I,I+1
*ENDDO
2 m+ P! D/ i3 WFINI
!加載與求解
7 n% k2 Z. `5 u" U0 Y  t/SOL
D,1,ALL
- E- ]2 {1 }8 R) P( q& @D,NN,ALL
ANTYPE,4    !瞬態(tài)分析
TRNOPT,FULL
! C0 z3 R4 _( U7 \, S- c, N! qTIMINT,OFF
OUTRES,BASIC,ALL    !輸出基本項，每一步都輸出
KBC,1       !指定載荷類型
# V; P# z& W( e) ]4 G* Q$ H" R/ rTIME,1E-10     !指定載荷步結束時間
- \" M; l, u9 b) Z8 dNSUB,4,8      !指定載荷子步數為4，最大子步數為8
+ `" [  V8 j9 k1 HSSTIF,ON      !打開應力剛化效應
" e7 G5 A$ r% `ACEL,0,9.8,0     !施加重力加速度
: \7 t# z  c0 h, d, |ALLS       !選中所有元素
SOLVE       !求解第一個載荷步
2 u" t, B( `9 t$ m$ @* k, o5 `TIMINT,ON
1 ^8 k8 Y4 }. v. r2 Y3 W*DO,II,1,NN,1
OUTRES,ALL,ALL    !輸出所有加載子步的所有結果
TIME,II*DT    !每個荷載步的求解時間
NSUBST,5    !每個荷載步包含五個子步
FDELE,ALL,ALL    !刪除所有的集中力荷載
% z" b5 a5 d; P* u( b" f9 BNSEL,S,LOC,X,(II-1)*DL !選擇集中力作用節(jié)點
F,ALL,FY,-FF    !施加集中荷載
ALLSEL    !選擇所有
SOLVE
*ENDDO
' ?+ l6 C8 N3 N+ Z0 l$ j& i! wSAVE
/ R0 _$ }( J- e2 cFINI
7 y2 E3 p& _" i7 V9 C2 \7 F/POST1
* ]' b# A: u# R0 _5 B) @( ySET,NE/2+1,LAST    !設定讀取的載荷步
# k5 ~- A% ^6 I, J2 T, n/ESHAPE,1    !顯示梁形狀
% M" n3 l3 C) \1 |PLNSOL,U,Y,0,1.0   !Y向位移
8 p" Z, o# w; f/ |PLNSOL,S,EQV,0,1.0   !等效應力
9 i" T" f& H; |6 lFINI
/POST26
) i% E; s: Z' W3 f/ \NSOL,2,NE/2+1,U,Y,UY_M       !取出梁中部節(jié)點Y方向的位移
( t$ w" ^7 h8 K6 L5 f$ e# g( BXVAR,1                   !定義時間變量為坐標橫軸
PLVAR,2                  !繪制位移隨時間變化的圖形
+ b0 M2 K- ~( [. a4 Y! _3 sDERIV,3,2,1,,,,,1        !對梁中部節(jié)點的Y向位移進行微分運算，得到Y向速度
4 A! J4 L' b+ mPLVAR,3                  !繪制速度隨時間變化的圖形
  s6 x5 d6 Q% p0 R, FDERIV,4,3,1,,,,,1        !對梁中部節(jié)點的Y向速度進行微分運算，得到Y向加速度
PLVAR,4                  !繪制加速度隨時間變化的圖形

小小實習生

為啥沒有高手回應呢？？？？

eddyzhang

自己先手算一下，看看有多少誤差

卡夏撒

你就是高手啊，自己解決吧，我不會

李曉輝

高人呀

samyfish

高手