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

小小实习生

属性：长度为20mm，直径为1mm，材料为45号钢，杨氏模量为2.1×1011pa，泊松比为0.3，密度为7.8×103kg/m³
6 R1 `3 U# s: F$ V2 H2 s力：大小1.37N，移动速度为：3mm/s.从一端移动到另一端
求梁中点处的位移相应

小小实习生

, U9 C. j0 c7 Y- G& t) Q/ Y
这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。
" k$ B; \" F  v5 y& S5 B% F; e# ?" GFINI
4 P5 r+ M* i$ @* C$ E0 j/CLE
7 n- L( F6 s& t  b8 N; F+ J!设置参数
: w) N* C: f- Y/CONFIG,NRES,2000
  e% T+ o% p8 e' r/ u0 y*SET,LONG,40    !梁长度
+ {, C" d/ n6 f$ H; X1 m*SET,NE,20    !单元数
*SET,DL,LONG/NE    !单元长度
9 U# {: O6 N0 W5 f( R* s, g*SET,NN,NE+1    !节点数
9 i* ?) u& P/ G2 @0 r" ~*SET,FF,5000    !移动的集中力
+ L$ W7 v* W& h+ F- @6 W*SET,VV,80    !移动速度
*SET,DT,DL/VV*3.6   !移动经过一个单元所需时间
!定义单元材料实常数
) ?( E. P5 o+ p/PREP7
ET,1,BEAM3    !梁单元
MP,EX,1,50E9    !弹性模量
MP,DENS,1,2500    !密度
) h( i: D( x, X+ f" w9 hMP,NUXY,1,0.2
R,1,0.5,0.041667,1   !几何常数
!建立模型
*DO,I,1,NN,1
9 l3 o% `( B$ ^N,I,(I-1)*DL    !建立
*ENDDO
- U9 q' _, F2 e' h# R*DO,I,1,NE
E,I,I+1
*ENDDO
1 X* _4 U4 q0 f+ s7 t9 Q3 sFINI
!加载与求解
& O, g- b/ R- o9 W/SOL
D,1,ALL
7 K2 ]) J% F2 P' R! D# H, L4 FD,NN,ALL
2 I+ U4 x2 e/ ^: b# ^6 PANTYPE,4    !瞬态分析
TRNOPT,FULL
& o/ B+ r% T3 o, [* J# I" k# cTIMINT,OFF
OUTRES,BASIC,ALL    !输出基本项，每一步都输出
! N" a' k1 H5 O# R$ l: k0 mKBC,1       !指定载荷类型
6 W7 o) X7 n; @4 `7 FTIME,1E-10     !指定载荷步结束时间
( a+ E9 K: q& l$ V. V( fNSUB,4,8      !指定载荷子步数为4，最大子步数为8
) D0 L/ ~0 F% n/ k4 d, ?0 FSSTIF,ON      !打开应力刚化效应
ACEL,0,9.8,0     !施加重力加速度
. h: d; i+ Q) C: LALLS       !选中所有元素
SOLVE       !求解第一个载荷步
TIMINT,ON
! N7 U1 c1 A- _) `# d) r8 L7 P$ I*DO,II,1,NN,1
OUTRES,ALL,ALL    !输出所有加载子步的所有结果
TIME,II*DT    !每个荷载步的求解时间
NSUBST,5    !每个荷载步包含五个子步
FDELE,ALL,ALL    !删除所有的集中力荷载
NSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点
F,ALL,FY,-FF    !施加集中荷载
ALLSEL    !选择所有
. ]# ?7 M  B9 L; [$ w& kSOLVE
*ENDDO
SAVE
$ ], `* h( o# D; b% C5 d8 ?FINI
$ j; i# ^6 u1 ~6 c/POST1
SET,NE/2+1,LAST    !设定读取的载荷步
$ b5 f4 R6 t; u- h  [/ D9 B+ X/ESHAPE,1    !显示梁形状
PLNSOL,U,Y,0,1.0   !Y向位移
PLNSOL,S,EQV,0,1.0   !等效应力
FINI
/POST26
6 D  C+ J: R' I1 t2 ANSOL,2,NE/2+1,U,Y,UY_M       !取出梁中部节点Y方向的位移
- y, r, `* a) A9 W* |  r5 ]5 p# f; sXVAR,1                   !定义时间变量为坐标横轴
PLVAR,2                  !绘制位移随时间变化的图形
: Q; _/ k+ W8 @6 s* BDERIV,3,2,1,,,,,1        !对梁中部节点的Y向位移进行微分运算，得到Y向速度
1 P# b/ ]' G: p3 b1 W0 T% q' M/ p4 ~PLVAR,3                  !绘制速度随时间变化的图形
DERIV,4,3,1,,,,,1        !对梁中部节点的Y向速度进行微分运算，得到Y向加速度
PLVAR,4                  !绘制加速度随时间变化的图形

小小实习生

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

eddyzhang

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

卡夏撒

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

李晓辉

高人呀

samyfish

高手