三版例子用五版写 不知错在哪

笨花part2

设定为Hertz contact,一个球以3.9的向下速度，瞬间碰撞三角形平面，求碰撞的持续时间T (The duration of the collision T)

圆球的参数如下:
半径: 0.0025
位置: (0.05, 0.04, 0.0025)
密度: 4000
初速度: z轴向下3.9
Local damping coefficient: 0
Cirtical damping ratio in the normal direction: 0
Cirtical damping ratio in the shear direction: 0
Shear moudulus: 1.545e+11
Poission ratio: 0.23
Friction: 0.092

墙的参数如下:
三角形座标 (0.0, 0.0, 0.0) (0.1, 0.0, 0.0) (0.05, 0.1, 0.0)
Friction: 0.092

解压包有三版与五版的程式,结果比较,解析解算法
比较结果，三版是6.10E-06与解析解相同，但五版的是8.10E-06，想请问是哪个命令弄错了，谢


//////////////////////////////////////////////////////////////////////
三版程式如下
FC3D 三版跑的
new
title 'Normal_impact_restitution_coefficient_1.0'
;
define simulated_time
       whilestepping
       currenttime=time
end
;
wall id 1 face (0.0, 0.0, 0.0) (0.1, 0.0, 0.0) (0.05, 0.1, 0.0)
wall property friction 0.092
;
damp default local 0.0
damp default viscous normal 0.0
damp default viscous shear 0.0
;
ball id 1 hertz rad 0.0025 x 0.05 y 0.04 z 0.0025  
property shear 1.545e+11 poiss 0.23 friction 0.092 density 4000 zvelocity -3.9
;
set dt max 1.0e-8
set hist_rep 10
simulated_time
;
history id 1 currenttime
history id 2 ball zp id 1
history id 3 ball zv id 1
history id 4 ball xs id 1
history id 5 ball ys id 1
cycle 1000
history write 1 2 3 4 5 file Results_1.0.HIS
return
;

//////////////////////////////////////////////////////////////////////
五版程式如下
;PFC3D 五版跑的
new
title 'Normal_impact_restitution_coefficient_1.0'
;
domain extent -5.0 5.0
;
cmat default type ball-facet   ...
             model hertz ...
             property hz_shear 1.545e+11 hz_poiss 0.23 fric 0.092 dp_nratio 0.0 dp_sratio 0.0
wall create id 1 vertices (0.0, 0.0, 0.0) (0.1, 0.0, 0.0) (0.05, 0.1, 0.0)
ball create id 1 x 0.05 y 0.04 z 0.0025 radius 0.0025
ball attribute density 4000 damp 0.0 zvelocity -3.9
;
set timestep maximum 1.0e-8
set hist_rep 10
history add id 1 mechanical solve time
ball history id 2 zposition id 1
ball history id 3 zvelocity id 1
ball history id 4 xspin id 1
ball history id 5 yspin id 1
cycle 1000
history write 1 2 3 4 5 file Results_damp_0.0.HIS
return
;

图档

我的形象

ffffffdddddd