応用分野5:多目的最適化

システムには複数の目的を満たすことが求められます。これらの問題には、互いに少しずつ異なるアプローチが必要です。

多目的最適化は、以下のように規定できます:
最小化 f =   [ ψ 01 ( x , b ) , ψ 02 ( x , b ) , , ψ 0 n ( x , b ) ] T MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabeqadiWa ceGabeqabeqadeqadeaakeaaieGaqaaaaaaaaaWdbiaa=zgacqGH9a qpcaWFGcWaamWaa8aabaWdbiaa=H8apaWaaSbaaSqaa8qacaaIWaGa aGymaaWdaeqaaOWdbmaabmaapaqaa8qacaWF4bGaaiilaiaa=jgaai aawIcacaGLPaaacaGGSaGaa8hOaiaa=H8apaWaaSbaaSqaa8qacaaI WaGaaGOmaaWdaeqaaOWdbmaabmaapaqaa8qacaWF4bGaaiilaiaa=j gaaiaawIcacaGLPaaacaGGSaGaa8hOaiabgAci8kaacYcacaWFGcGa a8hYd8aadaWgaaWcbaWdbiaaicdacaqGUbaapaqabaGcpeWaaeWaa8 aabaWdbiaa=HhacaGGSaGaa8NyaaGaayjkaiaawMcaaaGaay5waiaa w2faa8aadaahaaWcbeqaa8qacaWFubaaaaaa@55D3@ (目的関数)
適用条件 ψ i ( x , b )   0 ,   i = 1 ,   ,   p MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabeqadiWa ceGabeqabeqadeqadeaakeaaieGaqaaaaaaaaaWdbiaa=H8apaWaaS baaSqaa8qacaWFPbaapaqabaGcpeWaaeWaa8aabaWdbiaa=HhacaGG SaGaa8NyaaGaayjkaiaawMcaaiaabckacqGHLjYScaaIWaGaaiilai aabckacaWFPbGaeyypa0JaaGymaiaacYcacaqGGcGaeyOjGWRaaiil aiaabckacaWFWbaaaa@4745@ (不等式の制約条件)
  ψ i ( x , b ) = 0 ,   i = p + 1 ,   ,   m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabeqadiWa ceGabeqabeqadeqadeaakeaaieGaqaaaaaaaaaWdbiaa=H8apaWaaS baaSqaa8qacaWFPbaapaqabaGcpeWaaeWaa8aabaWdbiaa=HhacaGG SaGaa8NyaaGaayjkaiaawMcaaiabg2da9iaaicdacaGGSaGaaeiOai aa=LgacqGH9aqpcaWFWbGaey4kaSIaaGymaiaacYcacaqGGcGaeyOj GWRaaiilaiaabckacaWFTbaaaa@4732@

b L   b     b U MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabeqadiWa ceGabeqabeqadeqadeaakeaaieGaqaaaaaaaaaWdbiaa=jgapaWaaS baaSqaa8qacaWFmbaapaqabaGcpeGaaeiOaiabgsMiJkaa=jgacaqG GcGaeyizImQaaeiOaiaa=jgapaWaaSbaaSqaa8qacaWFvbaapaqaba aaaa@3E0E@

(等式の制約条件)

(設計の制限値)

