ボンドグラフ

単純な質量-バネ-ダンパーシステムと、その等価なボンドグラフ形式

ボンドグラフは、物理的な動的システムのグラフィカルな表現です。これにより、システムを状態空間表現に変換できます。ブロック線図シグナルフローグラフに似ていますが、ボンドグラフのアークは物理エネルギー双方向の交換を表すのに対し、ブロック線図やシグナルフローグラフのアークは一方向の情報の流れを表すという大きな違いがあります。ボンドグラフは、複数のエネルギー領域(例:機械、電気、油圧など)に対応し、領域に依存しません。つまり、ボンドグラフは複数の領域をシームレスに組み込むことができます

ボンドグラフは、「シングルポート」、「ダブルポート」、「マルチポート」の要素を結び付ける「ボンド」で構成されています(詳細は下記参照)。各ボンドは、エネルギーの瞬間的な流れ(dE / dt)または電力を表します。各ボンドのフローは、共役変数に似た電力変数と呼ばれる2つの変数で表され、その積はボンドの瞬間電力です。電力変数は、フロー努力の2つの部分に分けられます。たとえば、電気システムのボンドでは、フローは電流であり、努力は電圧です。この例では、電流と電圧を掛け合わせることで、ボンドの瞬間電力を得ることができます

結合には、ここで簡単に説明し、以下で詳しく説明する2つの特徴があります。1つは「半矢印」符号規則です。これは、正のエネルギーの流れの想定される方向を定義します。電気回路図や自由体図と同様に、正の方向の選択は任意ですが、解析者は選択した定義に一貫していなければならないという注意点があります。もう1つの特徴は「因果関係」です。これは、結合の片端にのみ配置される垂直バーです。これは任意ではありません。以下で説明するように、特定のポートに適切な因果関係を割り当てるための規則と、ポート間の優先順位に関する規則があります。因果関係は、エフォートとフローの間の数学的関係を説明します。因果関係の位置は、どの電力変数が従属し、どの電力変数が独立しているかを示します。

モデル化される物理システムのダイナミクスが大きく変化する時間スケールで動作する場合、ハイブリッド結合グラフを使用することで、高速連続時間の動作を瞬間現象としてモデル化できます。結合グラフはヘンリー・ペインターによって発明されました[1]

ボンドグラフのシステム

多くのシステムは、ボンドグラフで使用される用語で表現できます。これらの用語は以下の表に示されています。

以下の表の表記規則:

  • 有効電力です。
  • 行列オブジェクトです
  • はベクトルオブジェクトです。
  • はxエルミート共役です。x転置複素共役です。xスカラーの場合、エルミート共役は複素共役と同じです。
  • は微分のためのオイラー表記です。ここで:
  • 収束係数:
