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

小小实习生

属性：长度为20mm，直径为1mm，材料为45号钢，杨氏模量为2.1×1011pa，泊松比为0.3，密度为7.8×103kg/m³
力：大小1.37N，移动速度为：3mm/s.从一端移动到另一端
求梁中点处的位移相应
( m7 i: a/ ?! L0 d

小小实习生

& i* i+ `2 d$ F- {' P
% n; Y! f' l3 r! V& ?! `# }这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。
FINI
8 b2 |/ `2 N3 b* f. }6 p% a/CLE
!设置参数
/CONFIG,NRES,2000
+ e- [" {+ C; ?& o*SET,LONG,40    !梁长度
! X) F( n0 G; k' w*SET,NE,20    !单元数
*SET,DL,LONG/NE    !单元长度
4 F1 T2 a8 g( @% x/ C*SET,NN,NE+1    !节点数
*SET,FF,5000    !移动的集中力
! \' ~& U! e& Q8 P" ~*SET,VV,80    !移动速度
*SET,DT,DL/VV*3.6   !移动经过一个单元所需时间
!定义单元材料实常数
+ @/ K( v( |; K/PREP7
' n4 O) O7 Z3 p# |4 s2 w- \$ pET,1,BEAM3    !梁单元
  U/ y) C' F# l2 O  \3 CMP,EX,1,50E9    !弹性模量
MP,DENS,1,2500    !密度
1 h+ o- C& i# p9 Z5 P7 U9 h! dMP,NUXY,1,0.2
R,1,0.5,0.041667,1   !几何常数
# ^: S5 s6 Y" ~& Z!建立模型
*DO,I,1,NN,1
- k# @6 P+ O$ mN,I,(I-1)*DL    !建立
) p0 t* ?& i8 U, P3 F7 V" x*ENDDO
1 i1 K' }! ^1 k7 b2 a' `*DO,I,1,NE
$ x4 C; V7 O* e2 ]. _E,I,I+1
*ENDDO
' q0 q  N- D# B  v* dFINI
( r4 f/ C/ O" c% n6 J!加载与求解
/SOL
& N" _1 b/ d6 L( PD,1,ALL
4 T% C  R4 G7 d: N  dD,NN,ALL
: [6 M7 \1 d6 f$ HANTYPE,4    !瞬态分析
/ _& |. |- `, vTRNOPT,FULL
TIMINT,OFF
OUTRES,BASIC,ALL    !输出基本项，每一步都输出
KBC,1       !指定载荷类型
TIME,1E-10     !指定载荷步结束时间
  x) W9 T/ S# E" j# E) f, H* PNSUB,4,8      !指定载荷子步数为4，最大子步数为8
SSTIF,ON      !打开应力刚化效应
ACEL,0,9.8,0     !施加重力加速度
ALLS       !选中所有元素
# P, g2 |9 ^, D& i1 z/ xSOLVE       !求解第一个载荷步
3 s7 j+ K6 B1 R! J; h' uTIMINT,ON
*DO,II,1,NN,1
. X1 I$ k: w' ?4 k8 bOUTRES,ALL,ALL    !输出所有加载子步的所有结果
TIME,II*DT    !每个荷载步的求解时间
  K. H3 R  l: oNSUBST,5    !每个荷载步包含五个子步
FDELE,ALL,ALL    !删除所有的集中力荷载
NSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点
F,ALL,FY,-FF    !施加集中荷载
1 p4 R' U/ P% o* Z6 GALLSEL    !选择所有
! v  [* S) f; F# FSOLVE
*ENDDO
( |+ ^7 W1 C9 e! y. E$ ]SAVE
" L/ F' K& T: Y$ P) }0 nFINI
% \) \( i; v" [" |6 B/POST1
SET,NE/2+1,LAST    !设定读取的载荷步
/ESHAPE,1    !显示梁形状
PLNSOL,U,Y,0,1.0   !Y向位移
  n' {& n) O% V- VPLNSOL,S,EQV,0,1.0   !等效应力
FINI
7 \+ o7 E# h/ r, s/POST26
NSOL,2,NE/2+1,U,Y,UY_M       !取出梁中部节点Y方向的位移
$ F" g. K0 }, V) Y* NXVAR,1                   !定义时间变量为坐标横轴
( L9 K9 T6 m* a, @: l5 x! uPLVAR,2                  !绘制位移随时间变化的图形
DERIV,3,2,1,,,,,1        !对梁中部节点的Y向位移进行微分运算，得到Y向速度
4 i7 D: _) ]' Z0 P9 h- t3 V& ]0 w* N. WPLVAR,3                  !绘制速度随时间变化的图形
DERIV,4,3,1,,,,,1        !对梁中部节点的Y向速度进行微分运算，得到Y向加速度
# w. C6 E. R# k. U, b% HPLVAR,4                  !绘制加速度随时间变化的图形
2 t6 {$ I9 L, }

小小实习生

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

eddyzhang

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

卡夏撒

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

李晓辉

高人呀

samyfish

高手