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

小小实习生

属性：长度为20mm，直径为1mm，材料为45号钢，杨氏模量为2.1×1011pa，泊松比为0.3，密度为7.8×103kg/m³
力：大小1.37N，移动速度为：3mm/s.从一端移动到另一端
2 c. p2 b$ v4 U, A# q- u! o7 C, M求梁中点处的位移相应

小小实习生

) m6 h& ]3 ?" @- t
这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。
FINI
4 K" a, ^4 B5 x% [  h7 {" Y/CLE
+ X& L3 E) B- K8 `8 y!设置参数
) Y( z7 F$ }/ s" o8 ]+ Z# z, |7 {/CONFIG,NRES,2000
# u; t7 r- f% Y*SET,LONG,40    !梁长度
% P8 m( @$ H9 l  ~*SET,NE,20    !单元数
*SET,DL,LONG/NE    !单元长度
, E$ n0 a+ ^# [- y* t1 Y*SET,NN,NE+1    !节点数
" ]8 H0 |7 i" H% d; g, V*SET,FF,5000    !移动的集中力
*SET,VV,80    !移动速度
*SET,DT,DL/VV*3.6   !移动经过一个单元所需时间
!定义单元材料实常数
/PREP7
ET,1,BEAM3    !梁单元
% x3 n! {4 M% v( h. r# [% VMP,EX,1,50E9    !弹性模量
' k; E0 W* c$ a* P: ]' O- HMP,DENS,1,2500    !密度
7 Q5 Q- L, \3 t& B* JMP,NUXY,1,0.2
R,1,0.5,0.041667,1   !几何常数
$ r* g! _+ J& h& B' t( @!建立模型
2 t7 i# ?7 o* T- V+ U: Z4 O*DO,I,1,NN,1
N,I,(I-1)*DL    !建立
9 y' C5 q  i$ r$ }2 n  J& _*ENDDO
5 i3 ^9 D6 B/ h# z8 R*DO,I,1,NE
9 r, P- _# U: ~% ME,I,I+1
*ENDDO
0 D2 C. F: C- F" wFINI
!加载与求解
/SOL
4 F- U3 c  {6 B& HD,1,ALL
  d# e1 a+ U$ [D,NN,ALL
ANTYPE,4    !瞬态分析
TRNOPT,FULL
2 Q- H: D6 w; I# p  UTIMINT,OFF
OUTRES,BASIC,ALL    !输出基本项，每一步都输出
KBC,1       !指定载荷类型
TIME,1E-10     !指定载荷步结束时间
6 ~/ _( U8 P3 q$ A+ CNSUB,4,8      !指定载荷子步数为4，最大子步数为8
+ x. D  g0 l8 y  i' I# n9 qSSTIF,ON      !打开应力刚化效应
ACEL,0,9.8,0     !施加重力加速度
ALLS       !选中所有元素
! O9 ?; Y+ S+ YSOLVE       !求解第一个载荷步
+ a4 |, u# Q: L7 v5 n1 w/ d( dTIMINT,ON
( [% s% s( ], U" |*DO,II,1,NN,1
OUTRES,ALL,ALL    !输出所有加载子步的所有结果
6 a8 A% z# q/ J8 ATIME,II*DT    !每个荷载步的求解时间
! {) @# G7 k7 }# _2 Q  {NSUBST,5    !每个荷载步包含五个子步
0 \3 {$ P0 e' j7 N. j9 C' ?/ SFDELE,ALL,ALL    !删除所有的集中力荷载
6 O1 h% [; R! C- R, L# vNSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点
F,ALL,FY,-FF    !施加集中荷载
9 D7 P& z  D3 K# ~. ^2 UALLSEL    !选择所有
$ P: y- n' r; lSOLVE
*ENDDO
0 f$ V, y" O% U) A! x( F1 {% k' ZSAVE
2 u- ~$ F- ]1 l+ wFINI
+ R" H0 l* y; C2 P+ u/POST1
: h# ]) J" l9 A9 s8 @7 t+ H4 `. ySET,NE/2+1,LAST    !设定读取的载荷步
/ESHAPE,1    !显示梁形状
9 `% R3 I' `3 V7 b, UPLNSOL,U,Y,0,1.0   !Y向位移
5 `& d& h0 U6 l& x4 T! APLNSOL,S,EQV,0,1.0   !等效应力
  e" a) i% t! P* ?6 OFINI
# N3 A/ E2 t8 V# m& @( W/POST26
6 B3 x/ |0 N1 L  X2 \NSOL,2,NE/2+1,U,Y,UY_M       !取出梁中部节点Y方向的位移
XVAR,1                   !定义时间变量为坐标横轴
PLVAR,2                  !绘制位移随时间变化的图形
- [7 O- X0 s7 H0 _. uDERIV,3,2,1,,,,,1        !对梁中部节点的Y向位移进行微分运算，得到Y向速度
PLVAR,3                  !绘制速度随时间变化的图形
3 [+ H( u4 ^# G8 T1 m! wDERIV,4,3,1,,,,,1        !对梁中部节点的Y向速度进行微分运算，得到Y向加速度
4 r/ Y3 w2 u, H! g  N" B' m$ rPLVAR,4                  !绘制加速度随时间变化的图形

小小实习生

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

eddyzhang

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

卡夏撒

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

李晓辉

高人呀

samyfish

高手