一般化流量一般化変位一般化努力一般化運動量一般化電力(電力システムの場合はワット単位)一般化エネルギー(電力システムの場合はジュール単位)
名称
説明変位の時間微分静的挙動に関連する性質変位単位あたりのエネルギー作用力の時間積分ある形態から別の形態へのエネルギーの変換閉鎖系における保存量
要素
名称ハイパーランスハイパーリジタンス
コンプライアンスリジタンス
抵抗イナータンス(またはアブラハンスマグナンス
特性電力消費要素電荷貯蔵要素

(状態変数:変位)

(共状態変数:作用力)

電力消費要素運動量貯蔵要素

(状態変数:運動量)

(共状態変数:流れ)

電力消費要素電力消費要素
定量的挙動1次元系(線形)の場合:

1次元系の場合:

インピーダンス:

N次元系の位置エネルギー:

位置エネルギー:

位置コエネルギー

1次元系の場合:

インピーダンス:

1次元系(線形)の場合:

1次元非線形抵抗の電力(は要素によって発生する力):

レイリー電力:

非線形抵抗のレイリー電力:

レイリー力:

N次元系の場合:

1次元系の場合:

インピーダンス:

N次元系運動エネルギー:

運動エネルギー:

運動コエネルギー:

1次元系の場合:

インピーダンス:

1次元系(線形)の場合:

1次元系の場合:

インピーダンス

1次元系(線形)の場合

1次元系の場合

インピーダンス

一般化された挙動能動的な力源からのエネルギー:

ラグランジアン

ハミルトニアン

ハミルトニアン力:

ラグランジアン力:

受動的な力:

電力方程式:

力方程式:

ラグランジアン方程式:

ハミルトニアン方程式:

がコエネルギー、がエネルギー、が状態変数、が共状態変数である場合

線形要素の場合:

縦方向機械動力システム[2]
流れに関連する変数:ポンポン/ポップ:フラウンス/パチパチ:ガクガク/スナップジャーク
加速度 速度 (流れ) :変位(変位) :アブセーション
:アブセレーション:アブセレーション:アブセレーション
作用に関連する変数引っ張り (作用力) 線形運動量 (運動量)
受動要素
コンプライアンス(C)抵抗(R)イナータンス(I)アブハンス(A)マグナンス(M)
バネ :バネ剛性ダンパー ダンパーパラメータ質量 :質量アブラハム・ローレンツ力

ここで

  • 透磁率
  • は電荷
  • 光速
磁気放射反力

ここで

片持ち梁 サイクロトロン放射抵抗

ここで

広い水域における角柱状浮体

ここで

  • :液体の密度
  • :面積
  • :重力加速度
粘性摩擦 粘性摩擦パラメータ
弾性棒

ここで

  • ヤング率
  • :面積
  • :ロッド長さ
運動移動度の逆数。運動移動度はどこにありますか?
ニュートンの万有引力の法則

ここで

  • 万有引力定数
  • :物体1の質量
  • :物体2の質量
空気と相互作用する形状(例:抗力

ここで

  • :接触面積
  • :空気力学的形状係数
  • :流体密度
クーロンの法則

ここで

  • 誘電率
  • :物体1の電荷
  • :物体2の電荷
アブスクエアダンパー[説明が必要]アブスクエアダンパーパラメータはどこにありますか?
カシミール力

ここで

乾燥摩擦

ここで

  • :動摩擦係数
  • :垂直抗力
ビオ・サバールの法則

ここで

  • 浸透率
  • :導線1の電流
  • :導線2の電流
  • :導線の長さ
断熱室内の流体を押すピストン

ここで

  • :ピストンの面積
  • :内部の初期圧力
  • :ピストンの初期位置
  • :ポアソン係数
角動力システム
流れに関連する変数角ジャーク 角加速度
角速度 (流れ) 角変位 (変位)
作用に関連する変数回転力 :力(力)
角運動量 (運動量)
受動要素
コンプライアンス(C)抵抗(R)イナータンス(I)
角バネ定数の逆数。ここでは角バネ定数です。角減衰。ここでは減衰定数です。質量慣性モーメント

種類

質量慣性モーメント。ここでは

ロッドのねじり

ここで

  • :せん断弾性係数
  • :断面極モーメント
  • :ロッド長さ
調速機(例:オルゴールに使用)調速機定数はどこにありますか?
曲げモーメント(片持ち梁)

ここで

  • :ヤング率
  • :断面二次モーメント
  • :梁の長さ
平行力場

ここで

  • :極軸と力場との間の角度(反時計回りが正)
  • :物体(例:振り子)の現在の角度
電力システム
流れに関連する変数:電気慣性電流 (流れ) 電荷 (変位)
作用に関連する変数:膨張電圧 (作用) 鎖交磁束 (運動量)
要素
ハイパーランス(H)コンプライアンス(C)抵抗(R)イナータンス(I)アブハンス(A)
周波数依存負性抵抗器(FDNR)

線形コンデンサ

ここで

  • 誘電率
  • :収束係数
線形抵抗器

ここで

  • :抵抗率
  • :収束係数
  • :温度係数
  • :現在の温度
  • :基準温度
線形インダクタ(ソレノイド)

ここで

  • :透磁率
  • :巻数
  • :面積
  • :長さ
周波数依存負性コンダクタンス(FDNC)

種類

ダイオード

ここで

  • 理想係数
  • :熱電圧
  • :漏れ電流
トロイド

ここで

  • :透磁率
  • :巻数
  • :断面積
  • :中心線に対するトロイドの半径
イントラジャイレータインタージャイレータ
コンプライアントジャイレータレジスティブジャイレータイナータントジャイレータ
ホール効果素子

ここで

  • :ホール係数
  • :垂直磁気誘導場
  • :垂直厚さ
誘導電動機

ここで

  • :電気角
  • :回転子バーの電圧
  • :固定子バーの電圧
  • :固定子バーを流れる電流
  • :回転子バーを流れる電流
DC電動機

ここで

  • :電磁トルク
  • :軸角速度
  • 界磁電流
  • :電機子電流
ファラデージャイレータ

ここで

  • :力
  • :ロッド速度
  • :ロッド電圧
  • :磁気誘導場
  • :ロッド長さ
ファラデーディスク

ここで

  • :磁気誘導場
  • :ディスク半径
変圧器内変圧器間
変圧器(交流信号のみ)

油圧/空気圧動力システム
流れに関連する変数体積流量 (フロー) 容積 (変位)
作用に関連する変数圧力 (作用力) :流体運動量(モーメント)
要素
コンプライアンス(C)抵抗(R)イナータンス(I)
配管弾性

ここで

  • :公称配管半径
  • :配管容積(応力なし)
  • :壁厚
  • :ヤング率
ダルシースポンジ

ここで

  • :動粘性
  • 浸透率
  • :収束係数
配管内の流体慣性

ここで

  • :流体密度
  • :収束係数
圧縮性流体(近似値)

ここで

  • :配管容積(応力なし)
  • :体積弾性率
  • :基準圧力における気体密度
  • :音速
バルブ

  • :流量係数
  • :流体密度
  • :流体通路の最小面積
面積を持つタンク

ここで

  • :地面からの高さ
  • :入口高さ
  • :流体密度
  • :曲線のスケーリング(面積寸法)
  • :寸法補正
  • :曲線パラメータ
  • :重力加速度
  • ケース(角柱の場合):
  • ケース(電流ダイオードの場合):
円筒のポアズイユ抵抗

ここで

  • :動粘性
  • :パイプの長さ
  • :パイプの半径
等温チャンバー

チャンバー定数

乱流抵抗

経験的パラメータ

圧縮性流体

は体積弾性係数

ノズル

ここで

  • :流体密度
  • :出口面積
  • :入口面積
断熱ブラダー

ここで

  • :基準圧力
  • :基準容積
  • :ポアソン係数
逆止弁

ここで

  • :経験定数
  • :逆バイアス絶対流量
ジャイレータ・コンデンサ発電システム
流れに関連する変数磁束 (変位)
作用に関連する変数起磁力 (運動量)
要素
コンプライアンス(C)抵抗(R)イナータンス(I)
パーミアンス

ここで

  • :透磁率
  • :収束係数
磁気複素インピーダンス(

磁気複素インダクタンス(

重力発電システム
流れに関連する変数:重力電流ループ(流量) :重力電荷(変位)
作用に関連する変数:重力電圧(作用力) :重力運動量(運動量)
要素
コンプライアンス(C)抵抗(R)イナータンス(I)
重力容量

種類

重力軌道抵抗

種類

重力インダクタンス

種類

電場体積電力密度システム
流れに関連する変数電流密度 (流れ) 変位電場 (変位)
作用に関連する変数電場 (作用力) :磁気ポテンシャルベクトル(運動量)
要素
コンプライアンス(C)抵抗(R)体積電力密度イナータンス(
電気誘電率

種類

電気抵抗率

種類

透磁率
磁場体積電力密度システム
流れに関連する変数:磁束密度(変位)
作用に関連する変数:磁場強度(作用力)
要素
コンプライアンス(C)抵抗(R)体積電力密度イナータンス(
磁気回路の透磁率

種類

重力電気体積電力密度システム
流れに関連する変数:質量流束(流れ) :質量の累積流束(変位)
作用に関連する変数重力加速度 (作用力)
要素
コンプライアンス(C)抵抗(R)体積電力密度イナータンス(
重力誘電率

種類

重力透磁率
重力磁気体積電力密度システム
流れに関連する変数:重力磁場(変位)
作用に関連する変数:重力磁場の強さ(作用)
要素
コンプライアンス(C)抵抗(R)体積電力密度イナータンス(
重力透磁率

熱出力-温度システム
流れに関連する変数:熱速度(流量) :全熱量(変位)
作用に関連する変数:温度(作用)
要素
コンプライアンス(C)抵抗(R)イナータンス(I)
等圧熱

ここで

  • :物体の質量密度
  • :物体の体積
  • :圧力比熱容量
伝導抵抗

ここで

  • :熱伝導率
  • :収束係数
等体積熱

定積熱容量

対流抵抗

ここで

  • :対流係数
  • :界面積
等温熱

ここで

  • 普遍気体定数
  • :モル数
  • :最終体積
  • :初期体積
シュテファン・ボルツマンの法則

ここで

  • :放射率
  • :シュテファン・ボルツマン定数
  • :界面積
連続体力学における体積出力密度システム
流れに関連する変数:ひずみ速度(流量) :ひずみ(変位)
作用に関連する変数:応力(作用)
要素
コンプライアンス(C)抵抗(R)体積電力密度イナータンス(
剛性の逆数

種類

粘性

種類

電力密度慣性モーメント:物質の密度

その他のシステム:

  • 熱力学的発電システム(フローはエントロピー速度、努力値は温度)
  • 電気化学的発電システム(フローは化学活性、努力値は化学ポテンシャル)
  • 熱化学的な発電システム(フローは質量速度、努力値は質量比エンタルピー)
  • マクロ経済学の通貨レートシステム(変位は商品、努力値は商品あたりの価格)
  • ミクロ経済学の通貨レートシステム(変位は人口、努力値は一人当たりGDP)

状態の四面体

状態の四面体

状態の四面体は、努力値とフローの変換を図式的に示す四面体です。隣の図は、一般化された形の四面体を示しています。四面体はエネルギー領域に応じて変更できます

状態正四面体を用いることで、正四面体上の任意の変数間の数学的関係を見つけることができます。これは、図の周りの矢印をたどり、途中で定数を掛け合わせることで行われます。例えば、一般化流量と一般化変位の関係を求めたい場合は、f ( t )から始めて積分し、 q ( t )を求めます。その他の方程式の例は以下をご覧ください。

一般化変位と一般化流量の関係

一般化流量と一般化作用力の関係

一般化流量と一般化運動量との関係

一般化運動量と一般化作用力の関係

定数Cを含む一般化流量と一般化作用力の関係

エネルギー領域を切り替えても、すべての数学的関係は同じままで、記号のみが変わります。これは以下の例で確認できます。

変位と速度の関係

電流と電圧の関係。 これはオームの法則としても知られています

力と変位の関係。フックの法則としても知られています。この式では負の符号が省略されています。これは、結合グラフにおける矢印の向きに符号が考慮されるためです。

電力システムの場合、共振周波数の式は次のとおりです。

電力密度システムの場合、共振波の速度の式は次のとおりです。

構成要素

エンジンがシャフトを介して車輪に接続されている場合、動力は回転機械領域で伝達されます。つまり、力と流量はそれぞれトルク (τ) と角速度 (ω) です。単語結合グラフは、単語によって構成要素を定義する結合グラフへの第一歩です。単語結合グラフとして、このシステムは次のようになります。半矢印は符号規則を提供するために使用されます。つまり、τωが正のときにエンジンが仕事をしている場合、図は次のようになります。このシステムは、より一般的な方法で表現することもできます。これは、単語の使用から、同じ項目を表す記号に変更することです。これらの記号は、上記で説明した一般化された形式に基づいています。エンジンが車輪にトルクをかけているため、システムの力の源として表されます。車輪は、システム上のインピーダンスとして表すことができます。さらに、トルクと角速度の記号は削除され、力と流量を表す一般化された記号に置き換えられます。この例では必須ではありませんが、式で追跡できるように、結合に番号を付けることがよくあります。簡略化された図を以下に示します。

結合において、エフォートは常にフローより上にあるため、関連情報を失うことなく、エフォート記号とフロー記号を完全に省略することも可能です。ただし、結合番号は省略しないでください。例を以下に示します。

結合番号は、後でボンドグラフから状態空間方程式に変換するときに重要になります。

要素の関連付け

直列関連付け

要素が次の動作を持つと仮定します。ここで、は一般的な関数(入力を微分/積分することもできます)、は要素の定数です。次に、1-接合にこのタイプの要素が多数あると仮定します。すると、接合全体の電圧は次のようになります。

並列関連付け

要素が次の動作を持つと仮定します。ここで、は一般的な関数(入力を微分/積分することもできます)、は要素の定数です。次に、0-接合にこのタイプの要素が多数あると仮定します。すると、次の式は有効です。

単ポート要素

単ポート要素とは、ボンドグラフ内の1つのポートしか持てない要素です。

ソースとシンク

ソースは、システムへの入力を表す要素です。システムに力またはフローを入力します。力またはフローを表す小文字の「e」または「f」が付いた大文字の「S」で表されます。ソースの矢印は常に要素から外を向いています。ソースの例としては、モーター(力、トルクの源)、電圧源(力の源)、電流源(フローの源)などがあります。

ここで、Jはジャンクションを示します

シンクはシステムの出力を表す要素です。ソースと同じように表されますが、矢印は要素から離れる方向ではなく、要素に向かう方向を指しています。

慣性

慣性要素は大文字の「I」で表され、常に電力が流入します。慣性要素はエネルギーを蓄える要素です。最も一般的なのは、機械システムでは質量、電気システムではインダクタです。

抵抗

抵抗要素は大文字の「R」で表され、常に電力が流入します。抵抗要素はエネルギーを消散させる要素です。最も一般的なのは、機械システムではダンパー、電気システムでは抵抗器です。

コンプライアンス

コンプライアンス要素は大文字の「C」で表され、常に電力が流入します。コンプライアンス要素は位置エネルギーを蓄える要素です。最も一般的なのは、機械システムではバネ、電気システムではコンデンサです。

2ポート要素

これらの要素には2つのポートがあります。システム間またはシステム内の電力を変換するために使用されます。一方から他方へ変換する場合、転送中に電力が失われることはありません。要素には定数が与えられます。この定数は、使用されている要素に応じて、変圧器定数またはジャイレータ定数と呼ばれます。これらの定数は通常、要素の下に比率として表示されます。

変圧器

変圧器は、流入量と流出量、および力と流入量との関係を適用します。例としては、理想的な電気変圧器てこなどがあります。

r変圧器の係数を表します 。これは、そしてを 意味します。

ジャイレータ

ジャイレータ、流入量と流出量、および力と流入量との関係を適用します。ジャイレータの例としては、電圧(電気力)を角速度(角機械流)に変換するDCモーターがあります。

つまり、そしてを意味します。

マルチポート要素

ジャンクションは、他の要素とは異なり、入力または出力の任意の数のポートを持つことができます。ジャンクションはポート間で電力を分割します。0ジャンクションと1ジャンクションという2つの異なるジャンクションがあり、これらはエフォートとフローがどのように伝達されるかのみ異なります。同じジャンクションを直列に接続することはできますが、異なるジャンクションを直列に接続することはできません。

0ジャンクション

0ジャンクションは、すべての力値(およびその時間積分/微分)が結合全体で等しいように動作しますが、流入するフロー値の合計は流出するフロー値の合計に等しくなります。つまり、すべてのフローの合計はゼロになります。電気回路では、0ジャンクションはノードであり、そのノードのすべてのコンポーネントによって共有される電圧を表します。機械回路では、0ジャンクションはコンポーネント間のジョイントであり、それに接続されているすべてのコンポーネントによって共有される力を表します。

例を以下に示します。

結果として得られる方程式:

1ジャンクション

1-ジャンクションは0-ジャンクションとは逆の動作をします。1-ジャンクションは、すべてのフロー値(およびその時間積分/微分)が結合全体で等しくなりますが、入力の力の合計は出力の力の合計と等しくなります。つまり、すべての力の合計はゼロになります。電気回路では、1-ジャンクションはコンポーネント間の直列接続を表します。機械回路では、1-ジャンクションは接続されているすべてのコンポーネントで共有される速度を表します。

例を以下に示します。

結果として得られる方程式:

因果関係

ボンドグラフには因果関係の概念があり、ボンドのどちら側が瞬間的な努力を決定し、どちら側が瞬間的な流れを決定するかを示します。システムを記述する動的方程式を定式化する際に、因果関係は各モデリング要素について、どの変数が従属し、どの変数が独立しているかを定義します。因果関係をあるモデリング要素から別のモデリング要素にグラフィカルに伝播させることで、大規模モデルの分析が容易になります。ボンドグラフモデルで因果関係の割り当てを完了することで、代数ループが存在するモデリング状況、つまり変数が自身の関数として再帰的に定義されている状況を検出できます

因果関係の例として、バッテリーと直列に接続されたコンデンサを考えてみましょう。コンデンサを瞬時に充電することは物理的に不可能であるため、コンデンサと並列に接続されたものは、必然的にコンデンサの両端の電圧(効果変数)と同じになります。同様に、インダクタは磁束を瞬時に変化させることができないため、インダクタと直列に接続された部品は、必然的にインダクタと同じ電流量を持ちます。コンデンサとインダクタは受動素子であるため、それぞれの電圧と電流量を無期限に維持することはできません。接続されている部品は、それぞれの電圧と電流量に影響を与えますが、それぞれ電流と電圧に影響を与えることによって間接的にのみ影響を与えます。

注:因果関係は対称的な関係です。一方が効果を「引き起こす」場合、もう一方は電流量を「引き起こす」ことになります。

結合グラフ表記では、因果関係の線を電力​​結合の一方の端に追加して、こちら側が電流量を定義していることを示すことができます。したがって、因果関係の線とは反対側が効果量を制御します

流れの源()は流れを定義するため、因果ストロークを担います。力の源()は力を定義するため、もう一方の端は因果ストロークを担います。

一定トルクのモーターで車輪を駆動する、つまり力の源()を考えてみましょう。これは次のように描かれます。

対称的に、因果ストローク(この場合は車輪)側が結合の流れを定義します

因果関係は適合性制約をもたらします。明らかに、電力結合の一方の端だけが努力を定義できるため、結合の一方の端(もう一方の端)だけが因果ストロークを持つことができます。さらに、時間依存の挙動を示す2つの受動的なコンポーネントとには、1種類の因果関係しかありません。つまり、コンポーネントはフローを決定し、コンポーネントは努力を定義します。したがって、接合点から、好ましい因果関係の方向性は次のとおりです。

これらの要素に対してこれが好ましい方法である理由は、状態の四面体によって示される方程式を考慮すると、さらに分析できます。

結果として得られる方程式には、独立した電力変数の積分が含まれます。これは、因果関係を逆にして微分を得る結果よりも好ましいです。方程式は以下の通りです

結合グラフでは、これらの要素の1つに因果バーが非優先的な方法で存在する可能性があります。そのような場合、「因果衝突」がその結合で発生したと言われます。因果衝突の結果は、グラフの状態空間方程式を記述するときにのみ見られます。これについては、そのセクションで詳しく説明します。

抵抗器には時間依存の動作はありません。電圧を印加すると瞬時に流れが生じ、または流れを印加すると瞬時に電圧が生じるため、抵抗器は因果結合のどちらの端にも存在できます。

変圧器は受動的であり、エネルギーを消散も蓄積もしないため、因果関係は変圧器を通過します。

ジャイレータは流れを力に変換し、力は流れに変換します。したがって、片側で流れが発生すると、反対側で力が発生し、その逆も同様です。

接合

0ジャンクションでは、エフォートは等しく、1ジャンクションではフローは等しくなります。したがって、因果関係のある結合では、0ジャンクションでは1つの結合のみがエフォートを引き起こし、1ジャンクションでは1つの結合のみがフローを引き起こします。したがって、ジャンクションの1つの結合の因果関係が分かれば、他の結合の因果関係も分かります。その1つの結合は「強い結合」と呼ばれます。簡単に言うと、0ジャンクションには1つの因果バーが必要であり、1ジャンクションには1つを除くすべての因果バーが必要です。

因果関係の決定

結合グラフの因果関係を決定するには、特定の手順に従う必要があります。これらの手順は次のとおりです。

  1. ソース因果バーを描く
  2. C結合とI結合の優先因果関係を描く
  3. 0ジャンクションと1ジャンクション、トランスフォーマー、ジャイレータの因果バーを描く
  4. R結合の因果バーを描く
  5. 因果関係の矛盾が生じた場合は、CまたはI結合を微分化に変更してください。

手順の概要を以下に示します

最初のステップは、ソースの因果関係を描くことです。ソースには因果関係は1つしかありません。その結果、下のグラフが生まれます。

次のステップは、C結合の好ましい因果関係を描くことです。

次に、0ジャンクション、1ジャンクション、トランスフォーマー、ジャイレータに因果関係を適用します

しかし、左側の0ジャンクションに問題があります。0ジャンクションにはジャンクションに因果バーが2つありますが、0ジャンクションにはジャンクションに1つだけ必要です。これは、優先因果関係にあることが原因です。これを修正する唯一の方法は、因果バーを反転することです。これにより因果矛盾が発生します。修正されたグラフを以下に示します。因果矛盾を表しています。

他のシステムからの変換

ボンドグラフを使用する主な利点の1つは、ボンドグラフを作成すれば、元のエネルギー領域は関係ないことです。以下は、エネルギー領域からボンドグラフに変換するときに適用する手順の一部です。

電磁気

ボンドグラフとして電磁気の問題を解く手順は次のとおりです。

  1. 各ノードに0ジャンクションを配置する
  2. 1ジャンクションを持つ電源、R、I、C、TR、GYボンドを挿入する
  3. 接地(変圧器またはジャイレータがある場合は両側)
  4. 電力の流れの方向を割り当てる
  5. 簡素化する

これらの手順は、以下の例でより明確に示されています。

線形機械

線形機械の問題をボンドグラフとして解く手順は次のとおりです。

  1. 異なる速度ごとに1-接合を配置する(通常は質量)
  2. 作用する1-接合間の、それぞれの0-接合にR結合とC結合を挿入する
  3. 作用する1-接合にソース結合とI結合を挿入する
  4. 電力の流れの方向を割り当てる
  5. 簡素化する

これらの手順は、以下の例でより明確に示されています。

簡素化

簡素化の手順は、システムが電磁力か線形機械かに関係なく同じです。手順は次のとおりです。

  1. 電力がゼロの結合を削除する(接地または速度がゼロのため)
  2. 結合が3つ未満の0-接合と1-接合を削除する
  3. 並列電力を簡素化する
  4. 0-接合を直列に結合する
  5. 1-接合を直列に結合する

これらの手順は、以下の例でより明確に示されています。

並列電力

並列電力とは、ボンドグラフ内で電力が並列に流れることです。並列電力の例を以下に示します

並列パワーは、0-ジャンクションと1-ジャンクションにおけるエフォートとフローの関係を思い出すことで簡略化できます。並列パワーを解くには、まずジャンクションのすべての方程式を書き出す必要があります。提供されている例では、方程式は以下の通りです。(エフォート/フロー変数が表す番号結合に注意してください。)

これらの方程式を操作することで、並列パワーを表す同等の0-ジャンクションと1-ジャンクションのセットを見つけることができます

例えば、となので方程式の変数を置き換えると、となり、なので、となることがわかります。この2つの力の変数が等しい関係は、0-接合によって説明できます。他の方程式を操作することで、1-接合の関係を記述するを見つけることができます。必要な関係を決定したら、新しい接合を使用して並列電力セクションを再描画できます。例の結果を以下に示します。

簡単な電気システム

電圧源、抵抗器、コンデンサが直列に接続された簡単な電気回路

最初のステップは、すべてのノードに0ジャンクションを描くことです。

次のステップは、それぞれの1ジャンクションで作用するすべての要素を追加することです。

次のステップは、グランドを選択することです。グランドは、電圧がないと仮定される単純な0ジャンクションです。この場合、グランドは上記の下線で示されている左下の0ジャンクションとして選択されます。次のステップは、ボンドグラフのすべての矢印を描くことです。ジャンクションの矢印はグランドに向ける必要があります(電流と同様の経路をたどります)。抵抗、イナータンス、コンプライアンス要素の場合、矢印は常に要素に向けられます。矢印を描いた結果は以下に示されており、0ジャンクションはグランドとして星印でマークされています。

これでボンドグラフが作成されたので、簡素化のプロセスを開始できます。最初のステップは、すべてのグランドノードを削除することです。下部の0ジャンクションは両方とも接地されているため、削除できます。結果を以下に示します

次に、結合数が3未満の接合部を削除します。これは、流れと力はこれらの接合部を変更せずに通過するため、削除することで描画量を減らすことができるためです。結果は以下の通りです。

最後のステップは、結合グラフに因果関係を適用することです。因果関係の適用については上記で説明しました。最終的な結合グラフを以下に示します。

高度な電気システム

電流源、抵抗器、コンデンサ、変圧器を備えた、より高度な電気システム

この回路で手順に従うと、簡略化前の結合グラフは以下のようになります。星印でマークされたノードはグランドを表します。

結合グラフを簡略化すると、以下の画像になります。

最後に、因果関係を適用すると、以下の結合グラフが得られます。星印の結合は因果関係の矛盾を表します。

単純な線形機械

壁に取り付けられたバネ上の質量で構成される、単純な線形機械システム。質量には何らかの力が加えられています。システムの画像を以下に示します

機械システムの場合、最初のステップは、それぞれの異なる速度に1-接合を配置することです。この場合、質量と壁の2つの異なる速度があります。通常、参照のために1-接合にラベルを付けると便利です。結果は以下の通りです。

次のステップは、R結合とC結合を、それらが作用する1-接合間のそれぞれの0-接合に描くことです。この例では、これらの結合はバネのC結合の1つだけです。これは、質量を表す1-接合と壁を表す1-接合の間で作用します。結果は以下の通りです

次に、作用する1-接合部にソースとI結合を追加します。1つのソース(作用力の源)と1つのI結合(質量の質量)があり、どちらも質量の1-接合部に作用します。結果を以下に示します。

次に、電力の流れを割り当てます。電気の例と同様に、電力は地面(この場合は壁の1-接合部)に向かって流れる必要があります。例外はR、C、またはI結合で、これらは常に要素に向かっています。結果として得られる結合グラフを以下に示します

結合グラフが生成されたので、簡略化できます。壁は接地されている(速度がゼロ)ため、その接合を削除できます。したがって、C結合がある0-接合も削除できます。削除すると、結合が3つ未満になります。簡略化された結合グラフは以下の通りです。

最後のステップは因果関係を適用することです。最終的な結合グラフは以下の通りです。

高度な線形機械

より高度な線形機械システムは以下の通りです。

上記の例と同様に、最初のステップは、各遠方速度に1-接合を作成することです。この例では、質量1、質量2、そして壁の3つの遠方速度があります。次に、すべての結合を接続し、電力の流れを割り当てます。結合は以下の通りです

次に、壁の1接合部と3つ未満の結合部を削除することで、ボンドグラフを簡素化するプロセスを開始します。ボンドグラフを以下に示します。

ボンドグラフには並列べき乗があります。並列べき乗の解き方は上記で説明しました。解き方を以下に示します。

最後に因果関係を適用します。最終的なボンドグラフを以下に示します。

状態方程式

ボンドグラフが完成すると、それを利用してシステムの状態空間表現方程式を生成できます。状態空間表現は、複雑な多階微分システムを1階方程式の連立として解くことができるため、特に強力です。状態方程式の一般的な形は、状態変数、つまりシステムの未知数の列行列です。状態変数の時間微分です。システムの入力の列行列です。そして、と、はシステムに基づく定数の行列です。システムの状態変数は、因果関係の矛盾のない各C結合とI結合の値です各I結合にはが、各C結合にはが与えられます

例えば、次のボンドグラフがあるとします。

次の、、および行列 が得られます

の行列は、状態の四面体で説明したように、状態変数とそれぞれの要素の関係を決定することによって解かれます。状態方程式を解くための最初のステップは、結合グラフのすべての支配方程式をリストすることです。下の表は、結合とその支配方程式の関係を示しています。

要素の種類結合名因果関係のある結合支配方程式
シングルポート要素ソース/シンク、S
抵抗、R:

散逸エネルギー

イナータンス、I:

運動エネルギー

コンプライアンス、C:

位置エネルギー

2端子素子変圧器、TR

ジャイレータ、GY

多端子素子0端子1つだけの

端子における因果関係バー

1端子1つだけの因果関係バー

端子から離れた場所

「♦」は優先因果関係を示します。

提供された例では、

支配方程式は次のとおりです。

これらの方程式を操作して状態方程式を生成できます。この例では、、、およびに関してとを関連付ける方程式を見つけようとしています

まず、状態の四面体から、式2から始めて、式4を のように並べ替えることができることを思い出してください式4では、式3によりを に置き換えることができ、式5に置き換えることができます。同様に、式7を使用して を に置き換えることができ、式7では を に置き換えることができ、式10に置き換えることができます。これらの置き換えにより、以下に示す最初の状態方程式が得られます。

2番目の状態方程式も同様に を思い出すことで解くことができます。2番目の状態方程式を以下に示します。

両方の方程式をさらに行列形式に並べ替えることができます。その結果を以下に示します。

この時点で、方程式は他の状態空間表現問題と同様に扱うことができます

ボンドグラフシステム生物学

ボンドグラフは、生理学や生物学を含む生命科学に関連するシステムをモデル化するために使用されてきました。[3]特に、生物物理学的システムをモデル化するためのボンドグラフの使用は、1970年代初頭にAharon KatchalskyGeorge Oster 、Alan Perelsonによって導入されました。 [4] [5]最近では、これらのアイデアはシステム生物学 の文脈で使用され、細胞生物学の生化学反応システムをモデル化するためのエネルギーベースのアプローチが提供されました。 [ 6] [7] [8]およびフィジオーム全体をモデル化するためのアプローチが提供されました[9]

ボンドグラフアプローチには、フィジュームの大規模な計算モデルを構築するための優れた基盤となる多くの特徴があります

  • エネルギーベースである
    • モデルは物理的に妥当です[10]
    • 反応速度論の詳細なバランス(ヴェグシャイダーの条件)は自動的に満たされます[11]
    • エネルギーの使用量と散逸を直接考慮できます
  • モジュラー型です。ボンドグラフのコンポーネント自体がボンドグラフになることができます[12] [13]
  • 物理領域間のエネルギー変換は単純に表される[9] - 以下を参照
  • シミュレーションに使用できるシンボリックコードは、自動的に生成できる[14]

変数

生化学システムのボンドグラフ変数は以下のとおりです。

効果とフロー(μv)の積は、ボンドグラフの定式化では常に、電力(J/s)であることに注意してください。

構成要素

以下に詳述するように、生化学システムをモデル化するために使用されるコンポーネントの主な特徴は、RコンポーネントとCコンポーネントが非線形であること、Iコンポーネントは不要であること、Rコンポーネントが2ポートReコンポーネントに置き換えられていることです。[8]

ジャンクションコンポーネント

ボンドグラフの0(ゼロ)コンポーネントと1(イチ)コンポーネントは、この文脈では違いはありません。

Cコンポーネント

Cコンポーネントはフローを積分して化学種の量を与えます。

効果、化学ポテンシャルは次の式で与えられます。[8] [15]

ここで、はに対応する化学ポテンシャル気体定数はケルビン単位の絶対温度です。

の式は、簡略化して次のように書き直すことができます。

ここで

この特定のC成分は特殊な形状をしているため、特別なRe成分に類似してCeという特別な名前が付けられることがあります。

Re成分

Re(反応)成分には、化学反応の左側(順方向)と右側(逆方向)に対応する2つのエネルギーポートがあります。順方向親和力と逆方向親和力は、それぞれ反応の左側と右側の種による正味の化学ポテンシャルとして定義されます。Re成分は、反応流を次のように与えます。[8]

ここで、(mol/s)は反応速度定数です。

反応流は、順方向親和力と逆方向親和力の差ではなく、両方に依存するため、通常のR成分を1接合定式で使用することはできないことに注意してください[6]

単純な反応のモデル化

反応[6][8]

反応の結合グラフ。C:AとC:Bは種Aと種Bを表し、Re:r1は反応を表します。結合と結合点は3つの成分を接続します。各結合に関連する努力と流れが示されています。

単純な反応は3つの成分で表されます。

  • C:Aは化学ポテンシャルを持つ種Aを表します。流れはです
  • C:Bは化学ポテンシャルを持つ種Bを表します。流れはです
  • Re_r1は流れ、順親和力、逆親和力を持つ反応を表します
  • 結合と結合点は、示されている努力と流れの変数で化学エネルギーを伝達します。

上記の式を用いると、流れは次のように表されます 。

ここで、下付き文字は種に対応します。これは単純な質量作用方程式です。

ここで

酵素触媒反応

酵素Eと酵素複合体Cを介して、種Aを種Bに可逆的に変換する酵素触媒反応。4つのC成分は4つの種を表し、2つのRe成分は2つの反応r1とr2を表す。結合グラフは、酵素Eが反応r2の逆側から反応r1の順側へリサイクルされていることを示す。

Keener & Sneyd [15]のセクション1.4で議論されているように、酵素と酵素複合体を介して種を種に可逆的に変換する酵素触媒反応は、次の2つの反応として表すことができます。

酵素複合体は種から形成され、種に分解して 酵素を放出します。図に示されている結合グラフは、酵素がどのようにリサイクルされるかを示しています。

結合グラフは、一般化ミカエリス・メンテン形であるこれらの反応の特性を導き出すために使用できます[6]

エネルギー変換

生命科学におけるエネルギー変換。結合グラフのTF成分は、エネルギー領域内またはエネルギー領域間のエネルギー変換を表します。この図は、エネルギー(J/mol)と流量(mol/s)を持つ化学領域と、エネルギーと流量を持つ一般領域間の変換に焦点を当てています。変圧器の係数はmです

結合グラフTF(トランスフォーマー)コンポーネントは、エネルギー領域内またはエネルギー領域間のエネルギー変換を表します。[16](TFコンポーネントはTD(トランスダクション)コンポーネント[4] [5]とも呼ばれています。TFの方が広く使用されています。[16] )このセクションでは、エネルギー効率(J/mol)と流量(mol/s)を持つ化学領域と、エネルギー効率と流量を持つ一般領域間の変換に焦点を当てています

TFコンポーネントの主な特徴は、エネルギーを散逸させることなく伝達することです。[16]したがって、図を参照すると、

変成器は係数m(適切な単位付き)を持つため、

エネルギー式から次のことが分かります。

化学量論

TF成分を用いた化学反応を表す結合グラフ。ここで、mは正の整数です。

化学反応の化学量論は各化学種の発生数を決定します。したがって、たとえば、この反応は1モルの化学種をmモルの化学種に変換 します

の場合は、上記の最初の例の単純な反応に対応します。一般的な化学種についても同じアプローチを用いると、反応の流れは次のようになります。[6] [8]

化学電気変換

このセクションでは、一般的な領域が電気領域であり、努力が(電気)電圧)、流れが電流( )である場合を考察します。分子の電荷が(クーロン)である荷電イオン流れを考えます。ここで、はクーロンで測定された電子の電荷です。したがって、イオン1モルに関連する電荷は、ここでアボガドロ定数です。等価電流は

はファラデー定数です。したがって、対応するTF係数は次のようになります。

(C/mol)

繰り返しますが、

この文脈では、結合グラフのTF成分は、活動電位[17] 膜輸送体[18] 心臓活動電位[19] 、 およびミトコンドリア電子伝達系[17]に関連するエネルギーフローをモデル化するために使用できます

化学機械変換

アクチンのような長くて硬い分子を考えてみましょう。長さ(m)のサブユニットが速度(mol/sec)で加えられます。先端速度は次式で表されます。[20]

ここではアボガドロ定数です

したがって、弾性率(m/mol)と

ここでは先端における対応する力です。

これらの式は、アクチンフィラメントの力/速度曲線を生成するために使用されています[20]。このアプローチは、ブラウン運動ラチェット法[21]の有用な代替手段を提供します。なぜなら、ボンドグラフTFコンポーネントは、細胞システムのモジュラーボンドグラフモデルで使用できる可能性があるからです。[7]

ボンドグラフモデリングに関する国際会議(ECMSおよびICBGM)

ボンドグラフモデリングに関する参考文献は、以下の会議から抜粋できます。

  • ECMS-2013 第27回ヨーロッパモデリング・シミュレーション会議、2013年5月27~30日、ノルウェー、オーレスン
  • ECMS-2008 第22回ヨーロッパモデリング・シミュレーション会議、2008年6月3日~6日、キプロス、ニコシア
  • ICBGM-2007:第8回ボンドグラフモデリング・シミュレーション国際会議、2007年1月15日~17日、米国カリフォルニア州サンディエゴ
  • ECMS-2006 第20回ヨーロッパモデリング・シミュレーション会議、2006年5月28日~31日、ドイツ、ボン
  • IMAACA-2005 国際地中海モデリング・マルチカンファレンス
  • ICBGM-2005 ボンドグラフモデリングとシミュレーションに関する国際会議、2005年1月23日~27日、米国ルイジアナ州ニューオーリンズ – 論文
  • ICBGM-2003 ボンドグラフモデリングとシミュレーションに関する国際会議 (ICBGM'2003) 2003年1月19日~23日、米国フロリダ州オーランド – 論文
  • 第14回ヨーロッパシミュレーションシンポジウム 2002年10月23日~26日 ドイツ ドレスデン
  • ESS'2001 第13回ヨーロッパシミュレーションシンポジウム 2001年10月18日~20日 フランス マルセイユ
  • ICBGM-2001 ボンドグラフモデリングとシミュレーションに関する国際会議 (ICBGM 2001) 米国アリゾナ州フェニックス
  • ヨーロッパシミュレーションマルチカンファレンス 2000年5月23日~26日 ベルギー ゲント
  • 第11回ヨーロッパシミュレーションシンポジウム 1999年10月26日~28日 ドイツ エアランゲン=ニュルンベルク フリードリヒ=アレクサンダー大学キャッスル
  • ICBGM-1999 ボンドグラフモデリングとシミュレーションに関する国際会議 1999年1月17日~20日 カリフォルニア州サンフランシスコ
  • ESS-97 第9回ヨーロッパシミュレーションシンポジウムおよび展示会 産業におけるシミュレーション、ドイツ、パッサウ、1997年10月19日~22日
  • ICBGM-1997 第3回ボンドグラフモデリングおよびシミュレーション国際会議、1997年1月12日~15日、アリゾナ州フェニックス、シェラトン・クレセントホテル
  • 第11回ヨーロッパシミュレーションマルチカンファレンス トルコ、イスタンブール、1997年6月1日~4日
  • ESM-1996 第10回ヨーロッパシミュレーションマルチカンファレンス ハンガリー、ブダペスト、1996年6月2日~6日
  • ICBGM-1995 ボンドグラフモデリングとシミュレーションに関する国際会議 (ICBGM'95)、1995年1月15~18日、ネバダ州ラスベガス

参照

参考文献

  1. ^ Paynter, Henry M. (1961). Analysis and Design of Engineering Systems . The MIT Press. ISBN 0-262-16004-8 {{cite book}}: ISBN / Date incompatibility (help)
  2. ^ 「工学システムのボンドグラフモデリング」(PDF)
  3. ^ Akbarpour Ghazani, Mehran; Pan, Michael; Tran, Kenneth; Rampadarath, Anand; Nickerson, David P. (2024). 「生物学と生理学におけるボンドグラフの多様な応用のレビュー」Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences . 480 (2294) 20230807. Bibcode : 2024RSPSA.48030807A. doi : 10.1098/rspa.2023.0807 . ISSN  1471-2946
  4. ^ ab オスター、ジョージ、ペレルソン、アラン、カチャルスキー、アハロン (1971). 「ネットワーク熱力学」 . Nature . 234 (5329): 393– 399.書誌コード: 1971Natur.234..393O. doi : 10.1038/234393a0. ISSN  0028-0836
  5. ^ ab Oster, George F.; Perelson, Alan S.; Katchalsky, Aharon (1973). 「ネットワーク熱力学:生体物理システムの動的モデリング」 . Quarterly Reviews of Biophysics . 6 (1): 1– 134. doi :10.1017/S0033583500000081. ISSN  0033-5835.
  6. ^ abcde Gawthrop, Peter J.; Crampin, Edmund J. (2014). 「ボンドグラフを用いた生化学サイクルのエネルギーベース解析」Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences . 470 (2171). arXiv : 1406.2447 . Bibcode :2014RSPSA.47040459G. doi :10.1098/rspa.2014.0459. ISSN  1364-5021. PMC 4197480. PMID 25383030  . 
  7. ^ ab Rajagopal, Vijay; Arumugam, Senthil; Hunter, Peter J.; Khadangi, Afshin; Chung, Joshua; Pan, Michael (2022). 「細胞フィジオーム:単一細胞生物学を予測するための計算生理学フレームワークには何が必要か?」 Annual Review of Biomedical Data Science . 5 (1): 341– 366. doi :10.1146/annurev-biodatasci-072018-021246. ISSN  2574-3414. PMID  35576556
  8. ^ abcdef Gawthrop, Peter J.; Pan, Michael (2022). 「生物システムのネットワーク熱力学:ボンドグラフアプローチ」. Mathematical Biosciences . 352 108899. doi : 10.1016/j.mbs.2022.108899 . PMID  36057321
  9. ^ ab Hunter, P.; de Bono, B.; Brooks, D.; Christie, R.; Hussan, J.; Lin, M.; Nickerson, D. (2025). 「フィジオーム・プロジェクトとデジタルツイン」. IEEE Reviews in Biomedical Engineering . 18 : 300–315 . Bibcode :2025IRBE...18..300H. doi :10.1109/RBME.2024.3490455. ISSN  1937-3333. PMID  39504298
  10. ^ Gawthrop, Peter J.; Cudmore, Peter; Crampin, Edmund J. (2020). 「生体分子システムの物理的に妥当なモデリング:ミトコンドリア電子伝達系の簡略化されたエネルギーベースモデル」Journal of Theoretical Biology . 493 110223. arXiv : 1905.12958 . Bibcode :2020JThBi.49310223G. doi :10.1016/j.jtbi.2020.110223. PMID  32119969
  11. ^ Gawthrop, Peter J.; Cursons, Joseph; Crampin, Edmund J. (2015). 「生化学ネットワークの階層的ボンドグラフモデリング」Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences . 471 (2184) 20150642. arXiv : 1503.01814 . Bibcode :2015RSPSA.47150642G. doi :10.1098/rspa.2015.0642. ISSN  1364-5021.
  12. ^ パン、マイケル;ゴースロップ、ピーター・J.;カーソンズ、ジョセフ;クランプン、エドマンド・J. (2021年10月13日). ウォールクヴィスト、アンダース(編). 「システム生物学における動的モデルのモジュール式アセンブリ」. PLOS Computational Biology . 17 (10) e1009513.書誌コード:2021PLSCB..17E9513P. doi : 10.1371/journal.pcbi.1009513 . ISSN  1553-7358  . PMC 8544865. PMID 34644304 
  13. ^ ピーター・J・ゴースロップ、マイケル・パン、エドマンド・J・クランプン (2021). 「ボンドグラフを用いたモジュラー動的生体分子モデリング:化学量論、熱力学、反応速度論、およびデータの統合」Journal of the Royal Society Interface . 18 (181) 20210478. doi :10.1098/rsif.2021.0478. ISSN  1742-5662 . PMC 8385351. PMID  34428949 
  14. ^ Cudmore, Peter; Pan, Michael; Gawthrop, Peter J.; Crampin, Edmund J. (2021). 「BondGraphToolsを用いた生物学におけるエネルギーベースモデルの分析とシミュレーション」. The European Physical Journal E. 44 ( 12) 148. doi :10.1140/epje/s10189-021-00152-4. ISSN  1292-8941. PMID  34904197.
  15. ^ Keener, James; Sneyd, James 編 (2009). 数理生理学. 学際応用数学. 第8巻第1号. ニューヨーク州ニューヨーク: Springer New York. doi :10.1007/978-0-387-75847-3. ISBN 978-0-387-75846-6
  16. ^ abc カルノップ、ディーン、マーゴリス、ドナルド・L、ローゼンバーグ、ロナルド・C (2012).システムダイナミクス:メカトロニクスシステムのモデリングとシミュレーション(第5版). ホーボーケン、ニュージャージー州: Wiley. ISBN 978-0-470-88908-4
  17. ^ ab Gawthrop, Peter J. (2017). 「化学浸透圧生体分子エネルギー伝達のボンドグラフモデリング」. IEEE Transactions on NanoBioscience . 16 (3): 177– 188. arXiv : 1611.04264 . Bibcode :2017ITNB...16..177G. doi :10.1109/TNB.2017.2674683. ISSN  1536-1241.
  18. ^ Pan, Michael; Gawthrop, Peter J.; Tran, Kenneth; Cursons, Joseph; Crampin, Edmund J. (2019). 「膜トランスポーターモデリングのための熱力学的枠組み」. Journal of Theoretical Biology . 481 : 10–23 . arXiv : 1806.04341 . Bibcode :2019JThBi.481...10P. doi :10.1016/j.jtbi.2018.09.034. PMID:  30273576.
  19. ^ パン, マイケル; ゴースロップ, ピーター J.; トラン, ケネス; カーソンズ, ジョセフ; クランプン, エドマンド J. (2018). 「心臓活動電位のボンドグラフモデリング:ドリフトと非一意な定常状態への影響」. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences . 474 (2214) 20180106. arXiv : 1802.04548 . Bibcode :2018RSPSA.47480106P. doi :10.1098/rspa.2018.0106. ISSN  1364-5021 . PMC 6030650. PMID  29977132 
  20. ^ ab Gawthrop, Peter J.; Pan, Michael; Rajagopal, Vijay (2025). 「ボンドグラフを用いた単一アクチンフィラメント重合のエネルギーベースモデリング」Journal of the Royal Society Interface . 22 (222) 20240404. doi :10.1098/rsif.2024.0404. ISSN  1742-5662 . PMC 11780404. PMID  39881657 
  21. ^ Peskin, CS; Odell, GM; Oster, GF (1993). 「細胞運動と熱変動:ブラウン運動のラチェット」. Biophysical Journal . 65 (1): 316– 324. Bibcode :1993BpJ....65..316P. doi :10.1016/S0006-3495(93)81035-X. PMC 1225726. PMID 8369439.   

参考文献

  • Kypuros, Javier (2013).ボンドグラフモデリングによるシステムダイナミクスと制御. Boca Raton: Taylor&Francis. doi :10.1201/b14676. ISBN 978-1-4665-6075-8
  • ペインター、ヘンリー・M. (1960).工学システムの分析と設計. MITプレス. ISBN 0-262-16004-8 {{cite book}}: ISBN / Date incompatibility (help)
  • ディーン・C・カーノップ、ドナルド・L・マーゴリス、ロナルド・C・ローゼンバーグ (1990)。『システムダイナミクス:統一的アプローチ』。ニューヨーク:ジョン・ワイリー・アンド・サンズ。ISBN   0-471-62171-4
  • ジーン・ウルリッヒ・トーマ (1975)。『ボンドグラフ:入門と応用』 。オックスフォード:パーガモン・プレス。ISBN 0-08-018882-6
  • Gawthrop, Peter J.; Smith, Lorcan PS (1996).メタモデリング:ボンドグラフと動的システム. ロンドン:Prentice Hall. doi :10.5281/zenodo.6998395. ISBN  0-13-489824-9
  • フォーブス・T・ブラウン (2007)。『エンジニアリングシステムダイナミクス:統一的グラフ中心アプローチ』。ボカラトン:テイラー・アンド・。ISBN 978-0-8493-9648-9
  • Mukherjee, Amalendu; Karmakar, Ranjit (2000).ボンドグラフによる工学システムのモデリングとシミュレーション. ボカラトン: CRC Press. ISBN 978-0-8493-0982-3
  • Gawthrop, PJ; Ballance, DJ (1999). 「第2章:階層的ボンドグラフの操作のための記号計算」. Munro, N. (編).制御システム解析および設計における記号的手法. ロンドン: 電気技術者協会. pp. 23-52. ISBN 0-85296-943-0
  • Borutzky, Wolfgang (2010).ボンドグラフ手法. ロンドン: Springer. doi :10.1007/978-1-84882-882-7. ISBN 978-1-84882-881-0
  • http://www.site.uottawa.ca/~rhabash/ESSModelFluid.pdf 流体領域におけるボンドグラフのモデリングについて説明します
  • http://www.dartmouth.edu/~sullivan/22files/Fluid_sys_anal_w_chart.pdf 流体領域における結合グラフのモデリングについて説明します
  • Simscape 公式 MATLAB/Simulink アドオンライブラリ(グラフィカル結合グラフプログラミング用)
  • BG V.2.1 フリーウェア MATLAB/Simulink アドオンライブラリ(グラフィカル結合グラフプログラミング用)
Retrieved from "https://en.wikipedia.org/w/index.php?title=Bond_graph&oldid=1322683375"