物理マテリアルの追加パラメータ

AGXDynamicsプラグインを利用時には以下の物理マテリアル(物性)を利用することができます。

サンプル

AGXDynamicsPluginのマテリアルのサンプルが以下にあります。 パラメータ値によって動作結果が異なることを確認することができます。

  • choreonoid/samples/AGXDynamics/agxMaterialSample.cnoid

マテリアルの設定手順

AGXSimulatorでリンク間の摩擦係数、反発係数などは、以下の手順で調整することができます。

  1. マテリアルファイルにMaterial、ContactMaterialを記述

  2. ボディファイルにマテリアルファイルで定義したMaterialを設定

マテリアルファイル

マテリアルファイルは摩擦係数や反発係数などの物性を記述したリストファイルです。このファイルは材質(Material)について、同じまたは異なる材質の接触物性(ContactMaterial)を記述することができます。ここで定義した材質名をBodyファイルに記述することで、モデルに材質を設定することができます。

マテリアルファイルはワールドアイテムのプロパティに設定することで、読み込まれます。デフォルトでは choreonoid/share/default/materials.yaml が設定されており、自動的に読み込まれます。

materials:
  -
    name: Ground
    roughness: 0.5
    viscosity: 0.0
  -
    name: agxMat5
    density: 1.0

contact_materials:
  -
    materials: [ Ground, agxMat5 ]
    youngs_modulus: 1.0e5
    restitution: 0.1
    spook_damping: 0.08
    friction: 0.416667
    surface_viscosity: 1.0e-8
    adhesion_force: 100
    adhesiv_overlap: 0.2
    friction_model: [ iterative, direct ]
    contact_reduction_mode: reduceGeometry
    contact_reduction_bin_resolution: 3

Materialパラメータの説明

バルクマテリアル

パラメータ

デフォルト値

単位

型

意味

density

1000

kg/m3

double

密度。リンクの質量、慣性テンソル、重心の自動計算に利用されます。

youngs_modulus

4.0e8

Pa

double

ヤング率。リンク(剛体)の硬さを表します。値が小さいとリンク同士が侵入しやすくなります。

viscosity

0.5

-

double

反発粘性。反発粘性のペアが反発係数となります。

spook_damping

0.075

s

double

スプークダンパ。リンク同士の侵入の緩和(拘束条件を満たす)に利用します。

poissonRatio

0.3

-

double

ポアソン比（2.27.0.0より非推奨）

サーフェスマテリアル

パラメータ

デフォルト値

単位

型

意味

roughness

0.5

-

double

表面粗さ。表面粗さのペアが摩擦係数となります。

surface_viscosity

5e-09

-

double

表面粘性。接面方向に働く粘性です。表面粘性のペアがContactMaterialのsurface_viscosityとなります。オイルなど濡れを表現する時に利用します。

adhesion_force

0.0

N

double

粘着力。形状が接触している時に法線方向に粘着力が働きます。接着剤のような振る舞いをさせたい時に利用します。

adhesiv_overlap

0.0

m

double

粘着力有効距離。リンクの侵入量>有効距離となると粘着力が有効になります。

注釈

ContactMaterialが定義されているものについては、ContactMaterialのパラメータが利用されます。Materialのサーフェスマテリアルは利用されません。

ワイヤーマテリアル

パラメータ

デフォルト値

単位

型

意味

wire_youngs_modulus_stretch

6e10

Pa

double

引張ヤング率

wire_spook_damping_stretch

0.075

s

double

引張スプークダンパ

wire_youngs_modulus_bend

6e10

Pa

double

曲げヤング率

wire_spook_damping_bend

0.075

s

double

曲げスプークダンパ

ContactMaterialパラメータの説明

パラメータ

デフォルト値

単位

型

意味

youngs_modulus

2.0e8

Pa

double

ヤング率

restitution

0.0

-

double

反発係数。0:完全非弾性衝突、1:完全弾性衝突

spook_damping

0.075

s

double

スプークダンパ

friction

0.5

-

double

摩擦係数

secondary_friction

-1.0

-

double

副方向摩擦係数。摩擦モデルにoriented_box、oriented_scaled_box、constant_normal_force_oriented_box、oriented_iterativeのいずれかを指定した場合に、secondary_friction>=0で有効となります。

surface_viscosity

1.0e-8

-

double

表面粘性係数。摩擦拘束に対するコンプライアンスです。

secondary_surface_viscosity

-1.0

-

double

副方向表面粘性係数。摩擦モデルにoriented_box、oriented_scaled_box、constant_normal_force_oriented_box、oriented_iterativeのいずれかを指定した場合に、secondary_friction>=0で有効となります。

adhesion_force

0.0

N

double

粘着力

adhesiv_overlap

0.0

m

double

粘着力有効距離

friction_model

[ default, default ]

-
string
string
摩擦モデル: default(iterative), iterative, box, scaled_box, oriented_box, oriented_scaled_box, constant_normal_force_oriented_box, oriented_iterative
ソルバ : default(split), split, direct, iterative, direct_and_iterative

contact_reduction_mode

default

-

string

接触点削減方式: default(reduceGeometry), reduceGeometry, reduceALL, reduceNone

contact_reduction_bin_resolution

