JP5091331B2 - 設計支援装置およびプログラム - Google Patents
設計支援装置およびプログラム Download PDFInfo
- Publication number
- JP5091331B2 JP5091331B2 JP2011004710A JP2011004710A JP5091331B2 JP 5091331 B2 JP5091331 B2 JP 5091331B2 JP 2011004710 A JP2011004710 A JP 2011004710A JP 2011004710 A JP2011004710 A JP 2011004710A JP 5091331 B2 JP5091331 B2 JP 5091331B2
- Authority
- JP
- Japan
- Prior art keywords
- heat resistance
- ratio
- strength
- rigidity
- strength ratio
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/80—Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
- Y02T10/82—Elements for improving aerodynamics
Description
本発明は、製品や部品等の構造物の構造設計を支援する設計支援装置及びプログラムに関し、例えば、構造物の剛性及び耐熱性を向上させるために強化が必要な部位を求めるためのトポロジー最適化処理の入力値として用いる、剛性及び耐熱性の強度割合を算出する設計支援装置及びプログラムに関する。
従来から、コンピュータに有限要素法による構造解析を実行させるプログラムや、当該プログラムを実装した設計支援装置（コンピュータ）が提案されており、構造物の設計において、前記プログラムや前記設計支援装置が広く利用されるようになっている。
例えば、特許文献１には、インストルメントパネル（以下、「インパネ」という）等の構造物に対する構造解析を行い、構造物の弱体部位を抽出したり、構造物における剛性強化が必要な部位を抽出して表示したりする設計支援装置（コンピュータ）が開示されている。
また、例えば、特許文献２には、均質化法によるトポロジー最適化手法により構造物の剛性を最大にする構造物形状を検出する形状解析処理を行う設計支援装置が開示されている。
また、例えば、特許文献２には、均質化法によるトポロジー最適化手法により構造物の剛性を最大にする構造物形状を検出する形状解析処理を行う設計支援装置が開示されている。
そして、上記のような構造解析を実行する機能を備えた設計支援装置（コンピュータ）を利用して、「構造物の剛性及び耐熱性を向上させるために強化が必要な部位（以下、単に「強度必要部位」という）」を求めるトポロジー最適化処理が行われ、トポロジー最適化処理により求めた強度必要部位に板厚を追加したり、リブを配置したりすることが行われている。
また、前記強度必要部位を求めるトポロジー最適化処理では、入力値として「剛性及び耐熱性の強度割合」を与える必要があるため、設計者は、当該トポロジー最適化処理に先だって、「剛性及び耐熱性の強度割合」を求める作業（処理）を行っていた。
ここで、従来技術による「剛性及び耐熱性の強度割合（単に「強度割合」ということもある）」の算出処理の手順について説明する。なお、以下では、構造物がインパネである場合を例にする。
また、前記強度必要部位を求めるトポロジー最適化処理では、入力値として「剛性及び耐熱性の強度割合」を与える必要があるため、設計者は、当該トポロジー最適化処理に先だって、「剛性及び耐熱性の強度割合」を求める作業（処理）を行っていた。
ここで、従来技術による「剛性及び耐熱性の強度割合（単に「強度割合」ということもある）」の算出処理の手順について説明する。なお、以下では、構造物がインパネである場合を例にする。
先ず、設計者は、前記設計支援装置に、インパネ（構造物）の「設計情報」及び「物性値情報、拘束条件、荷重条件、体積密度等の解析条件」を入力する。
次に、前記設計支援装置に、「設計情報」及び「解析条件」を用いて、剛性割合を優位にした強度割合（例えば、剛性を１００％、耐熱性を０％とした強度割合）」で有限要素法によるトポロジー最適化処理を実行させ、インパネの剛性・耐熱性の強度必要部位を算出し、ディスプレイ等の表示手段に、図５（ａ）に示すような、インパネ２００の強度必要部位を示した解析結果２１０を表示させる。
なお、図中では、インパネ２００の各領域が、算出された密度の値毎に色分けされ示されており、所定以上の密度の領域が強度部位２０１であることを表している。
なお、図中では、インパネ２００の各領域が、算出された密度の値毎に色分けされ示されており、所定以上の密度の領域が強度部位２０１であることを表している。
次に、設計者は、上記の「剛性を優位にした強度割合」を用いたトポロジー最適化処理により求めた強度必要部位に追加する板厚を選定し（経験則などから選定し）、前記入力された設計情報に、前記選定した板厚のデータを入力する（強度必要部位の板厚を選定したものに変更する）。
また、上記板厚データの入力を終えると、前記設計支援装置に、有限要素法による剛性ＣＡＥ（Computer Aided Engineering）を実行させ、前記強度必要部位の「たわみ量」を求める。また、前記設計支援装置に、有限要素法による耐熱性ＣＡＥを実行させ、前記強度必要部位の「熱変形量」を求める。
また、上記板厚データの入力を終えると、前記設計支援装置に、有限要素法による剛性ＣＡＥ（Computer Aided Engineering）を実行させ、前記強度必要部位の「たわみ量」を求める。また、前記設計支援装置に、有限要素法による耐熱性ＣＡＥを実行させ、前記強度必要部位の「熱変形量」を求める。
次に、設計者は、図６に示すような「剛性及び耐熱性の強度割合」を横軸（ｘ軸）にとり、「たわみ量」及び「熱変形量」を縦軸（ｙ軸）にとった座標面Ｇを作成する。
また、設計者は、前記座標面Ｇに、「剛性目標の範囲（「たわみ量」がｅ以下）を示す剛性目標線（ｙ＝ｅ）」及び「耐熱性目標の範囲（「熱変形量」がｆ以下）を示す耐熱性目標線（ｙ＝ｆ）」を記載する。
次に、設計者は、前記座標面Ｇにおいて、剛性が１００（耐熱性が０）のｘ座標上に、前記求めた「たわみ量」の値をプロット（図示するＢ１点をプロット）し、剛性が１００（耐熱性が０）のｘ座標上に、前記求めた「熱変形量」の値をプロット（図示するＢ２点をプロット）する。
また、設計者は、前記座標面Ｇに、「剛性目標の範囲（「たわみ量」がｅ以下）を示す剛性目標線（ｙ＝ｅ）」及び「耐熱性目標の範囲（「熱変形量」がｆ以下）を示す耐熱性目標線（ｙ＝ｆ）」を記載する。
次に、設計者は、前記座標面Ｇにおいて、剛性が１００（耐熱性が０）のｘ座標上に、前記求めた「たわみ量」の値をプロット（図示するＢ１点をプロット）し、剛性が１００（耐熱性が０）のｘ座標上に、前記求めた「熱変形量」の値をプロット（図示するＢ２点をプロット）する。
次に、設計者は、前記設計支援装置に、「設計情報」及び「解析条件」を用いた「耐熱性を優位にした強度割合（例えば、耐熱性を１００％、剛性を０％とした強度割合）」を用いた有限要素法によるトポロジー最適化処理を実行させ、インパネの剛性・耐熱性の強度必要部位を算出し、ディスプレイ等の表示手段に、図５（ｂ）に示すような、インパネ２００の強度必要部位を示した解析結果２１１を表示させる。
なお、図中では、インパネ２００の各領域が、算出された密度の値毎に色分けされ示されており、所定以上の密度の領域が強度部位２０２であることを表している。
なお、図中では、インパネ２００の各領域が、算出された密度の値毎に色分けされ示されており、所定以上の密度の領域が強度部位２０２であることを表している。
次に、設計者は、上記「耐熱性を優位にした強度割合」を用いたトポロジー最適化処理で求めた強度必要部位の板厚を選定し（強度必要部位の板厚を選定した板厚に変更し）、上記同様の処理を行う。
具体的には、設計者は、前記設計支援装置に、有限要素法による剛性ＣＡＥを実行させ、前記強度必要部位の「たわみ量」を求める。また、前記設計支援装置に、有限要素法による耐熱性ＣＡＥを実行させ、前記強度必要部位の「熱変形量」を求める。
次に、設計者は、前記座標面Ｇにおいて、剛性が０（耐熱性が１００）のｘ座標上に、前記求めた「たわみ量」の値をプロット（図示するＢ３点をプロット）し、剛性が０（耐熱性が１００）のｘ座標上に、前記求めた「熱変形量」の値をプロット（図示するＢ４点をプロット）する。
そして、設計者は、前記座標面上のＢ１点及びＢ３点を結んだ剛性ライン１１０を作成すると共に、Ｂ２点及びＢ４点を結んだ耐熱性ライン１１２を作成する。
具体的には、設計者は、前記設計支援装置に、有限要素法による剛性ＣＡＥを実行させ、前記強度必要部位の「たわみ量」を求める。また、前記設計支援装置に、有限要素法による耐熱性ＣＡＥを実行させ、前記強度必要部位の「熱変形量」を求める。
次に、設計者は、前記座標面Ｇにおいて、剛性が０（耐熱性が１００）のｘ座標上に、前記求めた「たわみ量」の値をプロット（図示するＢ３点をプロット）し、剛性が０（耐熱性が１００）のｘ座標上に、前記求めた「熱変形量」の値をプロット（図示するＢ４点をプロット）する。
そして、設計者は、前記座標面上のＢ１点及びＢ３点を結んだ剛性ライン１１０を作成すると共に、Ｂ２点及びＢ４点を結んだ耐熱性ライン１１２を作成する。
次に、設計者は、上記座標面Ｇ上で作成した剛性ライン１１０及び剛性目標線（ｙ＝ｅ）に基づいて、剛性ライン１１０上においてｙ座標が「ｅ」以下となるｘ座標の範囲（剛性達成範囲ＲＡ）を特定する。また、耐熱性ライン１１２及び耐熱性目標線（ｙ＝ｆ）に基づいて、耐熱性ライン１１２上においてｙ座標が「ｆ」以下となるｘ座標の範囲（耐熱性達成範囲ＲＢ）を特定する（図６参照）。
そして、設計者は、「剛性達成範囲ＲＡ」及び「耐熱性達成範囲ＲＢ」の両者を満たす「剛性及び耐熱性の両性能達成範囲ＲＡＢ」を特定し、その特定した「両性能達成範囲ＲＡＢ」から強度割合を選定する。
なお、図６に示す例では、「両性能達成範囲ＲＡＢ」の領域が大き過ぎるため、設計者が最適な強度割合（剛性及び耐熱性の強度割合）を選択することが困難であった。
そのため、図示するような大きさの「両性能達成範囲ＲＡＢ」が得られたときには、板厚を変更した上で、再度、上記同様のＣＡＥ処理（剛性ＣＡＥ、耐熱性ＣＡＥ）を行い、「両性能達成範囲ＲＡＢ」を求めていた。すなわち、ある程度の大きさの「両性能達成範囲ＲＡＢ」が得られるまで、前記板厚変更及び前記ＣＡＥ処理を繰り返し行っていた。
なお、図６に示す例では、「両性能達成範囲ＲＡＢ」の領域が大き過ぎるため、設計者が最適な強度割合（剛性及び耐熱性の強度割合）を選択することが困難であった。
そのため、図示するような大きさの「両性能達成範囲ＲＡＢ」が得られたときには、板厚を変更した上で、再度、上記同様のＣＡＥ処理（剛性ＣＡＥ、耐熱性ＣＡＥ）を行い、「両性能達成範囲ＲＡＢ」を求めていた。すなわち、ある程度の大きさの「両性能達成範囲ＲＡＢ」が得られるまで、前記板厚変更及び前記ＣＡＥ処理を繰り返し行っていた。
また、図７に示すように、強度必要部位に設定する板厚によっては、「剛性達成範囲ＲＡ」及び「耐熱性達成範囲ＲＢ」の両者を満たす領域（剛性及び耐熱性の両性能達成範囲ＲＡＢ）が存在しないことがある。このような場合、設計者は、「両性能達成範囲ＲＡＢ」が得られるまで、前記板厚変更及び前記ＣＡＥ処理を繰り返し行っていた。
しかしながら、上述した従来技術による「剛性及び耐熱性の強度割合」の算出方法は、以下に示す技術的課題を有している。
具体的には、上述した従来技術の前記強度割合の算出方法では、設計者自身が「板厚変更及びＣＡＥ処理（剛性ＣＡＥ、耐熱性ＣＡＥ）」の操作を複数回繰り返して行わないと、ある程度の大きさに絞り込んだ「剛性及び耐熱性の両性能達成範囲ＲＡＢ（図６参照）」を得ることができなかった。そのため、前記強度割合の算出方法では、工数が大きく、設計作業に長時間かかっていた。
すなわち、前記強度割合の算出方法は、工数が大きく面倒であると共に、設計に長時間費やされるという技術的課題を有していた。
具体的には、上述した従来技術の前記強度割合の算出方法では、設計者自身が「板厚変更及びＣＡＥ処理（剛性ＣＡＥ、耐熱性ＣＡＥ）」の操作を複数回繰り返して行わないと、ある程度の大きさに絞り込んだ「剛性及び耐熱性の両性能達成範囲ＲＡＢ（図６参照）」を得ることができなかった。そのため、前記強度割合の算出方法では、工数が大きく、設計作業に長時間かかっていた。
すなわち、前記強度割合の算出方法は、工数が大きく面倒であると共に、設計に長時間費やされるという技術的課題を有していた。
本発明は、上記の技術的課題を解決するためになされたものであり、本発明の目的は、トポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合の算出工数を低減し、当該強度割合の算出時間を短縮させる設計支援装置を提供することにある。
また、本発明の目的は、コンピュータに、トポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合を算出する処理を実行させるプログラムであって、当該強度割合の算出工数を低減し、当該強度割合の算出時間を短縮させるプログラムを提供することにある。
また、本発明の目的は、コンピュータに、トポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合を算出する処理を実行させるプログラムであって、当該強度割合の算出工数を低減し、当該強度割合の算出時間を短縮させるプログラムを提供することにある。
上記課題を解決するために本発明は、構造物のトポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合を算出する設計支援装置に適用される。
そして、前記設計支援装置は、前記構造物の設計情報及び解析条件を用いて、剛性割合を優位にした強度割合で前記構造物の第１トポロジー最適化処理を行うと共に耐熱性割合を優位にした強度割合で前記構造物の第２トポロジー最適化処理を行う構造解析部を備え、前記構造解析部は、前記第１トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第１たわみ量及び第１熱変形量を算出すると共に、前記第２トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第２たわみ量及び第２熱変形量を算出する解析処理と、前記第１たわみ量及び前記第２たわみ量を用いて一次関数からなる剛性関数を算出し、該剛性関数及び所定剛性目標値を用いて、該剛性目標を満たし且つ剛性割合が最小となる第１強度割合を算出すると共に、前記第１熱変形量及び第２熱変形量を用いて一次関数からなる耐熱性関数を算出し、該耐熱性関数及び所定耐熱性目標値を用いて、耐熱性目標を満たし且つ耐熱性割合が最小となる第２強度割合を算出する強度割合算出処理とを実行し、さらに、剛性及び耐熱性の強度割合をｘ軸とし、たわみ量及び熱変形量をｙ軸としたｘｙ直交座標上の前記剛性関数及び前記耐熱性関数において、傾きをそれぞれ固定し、ｙ切片をそれぞれ変数として、当該両関数をｙ軸に沿ってシフトし、前記第１強度割合の剛性割合よりも前記第２強度割合が大きく、且つ、前記第１強度割合と前記第２強度割合との差分が所定値以下となる条件を満たす第２の剛性関数と第２の耐熱性関数とを算出し、前記第２の剛性関数と第２の耐熱性関数とから前記構造物の強度必要部位の第２の板厚を算出し、前記第２の板厚に基づき前記解析処理と前記強度割合算出処理とを再実行することを特徴としている。
そして、前記設計支援装置は、前記構造物の設計情報及び解析条件を用いて、剛性割合を優位にした強度割合で前記構造物の第１トポロジー最適化処理を行うと共に耐熱性割合を優位にした強度割合で前記構造物の第２トポロジー最適化処理を行う構造解析部を備え、前記構造解析部は、前記第１トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第１たわみ量及び第１熱変形量を算出すると共に、前記第２トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第２たわみ量及び第２熱変形量を算出する解析処理と、前記第１たわみ量及び前記第２たわみ量を用いて一次関数からなる剛性関数を算出し、該剛性関数及び所定剛性目標値を用いて、該剛性目標を満たし且つ剛性割合が最小となる第１強度割合を算出すると共に、前記第１熱変形量及び第２熱変形量を用いて一次関数からなる耐熱性関数を算出し、該耐熱性関数及び所定耐熱性目標値を用いて、耐熱性目標を満たし且つ耐熱性割合が最小となる第２強度割合を算出する強度割合算出処理とを実行し、さらに、剛性及び耐熱性の強度割合をｘ軸とし、たわみ量及び熱変形量をｙ軸としたｘｙ直交座標上の前記剛性関数及び前記耐熱性関数において、傾きをそれぞれ固定し、ｙ切片をそれぞれ変数として、当該両関数をｙ軸に沿ってシフトし、前記第１強度割合の剛性割合よりも前記第２強度割合が大きく、且つ、前記第１強度割合と前記第２強度割合との差分が所定値以下となる条件を満たす第２の剛性関数と第２の耐熱性関数とを算出し、前記第２の剛性関数と第２の耐熱性関数とから前記構造物の強度必要部位の第２の板厚を算出し、前記第２の板厚に基づき前記解析処理と前記強度割合算出処理とを再実行することを特徴としている。
上記のように構成したのは、本願発明者らが、剛性関数（上述した図６の剛性ライン１１０に相当）から求めた第１強度割合と、耐熱性関数（上述した図６の耐熱性ライン１１２に相当）から求めた第２強度割合との差分が所定値以下の場合の第１強度割合（又は第２強度割合）を入力値にしたトポロジー最適化処理を実行した場合に、剛性及び耐熱性を向上させるための強度必要部位が高精度に抽出されることを見出したためである。
そして、本発明によれば、設計支援装置の構造解析部が、前記トポロジー最適化処理の入力値として用いる好適な強度割合を算出するため、上述した従来技術と比べて、設計者の工数が削減されると共に、設計者の手間が大幅に軽減される。
特に、強度必要部位の板厚初期値に基づき得られた剛性関数と耐熱性関数の傾きを固定し、ｙ切片を変数として当該両関数をｙ軸に沿ってシフトすることにより、最終的に求めたい値に少なくとも近似する第２の板厚、及び剛性と耐熱性との強度割合の値を容易に得ることができる。したがって、強度割合（剛性割合）を算出するために費やされる設計時間を大幅に短縮することができる。
そして、本発明によれば、設計支援装置の構造解析部が、前記トポロジー最適化処理の入力値として用いる好適な強度割合を算出するため、上述した従来技術と比べて、設計者の工数が削減されると共に、設計者の手間が大幅に軽減される。
特に、強度必要部位の板厚初期値に基づき得られた剛性関数と耐熱性関数の傾きを固定し、ｙ切片を変数として当該両関数をｙ軸に沿ってシフトすることにより、最終的に求めたい値に少なくとも近似する第２の板厚、及び剛性と耐熱性との強度割合の値を容易に得ることができる。したがって、強度割合（剛性割合）を算出するために費やされる設計時間を大幅に短縮することができる。
また、前記第２の板厚に基づき実行された前記解析処理と前記強度割合算出処理の結果、前記構造解析部は、算出された前記第１強度割合と前記第２強度割合の剛性割合とを比較し、該第１強度割合の剛性割合の方が大きい値の場合には、前記第２強度割合の方が大きい値となるまで、前記第２の板厚を所定値増加させ、該板厚に基づく前記解析処理と、前記強度割合算出処理とを繰り返し実行することが望ましい。
また、前記第１強度割合と前記第２強度割合の剛性割合との比較処理の結果、該第２強度割合の剛性割合の方が大きい値の場合に、前記構造解析部は、前記第１強度割合と前記第２強度割合の差分を算出し、該差分が所定値より大きい場合には、該差分が所定値以下となるまで、前記第２の板厚を所定値減少させ、該板厚に基づく前記解析処理と、前記強度割合算出処理とを繰り返し実行することが望ましい。
この構成により、確実に剛性及び耐熱性の両性能達成範囲が存在し、且つ、その範囲が所定値以下となる強度割合を自動的に短時間で確実に算出することが可能となる。
また、前記第１強度割合と前記第２強度割合の剛性割合との比較処理の結果、該第２強度割合の剛性割合の方が大きい値の場合に、前記構造解析部は、前記第１強度割合と前記第２強度割合の差分を算出し、該差分が所定値より大きい場合には、該差分が所定値以下となるまで、前記第２の板厚を所定値減少させ、該板厚に基づく前記解析処理と、前記強度割合算出処理とを繰り返し実行することが望ましい。
この構成により、確実に剛性及び耐熱性の両性能達成範囲が存在し、且つ、その範囲が所定値以下となる強度割合を自動的に短時間で確実に算出することが可能となる。
また、上記課題を解決するためになされた本発明は、コンピュータに、構造物のトポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合を算出する処理を実行させるプログラムに適用される。
そして、前記プログラムは、前記構造物の設計情報及び解析条件を用いて、剛性割合を優位にした強度割合で前記構造物の第１トポロジー最適化処理を行うと共に耐熱性割合を優位にした強度割合で前記構造物の第２トポロジー最適化処理を行うステップと、前記第１トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第１たわみ量及び第１熱変形量を算出すると共に、前記第２トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第２たわみ量及び第２熱変形量を算出する解析処理を行うステップと、前記第１たわみ量及び前記第２たわみ量を用いて一次関数からなる剛性関数を算出し、該剛性関数及び所定剛性目標値を用いて、該剛性目標を満たし且つ剛性割合が最小となる第１強度割合を算出すると共に、前記第１熱変形量及び第２熱変形量を用いて一次関数からなる耐熱性関数を算出し、該耐熱性関数及び所定耐熱性目標値を用いて、耐熱性目標を満たし且つ耐熱性割合が最小となる第２強度割合を算出する強度割合算出処理を行うステップと、剛性及び耐熱性の強度割合をｘ軸とし、たわみ量及び熱変形量をｙ軸としたｘｙ直交座標上の前記剛性関数及び前記耐熱性関数において、傾きをそれぞれ固定し、ｙ切片をそれぞれ変数として、当該両関数をｙ軸に沿ってシフトし、前記第１強度割合の剛性割合よりも前記第２強度割合が大きく、且つ、前記第１強度割合と前記第２強度割合との差分が所定値以下となる条件を満たす第２の剛性関数と第２の耐熱性関数とを算出するステップと、前記第２の剛性関数と第２の耐熱性関数とから前記構造物の強度必要部位の第２の板厚を算出し、前記第２の板厚に基づき前記解析処理と前記強度割合算出処理とを再実行するステップとを有することを特徴としている。
そして、前記プログラムは、前記構造物の設計情報及び解析条件を用いて、剛性割合を優位にした強度割合で前記構造物の第１トポロジー最適化処理を行うと共に耐熱性割合を優位にした強度割合で前記構造物の第２トポロジー最適化処理を行うステップと、前記第１トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第１たわみ量及び第１熱変形量を算出すると共に、前記第２トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第２たわみ量及び第２熱変形量を算出する解析処理を行うステップと、前記第１たわみ量及び前記第２たわみ量を用いて一次関数からなる剛性関数を算出し、該剛性関数及び所定剛性目標値を用いて、該剛性目標を満たし且つ剛性割合が最小となる第１強度割合を算出すると共に、前記第１熱変形量及び第２熱変形量を用いて一次関数からなる耐熱性関数を算出し、該耐熱性関数及び所定耐熱性目標値を用いて、耐熱性目標を満たし且つ耐熱性割合が最小となる第２強度割合を算出する強度割合算出処理を行うステップと、剛性及び耐熱性の強度割合をｘ軸とし、たわみ量及び熱変形量をｙ軸としたｘｙ直交座標上の前記剛性関数及び前記耐熱性関数において、傾きをそれぞれ固定し、ｙ切片をそれぞれ変数として、当該両関数をｙ軸に沿ってシフトし、前記第１強度割合の剛性割合よりも前記第２強度割合が大きく、且つ、前記第１強度割合と前記第２強度割合との差分が所定値以下となる条件を満たす第２の剛性関数と第２の耐熱性関数とを算出するステップと、前記第２の剛性関数と第２の耐熱性関数とから前記構造物の強度必要部位の第２の板厚を算出し、前記第２の板厚に基づき前記解析処理と前記強度割合算出処理とを再実行するステップとを有することを特徴としている。
本発明によれば、トポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合の算出工数を低減し、当該強度割合の算出時間を短縮させる設計支援装置を提供することができる。
また、本発明によれば、コンピュータに、トポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合を算出する処理を実行させるプログラムであって、当該強度割合の算出工数を低減し、当該強度割合の算出時間を短縮させるプログラムを提供することができる。
また、本発明によれば、コンピュータに、トポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合を算出する処理を実行させるプログラムであって、当該強度割合の算出工数を低減し、当該強度割合の算出時間を短縮させるプログラムを提供することができる。
以下、本発明の実施形態について、図面に基づいて説明する。
先ず、本実施形態の設計支援装置の機能構成を図１に基づいて説明する。
図１は、本発明の実施形態の設計支援装置の機能ブロック図である。
先ず、本実施形態の設計支援装置の機能構成を図１に基づいて説明する。
図１は、本発明の実施形態の設計支援装置の機能ブロック図である。
図示するように、本実施形態の設計支援装置Ｗは、トポロジー最適化処理の入力値として用いる「剛性及び耐熱性の強度割合」を算出する情報処理装置１と、設計者からの各種要求を受け付ける入力装置２と、情報処理装置１が行った解析結果を出力する出力装置３とを備えている。また、情報処理装置１は、ＬＡＮ（Local Area Network）等のネットワークＮＷを介して、ＣＡＤ装置４に接続されている。
ここで、入力装置２は、キーボードやマウス等により構成され、設計者からの各種要求や解析条件（物性値情報、拘束条件、荷重条件、体積密度等）等を受け付け情報処理装置１に出力する。
出力装置３は、液晶ディスプレイ等により構成され、情報処理装置１が出力する画像情報を表示する。
また、ＣＡＤ装置４には、構造解析を行う対象の構造物のＣＡＤ情報（例えば、自動車の構成部品の設計情報）が格納されている。そして、ＣＡＤ装置４は、情報処理装置１からの要求にしたがい、情報処理装置１にＣＡＤ情報を出力する。なお、本実施形態のＣＡＤ装置４は、公知の技術により実現されるため、詳細な説明を省略する。
出力装置３は、液晶ディスプレイ等により構成され、情報処理装置１が出力する画像情報を表示する。
また、ＣＡＤ装置４には、構造解析を行う対象の構造物のＣＡＤ情報（例えば、自動車の構成部品の設計情報）が格納されている。そして、ＣＡＤ装置４は、情報処理装置１からの要求にしたがい、情報処理装置１にＣＡＤ情報を出力する。なお、本実施形態のＣＡＤ装置４は、公知の技術により実現されるため、詳細な説明を省略する。
また、情報処理装置１は、制御部１０、データ取得部２０、構造解析部３０および出力部４０を備えている。
制御部１０は、情報処理装置１の全体の動作を制御する。また、制御部１０は、入力装置２を介して、設計者が入力する各種要求を受け付ける。そして、制御部１０は、上記の受け付けた要求にしたがい、データ取得部２０、構造解析部３０、および出力部４０を制御して、設計者からの要求に応じた各種の処理を行う。
制御部１０は、情報処理装置１の全体の動作を制御する。また、制御部１０は、入力装置２を介して、設計者が入力する各種要求を受け付ける。そして、制御部１０は、上記の受け付けた要求にしたがい、データ取得部２０、構造解析部３０、および出力部４０を制御して、設計者からの要求に応じた各種の処理を行う。
また、データ取得部２０は、ネットワークＮＷに接続されている外部装置（例えば、ＣＡＤ装置４）と通信を行い、外部装置との間でデータの授受を行う。例えば、データ取得部２０は、ネットワークＮＷを介して、ＣＡＤ装置４にアクセスし、ＣＡＤ装置４に格納されている設計情報（ＣＡＤ情報）を取得する。
また、データ取得部２０は、入力装置２を介して、設計者が入力する解析対象の構造物の解析条件の入力を受け付ける。
また、データ取得部２０は、入力装置２を介して、設計者が入力する解析対象の構造物の解析条件の入力を受け付ける。
構造解析部３０は、ＣＡＤ装置４から取得した構造物の「設計情報」および「解析条件」を用いて、後述する図３に示す各処理ステップを実行し、構造物の剛性及び耐熱性を向上させるために強化が必要な部位（強度必要部位）を求めるためのトポロジー最適化処理の入力値として用いる、剛性及び耐熱性の強度割合を算出する。
また、出力部４０は、構造解析部３０から解析結果を取得し、その解析結果を示す画像情報を生成し、出力装置３に、その生成した画像情報を出力する。
また、出力部４０は、構造解析部３０から解析結果を取得し、その解析結果を示す画像情報を生成し、出力装置３に、その生成した画像情報を出力する。
つぎに、本実施形態の情報処理装置１のハードウェア構成を説明する。
図２は、本実施形態の情報処理装置のハードウェア構成図である。
図示するように、情報処理装置１は、ＣＰＵ（Central Processing Unit）５０と、ＲＡＭ（Random Access Memory）等により構成された主記憶装置５１と、Ｉ／Ｏインタフェース５２と、ハードディスク等により構成された補助記憶装置５３と、ネットワークＮＷに接続されている装置との間で行うデータ授受の制御を行うネットワークインタフェース５４とを有する。
また、補助記憶装置５３には、上述した各部（制御部１０、データ取得部２０、構造解析部３０、および出力部４０）の機能を実現するためのプログラム（設計支援プログラム５５）が格納されている。
図２は、本実施形態の情報処理装置のハードウェア構成図である。
図示するように、情報処理装置１は、ＣＰＵ（Central Processing Unit）５０と、ＲＡＭ（Random Access Memory）等により構成された主記憶装置５１と、Ｉ／Ｏインタフェース５２と、ハードディスク等により構成された補助記憶装置５３と、ネットワークＮＷに接続されている装置との間で行うデータ授受の制御を行うネットワークインタフェース５４とを有する。
また、補助記憶装置５３には、上述した各部（制御部１０、データ取得部２０、構造解析部３０、および出力部４０）の機能を実現するためのプログラム（設計支援プログラム５５）が格納されている。
そして、情報処理装置１の各部（制御部１０、データ取得部２０、構造解析部３０、および出力部４０）の機能は、ＣＰＵ５０が補助記憶装置５３に格納されている前記プログラムを主記憶装置５１にロードして実行することにより実現される。
続いて、本実施形態の設計支援装置Ｗが行う構造物の剛性及び耐熱性を向上させる強度必要部位を求めるトポロジー最適化処理の入力値として用いる、剛性及び耐熱性の強度割合を算出する処理について図３及び図４に基づいて説明する。
ここで、図３は、本発明の実施形態の設計支援装置が行うトポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合の算出処理の手順を示したフローチャートである。また、図４は、本発明の実施形態の剛性及び耐熱性の強度割合の算出処理において用いる剛性ラインを示す剛性関数、耐熱性ラインを示す耐熱性関数、剛性目標線、及び耐熱線目標線を示した模式図である。
なお、以下の説明では、解析対象の構造物が自動車のインパネ（合成樹脂製のインパネ）である場合を例にする。また、以下の処理ステップの中で行われるトポロジー最適化処理及びＣＡＥ処理（剛性ＣＡＥ、耐熱性ＣＡＥ）は、周知のものと同様であるため詳細な説明を省略する。
ここで、図３は、本発明の実施形態の設計支援装置が行うトポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合の算出処理の手順を示したフローチャートである。また、図４は、本発明の実施形態の剛性及び耐熱性の強度割合の算出処理において用いる剛性ラインを示す剛性関数、耐熱性ラインを示す耐熱性関数、剛性目標線、及び耐熱線目標線を示した模式図である。
なお、以下の説明では、解析対象の構造物が自動車のインパネ（合成樹脂製のインパネ）である場合を例にする。また、以下の処理ステップの中で行われるトポロジー最適化処理及びＣＡＥ処理（剛性ＣＡＥ、耐熱性ＣＡＥ）は、周知のものと同様であるため詳細な説明を省略する。
先ず、情報処理装置１のデータ取得部２０が、解析対象のデータの読み込みを行う（Ｓ１）。
具体的には、データ取得部２０は、ネットワークＮＷを介して、ＣＡＤ装置４にアクセスし、ＣＡＤ装置４に格納されているインパネの設計情報（ＣＡＤ情報）を取得し、情報処理装置１のメモリ（主記憶装置５１又は補助記憶装置５３）に、前記取得したインパネの設計情報を格納する。
また、データ取得部２０は、設計者が入力装置２を介して入力する「インパネの解析条件」を受け付け（受信し）、前記メモリ（主記憶装置５１又は補助記憶装置５３）に、前記解析条件を格納する。
具体的には、データ取得部２０は、ネットワークＮＷを介して、ＣＡＤ装置４にアクセスし、ＣＡＤ装置４に格納されているインパネの設計情報（ＣＡＤ情報）を取得し、情報処理装置１のメモリ（主記憶装置５１又は補助記憶装置５３）に、前記取得したインパネの設計情報を格納する。
また、データ取得部２０は、設計者が入力装置２を介して入力する「インパネの解析条件」を受け付け（受信し）、前記メモリ（主記憶装置５１又は補助記憶装置５３）に、前記解析条件を格納する。
そして、Ｓ１の処理が終わると、構造解析部３０は、「Ｓ２〜Ｓ４の処理」及び「Ｓ５〜Ｓ７の処理」を行う。
なお、構造解析部３０は、「Ｓ２〜Ｓ４の処理」及び「Ｓ５〜Ｓ７の処理」を並列処理で実行するように構成されていてもよいし、「Ｓ２〜Ｓ４（又はＳ５〜Ｓ７）の処理」を終えた後に「Ｓ５〜Ｓ７（又はＳ２〜Ｓ４）の処理」を実行するように構成されていてもよい。
なお、以下では、説明の便宜上、「Ｓ２〜Ｓ４の処理」を説明してから「Ｓ５〜Ｓ７の処理」を説明する。
なお、構造解析部３０は、「Ｓ２〜Ｓ４の処理」及び「Ｓ５〜Ｓ７の処理」を並列処理で実行するように構成されていてもよいし、「Ｓ２〜Ｓ４（又はＳ５〜Ｓ７）の処理」を終えた後に「Ｓ５〜Ｓ７（又はＳ２〜Ｓ４）の処理」を実行するように構成されていてもよい。
なお、以下では、説明の便宜上、「Ｓ２〜Ｓ４の処理」を説明してから「Ｓ５〜Ｓ７の処理」を説明する。
具体的には、Ｓ２では、構造解析部３０は、上述したＳ１により前記メモリに格納された「インパネの設計情報」および「インパネの解析条件」を用いて、剛性割合を優位にした強度割合で構造物（インパネ）のトポロジー最適化処理（第１トポロジー最適化処理）を行い、インパネ（構造物）の強度必要部位を抽出する。
なお、本実施形態では、剛性割合を優位にした強度割合に、「剛性を１００％とし、耐熱性を０％とした強度割合」を用いる。また、本実施形態では、１箇所の強度必要部位が抽出される場合を例にする。
なお、本実施形態では、剛性割合を優位にした強度割合に、「剛性を１００％とし、耐熱性を０％とした強度割合」を用いる。また、本実施形態では、１箇所の強度必要部位が抽出される場合を例にする。
また、本ステップ（Ｓ２）で用いられる「剛性割合を優位にした強度割合」は、デフォルト値として、構造解析部３０に予め設定されている。
なお、構造解析部３０は、デフォルト値として設定された「剛性割合を優位にした強度割合」を利用しないで、設計者が入力する任意の「剛性割合を優位にした強度割合」を用いたトポロジー最適化処理も実行できるようになっている。この場合、構造解析部３０は、設計者が入力する「剛性割合を優位にした強度割合」を示すデータを受け付け（データ取得部２０を介して受け付け）、その受け付けた「剛性割合を優位にした強度割合」を用いてトポロジー最適化処理を行う。
また、本ステップ（Ｓ２）において、ディスプレイ等の表示手段に、上述した図５に示すような、インパネ２００の強度必要部位を示した解析結果２１０を表示させるようにしてもよい。
なお、構造解析部３０は、デフォルト値として設定された「剛性割合を優位にした強度割合」を利用しないで、設計者が入力する任意の「剛性割合を優位にした強度割合」を用いたトポロジー最適化処理も実行できるようになっている。この場合、構造解析部３０は、設計者が入力する「剛性割合を優位にした強度割合」を示すデータを受け付け（データ取得部２０を介して受け付け）、その受け付けた「剛性割合を優位にした強度割合」を用いてトポロジー最適化処理を行う。
また、本ステップ（Ｓ２）において、ディスプレイ等の表示手段に、上述した図５に示すような、インパネ２００の強度必要部位を示した解析結果２１０を表示させるようにしてもよい。
次に、構造解析部３０は、Ｓ２において求めた強度必要部位に対する板厚初期値設定を行う（Ｓ３）。
具体的には、構造解析部３０は、Ｓ２において求めた強度必要部位に対して、密度の値に応じて定められた板厚を設定する（変更する）。すなわち、構造解析部３０は、前記設計情報に含まれるインパネの板厚のうち、強度必要部位の板厚を前記トポロジー最適化処理により求められた密度に応じて設定する。
例えば、構造解析部３０は、密度が「ｒ１」の範囲の強度必要部位に「ｔ＝Ｗ０（ｍｍ）」の厚さの板厚を初期値として設定する。
具体的には、構造解析部３０は、Ｓ２において求めた強度必要部位に対して、密度の値に応じて定められた板厚を設定する（変更する）。すなわち、構造解析部３０は、前記設計情報に含まれるインパネの板厚のうち、強度必要部位の板厚を前記トポロジー最適化処理により求められた密度に応じて設定する。
例えば、構造解析部３０は、密度が「ｒ１」の範囲の強度必要部位に「ｔ＝Ｗ０（ｍｍ）」の厚さの板厚を初期値として設定する。
次に、構造解析部３０は、剛性ＣＡＥ（剛性解析）及び耐熱性ＣＡＥ（耐熱性解析）を実行して、板厚設定した強度必要部位の「たわみ量」及び「熱変形量」を算出する（Ｓ４）。
具体的には、構造解析部３０は、Ｓ３において板厚設定した「インパネの設計情報」および「インパネの解析条件」を用いて、有限要素法による剛性ＣＡＥを実行し、前記強度必要部位の「たわみ量（説明の便宜上、「第１たわみ量」という）」を求める。
また、構造解析部３０は、Ｓ３において板厚設定した「インパネの設計情報」および「インパネの解析条件」を用いて、有限要素法による耐熱性ＣＡＥを実行し、前記強度必要部位の「熱変形量（説明の便宜上、「第１熱変形量」という）」を求める。
具体的には、構造解析部３０は、Ｓ３において板厚設定した「インパネの設計情報」および「インパネの解析条件」を用いて、有限要素法による剛性ＣＡＥを実行し、前記強度必要部位の「たわみ量（説明の便宜上、「第１たわみ量」という）」を求める。
また、構造解析部３０は、Ｓ３において板厚設定した「インパネの設計情報」および「インパネの解析条件」を用いて、有限要素法による耐熱性ＣＡＥを実行し、前記強度必要部位の「熱変形量（説明の便宜上、「第１熱変形量」という）」を求める。
次に、構造解析部３０は、Ｓ５〜Ｓ７の処理を行う。
具体的には、Ｓ５では、構造解析部３０は、Ｓ１で読み込んだ「インパネの設計情報」および「インパネの解析条件」を用いて、耐熱性割合を優位にした強度割合で構造物（インパネ）のトポロジー最適化処理（第２トポロジー最適化処理）を行い、インパネ（構造物）の強度必要部位を抽出する。
なお、本実施形態では、耐熱性割合を優位にした強度割合に、「剛性を０％とし、耐熱性を１００％とした強度割合」を用いる。
また、Ｓ５の処理は、上述したＳ２の処理と、強度割合が異なる以外は同じであるため、詳細な説明を省略する。
具体的には、Ｓ５では、構造解析部３０は、Ｓ１で読み込んだ「インパネの設計情報」および「インパネの解析条件」を用いて、耐熱性割合を優位にした強度割合で構造物（インパネ）のトポロジー最適化処理（第２トポロジー最適化処理）を行い、インパネ（構造物）の強度必要部位を抽出する。
なお、本実施形態では、耐熱性割合を優位にした強度割合に、「剛性を０％とし、耐熱性を１００％とした強度割合」を用いる。
また、Ｓ５の処理は、上述したＳ２の処理と、強度割合が異なる以外は同じであるため、詳細な説明を省略する。
次に、構造解析部３０は、上述したＳ３と同様の手順により、Ｓ７において求めた強度必要部位に対する板厚設定を行う（Ｓ６）。
次に、構造解析部３０は、上述したＳ４と同様の手順により、剛性ＣＡＥ及び耐熱性ＣＡＥを実行して、板厚変更した強度必要部位の「たわみ量」及び「熱変形量」を算出する（Ｓ７）。
具体的には、構造解析部３０は、Ｓ６において板厚設定した「インパネの設計情報」および「インパネの解析条件」を用いて、有限要素法による剛性ＣＡＥを実行し、前記強度必要部位の「たわみ量（説明の便宜上、「第２たわみ量」という）」を求める。
また、構造解析部３０は、Ｓ６において板厚設定した「インパネの設計情報」および「インパネの解析条件」を用いて、有限要素法による耐熱性ＣＡＥを実行し、前記強度必要部位の「熱変形量（説明の便宜上、「第２熱変形量」という）」を求める。
具体的には、構造解析部３０は、Ｓ６において板厚設定した「インパネの設計情報」および「インパネの解析条件」を用いて、有限要素法による剛性ＣＡＥを実行し、前記強度必要部位の「たわみ量（説明の便宜上、「第２たわみ量」という）」を求める。
また、構造解析部３０は、Ｓ６において板厚設定した「インパネの設計情報」および「インパネの解析条件」を用いて、有限要素法による耐熱性ＣＡＥを実行し、前記強度必要部位の「熱変形量（説明の便宜上、「第２熱変形量」という）」を求める。
次に、構造解析部３０は、Ｓ４及びＳ７で求めた「たわみ量（第１たわみ量、第２たわみ量）」及び「熱変形量（第２熱変形量、第２熱変形量）」を用いて、「剛性ラインを示す剛性関数（ｙ＝ａｘ＋ｂ）」及び「耐熱性ラインを示す耐熱性関数（ｙ＝ｃｘ＋ｄ）」を算出する（Ｓ８）。
具体的には、Ｓ８では、構造解析部３０は、図４に示すような「剛性及び耐熱性の強度割合（図中では剛性割合の値を用いている）」を横軸（ｘ軸）にとり、「たわみ量」及び「熱変形量」を縦軸（ｙ軸）にとったｘｙ直交座標面Ｇを作成する。
ここで、前記座標面Ｇは、左端ＧＬのｘ座標が「１００（剛性が１００％、耐熱性が０％）」を示し、右端ＧＲのｘ座標が「０（剛性が０％、耐熱性が１００％）」を示している。
すなわち、図示する座標面Ｇは、左端ＧＬから右端ＧＲに向けてｘ１方向に進むにつれて、ｘ軸の値が減少するようになっている（剛性の割合が減少して耐熱性の割合が増加するようになっている）。
ここで、前記座標面Ｇは、左端ＧＬのｘ座標が「１００（剛性が１００％、耐熱性が０％）」を示し、右端ＧＲのｘ座標が「０（剛性が０％、耐熱性が１００％）」を示している。
すなわち、図示する座標面Ｇは、左端ＧＬから右端ＧＲに向けてｘ１方向に進むにつれて、ｘ軸の値が減少するようになっている（剛性の割合が減少して耐熱性の割合が増加するようになっている）。
そして、Ｓ８では、構造解析部３０は、前記座標面Ｇにおいて、剛性が１００（耐熱性が０）のｘ座標上（左端ＧＬのｘ座標上）に、Ｓ４で求めた「第１たわみ量」の値をプロット（Ａ１）し、剛性が０（耐熱性が１００）のｘ座標上（右端ＧＲのｘ座標上）に、Ｓ６で求めた「第２たわみ量」の値をプロット（Ａ３）する。また、構造解析部３０は、前記プロットした２点（Ａ１、Ａ３）を用いて、当該２点（Ａ１、Ａ３）を通る剛性ラインを示す剛性関数（ｙ＝ａｘ＋ｂ）を算出する。尚、前記剛性ラインは、実際には図示するように僅かに曲線となるが、この剛性関数の算出においては、近似式として一次関数が算出される。
ここで、上記「Ａ１」のｘ座標の値が「１００」となっているのは、Ｓ４の剛性ＣＡＥにより得られた「第１たわみ量」が、Ｓ２のトポロジー最適化処理により得られた強度必要部位に対して実施されたものであるためである。
また、上記「Ａ３」のｘ座標の値が「０」となっているのは、Ｓ６の剛性ＣＡＥにより得られた「第２たわみ量」が、Ｓ５のトポロジー最適化処理により得られた強度必要部位に対して実施されたものであるためである。
すなわち、Ｓ２及びＳ５において実施するトポロジー最適化処理における強度割合により、「Ａ１」及び「Ａ３」のｘ座標の値は変わってくる。
また、上記「Ａ３」のｘ座標の値が「０」となっているのは、Ｓ６の剛性ＣＡＥにより得られた「第２たわみ量」が、Ｓ５のトポロジー最適化処理により得られた強度必要部位に対して実施されたものであるためである。
すなわち、Ｓ２及びＳ５において実施するトポロジー最適化処理における強度割合により、「Ａ１」及び「Ａ３」のｘ座標の値は変わってくる。
さらに、構造解析部３０は、Ｓ８において、前記座標面Ｇにおいて、剛性が１００（耐熱性が０）のｘ座標上（左端ＧＬのｘ座標上）に、Ｓ４で求めた「第１熱変形量」の値をプロット（Ａ２）し、剛性が０（耐熱性が１００）のｘ座標上（右端ＧＲのｘ座標上）に、Ｓ６で求めた「第２熱変形量」の値をプロット（Ａ４）する。また、構造解析部３０は、前記プロットした２点（Ａ２、Ａ４）を用いて、当該２点（Ａ２、Ａ４）を通る耐熱性ラインを示す耐熱性関数（ｙ＝ｃｘ＋ｄ）を算出する。尚、前記剛性ラインは、実際には図示するように僅かに曲線となるが、この剛性関数の算出においては、近似式として一次関数が算出される。
ここで、上記「Ａ２」のｘ座標の値が「１００」となっているのは、Ｓ４の耐熱性ＣＡＥにより得られた「第１熱変形量」が、Ｓ２のトポロジー最適化処理により得られた強度必要部位に対して実施されたものであるためである。
また、上記「Ａ４」のｘ座標の値が「０」となっているのは、Ｓ６の耐熱性ＣＡＥにより得られた「第２熱変形量」が、Ｓ５のトポロジー最適化処理により得られた強度必要部位に対して実施されたものであるためである。
ここで、上記「Ａ２」のｘ座標の値が「１００」となっているのは、Ｓ４の耐熱性ＣＡＥにより得られた「第１熱変形量」が、Ｓ２のトポロジー最適化処理により得られた強度必要部位に対して実施されたものであるためである。
また、上記「Ａ４」のｘ座標の値が「０」となっているのは、Ｓ６の耐熱性ＣＡＥにより得られた「第２熱変形量」が、Ｓ５のトポロジー最適化処理により得られた強度必要部位に対して実施されたものであるためである。
次に、構造解析部３０は、上記剛性関数（ｙ＝ａｘ＋ｂ）を用いて第１強度割合ｍを算出すると共に、上記耐熱性関数（ｙ＝ｃｘ＋ｄ）を用いて第２強度割合ｎを算出する（Ｓ９）。
具体的には、Ｓ９では、構造解析部３０は、「剛性目標の範囲（「たわみ量」がｅ以下）を示す剛性目標線（ｙ＝ｅ）と、剛性ラインを示す剛性関数（ｙ＝ａｘ＋ｂ）との交点（（ｅ−ｂ）／ａ，ｅ）を算出して、その交点のｘ座標（ｘ＝（ｅ−ｂ）／ａ）の値を「第１強度割合（ｍ）」とする。
すなわち、ここでは、剛性関数（ｙ＝ａｘ＋ｂ）上において、上記の剛性目標の要件を満たし且つ剛性割合が最小となるＸ座標の値である「第１強度割合（ｍ）」が算出される。
なお、「ｍ（第１強度割合）」は、剛性割合が「ｍ」であり、耐熱性割合が「１００−ｍ」であることを示している。
すなわち、ここでは、剛性関数（ｙ＝ａｘ＋ｂ）上において、上記の剛性目標の要件を満たし且つ剛性割合が最小となるＸ座標の値である「第１強度割合（ｍ）」が算出される。
なお、「ｍ（第１強度割合）」は、剛性割合が「ｍ」であり、耐熱性割合が「１００−ｍ」であることを示している。
さらに、構造解析部３０は、Ｓ９において、耐熱性目標の範囲（「熱変形量」がｆ以下）を示す耐熱性目標線（ｙ＝ｆ）と、耐熱性関数（ｙ＝ｃｘ＋ｄ）との交点（（ｆ−ｄ）／ｃ，ｆ）を算出して、その交点のｘ座標（ｘ＝（ｆ−ｄ）／ｃ）の値を「第２強度割合（ｎ）」とする。
すなわち、耐熱性関数（ｙ＝ｃｘ＋ｄ）上において、上記の耐熱性目標の要件を満たし且つ耐熱性割合が最小となるＸ座標の値である「第２強度割合（ｎ）」が算出される。
なお、「ｎ（第２強度割合）」は、剛性割合が「ｎ」であり、耐熱性割合が「１００−ｎ」であることを示している。
すなわち、耐熱性関数（ｙ＝ｃｘ＋ｄ）上において、上記の耐熱性目標の要件を満たし且つ耐熱性割合が最小となるＸ座標の値である「第２強度割合（ｎ）」が算出される。
なお、「ｎ（第２強度割合）」は、剛性割合が「ｎ」であり、耐熱性割合が「１００−ｎ」であることを示している。
次に、構造解析部３０は、Ｓ９で求めた「ｍ」と「ｎ」とを比較し、「ｎ＞ｍ」、且つ「ｎ−ｍ≦閾値」の条件を満たす場合には、出力部３０を介して、出力装置３に「ｎ」及び「ｍ」の値を示す情報（例えば、画像情報）を出力して処理を終了する。また、いずれか一方の条件でも満たさない場合、Ｓ１１の処理に進む（ステップＳ１０）。
なお、「ｎ＞ｍ」の要件を満たす場合とは、図４に示すように「剛性及び耐熱性の両性能達成範囲ＲＡＢ」が存在する場合を示し、「ｎ＞ｍ」の要件を満たさない場合とは、上述した図７に示すような「剛性及び耐熱性の両性能達成範囲ＲＡＢ」が存在しないケースを示している。
また、ｎ−ｍ≦閾値の場合とは、「ｎ−ｍ」すなわち「剛性及び耐熱性の両性能達成範囲ＲＡＢ」の大きさが、「閾値（所定値）」以下であるケースを示している。
なお、「閾値（所定値）」の値は、ある程度小さい値（例えば、１未満の値）に設定されていることが望ましい。
このようにＳ１０では、剛性及び耐熱性の両性能達成範囲ＲＡＢの有無の判定を行い、さらに前記両性能達成範囲ＲＡＢを有する場合、「ｎ−ｍ」の値が閾値以下であるか否かの判定を行っている。
なお、「ｎ＞ｍ」の要件を満たす場合とは、図４に示すように「剛性及び耐熱性の両性能達成範囲ＲＡＢ」が存在する場合を示し、「ｎ＞ｍ」の要件を満たさない場合とは、上述した図７に示すような「剛性及び耐熱性の両性能達成範囲ＲＡＢ」が存在しないケースを示している。
また、ｎ−ｍ≦閾値の場合とは、「ｎ−ｍ」すなわち「剛性及び耐熱性の両性能達成範囲ＲＡＢ」の大きさが、「閾値（所定値）」以下であるケースを示している。
なお、「閾値（所定値）」の値は、ある程度小さい値（例えば、１未満の値）に設定されていることが望ましい。
このようにＳ１０では、剛性及び耐熱性の両性能達成範囲ＲＡＢの有無の判定を行い、さらに前記両性能達成範囲ＲＡＢを有する場合、「ｎ−ｍ」の値が閾値以下であるか否かの判定を行っている。
また、このように構成するのは、本願発明者らが研究を重ねた結果、本処理ステップにより算出される差分Δ（ｎ−ｍ）が小さい場合の「ｎ（或いは「ｍ」）」の値が、前記トポロジー最適化処理の入力値として用いる強度割合（剛性割合）」に好適であることを見出したことによる。
すなわち、前記差分Δ（ｎ−ｍ）が小さい場合の「ｎ（或いは「ｍ」）」の値を剛性割合にした強度割合を入力値にした前記トポロジー最適化処理を実行した場合に、剛性及び耐熱性を向上させるための強度必要部位を正確に抽出できることを見出したためである。
すなわち、前記差分Δ（ｎ−ｍ）が小さい場合の「ｎ（或いは「ｍ」）」の値を剛性割合にした強度割合を入力値にした前記トポロジー最適化処理を実行した場合に、剛性及び耐熱性を向上させるための強度必要部位を正確に抽出できることを見出したためである。
具体的には、構造解析部３０は、「ｎ−ｍ≦閾値」であれば、上術したように、出力装置３に「ｎ」及び「ｍ」の値を示す情報（例えば、画像情報）を出力する。
設計者は、出力装置３に出力された「ｎ」及び「ｍ」の値により、剛性及び耐熱性の強度割合を簡単に選定することができる。例えば、「ｎ」が「４５．８」であり、「ｍ」が「４５．１」である場合には、設計者は、上記「ｎ」及び「ｍ」より、剛性割合を「４５％」と選定する。この場合、耐熱性割合は、「５５％」に定まる。
なお、本ステップにおいて、構造解析部３０が、「ｎ」及び「ｍ」の値から剛性割合を算出し（例えば、ｎ及びｍの平均値を求めて小数点を四捨五入した値を剛性割合とする）、その剛性割合に対応する耐熱性割合を算出し、当該算出した剛性割合及び耐熱性割合を出力するように構成されていてもよい。
設計者は、出力装置３に出力された「ｎ」及び「ｍ」の値により、剛性及び耐熱性の強度割合を簡単に選定することができる。例えば、「ｎ」が「４５．８」であり、「ｍ」が「４５．１」である場合には、設計者は、上記「ｎ」及び「ｍ」より、剛性割合を「４５％」と選定する。この場合、耐熱性割合は、「５５％」に定まる。
なお、本ステップにおいて、構造解析部３０が、「ｎ」及び「ｍ」の値から剛性割合を算出し（例えば、ｎ及びｍの平均値を求めて小数点を四捨五入した値を剛性割合とする）、その剛性割合に対応する耐熱性割合を算出し、当該算出した剛性割合及び耐熱性割合を出力するように構成されていてもよい。
次に、Ｓ１０の条件を満たさないと判定された場合に進むＳ１１以降の処理を説明する。
Ｓ１０の条件を満たさないとする判定が初回の場合、Ｓ２、Ｓ５からの処理に戻るが、新たに板厚を設定するためにステップＳ１２へ進む（ステップＳ１１）。一方、Ｓ１０の条件を満たさないとする判定が２回目以上の場合はＳ１４の処理に進む（ステップＳ１１）。
具体的には、Ｓ１０の条件を満たさないとする判定が初回の場合、構造解析部３０は、上記Ｓ８で求められた剛性関数（ｙ＝ａｘ＋ｂ）、及び耐熱性関数（ｙ＝ｃｘ＋ｄ）において、それぞれｂ、ｄを変数として、Ｓ１０の条件、即ち「ｎ＞ｍ」且つ「ｎ−ｍ≦閾値」を満たすｂ＝ｂ１、ｄ＝ｄ１を設定する（ステップＳ１２）。
Ｓ１０の条件を満たさないとする判定が初回の場合、Ｓ２、Ｓ５からの処理に戻るが、新たに板厚を設定するためにステップＳ１２へ進む（ステップＳ１１）。一方、Ｓ１０の条件を満たさないとする判定が２回目以上の場合はＳ１４の処理に進む（ステップＳ１１）。
具体的には、Ｓ１０の条件を満たさないとする判定が初回の場合、構造解析部３０は、上記Ｓ８で求められた剛性関数（ｙ＝ａｘ＋ｂ）、及び耐熱性関数（ｙ＝ｃｘ＋ｄ）において、それぞれｂ、ｄを変数として、Ｓ１０の条件、即ち「ｎ＞ｍ」且つ「ｎ−ｍ≦閾値」を満たすｂ＝ｂ１、ｄ＝ｄ１を設定する（ステップＳ１２）。
例えば、Ｓ１０の条件のうち、「ｎ＞ｍ」の条件を満たさない場合、すなわち、図７に示したように剛性及び耐熱性の両性能達成範囲ＲＡＢが存在しない場合、剛性関数（ｙ＝ａｘ＋ｂ）と耐熱性関数（ｙ＝ｃｘ＋ｄ）の傾きａ、ｃをそれぞれ固定し、剛性ライン及び耐熱性ラインをＹ軸に沿って下方にシフトする。これにより、ｙ切片（ｂ，ｄ）の変化に伴い上記ｎ，ｍの値が変化し、「ｎ＞ｍ」、且つ「ｎ−ｍ≦閾値」を満たすｂ＝ｂ１，ｄ＝ｄ１を求めることができる。
尚、「ｎ−ｍ≦閾値」の条件を満たさない場合、すなわち図６に示したように剛性及び耐熱性の両性能達成範囲ＲＡＢが広すぎる場合にも、同様に剛性関数（ｙ＝ａｘ＋ｂ）と耐熱性関数（ｙ＝ｃｘ＋ｄ）の傾きａ、ｃをそれぞれ固定し、剛性ライン及び耐熱性ラインをＹ軸に沿って上方にシフトすることによって、Ｓ１０の条件を満たすｂ＝ｂ１，ｄ＝ｄ１を求めることができる。
尚、「ｎ−ｍ≦閾値」の条件を満たさない場合、すなわち図６に示したように剛性及び耐熱性の両性能達成範囲ＲＡＢが広すぎる場合にも、同様に剛性関数（ｙ＝ａｘ＋ｂ）と耐熱性関数（ｙ＝ｃｘ＋ｄ）の傾きａ、ｃをそれぞれ固定し、剛性ライン及び耐熱性ラインをＹ軸に沿って上方にシフトすることによって、Ｓ１０の条件を満たすｂ＝ｂ１，ｄ＝ｄ１を求めることができる。
さらに、構造解析部３０は、Ｓ４，Ｓ７での剛性ＣＡＥ及び耐熱性ＣＡＥに基づき、第２の構造関数（ｙ＝ａｘ＋ｂ１）、及び第２の耐熱性関数（ｙ＝ｃｘ＋ｄ１）から、第２の板厚Ｗ１を逆算処理により求める（ステップＳ１３）。
構造解析部３０は、上記求めた第２の板厚Ｗ１を新たな板厚として設定し、上述したＳ４、Ｓ７の処理に戻る。
このようにＳ１０での条件を満たさない初回時に、Ｓ１２及びＳ１３の処理を実行するのは、これらの処理により、最終的に求めたい値に少なくとも近似する（剛性と耐熱性との）強度割合の値を得ることができるためである。
詳しく説明すると、本願発明者らが研究を重ねた結果、板厚を変化させて剛性ＣＡＥ及び耐熱性ＣＡＥをそれぞれ実施した場合、剛性ライン及び耐熱性ラインの傾きの変化よりも、たわみ量或いは熱変形量が大きく変化することを見出した。そのため、剛性ライン及び耐熱性ラインにおいて、変形量（すなわちｙ切片「ｂ」「ｄ」）を変化させ、所望のｘ＝ｍ、ｘ＝ｎとしたときの板厚を求めれば、それが設計上最も好ましい板厚に近似することを見出すに至った。
構造解析部３０は、上記求めた第２の板厚Ｗ１を新たな板厚として設定し、上述したＳ４、Ｓ７の処理に戻る。
このようにＳ１０での条件を満たさない初回時に、Ｓ１２及びＳ１３の処理を実行するのは、これらの処理により、最終的に求めたい値に少なくとも近似する（剛性と耐熱性との）強度割合の値を得ることができるためである。
詳しく説明すると、本願発明者らが研究を重ねた結果、板厚を変化させて剛性ＣＡＥ及び耐熱性ＣＡＥをそれぞれ実施した場合、剛性ライン及び耐熱性ラインの傾きの変化よりも、たわみ量或いは熱変形量が大きく変化することを見出した。そのため、剛性ライン及び耐熱性ラインにおいて、変形量（すなわちｙ切片「ｂ」「ｄ」）を変化させ、所望のｘ＝ｍ、ｘ＝ｎとしたときの板厚を求めれば、それが設計上最も好ましい板厚に近似することを見出すに至った。
上記のように新たな板厚Ｗ１によりＳ４，Ｓ７でのＣＡＥが実施されると、改めて構造関数（ｙ＝ａｘ＋ｂ）、及び耐熱性関数（ｙ＝ｃｘ＋ｄ）が算出される（ステップＳ８）。
また、Ｓ９において第１強度割合（ｍ）及び第２強度割合（ｎ）が算出され、再びＳ１０において、「ｎ＞ｍ」且つ「ｎ−ｍ≦閾値」の条件を満たしているか否かが判定される。
ここで、上記Ｓ１０の条件を満たす場合には、出力部３０を介して、出力装置３に「ｎ」及び「ｍ」の値を示す情報（例えば、画像情報）を出力して処理を終了する。また、いずれか一方の条件でも満たさない場合、Ｓ１１の処理に進む。
また、Ｓ９において第１強度割合（ｍ）及び第２強度割合（ｎ）が算出され、再びＳ１０において、「ｎ＞ｍ」且つ「ｎ−ｍ≦閾値」の条件を満たしているか否かが判定される。
ここで、上記Ｓ１０の条件を満たす場合には、出力部３０を介して、出力装置３に「ｎ」及び「ｍ」の値を示す情報（例えば、画像情報）を出力して処理を終了する。また、いずれか一方の条件でも満たさない場合、Ｓ１１の処理に進む。
ここで、Ｓ１０の条件を満たさないとする判定は２回目であるため、Ｓ１４での処理に進む（ステップＳ１４）。
Ｓ１４では、「ｎ＞ｍ」であるか、すなわち「剛性及び耐熱性の両性能達成範囲ＲＡＢ」が存在するか否かが判定される。
ここで、「ｎ＞ｍ」、すなわち「剛性及び耐熱性の両性能達成範囲ＲＡＢ」が存在しない場合、Ｓ１３において板厚変更した「インパネの設計情報」を用いて、当該板厚変更したインパネの強度必要部位に対して所定の板厚（Ｔ１）を加算する（ステップＳ１６）。そして、Ｓ１６の処理を終えると上述したＳ４、Ｓ７の処理に戻る。
Ｓ１４では、「ｎ＞ｍ」であるか、すなわち「剛性及び耐熱性の両性能達成範囲ＲＡＢ」が存在するか否かが判定される。
ここで、「ｎ＞ｍ」、すなわち「剛性及び耐熱性の両性能達成範囲ＲＡＢ」が存在しない場合、Ｓ１３において板厚変更した「インパネの設計情報」を用いて、当該板厚変更したインパネの強度必要部位に対して所定の板厚（Ｔ１）を加算する（ステップＳ１６）。そして、Ｓ１６の処理を終えると上述したＳ４、Ｓ７の処理に戻る。
そして、Ｓ４，Ｓ７でのＣＡＥの処理が終わると、再び、Ｓ８〜Ｓ１０の処理を行い、「ｎ＞ｍ」且つ「ｎ−ｍ≦閾値」の条件を満たしているか否かが判定される。
ここで、「ｎ＞ｍ」の条件を満たさない場合には、当該条件を満たすまで、Ｓ１１，Ｓ１４，Ｓ１６，Ｓ４，Ｓ７，及びＳ８〜Ｓ１０の処理を繰り返す。
ここで、「ｎ＞ｍ」の条件を満たさない場合には、当該条件を満たすまで、Ｓ１１，Ｓ１４，Ｓ１６，Ｓ４，Ｓ７，及びＳ８〜Ｓ１０の処理を繰り返す。
次に、Ｓ１４において「剛性及び耐熱性の両性能達成範囲ＲＡＢ」があると判定された場合（「ｎ＞ｍ」を満たす場合）に進むＳ１５の処理を説明する。
Ｓ１５では、構造解析部３０は、Ｓ９で改めて算出された「ｍ」及び「ｎ」から、差分Δ（ｎ−ｍ）が「閾値（所定値）」以下であるか否かを判定する（ｎ−ｍ≦閾値）。
Ｓ１５では、構造解析部３０は、Ｓ９で改めて算出された「ｍ」及び「ｎ」から、差分Δ（ｎ−ｍ）が「閾値（所定値）」以下であるか否かを判定する（ｎ−ｍ≦閾値）。
そして、構造解析部３０は、「ｎ−ｍ≦閾値」であれば、出力部３０を介して、出力装置３に「ｎ」及び「ｍ」の値を示す情報（例えば、画像情報）を出力して処理を終了し、「ｎ−ｍ≦閾値」を満たさなければ、Ｓ１７の処理に進む。
Ｓ１７において、「ｎ−ｍ≦閾値」を満たしていないと判定した場合、構造解析部３０は、前段のＳ１４の判定処理時の所定板厚を所定の厚さＴ２だけ減算した上で、上述したＳ４、Ｓ７の処理に戻る。
そして、Ｓ４，Ｓ７でのＣＡＥの処理が終わると、再び、Ｓ８〜Ｓ１０の処理を行い、「ｎ＞ｍ」且つ「ｎ−ｍ≦閾値」の条件を満たしているか否かが判定される。
ここで、「ｎ−ｍ≦閾値」の条件を満たさない場合には、当該条件を満たすまで、Ｓ１１，Ｓ１４，Ｓ１５，Ｓ１７，Ｓ４，Ｓ７，及びＳ８−Ｓ１０の処理を繰り返す。
以上の構成により、構造解析部３０が、トポロジー最適化処理の入力値として好適な強度割合（剛性割合）を自動的に比較的短時間で導き出すことが可能になる。
Ｓ１７において、「ｎ−ｍ≦閾値」を満たしていないと判定した場合、構造解析部３０は、前段のＳ１４の判定処理時の所定板厚を所定の厚さＴ２だけ減算した上で、上述したＳ４、Ｓ７の処理に戻る。
そして、Ｓ４，Ｓ７でのＣＡＥの処理が終わると、再び、Ｓ８〜Ｓ１０の処理を行い、「ｎ＞ｍ」且つ「ｎ−ｍ≦閾値」の条件を満たしているか否かが判定される。
ここで、「ｎ−ｍ≦閾値」の条件を満たさない場合には、当該条件を満たすまで、Ｓ１１，Ｓ１４，Ｓ１５，Ｓ１７，Ｓ４，Ｓ７，及びＳ８−Ｓ１０の処理を繰り返す。
以上の構成により、構造解析部３０が、トポロジー最適化処理の入力値として好適な強度割合（剛性割合）を自動的に比較的短時間で導き出すことが可能になる。
このように、本実施形態によれば、設計支援装置Ｗの構造解析部３０が、トポロジー最適化処理の入力値として用いる好適な強度割合（剛性割合）を自動的に算出するため、上述した従来技術と比べて、設計者の工数が削減されると共に、設計者の手間が大幅に軽減される。
特に、強度必要部位の板厚初期値に基づき得られた剛性関数（剛性ライン）と耐熱性関数（耐熱性ライン）の傾きを固定し、ｙ切片を変数として当該両関数をｙ軸に沿ってシフトすることにより、最終的に求めたい値に少なくとも近似する板厚、及び剛性と耐熱性との強度割合の値を容易に得ることができる。したがって、強度割合を算出するために費やされる設計時間を大幅に短縮することができる。
特に、強度必要部位の板厚初期値に基づき得られた剛性関数（剛性ライン）と耐熱性関数（耐熱性ライン）の傾きを固定し、ｙ切片を変数として当該両関数をｙ軸に沿ってシフトすることにより、最終的に求めたい値に少なくとも近似する板厚、及び剛性と耐熱性との強度割合の値を容易に得ることができる。したがって、強度割合を算出するために費やされる設計時間を大幅に短縮することができる。
また、本実施形態によれば、設計支援装置Ｗの構造解析部３０が、前記トポロジー最適化処理の入力値として用いる強度割合（剛性割合）を算出するため、設計者の習熟度に関係なく、効果的なトポロジー最適化処理が実現される。
なお、本発明は、上述した実施形態に限定されるものではなく、その要旨の範囲内において種々の変更が可能である。
なお、本発明は、上述した実施形態に限定されるものではなく、その要旨の範囲内において種々の変更が可能である。
例えば、上述した実施形態では、構造物として自動車のインパネを例にしたがあくまでもこれは一例に過ぎない。樹脂製の製品や部品の設計であれば、どのようなものにも適用することができる。
Ｗ 設計支援装置
１ 情報処理装置
１０ 制御部（情報処理装置）
２０ データ取得部（情報処理装置）
３０ 構造解析部（情報処理装置）
４０ 出力部（情報処理装置）
２ 入力装置
３ 出力装置
４ ＣＡＤ装置
５０ ＣＰＵ
５１ 主記憶装置
５２ Ｉ／Ｏインタフェース
５３ 補助記憶装置
５４ ＮＷインタフェース
５５ 設計支援プログラム
２００ インストルメントパネル
２０１、２０２ 強度必要部位（インストルメントパネル）
１ 情報処理装置
１０ 制御部（情報処理装置）
２０ データ取得部（情報処理装置）
３０ 構造解析部（情報処理装置）
４０ 出力部（情報処理装置）
２ 入力装置
３ 出力装置
４ ＣＡＤ装置
５０ ＣＰＵ
５１ 主記憶装置
５２ Ｉ／Ｏインタフェース
５３ 補助記憶装置
５４ ＮＷインタフェース
５５ 設計支援プログラム
２００ インストルメントパネル
２０１、２０２ 強度必要部位（インストルメントパネル）
Claims (4)
- 構造物のトポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合を算出する設計支援装置であって、
前記構造物の設計情報及び解析条件を用いて、剛性割合を優位にした強度割合で前記構造物の第１トポロジー最適化処理を行うと共に耐熱性割合を優位にした強度割合で前記構造物の第２トポロジー最適化処理を行う構造解析部を備え、
前記構造解析部は、
前記第１トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第１たわみ量及び第１熱変形量を算出すると共に、前記第２トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第２たわみ量及び第２熱変形量を算出する解析処理と、
前記第１たわみ量及び前記第２たわみ量を用いて一次関数からなる剛性関数を算出し、該剛性関数及び所定剛性目標値を用いて、該剛性目標を満たし且つ剛性割合が最小となる第１強度割合を算出すると共に、前記第１熱変形量及び第２熱変形量を用いて一次関数からなる耐熱性関数を算出し、該耐熱性関数及び所定耐熱性目標値を用いて、耐熱性目標を満たし且つ耐熱性割合が最小となる第２強度割合を算出する強度割合算出処理とを実行し、
さらに、剛性及び耐熱性の強度割合をｘ軸とし、たわみ量及び熱変形量をｙ軸としたｘｙ直交座標上の前記剛性関数及び前記耐熱性関数において、傾きをそれぞれ固定し、ｙ切片をそれぞれ変数として、当該両関数をｙ軸に沿ってシフトし、前記第１強度割合の剛性割合よりも前記第２強度割合が大きく、且つ、前記第１強度割合と前記第２強度割合との差分が所定値以下となる条件を満たす第２の剛性関数と第２の耐熱性関数とを算出し、
前記第２の剛性関数と第２の耐熱性関数とから前記構造物の強度必要部位の第２の板厚を算出し、前記第２の板厚に基づき前記解析処理と前記強度割合算出処理とを再実行することを特徴とする設計支援装置。 - 前記第２の板厚に基づき実行された前記解析処理と前記強度割合算出処理の結果、
前記構造解析部は、算出された前記第１強度割合と前記第２強度割合の剛性割合とを比較し、該第１強度割合の剛性割合の方が大きい値の場合には、前記第２強度割合の方が大きい値となるまで、前記第２の板厚を所定値加算し、該板厚に基づく前記解析処理と、前記強度割合算出処理とを繰り返し実行することを特徴とする請求項１に記載された設計支援装置。 - 前記第１強度割合と前記第２強度割合の剛性割合との比較処理の結果、該第２強度割合の剛性割合の方が大きい値の場合に、
前記構造解析部は、前記第１強度割合と前記第２強度割合の差分を算出し、該差分が所定値より大きい場合には、該差分が所定値以下となるまで、前記板厚を所定値減算し、該板厚に基づく前記解析処理と、前記強度割合算出処理とを繰り返し実行することを特徴とする請求項１または請求項２に記載された設計支援装置。 - コンピュータに、構造物のトポロジー最適化処理の入力値として用いる剛性及び耐熱性の強度割合を算出する処理を実行させるプログラムであって、
前記構造物の設計情報及び解析条件を用いて、剛性割合を優位にした強度割合で前記構造物の第１トポロジー最適化処理を行うと共に耐熱性割合を優位にした強度割合で前記構造物の第２トポロジー最適化処理を行うステップと、
前記第１トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第１たわみ量及び第１熱変形量を算出すると共に、前記第２トポロジー最適化処理で抽出された前記構造物の強度必要部位の板厚を変更した上で剛性解析および耐熱性解析を実施して該強度必要部位の第２たわみ量及び第２熱変形量を算出する解析処理を行うステップと、
前記第１たわみ量及び前記第２たわみ量を用いて一次関数からなる剛性関数を算出し、該剛性関数及び所定剛性目標値を用いて、該剛性目標を満たし且つ剛性割合が最小となる第１強度割合を算出すると共に、前記第１熱変形量及び第２熱変形量を用いて一次関数からなる耐熱性関数を算出し、該耐熱性関数及び所定耐熱性目標値を用いて、耐熱性目標を満たし且つ耐熱性割合が最小となる第２強度割合を算出する強度割合算出処理を行うステップと、
剛性及び耐熱性の強度割合をｘ軸とし、たわみ量及び熱変形量をｙ軸としたｘｙ直交座標上の前記剛性関数及び前記耐熱性関数において、傾きをそれぞれ固定し、ｙ切片をそれぞれ変数として、当該両関数をｙ軸に沿ってシフトし、前記第１強度割合の剛性割合よりも前記第２強度割合が大きく、且つ、前記第１強度割合と前記第２強度割合との差分が所定値以下となる条件を満たす第２の剛性関数と第２の耐熱性関数とを算出するステップと、
前記第２の剛性関数と第２の耐熱性関数とから前記構造物の強度必要部位の第２の板厚を算出し、前記第２の板厚に基づき前記解析処理と前記強度割合算出処理とを再実行するステップとを有することを特徴とするプログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011004710A JP5091331B2 (ja) | 2011-01-13 | 2011-01-13 | 設計支援装置およびプログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011004710A JP5091331B2 (ja) | 2011-01-13 | 2011-01-13 | 設計支援装置およびプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2012146172A JP2012146172A (ja) | 2012-08-02 |
| JP5091331B2 true JP5091331B2 (ja) | 2012-12-05 |
Family
ID=46789669
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2011004710A Expired - Fee Related JP5091331B2 (ja) | 2011-01-13 | 2011-01-13 | 設計支援装置およびプログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5091331B2 (ja) |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003320530A (ja) * | 2002-05-08 | 2003-11-11 | Toyota Motor Corp | 複層樹脂材の熱変形予測方法およびそのプログラム |
| JP4590377B2 (ja) * | 2006-07-14 | 2010-12-01 | 関東自動車工業株式会社 | 樹脂成形品の剛性最適化解析プログラム及びその記録媒体 |
| JP4744541B2 (ja) * | 2008-03-11 | 2011-08-10 | 関東自動車工業株式会社 | 設計支援装置 |
-
2011
- 2011-01-13 JP JP2011004710A patent/JP5091331B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2012146172A (ja) | 2012-08-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Amir et al. | Achieving stress-constrained topological design via length scale control | |
| US10346554B2 (en) | System and method for isogeometric analysis | |
| JP7446579B2 (ja) | 曲面のフィッティング処理方法、フィッティング処理装置およびフィッティング処理プログラム、並びに、該フィッティング処理プログラムを記憶したコンピュータ読取可能な記憶媒体 | |
| Sehmi et al. | Review of topology optimisation refinement processes for sheet metal manufacturing in the automotive industry | |
| JP2018106521A (ja) | 解析装置、解析方法、及びコンピュータプログラム | |
| JP5403194B1 (ja) | 解析装置、解析方法、及びコンピュータプログラム | |
| Davies et al. | On the utility of asymptotic bandwidth selectors for spatially adaptive kernel density estimation | |
| JP4744541B2 (ja) | 設計支援装置 | |
| JP4994475B2 (ja) | 設計支援装置およびプログラム | |
| JP5091331B2 (ja) | 設計支援装置およびプログラム | |
| JP4879244B2 (ja) | 部品剛性の最適構造作製プログラムおよび情報処理装置 | |
| JP2016198997A (ja) | コンピュータ支援の樹脂挙動解析装置 | |
| JP5059153B2 (ja) | 設計支援装置およびプログラム | |
| Debray et al. | Parametric design and optimization of addendum surfaces for sheet metal forming process | |
| JP2019003324A (ja) | Ｃａｅ解析支援装置 | |
| JP4528962B2 (ja) | 設計支援方法 | |
| JP2005301349A (ja) | 解析モデルの作成支援装置 | |
| Lee et al. | Shape optimization of mobile phone folder module for structural strength | |
| JP2001229407A (ja) | 数値解析用モデル作成装置、数値解析用モデル作成方法および記憶媒体 | |
| US8325183B2 (en) | System and method for determining a position for an addendum mesh node | |
| JP3669490B2 (ja) | 圧延解析システム及びプログラム | |
| JP4744543B2 (ja) | 設計支援装置 | |
| JP2009142828A (ja) | 曲げ加工方法及び曲げ加工支援装置 | |
| JP2011107927A (ja) | 設計支援プログラムおよび設計支援装置 | |
| Barros et al. | Trimming of 3D solid finite element meshes: sheet metal forming tests and applications |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120501 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120626 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20120904 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120913 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150921 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5091331 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |