次期 Siv3D 用 3D 物理演算機能の紹介

Siv3D AdventCalendar 2017 2 日目の記事です．

概要

本記事では，次期 Siv3D である OpenSiv3D の 3D 機能に向けた，物理演算機能のプロトタイプを紹介します．

Main.cpp

# include <Siv3D.hpp>
# include <HamFramework.hpp>

void Main()
{
    Physics3DWorld world;
    Physics3DBody ground = world.createBox(Vec3(0, -50, 0), Vec3(50, 50, 50), 0.0);

    Array<Physics3DBody> bodies;

    for (int y = 0; y < 10; y++)
    {
        for (int z = 0; z < 10; z++)
        {
            for (int x = 0; x < 10; x++)
            {
                const Vec3 pos(x, 10 + y, z);

                bodies.push_back(world.createBox(pos, Vec3::One * 0.5, 1.0));
            }
        }
    }


    while (System::Update())
    {
        world.update();

        ground.draw();

        for (const auto& body : bodies)
        {
            body.draw();
        }
    }
}

3D 物体の扱いと物理演算

C++で 3D のゲームやアプリケーションなどを開発する際に，3D 物体の扱いには苦労した経験はありませんか．3D の物体同士の相互作用には，面倒な衝突判定・幾何学計算・数値計算を行う必要があります．この面倒さを解決する手段として，3D 物理演算ライブラリが使用される場合があります．3D 物理演算ライブラリとは，3D の形状同士の衝突判定や，物理的な挙動の計算をサポートする機能群であり，Open Dynamics Engine(ODE)，Havok Physics，Bullet Physics などが知られています．現在は次期 Siv3D である OpenSiv3D に向けて，短い C++ ソースコードで 3D 物理演算 (Bullet Physics) を利用できる機能を開発しています．本記事では物理演算機能のプロトタイプを紹介します．

• OpenSiv3Dの3D機能に搭載予定で現時点では使用することはできません．
• 本プロトタイプは Siv3D August 2016 で開発しています．

物体の生成

Main.cpp

# include <Siv3D.hpp>
# include <HamFramework.hpp>

void Main()
{
    Physics3DWorld world;
    Physics3DBody ground = world.createBox(Vec3(0, -50, 0), Vec3(50, 50, 50), 0.0);

    Physics3DBody body = world.createBox({0.0, 5.0, 0.0}, Vec3::One * 0.5, 1.0);

    while (System::Update())
    {
        world.update();

        ground.draw();

        body.draw();
    }
}

マウスとのインタラクション

Main.cpp

# include <Siv3D.hpp>
# include <HamFramework.hpp>

void Main()
{
    Physics3DWorld world;
    Physics3DBody ground = world.createBox(Vec3(0, -50, 0), Vec3(50, 50, 50), 0.0);

    Array<Physics3DBody> bodies;

    bodies.push_back(world.createBox(RandomVec3({ -1.0, 1.0 }, { 0.0, 1.0 }, { -1.0, 1.0 }), Vec3::One * 0.5, 1.0));
    bodies.push_back(world.createSphere(RandomVec3({ -1.0, 1.0 }, { 0.0, 1.0 }, { -1.0, 1.0 }), 0.5, 1.0));
    bodies.push_back(world.createCylinder(RandomVec3({ -1.0, 1.0 }, { 0.0, 1.0 }, { -1.0, 1.0 }), 0.5, 1.0, 1.0));
    bodies.push_back(world.createCapsule(RandomVec3({ -1.0, 1.0 }, { 0.0, 1.0 }, { -1.0, 1.0 }), 0.5, 1.0, 1.0));

    Physics3DCompoundShape compound;

    compound.addShape(Physics3DCapsuleShape(0.3, 2.0));
    compound.addShape(Physics3DCapsuleShape(0.3, 2.0), Vec3(0.0, 1.0, 0.0), Quaternion(Vec3(0.0, 0.0, 1.0), Math::HalfPi));
    compound.addShape(Physics3DCapsuleShape(0.3, 2.0), Vec3(0.0, -1.0, 0.0), Quaternion(Vec3(1.0, 0.0, 0.0), Math::HalfPi));

    bodies.push_back(world.createBody(compound, RandomVec3({ -1.0, 1.0 }, { 0.0, 1.0 }, { -1.0, 1.0 }), 1.0));

    Physics3DMousePickManager mousePick;

    while (System::Update())
    {
        mousePick.update(world);

        world.update();

        ground.draw();

        for (const auto& body : bodies)
        {
            body.draw();
        }

        mousePick.draw();
    }
}

拘束条件

Main.cpp

# include <Siv3D.hpp>
# include <HamFramework.hpp>

void Main()
{
    Physics3DWorld world;
    Physics3DBody ground = world.createBox(Vec3(0, -50, 0), Vec3(50, 50, 50), 0.0);

    Array<Physics3DBody> bodies;
    Array<Physics3DPointToPointConstraint> constraints;

    const double plankWidth = 0.4;
    const double plankHeight = 0.2;
    const double plankBreadth = 1.0;
    const double plankOffset = plankWidth;
    const double bridgeWidth = plankWidth * 10.0 + plankOffset*(10.0 - 1.0);
    const double bridgeHeight = 5.0;
    const double halfBridgeWidth = bridgeWidth*0.5;
    const double mass = 1.0;
    const int lastBoxIndex = 10 - 1;

    for (int i = 0; i < 10; ++i)
    {
        double t = static_cast<double>(i) / lastBoxIndex;
        t = -(t * 2.0 - 1.0) * halfBridgeWidth;
        bodies.push_back(world.createBox(Vec3(t, bridgeHeight, 0), Vec3(plankWidth, plankHeight, plankBreadth), (i == 0 || i == lastBoxIndex) ? 0.0 : mass));
    }

    for (int i = 0; i<10 - 1; ++i)
    {
        constraints.push_back(world.createPointToPointConstraint(bodies[i], bodies[i + 1], Vec3(-0.5, 0, -0.5), Vec3(0.5, 0, -0.5)));
        constraints.push_back(world.createPointToPointConstraint(bodies[i], bodies[i + 1], Vec3(-0.5, 0, 0.5), Vec3(0.5, 0, 0.5)));
    }

    while (System::Update())
    {
        if (Input::MouseL.clicked)
        {
            bodies.push_back(world.createSphere(Vec3(Random(-1.0, 1.0), 10.0 + Random(-1.0, 1.0), Random(-1.0, 1.0)), 0.5, 1.0));
        }

        world.update();

        ground.draw();

        for (const auto& body : bodies)
        {
            body.draw();
        }
    }
}

応用例

ギア

車

おわりに

本記事では，次期 Siv3D である OpenSiv3D の 3D 機能に向けた，物理演算機能のプロトタイプを紹介しました．本機能により，より多くの 3D ゲームや 3D アプリケーションなどが開発されることを期待しています．