0

-

uint8_t

接触点削減ビン解像度。0の場合はAGXSimulatorアイテムのパラメータを利用します。

primary_direction

[ 0, 0, 0 ]

Unit vector

Vec3

摩擦モデルorientedBox指定時の主要方向ベクトル

reference_body

-

-

string

摩擦モデルorientedBox指定時の参照Body名

reference_link

-

-

string

摩擦モデルorientedBox指定時の参照Link名

注釈

AGXDynamicsは動摩擦係数、静止摩擦係数の区別がありません。実際、値の差は10-20%程度であり、ほとんどの状況では気にしなくて良いとの考えです。

注釈

摩擦モデルについてはChoreonoid 1.7で利用可能なものから追加されています。また、iterative, constant_normal_force_oriented_box については、それぞれ1.7までの cone, orientedBox に対応します。

ContactMaterialが定義されていない場合

全てのMaterialのペアの物性がContactMaterialに記述されているのが望ましいのですが、難しいと思います。
ContactMaterialが設定されていない場合にはMaterialに記述されているパラメータついて以下の式に従って値を算出します。
Materialにもパラメータが設定されていない場合にはデフォルト値が適用されます。

  • youngs_modulus = (m1.youngs_modulus * m2.youngs_modulus)/(m1.youngs_modulus + m2.youngs_modulus)

  • restitution = sqrt((1-m1.viscosity) * (1-m2.viscosity))

  • spook_damping = max(m1.spook_damping, m2.spook_damping)

  • friction = sqrt(m1.roughness * m2.roughness)

  • surface_viscosity = m1.surface_viscosity + m2.surface_viscosity

  • adhesion_force = m1.adhesion_force + m2.adhesion_force

ボディファイルのマテリアルの記述方法

ボディファイルのマテリアルの記述方法について説明します。
重心、質量、慣性はmassTypeで直接指定か密度を使った自動計算を選択することができます。
デフォルトはmassです。
massType: mass             # 直接指定
massType: density          # 密度を使った自動計算
また、材質はmaterialでマテリアルファイルに定義されているマテリアルか直接指定を選択することができます。
デフォルトはマテリアルファイルに定義されているDefaultまたはdefualtです。
material: Default          # デフォルトマテリアル
material: Ground           # マテリアル
material: useLinkInfo      # 直接指定

以下は記述例です。

注釈

現在のところ、densityを使った重心、質量、慣性テンソルの計算結果はAGXDynamics内部で保持しており、ChorenoidのリンクやGUIから取得、確認することはできません。

従来記法

  • 従来のChoreonoidの記法です

  • 記載されいているcenter_of_mass, mass, inertiaを利用します

  • Materialはdensityを除いて、defaultとなります

  • ContactMaterialはdefault vs xxxxx となります

links:
  -
    name: box1
    center_of_mass: [ 0, 0, 0 ]
    mass: 1.0
    inertia: [
      0.02, 0,    0,
      0,    0.02, 0,
      0,    0,    0.02 ]

マテリアルファイルの利用(推奨)

  • densityを含むマテリアルファイルに記述されたパラメータを使います

links:
  -
    name: box1
    massType: density     # 密度を利用して重心、質量、慣性テンソルを自動計算する
    material: steel       # マテリアルファイルのsteelを利用
    density: 1.0          # densityが記述されている場合はsteelのdensityを
                          # オーバライドして、直接記述されているものを利用します

従来記法+マテリアルリストの利用(推奨)

  • massType: massで直接記述されている重心、質量、慣性テンソルを利用します

  • その他のマテリアルパラメータはマテリアルファイルのsteelを利用します

links:
  -
    name: box1
    massType: mass      # 直接記述された重心、質量、慣性テンソルを利用する
    center_of_mass: [ 0, 0, 0 ]
    mass: 1.0
    inertia: [
      0.02, 0,    0,
      0,    0.02, 0,
      0,    0,    0.02 ]
    material: steel     # マテリアルファイルのsteelを利用

直接記述(非推奨)

  • material: useLinkInfoとするとボディファイルに記述されたMaterialのパラメータを利用することができます

  • ContactMaterialが定義されていない場合 に従ってContactMaterialの値が計算されます

links:
  -
    name: box1
    massType: density
    material: useLinkInfo
    density: 1.0
    youngs_modulus:
    poissonRatio:
    viscosity:
    spook_damping:
    roughness:
    surface_viscosity:
    adhesion_force:
    adhesiv_overlap:

全記述(非推奨)

  • すべてが記述されている場合です

  • どのパラメータが利用されているのか判別がしずらいのでおすすめしません

links:
  -
    name: box1
    massType: density               # 密度を利用して重心、質量、慣性テンソルを自動計算する
    center_of_mass: [ 0, 0, 0 ]
    mass: 1.0
    inertia: [
      0.02, 0,    0,
      0,    0.02, 0,
      0,    0,    0.02 ]
    material: steel                 # materialリストを利用
    density: 1.0                    # 記述されたdensityを利用
    youngs_modulus:                  # 以下は使用されない
    poissonRatio:
    viscosity:
    spook_damping:
    roughness:
    surface_viscosity:
    adhesion_force:
    adhesiv_overlap: