求助：悬臂梁移动载荷分析。

小小实习生

属性：长度为20mm，直径为1mm，材料为45号钢，杨氏模量为2.1×1011pa，泊松比为0.3，密度为7.8×103kg/m³
2 g  o) ^- C6 \. Z( u力：大小1.37N，移动速度为：3mm/s.从一端移动到另一端
+ d6 q4 }6 G) i, J2 W2 ?6 }$ V; x求梁中点处的位移相应
8 ~8 j# x- D1 b$ S

小小实习生

$ C0 w1 C2 O$ T3 k1 F7 r% i9 W- P5 M
这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。
FINI
/CLE
!设置参数
$ @  x# Y0 z9 d/CONFIG,NRES,2000
*SET,LONG,40    !梁长度
*SET,NE,20    !单元数
3 c& a+ [  r/ g*SET,DL,LONG/NE    !单元长度
" v/ d" T6 `& Q* x*SET,NN,NE+1    !节点数
  I; U9 o6 O8 q- N5 p5 C*SET,FF,5000    !移动的集中力
*SET,VV,80    !移动速度
- P; h4 X& G# e3 V*SET,DT,DL/VV*3.6   !移动经过一个单元所需时间
!定义单元材料实常数
6 q3 T  a2 t) c% ]  k* z- Y; G" N/PREP7
ET,1,BEAM3    !梁单元
MP,EX,1,50E9    !弹性模量
! i8 J% M; _: aMP,DENS,1,2500    !密度
5 d2 _) z. {# Z$ xMP,NUXY,1,0.2
& t( b! [; T4 pR,1,0.5,0.041667,1   !几何常数
!建立模型
*DO,I,1,NN,1
N,I,(I-1)*DL    !建立
*ENDDO
*DO,I,1,NE
E,I,I+1
6 j# D9 J  f# M! @" X9 _*ENDDO
FINI
/ P+ |% ~" \  y. Z% @2 K- k!加载与求解
/SOL
D,1,ALL
: x  m9 ], }1 T+ pD,NN,ALL
ANTYPE,4    !瞬态分析
) J# z; i3 t  n0 R) ^- p8 f! B1 |$ rTRNOPT,FULL
TIMINT,OFF
0 ?. k; `7 K; X4 k  j+ M; xOUTRES,BASIC,ALL    !输出基本项，每一步都输出
KBC,1       !指定载荷类型
* N" `2 R% z+ a  NTIME,1E-10     !指定载荷步结束时间
NSUB,4,8      !指定载荷子步数为4，最大子步数为8
SSTIF,ON      !打开应力刚化效应
ACEL,0,9.8,0     !施加重力加速度
0 B# C" O, n' V% k- a5 ~0 S& TALLS       !选中所有元素
SOLVE       !求解第一个载荷步
- C* H" u* ~" a8 o! kTIMINT,ON
% I# a! r/ I; D# j6 H  U; \*DO,II,1,NN,1
" s( k$ ?) h' O. G8 POUTRES,ALL,ALL    !输出所有加载子步的所有结果
+ w$ D/ [* a: ]; y# L8 }TIME,II*DT    !每个荷载步的求解时间
* G; [) u: ^7 WNSUBST,5    !每个荷载步包含五个子步
( p8 n% t. ^5 i# T! `FDELE,ALL,ALL    !删除所有的集中力荷载
NSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点
F,ALL,FY,-FF    !施加集中荷载
7 N0 v+ e2 L; L; }+ O* G6 PALLSEL    !选择所有
+ A* m$ e/ K$ q% X0 u6 lSOLVE
" x+ S0 F! }0 R5 Z*ENDDO
. O- o. l* N: z) H/ z( wSAVE
FINI
/POST1
: g' x+ w1 Z9 q- O# p; vSET,NE/2+1,LAST    !设定读取的载荷步
/ESHAPE,1    !显示梁形状
PLNSOL,U,Y,0,1.0   !Y向位移
1 v" C7 ~3 K, u- V4 zPLNSOL,S,EQV,0,1.0   !等效应力
FINI
5 C, n1 C* V. Q6 d% x) t/POST26
5 V8 X# \, R) v( v. Z2 B# WNSOL,2,NE/2+1,U,Y,UY_M       !取出梁中部节点Y方向的位移
, A6 |' U0 K$ VXVAR,1                   !定义时间变量为坐标横轴
0 S5 d  ]5 j( ~8 x2 ~8 V5 N4 l* XPLVAR,2                  !绘制位移随时间变化的图形
DERIV,3,2,1,,,,,1        !对梁中部节点的Y向位移进行微分运算，得到Y向速度
. h& c) I9 B6 \6 aPLVAR,3                  !绘制速度随时间变化的图形
DERIV,4,3,1,,,,,1        !对梁中部节点的Y向速度进行微分运算，得到Y向加速度
PLVAR,4                  !绘制加速度随时间变化的图形
8 g! Q# w' m' T' i5 S
4 `) z! Q5 H% z" Y! r' U) T0 u

小小实习生

为啥没有高手回应呢？？？？

eddyzhang

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

卡夏撒

你就是高手啊，自己解决吧，我不会

李晓辉

高人呀

samyfish

高手