用瞬态分析计算简谐载荷加载问题

xinxinxiangsun

之前搜了一下，发现这个问题貌似已经讨论过很多了，学习了一下，编了如下的一段命令流，但是求解后好像还是不行的。恳请各位前辈指点一二
*afun,deg
/solu                  
anty,trans
trnopt,full
d,52000,all,,,52310,1
k1=52502
k2=52402
k3=52501
k4=52401
dt=60/740/4
*do,j,1,20
fk,k1,fx,135864*cos(4440*dt*j+95.4+42.3)
fk,k1,fy,135864*sin(4440*dt*j+95.4+42.3)
fk,k2,fx,135864*cos(4440*dt*j+95.4+42.3)
fk,k2,fy,135864*sin(4440*dt*j+95.4+42.3)
fk,k3,fx,67680*cos(4440*dt*j+42.3)
fk,k3,fy,67680*sin(4440*dt*j+42.3)
fk,k4,fx,67680*cos(4440*dt*j+42.3)
fk,k4,fy,67680*sin(4440*dt*j+42.3)
time,dt*j
nsubst,1,,,1
kbc,0
outres,all,all
solve
*enddo

下面是我用这个命令流求解的一些结果，好像是有问题的，应该是五个循环才对。

[ 本帖最后由 xinxinxiangsun 于 2008-5-25 12:19 编辑 ]

16443

问题没有说清楚，你认为哪里不正确。
检查一下的后处理。

xinxinxiangsun

是这样的，我想用瞬态分析计算在简谐力作用下结构某点的位移——时间曲线。
尝试在加载了五个周期：dt=60/740/4，其中740是转速，转/分；dt是想把一个周期分为四个子部。do语句想计算5个周期。

得到的位移——时间结果应该是有五个周期啊，为什么就是一条直线呢？
如果不是载荷加载有问题，会不会是我后处理有错误呢？

这是在post26里列表显示的结果：
***** ANSYS POST26 VARIABLE LISTING *****


    TIME       2294 UX      
                  UX_2   
  0.20270E-01     0.00000   
  0.40541E-01     0.00000   
  0.60811E-01     0.00000   
  0.81081E-01     0.00000   
  0.10135         0.00000   
  0.12162         0.00000   
  0.14189         0.00000   
  0.16216         0.00000   
  0.18243         0.00000   
  0.20270         0.00000   
  0.22297         0.00000   
  0.24324         0.00000   
  0.26351         0.00000   
  0.28378         0.00000   
  0.30405         0.00000   
  0.32432         0.00000   
  0.34459         0.00000   
  0.36486         0.00000   
  0.38514         0.00000   
  0.40541         0.00000   
为什么都是0啊？

下面是我后处理部分的截图：我是计算完成后直接进入TimeHist Post,然后选中我要观察的节点。

[ 本帖最后由 xinxinxiangsun 于 2008-5-25 14:23 编辑 ]

16443

可能是你加载没有成功。看看err文件的信息。

xinxinxiangsun

上次可能是我加载有问题，这次我又改了一下，直接用节点号添加力；还加入了重力加速度和阻尼，下面是我的命令流：
*afun,deg
/solu

d,52000,all,,,52310,1                    
acel,0,9.8
alph,0.1
anty,trans
trnopt,full
dt=60/740/4
*do,j,1,20
f,52502,fx,135864*cos(4440*dt*j+95.4+42.3)
f,52502,fy,135864*sin(4440*dt*j+95.4+42.3)
f,52402,fx,135864*cos(4440*dt*j+95.4+42.3)
f,52402,fy,135864*sin(4440*dt*j+95.4+42.3)
f,52501,fx,67680*cos(4440*dt*j+42.3)
f,52501,fy,67680*sin(4440*dt*j+42.3)
f,52401,fx,67680*cos(4440*dt*j+42.3)
f,52401,fy,67680*sin(4440*dt*j+42.3)
time,dt*j
nsubst,1,,,1
kbc,0
autots,on
outres,all,all
lswrite
*enddo
lssolve,1,20,1

但是改过之后好像还是不行，得到的结果也不是正弦曲线：
会不会是我的重力或者阻尼加载的不对呢？

xinxinxiangsun

又重新改了一下命令流，去掉了重力加速度，将子部划分的更细
*afun,deg
/solu

d,52000,all,,,52310,1                    

alph,0.1
anty,trans
trnopt,full
dt=60/740/16
*do,j,1,160
f,52502,fx,135864*cos(4440*dt*j+95.4+42.3)
f,52502,fy,135864*sin(4440*dt*j+95.4+42.3)
f,52402,fx,135864*cos(4440*dt*j+95.4+42.3)
f,52402,fy,135864*sin(4440*dt*j+95.4+42.3)
f,52501,fx,67680*cos(4440*dt*j+42.3)
f,52501,fy,67680*sin(4440*dt*j+42.3)
f,52401,fx,67680*cos(4440*dt*j+42.3)
f,52401,fy,67680*sin(4440*dt*j+42.3)
time,dt*j
nsubst,1,,,1
kbc,0
autots,on
outres,all,all
lswrite
*enddo
lssolve,1,160,1

但是好像出来的结果还是有问题的,我是想吧一个正弦周期分成16份（dt）然后计算10个周期，但是得到的结果好像不是十个周期。当使dt=60/740/4，j是1-20，也就是把正弦周期分成4份，计算5个周期时结果就又变成5个周期了。

不知道那位大侠能帮忙解释一下


下面的图片是分16份，计算10个周期的结果截图。

16443

1、dt=60/740/16这个值不确定，带来新的数值误差。
2、考虑这个函数135864*cos(4440*dt*j+95.4+42.3)中个参数的正确性。

xinxinxiangsun

楼上的意思是说dt的设定不是一个整数，在计算的时候会产生累计误差是吧。
你的第二条意见我没太理解，是说这里的dt*j设定的有问题么？

xinxinxiangsun

我又改了一下，
dt=0.02
*do,j,1,50
f,52502,fx,135864*cos(4440*0.02*j+95.4+42.3)
…………
f,52401,fy,67680*sin(4440*0.02*j+42.3)
time,0.02*j
…………
lssolve,1,50,1
这回计算出来的时间范围是对的，但是得到的结果仍然不是想像中的那种正弦分布的。不知道是那里设的有问题哦。

ps:用这个命令流计算，耗费时间很长，是不是我的载荷步读写计算方法上面有问题呢？请大侠不吝赐教:@)

shr20604043

可以定义一个函数加载，这样比较准确。载荷步态多很费时间。

gdayang

我也正在做这方面的东西，想在流体中加一个振动源，是给定节点的位移变化载荷还是压力载荷，请指教！