核爆冲击波效果 粒子系统

hjz000

本人想做一个核爆冲击波的效果，用粒子系统实现，各位大侠有没有好的建议，或例子供参考一下。

kodong122

好恐怖，怎么样子的

liuzhiyu123

osg现有的估计没有这个功能，考虑其他的方案或者成熟的粒子引擎

liushui

我们来模拟，随便想个简单的算法，控制粒子移动就行了，

liushui

我们来模拟，随便想个简单的算法，控制粒子移动就行了，

liushui

不专业的模拟，随便想想就能数学建模吧

hjz000

不专业也行啊，大概要个样子就行，该如何建模呢

array

核爆炸的数学模型肯定是需要您自己去建立的，然后用osgParticle或者Spark这样的库去实现都是可以的

hjz000

d:\1.jpg

hjz000

hjz000

代码如下，但是问题是，好像所有的粒子的运行轨迹一致导致，出现一个个很大的亮点，并没有预想的分散开形成预想的效果，问题到底在哪如何修改呢？项目难点，下周四要完成，各位大侠帮忙瞧一瞧，给个建议。
// apply a vortex-like acceleration. This code is not optimized,
    // it's here only for demonstration purposes.

    void operate(osgParticle:article *P, double dt)
    {
                static osgParticle::Particle *gP = P;
                double existTime = dt;

                //平面角
                static double circulTheta=360;

                //俯仰角
                static double circulPitch = 0;
                //球半径a
                static double circurRadius=30;

                //粒子运行1.5秒后进行球体轨迹运动
                if (P->getAge()>1.5)
                {
                        circulTheta -=0.0045;
                        circulPitch -=0.0045;
                        if (circulTheta>90)
                        {
                                circulTheta = 0;
                        }
                        if (circulPitch>360)
                        {
                                circulPitch = 0;
                        }
                        circulPitch = P->getCurrentAlpha();
                        circulPitch = P->ge;
                        osg::Vec3 tempPos = P->getPosition();               
                       
                        //计算粒子轨迹
                        y = circurRadius*cos((360-circulTheta)*osg::PI/180);
                        z = -(circurRadius*sin((360-circulTheta)*osg::PI/180));
                        x = circurRadius*sin(circulPitch*osg::PI/180);       
               
                        osg::Vec3 newpos;// = P->getPosition() + v * dt;
                        newpos[0] = x;
                        newpos[1] = y;
                        newpos[2] = z;


                        P->setPosition(newpos);  
                }               
        
       
                static double theta=360;

                //俯仰角
                static double pitch = 360;
                //双钮线最大半径a
                static double maxRadius=100;
                if (P->getAge()>5)
                {
                       
                        theta -=0.0045;
                        pitch -=0.0045;
                        if (theta<270)
                        {
                                theta = 360;
                        }
                        if (pitch<0)
                        {
                                pitch = 360;
                        }
                        osg::Vec3 tempPos = P->getPosition();                       

                        ////双钮线半径
                        lou = sqrt(maxRadius*maxRadius*cos(2*theta*osg::PI/180));               

                        y = lou*cos((360-theta)*osg::PI/180);
                        z = -(lou*sin((360-theta)*osg::PI/180));
                        x = lou*sin(pitch*osg::PI/180);                       

                        //  compute new position
                        osg::Vec3 newpos;// = P->getPosition() + v * dt;
                        newpos[0] = x;
                        newpos[1] = y;
                        newpos[2] = z;

                        P->setPosition(newpos);   

                }   
          
    }

array

这可能是您的算法存在问题吧？您可以尝试打印一下Newpos来看看自己的计算是否正确？