コスト関数と制約関数の各関数 ψ ( x , b ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabeqadiWa ceGabeqabeqadeqadeaakeaaieGaqaaaaaaaaaWdbiaa=H8adaqada WdaeaapeGaa8hEaiaacYcacaWFIbaacaGLOaGaayzkaaaaaa@37A1@ は、次の形式をとると想定されます:

ψ ( x , b ) =   ψ 0 ( x , b ) + T 0 T f L   ( x , b )   d t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeGabeqadiWa ceGabeqabeqadeqadmaakeaaieGaqaaaaaaaaaWdbiaa=H8adaqada WdaeaapeGaa8hEaiaacYcacaWFIbaacaGLOaGaayzkaaGaeyypa0Ja a8hOaiaa=H8apaWaaSbaaSqaa8qacaaIWaaapaqabaGcpeWaaeWaa8 aabaWdbiaa=HhacaGGSaGaa8NyaaGaayjkaiaawMcaaiabgUcaRiaa =bkadaGfWbqabSWdaeaapeGaa8hva8aadaWgaaadbaWdbiaaicdaa8 aabeaaaSqaa8qacaWFubWdamaaBaaameaapeGaa8NzaaWdaeqaaaqd baWdbiabgUIiYdaakiaa=XeacaWFGcWaaeWaa8aabaWdbiaa=Hhaca GGSaGaa8NyaaGaayjkaiaawMcaaiaa=bkacaWFKbGaa8hDaaaa@5260@

MotionSolveでは、コスト関数を次のように再定義することで、多目的問題が単一目的の問題に変換されます:
最小化: f =   j = 0 n w j ψ 0 j ( x , b ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeGabeqadiWa ceGabeqabeqadeqadmaakeaaieGaqaaaaaaaaaWdbiaa=zgacqGH9a qpcaWFGcWaaybCaeqal8aabaWdbiaa=PgacqGH9aqpcaaIWaaapaqa a8qacaWFUbaan8aabaWdbiabggHiLdaakiaa=DhapaWaaSbaaSqaa8 qacaWFQbaapaqabaGcpeGaa8hYd8aadaWgaaWcbaWdbiaaicdaieaa caGFQbaapaqabaGcpeWaaeWaa8aabaWdbiaa=HhacaGGSaGaa8Nyaa GaayjkaiaawMcaaaaa@4576@ (目的関数)
適用条件: ψ i ( x , b ) 0 , i = 1 , , p MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeGabeqadiWa ceGabeqabeqadeqadmaakeaaieGaqaaaaaaaaaWdbiaa=H8apaWaaS baaSqaa8qacaWFPbaapaqabaGcpeWaaeWaa8aabaWdbiaa=HhacaGG SaGaa8NyaaGaayjkaiaawMcaaGqaaiaa+bkacqGHLjYScaaIWaGaai ilaiaa+bkacaWFPbGaeyypa0JaaGymaiaacYcacaGFGcGaeyOjGWRa aiilaiaa+bkacaWFWbaaaa@4746@

ψ i ( x , b ) = 0 , i = p + 1 , , m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeGabeqadiWa ceGabeqabeqadeqadmaakeaaieGaqaaaaaaaaaWdbiaa=H8apaWaaS baaSqaa8qacaWFPbaapaqabaGcpeWaaeWaa8aabaWdbiaa=HhacaGG SaGaa8NyaaGaayjkaiaawMcaaiabg2da9iaaicdacaGGSaacbaGaa4 hOaiaa=LgacqGH9aqpcaWFWbGaey4kaSIaaGymaiaacYcacaGFGcGa eyOjGWRaaiilaiaa+bkacaWFTbaaaa@4736@

b L   b     b U MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeGabeqadiWa ceGabeqabeqadeqadmaakeaaieGaqaaaaaaaaaWdbiaa=jgapaWaaS baaSqaa8qacaWFmbaapaqabaacbaGcpeGaa4hOaiabgsMiJkaa=jga caGFGcGaeyizImQaa4hOaiaa=jgapaWaaSbaaSqaa8qacaWFvbaapa qabaaaaa@3E12@

(不等式の制約条件)

(等式の制約条件)

(設計の制限値)

w j > 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeGabeqadiWa ceGabeqabeqadeqadmaakeaaieGaqaaaaaaaaaWdbiaa=DhapaWaaS baaSqaa8qacaWFQbaapaqabaGcpeGaeyOpa4JaaGimaaaa@363E@ は、さまざまな目的の相対的な重要度を反映した重み付けのセットです。

例1

設計目標
設計目標は、運動学的SLAサスペンションを考慮し、サスペンションのハードポイントを変更して目的のライド挙動を得ることです。


図 1.
ライド挙動の特性は、次の6つの目的で記述できます:
  • 設計位置のライドステア = -4°/m
  • フルジャウンス位置のライドステア = -4°/m
  • 設計位置のライドキャンバー = +5°/m
  • フルジャウンス位置のライドキャンバー = +5°/m
  • 設計位置のライドキャスター = -8°/m
  • フルジャウンス位置のライドキャスター = -8°/m

ライドステア =   s t e e r A n g l e     r i d e H e i g h t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeGabeqadiWa ceGabeqabeqadeqadmaakeaaqaaaaaaaaaWdbmaalaaapaqaa8qacq GHciITieaacaWFGcacbiGaa43Caiaa+rhacaGFLbGaa4xzaiaa+jha caGFbbGaa4NBaiaa+DgacaGFSbGaa4xzaaWdaeaapeGaa8hOaiabgk Gi2kaa=bkacaGFYbGaa4xAaiaa+rgacaGFLbGaa4hsaiaa+vgacaGF PbGaa43zaiaa+HgacaGF0baaaaaa@4AB2@

ライドキャンバー =   s t e e r A n g l e     r i d e H e i g h t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeGabeqadiWa ceGabeqabeqadeqadmaakeaaqaaaaaaaaaWdbmaalaaapaqaa8qacq GHciITieaacaWFGcacbiGaa43Caiaa+rhacaGFLbGaa4xzaiaa+jha caGFbbGaa4NBaiaa+DgacaGFSbGaa4xzaaWdaeaapeGaa8hOaiabgk Gi2kaa=bkacaGFYbGaa4xAaiaa+rgacaGFLbGaa4hsaiaa+vgacaGF PbGaa43zaiaa+HgacaGF0baaaaaa@4AB2@

ライドキャスター =   c a s t e r A n g l e   r i d e H e i g h t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qiqrFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeGabeqadiWa ceGabeqabeqadeqadmaakeaaqaaaaaaaaaWdbmaalaaapaqaa8qacq GHciITieaacaWFGcacbiGaa43yaiaa+fgacaGFZbGaa4hDaiaa+vga caGFYbGaa4xqaiaa+5gacaGFNbGaa4hBaiaa+vgaa8aabaWdbiabgk Gi2kaa=bkacaGFYbGaa4xAaiaa+rgacaGFLbGaa4hsaiaa+vgacaGF PbGaa43zaiaa+HgacaGF0baaaaaa@4A70@

設計変数
サスペンションの中でさまざまなハードポイントの位置を設計変数として使用できます。この問題には、26個の設計変数があります。
アッパーコントロールアームのインボードのハードポイント(前後)XYZ (6個の設計変数)
ロアーコントロールアームのインボードのハードポイント(前後)XYZ (6個の設計変数)
アッパーボールジョイントのハードポイントXYZ (3個の設計変数)
ロアーボールジョイントのハードポイントXYZ (3個の設計変数)
インナータイロッドとアウタータイロッドXYZ (6個の設計変数)
スピンドルのアラインメントポイントの位置XZ (2個の設計変数)
結果
この例では、この最適化方法がきわめて良好に機能します。
応答 ターゲット値 実際の値
設計位置のライドステア -4 °/m -4.107 °/m
フルジャウンス位置のライドステア -4 °/m -3.898 °/m
設計位置のライドキャンバー +5 °/m +4.937 °/m
フルジャウンス位置のライドキャンバー +5 °/m +5.024 °/m
設計位置のライドキャスター -8 °/m -8.402 °/m
フルジャウンス位置のライドキャスター -8 °/m -7.594 °/m

例2

これはバス用とした別のサスペンションの例です。ここで実行するタスクは、必須の基準に合致するようにステアリングシステムを設計することです。適切な設計のステアリングシステムは以下のとおりです:
  • ステアリングホイール角がゼロの場合、左右のホイールのステアもゼロになります。
  • 最大ステアおよび最小ステアでは、左右のホイールのステアリング角の差異がゼロになる必要があります。
  • ステアリング比とは、ステアリングホイールの回転数とホイールの回転数の比率を指します。ステアリングホイール角がゼロの状態のステアリング比は9.27になる必要があります。


図 2.
この問題には次の4つの目的があります:
  • ステアリングホイール角がゼロの状態の左ホイールのステア角 = 0
  • ステアリングホイール角がゼロの状態の右ホイールのステア角 = 0
  • 左右のホイールのステアリング角の差 = 0
  • ステアリングホイール角がゼロの状態のステアリング比 = 9.27
設計変数
  • ドラッグリンクボール位置のX-Y-Z座標
  • ピットマンアームボール位置のX-Y-Z座標
結果
この例では、この最適化方法がきわめて良好に機能します。取得される結果を以下に示します:
応答 ターゲット値 実際の値
SWAゼロの左ステア +0.001°
SWAゼロの右ステア +0.000°
ステアリングの差(左対右) -0.003°
SWAゼロのステアリング比 9.27 9.271