10 月 19
時が経つのは本当に早いもので，この間暑いと思っていたらもうすっかり秋ですね．１０月１６～１８日と妙高高原にある大学のセミナーハウスに研究室の学生と行き，まじめに輪講やテニスを楽しみました．山はもうすっかり紅葉が色づき美しかったのが印象的でした．
さて，Open Dynamics Engine プチキャンプもデビューしてから２ヶ月がすぎ，新しいロボコンマガジン（No.48）が１０月１４日（土）に発売されました．
今回は連載２回目となります．連載が全４回しかないのにも関わらず，グリーディ（greedy)にも最終回は２脚歩行ロボットのプチシミュレータを作ることを目標としているため，今回は車輪型ロボット，次回はロボットアームのプチシミュレータとドンドン作っていきます．
ロボコンマガジンでは紙面が限られているため詳しく説明できません．このページではそれを補うため，サンプルプログラムの全ソースコードと豪華付録？として連載に登場したAPI集を掲載しています．
- 君もサンプルプログラムをゲットしよう！
- プログラム１（自由落下：衝突検出なし）
ゆっくり落下していくと．．．
- プログラム2（自由落下：衝突検出あり）
地面に着陸
- プログラム3（プチシミュレータ）
ロングシュート！
- 注1：これらのプログラムはODE教育用のサンプルプログラムとして書いたもので，ロボコンマガジンの記事で説明している全サンプルプログラムです．上のプログラムを自由にダウンロードして遊んでください．これらのプログラムは無保証です．なお，再配布及び無断転載等はお止めください．
- 注2：mingw, linux, cygwin. OS Xに対応しています．ご自分の環境のmakefile（例えばmakefile.linux）をmakefileにコピーしてからmakeしてください．なお，オリジナルのmakefileはmingw用になっていますのでmingwを利用されている方はコピーする必要はありません．
- makeを実行したとき[cannot found -lode]などとコンピュータに怒られた方はodeのライブラリが所定の位置にインストールされていません。mingwやcygwinターミナルを開き以下のコマンドを実行してください。
- cd /home/ユーザ名/src/ode-0.6
- make install
今回の目玉は上にあるロボカップ中型ロボットリーグのプチシミュレータです．画面の中央より少し左に見える赤い物体はボールで，中央にある黒い物体はロボットで，Fu Fightersのロングシュートをシミュレートしています．遊び方やダウンロードの方法はロボコンマガジンNo.48の記事を読んでください．
- サンプルプログラム
- プログラム１
- プログラム2
- プログラム3 実行開始時の設定ではカメラを遠方に設定しているためロボット等がよく見えません． ウインドウの隅をドラッグしてウインドウを大きくするか，マウスボタン（左，右）を押しなら画面をドラッグすると視線，大きさなどを変更できます．
- makeを実行したとき[-lode cannot found]などと怒られた方はodeのライブラリが所定の位置にインストールされていません。mingwやcygwinターミナルを開き以下のコマンドを実行してください。
- cd /home/ユーザ名/src/ode-0.6
- dBodyID dBodyCreate(dWorldID world) 剛体を作り、そのID番号を戻り値として返す
- const dReal *dBodyGetPosition(dBodyID body) 剛体bodyの重心位置を取得しそのポインタを返す
- const dReal *dBodyGetRotation(dBodyID body) 回転行列を取得しそのポインタを返す
- void dBodySetMass(dBodyID body, const dMass *mass) 剛体bodyに質量パラメータmassを設定する
- void dBodySetPosition(dBodyID body, dReal x, dReal y, dReal z) 剛体bodyの重心を絶対座標系(x,y,z)の位置に設定する
- void dBodySetRotation(dBodyID body, const dMatrix3 R) 剛体bodyを回転行列Rで表す姿勢に設定する
- dGeomID dCreateBox(dSpaceID space, dReal lx, dReal ly, dReal lz) spaceに各辺の長さ(lx,ly,lz)の直方体ジオメトリを生成する
- dGeomID dCreateCylinder(dSpaceID space, dReal r, dReal l) 半径r、長さlの円柱ジオメトリを生成し、そのID番号を返す
- dGeomID dCreatePlane(dSpaceID space ,dReal a, dReal b, dReal c, dReal d) spaceにax+by+cz=dの平面ジオメトリを生成する
- void dGeomSetBody(dGeomID geom, dBodyID body) 物体の2つの属性であるジオメトリgeomと剛体bodyを関連づける
- dSpaceID dHashSpaceCreate(0) 衝突計算用スペースを生成し、そのID番号を返す
- void dJointGroupEmpty(dJointGroupID) 接触点が格納されているジョイントグループを空にする
- void dJointAttach(dJointID joint, dBodyID body1, dBodyID dBody2) 剛体body1とbody2にジョイントjointを取り付ける
- dJointGroupID dJointGroupCreate(0) 接触点のグループを格納するジョイントグループを生成し、そのID番号を返す
- void dJointSetHingeAnchor(dJointID joint, dReal x, dReal y, dReal z) ヒンジジョイントの回転軸の中心点(x,y,z)を設定する
- void dJointSetHingeAxis(dJointID joint, dReal x, dReal y, dReal z) ヒンジジョイントの回転軸ベクトル(x,y,z)を設定する
- void dMassSetBoxTotal(dMass *mass, dReal m, dReal lx, dReal ly, dReal lz) サイズ(lx,ly,lz)、全質mの直方体の質量パラメータmassを計算する
- void dMassSetCylinderTotal(dMass *mass, dReal m, int dir, dReal r, dReal l) 方向dir(1=x,2=y,3=z）、半径r、長さl、全質量mの円柱の質量パラメータmassを計算する
- void dMassSetZero(dMass *mass) 質量パラメータmassの要素を0に初期化する
- void dRFromAxisAndAngle(dMatrix3 R, dReal ax, dReal ay, dReal az, dReal angle)軸ベクトル(ax,ay,az)周りにangle[rad]回転させた時の回転行列Rを求める
- void dSpaceCollide(dSpaceID space, void *data, dNearCallback *callback) spaceで衝突する可能性のある2つのジオメトリを選び出し、引数のcallback関数を呼び出す。dataはcallback関数へ渡すデータへのポインタ
- dWorldID dWorldCreate() 物理計算用ワールドを生成し、そのID(識別)番号を戻り値として返す
- void dWorldDestroy(dWorldID world) 物理計算用ワールドworldを破壊する
- void dWorldSetGravity(dWorldID world, dReal x, dReal y, dReal z) ワールドworldに(x,y,z)[m/s2]の重力加速度を設定する
- void dWorldStep(dWorldID world, dReal step_size) ワールドworldのシミュレーションを1ステップ進める。2番目の引数step_sizeは積分のステップ時間[s]
- void dsDrawBoxD(const double pos[3], const double R[12], const double sides[3])位置posに回転行列Rを取るサイズsidesの直方体を描画する
- void dsSetViewPoint(float xyz[3], float hpr[3]) 視点xyzと視線hpr[°]を設定する
- void dsSimulationLoop(int argc, char **argv,int w, int h, dsfunction *fn) シミュレーションループを呼び出す。1、2番目の引数はmain関数の引数で、3,4番目は描画ウインドウの大きさ(横w,縦h)で、5番目の引数は描画関数クラス
投稿：06年10月19日
shanさん、
新年あけましておめでとうございます。今年もdemura.netをよろしくお願いします。
妙高高原にスキーにいって返事が遅くなりました。今年は雪が少ないです。
さて、simLoop関数はサンプルプログラムでは見えませんがwhileループの中にあり、シミュレーションの各ステップ毎によびだされています。
simLoop関数の中に各ステップ毎のロボットの動作を記述してください。
このサンプルを改造してロボットにあらかじめプログラムした動きをさせたいのですが、どこを変えればいいのでしょうか？
最終的にはstatic void simLoopのところにleft_velとright_velの値を書き出すようにすればいいかと思いますが、どこにループを作ってロボットの動きを計算させたらいいかわかりません。
よろしくお願いします。
片岡さん、
いつも貴重なコメントありがとうございます。大変助かります。
ソースコードの該当箇所を2から2.0に変更しました。
VCでのコンパイルエラーの報告です。
サンプルのプログラム３で２８２行目の部分ですが
dReal vmax45 = POWER * 0.01 * 8.0 / sqrt(2);
関数sqrt()の引数を2.0か2.0fに変えないとオーバーロードが解決出来ないというエラーが発生してしまいます。