[Havok] - Havok Physics物理引擎的学习入门
来源：http://www.cnblogs.com/hellohuan/archive/2008/09/03/1283410.html
这几天在看Intel公司刚刚收购开放下载的havok物理引擎，有几点感受。为了研究havok，多么猥琐阴险狡诈的招都被我使了出来，在google上搜过n个正则关键字，去过专业的物理研究论坛，去了havok的官方网站，去了Inter的havok论坛。最后没找到多少资料，没找到一个小例子。所以只得再havok SDKs的Doc和Demo下面逛，看英文，猜代码意图。浑浑噩噩的过了两天。
Havok的demo看过的肯定知道，框架做得相当复杂和诡异，其实并不适合去学习。The framework of the demos is very complex and esoteric. I do not want to understand these complex framework, Because I have a practical application in my framework, companies have the framework of the company, no one will want to make their own framework replaced for a physical engine.这是我发在havok论坛上的帖子。
这也难怪，havok引擎刚刚开放下载，国内应用的都少，教程更是别提。我这几天算是过得比较痛苦，写个简单的入门文档，希望抛砖引玉，帮大家简单入入门，完全没有误人子弟的意思，大家别拍我~ ~我做的主要是havok的刚体碰撞。简单的20 * 20 的正方体自由下落在一平面上。
首先，关于Havok的几个简单的物理概念：
hkpWorld： havok的物理世界场景 场景有其自己的属性，如万有引力等，每个havok的仿真场景可以有一个或者多个hkpWorld的实例。
hkpRigidBody： havok的刚体，刚体可以有自己的运动属性，并能通过setMotionType来设置修改。
hkVisualDebugger： 关于havok的一个很好的调试工具
下面是整个程序的简单代码，没有用图形库输出图形，可以通过havok自带的工具VisualDebugger工具来查看刚体。省略了部分内存方面的代码，我也不是很明白。
int main( int argc, const char ** argv)
{
// 创建physics world 通过填充hkpWorldCinfo
hkpWorld * physicsWorld;
{
// 创建物理世界：
hkpWorldCinfo worldInfo;
worldInfo.m_simulationType = hkpWorldCinfo::SIMULATION_TYPE_MULTITHREADED;
worldInfo.m_gravity. set ( 0.0f , - 9.8f , 0.0f );
physicsWorld = new hkpWorld(worldInfo);
}
physicsWorld -> markForWrite();
{
hkpAgentRegisterUtil::registerAllAgents( physicsWorld -> getCollisionDispatcher() );
}
// 在这里创建我们的刚体模拟
setupPhysics( physicsWorld );
// 初始化VisualDebugger工具
hkVisualDebugger * vdb = setupVisualDebugger(context);
hkReal timestep = 1 .f / 60 .f;
int numSteps = int ( 60 .f / timestep);
for ( int i = 0 ; i < numSteps; ++ i )
{
// 可以得到刚体的一些实时信息，借助图形库进行渲染
hkVector4 pos = g_ball -> getPosition();
printf( " [%f,%f,%f] " n " , pos(0), pos(1), pos(2));
}
hkBaseSystem::quit();
return 0 ;
}
void setupPhysics(hkpWorld * physicsWorld)
{
// 创建地面的 静态刚体
{
hkpRigidBodyCinfo boxInfo;
hkVector4 boxSize( 5.0f , 0.5f , 5.0f );
hkpBoxShape * boxShape = new hkpBoxShape(boxSize);
boxInfo.m_shape = boxShape;
boxInfo.m_motionType = hkpMotion::MOTION_FIXED;
boxInfo.m_position. set ( 0.0f , 0.0f , 0.0f );
boxInfo.m_restitution = 0.9f ;
hkpRigidBody * floor = new hkpRigidBody(boxInfo);
boxShape -> removeReference();
g_ball = floor;
physicsWorld -> addEntity(floor);
floor -> removeReference();
}
// 创建 20 * 20 的盒子
const hkReal boxDim = 0.2f ; // This is the size of the cube side of the box
const hkReal boxRadius = boxDim * 0.01f ; // 盒子的角度 这个有问题 原来这个是大小
const hkReal gapx = boxDim * 0.05f ; // This is the gap betwen boxes 间距是盒子大小的多少比率
const hkReal gapy = boxRadius;
const hkReal heightOffGround = 10.0f ; // This is the height of the pyramid off the gound
int pyramidSize = 20 ; // This is the number of boxes along the base (also vertically)
// 长方体的长宽
hkReal extendedBoxDimX = boxDim + gapx;
// hkReal extendedBoxDimY = boxDim + gapy;
hkReal extendedBoxDimZ = boxDim + gapx;
hkVector4 startPos( 0.0f , heightOffGround + gapy + boxDim * 0.5f , 0.0f );
// Build pyramid 构建
{
hkVector4 boxRadii(boxDim * .5f, boxDim * .5f, boxDim * .5f);
hkpShape * boxShape = new hkpBoxShape( boxRadii , boxRadius );
// 变成横向的一块下来 20 * 20的盒子
for ( int i = 0 ; i < pyramidSize; i ++ )
{
// This constructs a row, from left to right 从左到右
hkVector4 start( 0 , 0 , extendedBoxDimZ * i);
for ( int j = 0 ; j < pyramidSize; j ++ )
{
hkVector4 boxPos(start);
hkVector4 shift(j * extendedBoxDimX, 0.0f , 0.0f );
boxPos.setAdd4(boxPos, shift);
boxPos.setAdd4(boxPos, startPos);
hkpRigidBodyCinfo boxInfo;
// 盒子的质量
boxInfo.m_mass = 100.0f ;
// calculate the correct inertia 惯性与质量有关
hkReal d = boxInfo.m_mass * boxDim * boxDim / 6.0f ;
// for small boxes increase inertia slightly
if ( boxRadius < 0.1f )
{
d *= 2.0f ;
if ( boxRadius < 0.03f )
{
d *= 2.0f ;
}
}
boxInfo.m_inertiaTensor.setDiagonal(d,d,d);
boxInfo.m_shape = boxShape;
boxInfo.m_motionType = hkpMotion::MOTION_DYNAMIC;
boxInfo.m_position = boxPos;
boxInfo.m_restitution = 0.5f ; // 表面的弹力
boxInfo.m_friction = 0.6f ; // 表面的摩擦
boxInfo.m_solverDeactivation = hkpRigidBodyCinfo::SOLVER_DEACTIVATION_MAX;
hkpRigidBody * boxRigidBody = new hkpRigidBody(boxInfo);
// 将所有的刚体添加到世界中
physicsWorld -> addEntity( boxRigidBody );
boxRigidBody -> removeReference();
}
}
boxShape -> removeReference();
}
}
我们只要在自己的代码中调用Havok API构建抽象的物理世界,抽象的物体,即给出物理世界的数据,就能直接通过hkVisualDebugger接口的socket连接发送世界数据到HavokVisualDebugger宿主程序,并在其可视化物理场景视图中显示出来。以下是在工具下的两张截图：
![[Havok] - Havok Physics物理引擎的学习入门_第1张图片](http://img.e-com-net.com/image/product/a385f1e19efc40a2b0672230923e3e01.jpg)
![[Havok] - Havok Physics物理引擎的学习入门_第2张图片](http://img.e-com-net.com/image/product/bbde3ec55bca4ee1ab3d0c933ea27ed1.jpg)
在demo框架中：
![[Havok] - Havok Physics物理引擎的学习入门_第3张图片](http://img.e-com-net.com/image/product/18504a44ce514efd8053e49d1db0a0ea.jpg)
程序的可运行代码：