【求助】加载正玄函数-节点响应成直线？

miranda

我想在悬臂梁的自由端施加载荷：P(t)=200*sin(628*t+30)，搜了一下论坛，我试了一下大家建议的方法，可是总是出错了啦。
（1）函数加载。步骤如下：
a) parameters-->function-->define/edit  
    定义了一正玄函数 200*sin(628*t+30),并命名为“fsint”;
b) parameters-->function-->read from file
    读取文件“fsint.func”,并定义“TABLE”名为“fsint”;
c) 求解。
  加载时给自由端横截面沿Y方向（截面很小）加载，选择“exsiting table”；
  求解。
d) 时间历程后处理查看。
查看节点随时间位移时是直线，怎么回事嘛？请大侠帮忙哦


（2）多步载荷求解。命令流如下：
*afun,deg
/solu
anty,trans
trnopt,full

*do,j,1,100
dt=0.005
f,12,fy,200*sin(628*dt*j)
time,dt*j
nsubst,1,,,1
kbc,0
autots,on
outres,all,all
lswrite
*enddo
lssolve,1,100,1

save
finish

显示节点位移时，正玄曲线不光滑？哪位大侠帮忙指点一下哦

16443

增加载荷步

miranda

首先谢谢楼上的指点:handshake

我又搜了一下了，重新试了table施加简谐载荷，p(t)=200*sin(628*t),施加位置为自由端的12节点处。
取自由端的16节点随时间位移响应，我觉得稳定状态该是周期性的正弦（余弦）曲线才对嘛，可仿真显示不是的。请大虾们帮忙指点一下了。
整体命令流如下：

/clear
/filname, ndisp
/title, nodal displacement of straigh cantilever beam
/prep7
！建立模型
ET,1,SOLiD45
  MP,EX,1,2.1e11
  MP,PRXY,1,0.25
  MP,DENS,1,7800
  blc4,,,0.06,0.02,1      
  lesize,11,0.1
  lesize,6,0.02
  lesize,7,0.03
  vmesh,all
  allsel

save
  finish

！模态分析
/solu
  antype,modal
  modopt,subsp,8，，，，on
  lumpm,1                        !集中质量矩阵
  nsel,s,loc,z,0
  d,all,all
  allsel
  nsel,all
  allsel
  save
  solve
  finish
/solu
  expass,on                 !激活扩展模块
  mxpand,3                 !指定模块展开的阶数
  save
  solve
  finish
                                        ！自由右端施加简谐载荷，并测节点位移响应
NT=68
DT=0.0005

*dim,p,array,NT

  /INPUT,modal,txt   
  /SOLU
  ANTYPE,TRANS               
  TRNOPT,FULL
                                             !ALPHAD,0.1           
                                             !BETAD,0.03            
  KBC,1                 
  NSUBST,1, , ,1         
  OUTRES,ALL,ALL         
  *do,i,1,NT
  TIME,i*DT
  AUTOTS,on               
  f,12,fy,p(i)

  LSWRITE, i
  *enddo
  LSSOLVE,1,NT,1
  
  SAVE
  finish
  
  /post26
  nsol,2,16,u,y,16uy
  plvar,2
  
  save
  finish

miranda

怎么没人响应嘛？

terentao

200*sin(628*t+30)  改为 200*sin(628*Time+30),

terentao

你的时间取得太长了，你的振动周期为1 /628.
参考下图

miranda

谢谢，指点
我按大虾讲的做了如下改动：
NT=33
DT=0.0005
*dim,p,array,NT
/INPUT,ndisp1,txt            !表名 ndisp1.txt

图22 是自由端中间节点施加正弦载荷的时间历程，
图23怎么还是不对呐，:@L

kathy810

如果定义的函数载荷，是正弦变化的加速度，怎么加载啊？

miranda

我想，这个图形的平衡线是否是倾斜的？
如果是的话，如何调整曲线平衡点为横坐标？？？