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

小小實習(xí)生

屬性：長度為20mm,，直徑為1mm,，材料為45號鋼，楊氏模量為2.1×1011pa,，泊松比為0.3,，密度為7.8×103kg/m³
0 ^8 k7 u- h. y" I) V" c. i力：大小1.37N，移動速度為：3mm/s.從一端移動到另一端
求梁中點處的位移相應(yīng)
0 L! @; f7 S7 M) c* S% u

小小實習(xí)生

( `, \2 D0 o/ y1 U9 }6 ]這個是經(jīng)典的簡支梁的分析,，移動載荷,，但我做出來的結(jié)果不是怎么的理想,。
FINI
% n8 j  k$ f9 j& S3 Y* z3 C/CLE
!設(shè)置參數(shù)
/CONFIG,NRES,2000
/ j$ l. Z, e4 w5 {5 H: X*SET,LONG,40    !梁長度
*SET,NE,20    !單元數(shù)
*SET,DL,LONG/NE    !單元長度
7 W4 F5 T! z" ^5 p# {3 j, M*SET,NN,NE+1    !節(jié)點數(shù)
*SET,FF,5000    !移動的集中力
& g, G( ?% r! D: c# t*SET,VV,80    !移動速度
*SET,DT,DL/VV*3.6   !移動經(jīng)過一個單元所需時間
: C1 o$ X3 Y, Y9 c0 H8 K1 a!定義單元材料實常數(shù)
: y2 H* {5 b. I3 v) t! _/PREP7
+ ]8 w8 v9 p$ H5 M2 VET,1,BEAM3    !梁單元
; M) i3 R5 g7 r, JMP,EX,1,50E9    !彈性模量
4 \- Y8 v( S; v. jMP,DENS,1,2500    !密度
MP,NUXY,1,0.2
. N+ F8 l& Y$ w0 F, ~8 G4 P5 SR,1,0.5,0.041667,1   !幾何常數(shù)
, o- N5 V7 |* B+ b' x!建立模型
*DO,I,1,NN,1
N,I,(I-1)*DL    !建立
*ENDDO
2 i6 U! @* Z9 t9 Y  Z( ~/ `: X*DO,I,1,NE
E,I,I+1
*ENDDO
2 k1 Y2 r0 @% U1 @" qFINI
8 R) |9 V' e- G2 R!加載與求解
/SOL
D,1,ALL
D,NN,ALL
4 S& q7 k$ s, D* IANTYPE,4    !瞬態(tài)分析
TRNOPT,FULL
* E2 b1 k: H  b% O; jTIMINT,OFF
0 a/ s9 u) `* R6 ^3 f, `OUTRES,BASIC,ALL    !輸出基本項，每一步都輸出
  T' C& {0 O  C- m  B4 j* ~# VKBC,1       !指定載荷類型
. c) q9 s+ S  k1 ]TIME,1E-10     !指定載荷步結(jié)束時間
, k- u" \( `5 s0 \NSUB,4,8      !指定載荷子步數(shù)為4,，最大子步數(shù)為8
0 F  o% Y" C8 wSSTIF,ON      !打開應(yīng)力剛化效應(yīng)
ACEL,0,9.8,0     !施加重力加速度
0 A+ ^! w; h! Q. Z) p5 J# Y# yALLS       !選中所有元素
SOLVE       !求解第一個載荷步
TIMINT,ON
*DO,II,1,NN,1
) K1 v& z0 q8 [6 H( m6 Y* COUTRES,ALL,ALL    !輸出所有加載子步的所有結(jié)果
, d5 M( I* O! _/ x* ]- j) [9 \TIME,II*DT    !每個荷載步的求解時間
NSUBST,5    !每個荷載步包含五個子步
FDELE,ALL,ALL    !刪除所有的集中力荷載
NSEL,S,LOC,X,(II-1)*DL !選擇集中力作用節(jié)點
2 b  a9 T% o" ]. tF,ALL,FY,-FF    !施加集中荷載
4 v' l3 h) a  }. z/ VALLSEL    !選擇所有
SOLVE
*ENDDO
SAVE
- U8 N8 d- Q/ @( [0 AFINI
, ]4 ]$ d4 y( J& p7 h3 U# x' I* o/POST1
3 W4 D0 M% k) Y; ASET,NE/2+1,LAST    !設(shè)定讀取的載荷步
/ESHAPE,1    !顯示梁形狀
PLNSOL,U,Y,0,1.0   !Y向位移
3 w: Y' ^, a1 dPLNSOL,S,EQV,0,1.0   !等效應(yīng)力
FINI
; O6 K) P9 z5 f7 L/POST26
NSOL,2,NE/2+1,U,Y,UY_M       !取出梁中部節(jié)點Y方向的位移
XVAR,1                   !定義時間變量為坐標(biāo)橫軸
PLVAR,2                  !繪制位移隨時間變化的圖形
DERIV,3,2,1,,,,,1        !對梁中部節(jié)點的Y向位移進(jìn)行微分運(yùn)算，得到Y(jié)向速度
PLVAR,3                  !繪制速度隨時間變化的圖形
DERIV,4,3,1,,,,,1        !對梁中部節(jié)點的Y向速度進(jìn)行微分運(yùn)算,，得到Y(jié)向加速度
9 {+ v  o% B* F, gPLVAR,4                  !繪制加速度隨時間變化的圖形

小小實習(xí)生

為啥沒有高手回應(yīng)呢,？？,？,？

eddyzhang

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

卡夏撒

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

李曉輝

高人呀

samyfish

高手