想不通啊，求助各位前辈，关于这个自由落体的小程序

72骑士

new
ball id=1 rad=1 x=5 y=10
prop kn=1e100 ks=1e100 dens=6000
set grav 0 -9.8
his ball yvel id 1
his ball ypos id 1
cyc 1
set dt 1e-3
cyc 2999
pl his 1 2
看fish ref的那个一个球下落变2个球的时候自己突发奇想，编了这么一个小程序，应该模拟小球自由落体的程序，追踪记录小球y方向速度和y方向位移。为了避免刚度造成的时步差异，我规定了1time-step为实际运行的1^(-3)秒。以下为刚度分别为1e6和1e100的图像。


1.发现刚度和密度均会影响速度以及位移曲线，刚度越小，密度越大，3秒（3000time-step）内下降高度越多，初速度越大，3s末速度越大。这是为什么？

2.设置刚度为1e100（可看作无穷大），能不能近似看作自由落体，程序规定的g=9.8m/平方s，可是从图像可以看出，3s内物体下降13.226m，速度达到8.84m/s，如此逆推出来，g=2.94m/平方s，这是怎么回事？

如果说的不清楚还可以再问我，谁能帮帮我解决一下心中的疑惑，万分感谢。

ccb137731

这个主要和参数选择有关。一般刚度系数的量级在1e5到1e10左右（具体数值记得不太清楚），得出的结果才会相对合理。这在魏群写的离散元一书中有介绍。你的问题主要原因是set dt 1e-3。颗粒接触只有在允许的时步内，才会得出合理的结果。你这样设置后，导致时步太大，因而造成这些现象。
像这些问题，没必要进行悬赏，我想只要知道的，都会帮你的。

ccb137731

这个主要和参数选择有关。一般刚度系数的量级在1e5到1e10左右（具体数值记得不太清楚），得出的结果才会相对合理。这在魏群写的离散元一书中有介绍。你的问题主要原因是set dt 1e-3。颗粒接触只有在允许的时步内，才会得出合理的结果。你这样设置后，导致时步太大，因而造成这些现象。

72骑士

首先谢谢版主的回答，可是感觉问题仍然没有解决，照版主提示做了以下修改。
new
ball id=1 rad=1 x=5 y=10
prop kn=1e8 ks=1e8 dens=6000
set grav 0 -9.8
fix x spin range id=1
his ball yvel id 1
his ball ypos id 1
cyc 300
pl his 1 2
请问这个程序是自由落体的模拟么？为什么结果仍然错误。

如图所示，物体下落实际时间2.866s（time），下落位移-2.535-9.963=-12.498m，物体速度变化：-8.571-（-0.4599）=-8.1111m/s
问题一：规定求初始位置10m ，为什么图上变成9.963m，而且具有初速度-0.4599m/s

问题二，如果由运动学公式忽略初速度的情况下：s=1/2gt^2，得出s=0.5*9.8*2.866^2=40.24m,而图上显示却只有12.498m的位移？

72骑士

该问题问过大神已经完全解决，阻尼问题，设置prop damp 0即可得到理论值，调整set hist_rep 1即可还原初速度0和初位移0，存在初位移的原因是一个cyc内循环历史记录过多，导致得到的点较离散

ccb137731

72骑士 发表于 2011-12-1 20:30
该问题问过大神已经完全解决，阻尼问题，设置prop damp 0即可得到理论值，调整set hist_rep 1即可还原初速 ...

呵呵，解决了就好！