振动函数（正余弦）的实现(页 1) - D03：PFC2D/PFC3D - SimWe仿真论坛 |Simwe.com|仿真|设计|有限元|虚拟仪器

zj6962369 发表于 2008-4-19 20:11

振动函数（正余弦）的实现

现在做振动平面的问题：平面位置按照正余弦的规律变化，粒子在平面上随之振动，设定频率与振幅实现不同的振动规律。%T2b*eop%eXv
以下程序调试不成功：望高手指点一二：
qq：108739875
e-mail: [email=zj6962369@163.com]zj6962369@163.com[/email](T\SQ)Y sp-f2FC
;---------------------shake------------------------------------------
DEF startup
   old_time = time 5T ^_$Qt+T&i
   omega = freq * 2.0 * pi ,I;C,m(M;yu
END
DEF shake (q9x^.Zmw%_&F-c(\
   while_stepping H_N{ D\*W
   real_time = time - old_time
   if real_time < tlength
      vel = ampl * sin(omega*real_time) !B:smy+^ `f
   else +c(k @+uY5J
      vel = 0.0
   end_if 1Huxa{$uK3M
END w,qm_%Ts
;-----------set shake parameters------------------------------------------------------
DEF set_shake%j^UFo-c
freq=25
ampl=0.003
tlength=1.0 'JTY#tp"gr;Y
beta=72
x_vel=vel*0.309  ;cos(beta)=0.309
z_vel=vel*0.9511  ;sin(beta)=0.9511
command
wall  id  1  xv=x_vel yv=0.0 zv=z_vel xs=0.0 ys=0.0 zs=0.0
endcommand
END^;`P6dnz^jNE
;-----------------------------------------------------------------------------1E1G`{.R+m[tY)^
;---------------------运行1s准备计量-------------------------------------
  hist nstep 10,R`-` @q2h5G t6{
  set dt 1e-6
  startup4S(Z#ib0Z+O
  set_shake N-QjSO0M3n/s
  set fishcall(0) shakei y@Mv;_a
  cycle 1000000

zj6962369 发表于 2008-4-23 09:01

:( 没人理，自己先顶一下吧……

hbt2008 发表于 2008-4-24 11:43

兄弟我加你了进群了支持你

luzheng111 发表于 2008-4-27 07:01

楼主，你对wall的赋值，没有进入loop的循环，应该是不能运行的。你把赋值放在loop循环里面，估计就可以了。参看shake那个例子。

zj6962369 发表于 2008-4-28 14:30

感谢各位！我试一下……

zj6962369 发表于 2008-4-28 16:27

程序没问题可以运行，只是参数设置不当不容易观察到，麻烦各位了，我对其中的参数设置，参数的单位还不太熟悉，熟悉得说一下吧，频率，振幅等的默认单位是多少啊？

luzheng111 发表于 2008-4-29 14:26

你其他的单位用的是什么，你的频率，振幅的单位就是什么，各个量的单位保持一致就可以。对了，set fishcall(0) shake这个怎么理解呢？fishcall我一直不太懂咯。

zj6962369 发表于 2008-5-1 10:12

不好意思，我也不太懂……试验出来的，可以实时调用后面的函数，:victory:

页: [1]

SimWe仿真论坛's Archiver

振动函数（正余弦）的实现