3维空间的物理运动

/**
     * 小球的物理运动
     */
    function SpherePhysicalRun() {


        //小球属性
        var u = 0.25, /*物体摩擦系数*/
            t = 2.0, /*时间 单位second*/
            g = -9.8;
        var w = 50 * lbstoN, /*物体重力*/
            n = 400, /*推力 单位'牛'*/
            n_agl = 30 * degTorad, /*推物体的角度 方便计算'角->弧'*/
            k_agl = 210 * degTorad;


        //计算滑动摩擦力 计算时 反作用力都以负值
        var fk = -u * w;
        print("fk=" + fk);


        //计算合力
        //拆分分量形式计算
        var fpVec = new SFVec3f(n * cos(n_agl), n * sin(n_agl), 0);
        print("fpVec=" + fpVec);
        var fkVec = new SFVec3f(-fk * cos(k_agl), -fk * sin(k_agl), 0);
        print("fkVec=" + fkVec);
        var fnetVec = fpVec.add(fkVec);
        print("FnetVec=" + fnetVec);
        //计算物体质量 w=mg
        var m = -w / g;


        //计算加速度 fnet=ma
        var aVec = fnetVec.multiply(1 / m);
        print("aVec=" + aVec);


        var oTrans = _person.translation;
        oTrans.setValue(-3, 0, 0);
        var curTrans = oTrans.clone();




        var timerId = timer.bindFraction(t, function (f) {
            var tt = f * t;
            //计算物体位移 方程 R=Vi*t+1/2*a*(t*t)
            var disVec = aVec.multiply(0.5 * tt * tt);
            f == 1 && print("disVec==" + disVec);
            oTrans.setValue(curTrans.add(disVec));
            _person.rotation.setValue(0, 0, 1, 9 * f * 360 * degTorad);
        }, function () {
//            oTrans.setValue(curTrans);
        });


    }