(57)【要約】【課題】より高い剛性、減衰性を実現することができる複合構造体の設計方法を提供する。【解決手段】複合構造体の設計方法は、複合構造体のモデルと性能を解析する環境とを構築する予備構築ステップＳ１０と、目標性能を達成する剛性材料に対する減衰材料の最適な３次元構造及び３次元配置を計算するための最適化条件を設定する最適化条件設定ステップＳ２０，Ｓ３０と、減衰材料の最適な３次元構造及び３次元配置を計算するための変数を設定する変数設定ステップＳ４０と、変数設定ステップＳ４０で設定された変数に基づいて計算された、減衰材料が配置された複合構造体の性能を解析する解析ステップＳ５０と、解析ステップＳ５０の結果が目標性能を達成したか否かを判定する判定ステップＳ６０と、を備えている。そして、判定ステップＳ６０により目標性能を達成していないと判定された場合には、変数設定ステップＳ４０にリターンする。【選択図】図３
2018206090.tif 000002
【特許請求の範囲】
【請求項１】
所定の剛性を有する剛性材料と、該剛性材料と比較して振動の減衰性の高い減衰材料とからなる複合構造体の設計方法であって、
前記複合構造体のモデルと、該複合構造体の性能を解析する環境とを構築する予備構築ステップと、
前記複合構造体の解析制約条件を与えて、該複合構造体が目標性能を達成する前記減衰材料の最適な３次元構造及び３次元配置を計算するための最適化条件を設定する最適化条件設定ステップと、
前記減衰材料の最適な３次元構造及び３次元配置を計算するための変数を設定する変数設定ステップと、
前記変数設定ステップで設定された変数に基づいて計算された、前記減衰材料が配置された前記複合構造体の性能を解析する解析ステップと、
前記解析ステップの結果が前記目標性能を達成したか否かを判定する判定ステップと、を備え、
前記判定ステップにより前記目標性能を達成していないと判定された場合には、前記変数設定ステップにリターンすることを特徴とする複合構造体の設計方法。
【請求項２】
前記解析ステップでは、前記複合構造体に外力を加えた場合の該複合構造体の所定箇所の変位である静剛性と、該複合構造体の固有振動数とを解析することを特徴とする請求項１に記載の複合構造体の設計方法。
【請求項３】
前記解析ステップでは、前記複合構造体に外力を加えた場合の振動の減衰速度である損失係数を解析することを特徴とする請求項１又は２に記載の複合構造体の設計方法。
【請求項４】
前記減衰材料の最適な３次元構造及び３次元配置の計算に有限要素法を用いることを特徴とする請求項１〜３の何れか１項に記載の複合構造体の設計方法。
【請求項５】
前記複合構造体中の前記剛性材料と前記減衰材料との境界を決定するためにレベルセット法を用い、
前記判定ステップから前記変数設定ステップにリターンして前記変数が変更される度に、前記レベルセット法で使用される関数が変更されることを特徴とする請求項１〜４の何れか１項に記載の複合構造体の設計方法。
【請求項６】
所定の剛性を有する剛性材料と、該剛性材料と比較して振動の減衰性の高い減衰材料とからなる複合構造体であって、
前記減衰材料は、有限要素法を用いて計算した最適な３次元構造及び３次元配置に基づいて、前記剛性材料の間に充填されていることを特徴とする複合構造体。
【請求項７】
前記剛性材料は、炭素繊維強化プラスチックであり、
前記減衰材料は、コルクラバー又は発泡体であることを特徴とする請求項６に記載の複合構造体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複合構造体の設計方法、及びその方法により設計された複合構造体に関する。
【背景技術】
【０００２】
現在、工作機械等に用いられる材料には、軽量、高剛性、高減衰性等の性能が求められているが、単一の材料で全ての要求性能を満足することは困難である。このため、炭素鋼と比較して剛性や減衰性の高い炭素繊維強化プラスチック（ＣＦＲＰ: Carbon Fiber Reinforced Plastics）と、高粘弾性を有する減衰材料とを組み合わせた複合材料が注目されている。
【０００３】
例えば、特許文献１には、カーボン繊維強化エポキシ樹脂やカーボン繊維強化プラスチック材料等の剛性材料の層の間に、粘弾性材料を含侵させた中間層を挿入した複合積層板の構造が開示されている。中間層は剛性材料と比べて柔らかい材料でできているため、複合積層板は機械特性に優れている（段落００１１，００１２，図１）。
【０００４】
このように、これまで、複合材料の減衰性を高めるため、剛性材料の表面や層間に減衰材料を付加する試み、また、減衰材料の配置の最適化に関する試みが多くなされてきた。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特表２０１０−５１６５１２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、剛性材料に対する減衰材料の最適化配置の手法は２次元的なものに止まっており、特許文献１に示されているような単純サンドウィッチ構造の積層構造体があるに過ぎなかった。
【０００７】
本発明は、このような事情に鑑みてなされたものであり、より高い剛性、減衰性を実現することができる複合構造体の設計方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
第１発明は、所定の剛性を有する剛性材料と、該剛性材料と比較して振動の減衰性の高い減衰材料とからなる複合構造体の設計方法であって、前記複合構造体のモデルと、該複合構造体の性能を解析する環境とを構築する予備構築ステップと、前記複合構造体の解析制約条件を与えて、該複合構造体が目標性能を達成する前記減衰材料の最適な３次元構造及び３次元配置を計算するための最適化条件を設定する最適化条件設定ステップと、前記減衰材料の最適な３次元構造及び３次元配置を計算するための変数を設定する変数設定ステップと、前記変数設定ステップで設定された変数に基づいて計算された、前記減衰材料が配置された前記複合構造体の性能を解析する解析ステップと、前記解析ステップの結果が前記目標性能を達成したか否かを判定する判定ステップと、を備え、前記判定ステップにより前記目標性能を達成していないと判定された場合には、前記変数設定ステップにリターンすることを特徴とする。
【０００９】
本発明の複合構造体の設計方法では、まず、予備構築ステップにおいて、設計する複合構造体のモデル（解析のための要素数や節点座標等）と、その性能を解析する環境（静荷重や拘束面の定義等）を構築する。次に、最適化条件設定ステップでは、複合構造体の解析制約条件（変位の限界値や固有振動数等）を与えて、複合構造体が目標性能を達成する減衰材料の３次元的な構造及び配置を最適化する最適化条件を設定する。
【００１０】
さらに、変数設定ステップで、減衰材料の最適な構造及び配置を計算するための変数を設定し、解析ステップで、設定された変数に基づいて計算された複合構造体の性能を解析する。最後に、判定ステップにて解析ステップの結果が目標性能を達成したか否かが判定され、達成していないと判定された場合には、変数設定ステップにリターンする。
【００１１】
これにより、目標性能が達成されるまで変数設定ステップにて新たな変数が設定されて、減衰材料の最適な構造及び配置が計算されるまでループする。すなわち、所定の剛性や減衰性を有する複合構造体を効率よく設計することができる。
【００１２】
第１発明の複合構造体の設計方法において、前記解析ステップでは、前記複合構造体に外力を加えた場合の該複合構造体の所定箇所の変位である静剛性と、該複合構造体の固有振動数とを解析することが好ましい。
【００１３】
この手法によれば、複合構造体の静剛性と固有振動数（共振周波数）とを解析することで、ある変数により設計された複合構造体の性能と目標性能とを比較して、減衰材料の最適な構造及び配置を評価することができる。
【００１４】
また、第１発明の複合構造体の設計方法において、前記解析ステップでは、前記複合構造体に外力を加えた場合の振動の減衰速度である損失係数を解析することが好ましい。
【００１５】
複合構造体の損失係数を解析することでも、設計された複合構造体の性能と目標性能とを比較して、減衰材料の最適な構造及び配置を評価することができる。なお、解析ステップにおいて、複合構造体の静剛性や動剛性に加えて、損失係数を解析するようにしてもよい。
【００１６】
また、第１発明の複合構造体の設計方法において、前記減衰材料の最適な３次元構造及び３次元配置の計算に有限要素法を用いることが好ましい。
【００１７】
減衰材料の最適な３次元構造及び３次元配置の計算する際には、有限要素法を用いることができる。有限要素法は、設計対象がどのような形状であっても、例えば、専用の解析ソフトウェアを用いて、迅速に目的の複合構造体を設計することができる。
【００１８】
また、第１発明の複合構造体の設計方法において、前記複合構造体中の前記剛性材料と前記減衰材料との境界を決定するためにレベルセット法を用い、前記判定ステップから前記変数設定ステップにリターンして前記変数が変更される度に、前記レベルセット法で使用される関数が変更されることが好ましい。
【００１９】
複合構造体中の剛性材料と減衰材料との境界を決定するための手法であるレベルセット法を用いることができる。判定ステップから変数設定ステップにリターンして変数が変更される度に、レベルセット法で使用される関数（レベルセット関数）が変更されるようにして、減衰材料の最適な構造及び配置を計算することができる。
【００２０】
第２発明は、所定の剛性を有する剛性材料と、該剛性材料と比較して振動の減衰性の高い減衰材料とからなる複合構造体であって、前記減衰材料は、有限要素法を用いて計算した最適な３次元構造及び３次元配置に基づいて前記剛性材料の間に充填されていることを特徴とする。
【００２１】
本発明の複合構造体は、有限要素法を用いることで、減衰材料が剛性材料に対して最適な３次元構造及び３次元配置で充填されている。これにより、設計された複合構造体は、目標性能（剛性や損失係数）を確実に満たすことができる。
【００２２】
第２発明の複合構造体において、前記剛性材料は、炭素繊維強化プラスチックであり、前記減衰材料は、コルクラバー又は発泡体であることが好ましい。
【００２３】
この構成によれば、炭素繊維強化プラスチックは、炭素鋼と比較して高い剛性、減衰性を有するので、剛性材料として適切である。また、コルクラバーやポリウレタン等の発泡体は、損失係数が高く、加工もし易いので、減衰材料として適切である。そして、これらを組み合わせることで、高剛性かつ高減衰性の複合構造体とすることができる。
【図面の簡単な説明】
【００２４】
【図１】設計する片持ち梁のモデルを説明する図。
【図２】(ａ)ＣＦＲＰの材料物性値テーブル。（ｂ）コルクラバーの材料物性値テーブル。
【図３】ＣＦＲＰに対する減衰材料の最適化配置のフローチャート。
【図４Ａ】最適化処理前の減衰材料の配置（単純サンドウィッチ構造）を説明する図。
【図４Ｂ】最適化処理後の減衰材料の配置（修正最適化構造）を説明する図。
【図５Ａ】最適化処理前後の片持ち梁の１次固有振動数（解析値）を比較したグラフ。
【図５Ｂ】最適化処理前後の片持ち梁の静的変位（解析値）を比較したグラフ。
【図５Ｃ】最適化処理前後の片持ち梁の損失係数（解析値）を比較したグラフ。
【図６Ａ】各構造の片持ち梁による１次固有振動数の実験値と解析値とを比較したグラフ。
【図６Ｂ】各構造の片持ち梁による損失係数の実験値と解析値とを比較したグラフ。
【発明を実施するための形態】
【００２５】
以下、本発明の複合構造体の設計方法の実施形態について説明する。
【００２６】
図１は、本発明の設計方法により設計する片持ち梁１の形状（設計領域）を示した図である。片持ち梁１は、長手方向が３００ｍｍ、幅方向が５０ｍｍ、厚み方向が２０ｍｍであって、図中の壁Ｗと接した面が固定面（固定端）となっている。
【００２７】
片持ち梁１は、基本的に炭素繊維強化プラスチック（以下、ＣＦＲＰ）を厚み方向に積層した構造である。ＣＦＲＰは、縦糸と横糸を交差させることで、両方向に均一な引張り強度を有する平織構造のものを使用する。なお、ＣＦＲＰの繊維方向は、片持ち梁１の長手方向及び幅方向に一致させている。
【００２８】
本発明の設計方法では、片持ち梁１を、後述する静剛性と動剛性の制約条件の下で、減衰性能が最大となるようにする。ＣＦＲＰは炭素鋼と比較して剛性や減衰性が高い材料であるが、片持ち梁１を工作機械等に適用するためには、より高い減衰性が必要となる。このため、片持ち梁１では、ＣＦＲＰの間に、ＣＦＲＰと比較して減衰性の高い減衰材料を充填する。
【００２９】
減衰材料としては、損失係数が高く、加工し易いコルクラバー（ＡＭＯＲＩＭ社、ＡＣＭ８７）を採用する。今回用いたＣＦＲＰとコルクラバーの材料物性値は、図２に示す通りである。
【００３０】
図２（ａ）には、ＣＦＲＰの密度［ｋｇ／ｍ3］、ヤング率［ＧＰａ］、せん断弾性係数［ＧＰａ］、ポアソン比及び損失係数が示されている。ここで、ヤング率とは、物体の引張り応力と単位長さ当りの伸び（ひずみ）との比であり、この値が大きくなるほど硬い材料といえる。
【００３１】
また、せん断弾性係数（剛性率）とは、せん断応力とせん断ひずみの比で定義され、せん断力による変形のし難さを示す数値である。また、ポアソン比とは、物体に弾性限界内で応力を加えたとき、応力に直角方向に発生するひずみと、応力方向に沿って発生するひずみとの比である。
【００３２】
なお、ＣＦＲＰは直交する３つの軸で弾性定数や線膨張係数が異なる直交異方弾性体であるため、これら数値が３種類ある。また、損失係数とは、物体を振動させた場合の減衰速度を示す数値である。
【００３３】
図２（ｂ）は、コルクラバーの密度［ｋｇ／ｍ3］、ヤング率［ＧＰａ］、ポアソン比及び損失係数が示されている。特に、コルクラバーの損失係数は、ＣＦＲＰの２７倍も大きい数値となっている。減衰材料として、コルクやポリウレタン、ポリスチレン等の発泡体を用いてもよい。
【００３４】
剛性の制約条件は、片持ち梁１の自由端である先端部（参照点）に１００Ｎを荷重したときの先端部のたわみが１ｍｍ以下（静剛性の制約）、かつ１次固有振動数が８０Ｈｚ以上（動剛性の制約）となるように設定した。
【００３５】
また、今回の設計方法では、片持ち梁１の最適構造を求めるために有限要素法を用いる。片持ち梁１は、長手方向を５０分割、幅方向を８分割、厚み方向を１５分割して、合計６０００の要素（メッシュ）に分割した。各要素は、８節点の６面体で、中間節点のない１次要素を用いた。そして、片持ち梁１の厚み方向の中央部分に減衰材料の層を挿入した状態（図４Ａの単純サンドウィッチ構造）を初期状態として最適化を行った。
【００３６】
次に、図３を参照して、複合構造体の設計方法、すなわち、ＣＦＲＰに対する減衰材料の最適な３次元構造及び３次元配置を行うためのフローチャートについて説明する。今回、有限要素法の解析ソフトウェアであるＡｂａｑｕｓ（登録商標）と数値計算ソフトウェアであるＭａｔｌａｂ（登録商標）とを連動させた。
【００３７】
まず、ステップＳ１０では、Ａｂａｑｕｓを用いて解析モデルと解析環境の構築を行う。具体的には、有限要素法の解析に必要な要素数や節点数を決定するとともに、物体に加わる時間的に変化しない荷重である静荷重や各節点の自由度を拘束する拘束条件の定義を行う。その後、ステップＳ２０に進む。
【００３８】
ステップＳ２０もＡｂａｑｕｓで行われる処理であり、最適化に必要なＣＦＲＰと減衰材料のデータを出力する。具体的には、ＣＦＲＰ及び減衰材料の要素マトリクス、荷重マトリクス、要素数及び節点座標を出力する。その後、ステップＳ３０に進む。
【００３９】
ステップＳ３０は、Ｍａｔｌａｂで行われる処理である。ここでは、最適化条件の設定が行われる。具体的には、節点番号とその座標、拘束される節点集合、静荷重のかかる節点集合の定義を行う。また、ステップＳ２０で出力した要素マトリクス及び荷重マトリクスの読込みと、静解析と動解析の定義を行う。その後、ステップＳ４０に進む。
【００４０】
ステップＳ４０もＭａｔｌａｂで行われる処理であり、設計変数α（初期値α1）の設定、更新を行う。ここで、設計変数αは、以下の（１）式に示すレベルセット関数φの膨張係数である。
【００４１】
【数１】
2018206090.tif 000003
なお、ξiはｉ要素の形状関数であり、Ｘは要素座標、ｔは最適化計算の繰り返し数である。
【００４２】
レベルセット関数φを用いるレベルセット法は、材料境界を定義する数学的手法の１つであり、ここでは、ＣＦＲＰと減衰材料の境界を定義するのに利用することができる。
【００４３】
また、ここでの数値計算は、逐次凸関数近似法のＧＣＭＭＡ（Globally Convergent version of the Method of Moving Asymptotes）法を用いた。ＧＣＭＭＡ法は、ステップ毎に目的関数及び制約関数をテイラー展開で凸関数に近似することで、最適化問題を簡易的に解いていく手法である。その後、ステップＳ５０に進む。
【００４４】
ステップＳ５０は、Ａｂａｑｕｓで行われる処理である。ここでは、ステップＳ４０で設定された設計変数αにより、静剛性と動剛性の解析を行う。具体的には、変位ベクトル、１次固有振動数及び固有ベクトルを出力する。その後、ステップＳ６０に進む。
【００４５】
ステップＳ６０は、Ｍａｔｌａｂで行われる処理である。ここでは、ステップ５０で出力された固有ベクトルを用いて損失係数や参照点の静的変位を算出し、収束条件を満たしているか否かを判定する。
【００４６】
そして、収束条件を満たしている場合（ステップＳ６０：ＹＥＳ）には最適化処理を終了し、収束条件を満たしていない場合（ステップＳ６０：ＮＯ）にはステップＳ４０にリターンする。
【００４７】
ステップＳ４０にリターンした場合には、Ｍａｔｌａｂにて設計変数αが更新される。そして、新たな設計変数α（＝α2）によりステップＳ５０以降の各処理が行われる。ＧＣＭＭＡ法では、レベルセット関数φを変化させることで境界を動かし、トポロジーの変化を表現することができる。以上の処理を複数回繰り返すことにより、発明者らは、片持ち梁１の最適化構造を求めることに成功した。
【００４８】
次に、図４Ａ、図４Ｂを参照して、最適化処理前後の片持ち梁の構造（減衰材料の３次元構造及び３次元配置）について説明する。
【００４９】
まず、図４Ａに、最適化処理を行う前の初期状態の片持ち梁１を示す。片持ち梁１は単純サンドウィッチ構造であり、厚み方向の中央の黒色部（厚み７．２５ｍｍ）が減衰材料２、その上面及び下面がＣＦＲＰ３である。なお、ここでも、図中の壁Ｗと接した面が固定面である。
【００５０】
次に、図４Ｂに、最適化処理を行った後の最適化構造の片持ち梁１’を示す。ここで、上記の最適化処理では、レベルセット関数φの値に基づいて、φ≧１の領域に減衰材料、φ≦−１の領域にＣＦＲＰ、−１＜φ＜１の領域にＣＦＲＰと減衰材料の中間的な性質を示す仮想的な物質（中間物質）を配置する計算を行った。
【００５１】
しかし、実際に片持ち梁を製造する場合には、中間物質を無くし、ＣＦＲＰと減衰材料のみで構成する必要がある。そこで、φ≧０．６の領域に減衰材料、φ＜０．６の領域にＣＦＲＰが配置されるように、新たなレベルセット関数φの閾値を設定した。そして、修正された最適化構造が、図示する片持ち梁１’である。すなわち、片持ち梁１’は、最初の最適化構造の性質を保持しつつ、ＣＦＲＰと減衰材料とを３次元的に再配置したもの（修正最適化構造）となっている。
【００５２】
片持ち梁１’の厚み方向の中央に略直線状に配置された黒色部が減衰材料２’、その上面及び下面がＣＦＲＰ３’である。これによれば、高ひずみエネルギー領域の近傍に減衰材料２’が配置されるように、最適化が進んだと考えられる。
【００５３】
次に、図５Ａ〜図５Ｃに、単純ＳＷ（サンドウィッチ）構造と修正最適化構造との各種性質を比較したグラフを示す。
【００５４】
まず、図５Ａは、両構造の１次固有振動数の解析値を比較したグラフである。図示するように、修正最適化構造の１次固有振動数は１０５．５［Ｈｚ］、単純ＳＷ構造の１次固有振動数は１０９．３［Ｈｚ］であった。単純ＳＷ構造の方が高い数値を示したが、修正最適化構造についても動剛性の制約条件である１次固有振動数（８０［Ｈｚ］）を十分満たす結果となった。
【００５５】
次に、図５Ｂは、両構造の静的変位の解析値を比較したグラフである。図示するように、修正最適化構造の静的変位は０．９９［ｍｍ］、単純ＳＷ構造の静的変位は１．０７［ｍｍ］であった。ここで、静剛性の制約条件は、参照点のたわみが１．０［ｍｍ］以下であったので、修正最適化構造は制約条件を満たすが、単純ＳＷ構造は制約条件を満たさないという結果となった。
【００５６】
次に、図５Ｃは、両構造の損失係数の解析値を比較したグラフである。図示するように、修正最適化構造の損失係数は０．１３２、単純ＳＷ構造の損失係数は０．１１６であった。修正最適化構造の損失係数は、単純ＳＷ構造の損失係数と比較して約１４％も高い値となり、目標である１０％以上の結果が得られた。
【００５７】
最後に、図６Ａ、図６Ｂを参照して、単純ＳＷ（サンドウィッチ）構造及び修正最適化構造の実験値と解析値とを比較したグラフを示す。なお、実験値とは、実際に、コルクラバーを加工してＣＦＲＰの板材の間に埋設させた複合構造体（単純ＳＷ構造と修正最適化構造）とを製造し、ハンマーでその中心を加振して計測した値である。
【００５８】
まず、図６Ａは、各構造の１次固有振動数の実験値及び解析値を比較したグラフである。単純ＳＷ構造と修正最適化構造の１次固有振動数の解析値は、図５Ａで示した数値（それぞれ１０９．３［Ｈｚ］、１０５．５［Ｈｚ］）であった。
【００５９】
一方、実験値は、それぞれ１０９．６［Ｈｚ］、１０６．４［Ｈｚ］となり、解析値に非常に近い数値が得られた。従って、図３の処理によって求めたＣＦＲＰと減衰材料とによる最適化構造の妥当性が示された。
【００６０】
次に、図６Ｂは、各構造の損失係数の実験値及び解析値を比較したグラフである。単純ＳＷ構造と修正最適化構造の損失係数の解析値は、図５Ｃで示した数値（それぞれ０．１１６、０．１３２）であった。
【００６１】
一方、実験値は、それぞれ０．１０１、０．１２４となり、実験値の方が１割程度低い結果が得られた。これは、試験片の把持が完全固定でないことや、解析においてはＣＦＲＰと減衰材料を貼りあわせる接着剤の影響を考慮していないこと等が原因と考えられる。
【００６２】
以上で説明したように、本発明の設計方法（図３の処理）によりＣＦＲＰと減衰材料（コルクラバー）の最適な３次元構造及び３次元配置の構造を求め、解析値を調べた。また、実際にＣＦＲＰとコルクラバーを加工して最適化構造を製造し、加振する実験を行った。この結果、特に、１次固有振動数は解析値と実験値とが非常に近い値を示したので、一連の最適化処理の妥当性が証明された。
【００６３】
今回、片持ち梁のモデルで最適化処理を行ったが、本発明の設計方法は、様々な形状の部品設計に適用することができる。例えば、工作機械のスピンドルを保持するクイル（円筒状部材）の設計に適用可能である。この部材の場合、円筒周囲に配置されたリニアモータによる吸引力等を考慮して制約条件を設定することで、ＣＦＲＰと減衰材料による最適な３次元構造及び３次元配置の構造を求めることができる。
【符号の説明】
【００６４】
１  片持ち梁（単純サンドウィッチ構造）
１’  片持ち梁（修正最適化構造）
２，２’  減衰材料
３，３’ＣＦＲＰ（剛性材料）
【図１】
2018206090.tif 000004
【図２】
2018206090.tif 000005
【図３】
2018206090.tif 000006
【図４Ａ】
2018206090.tif 000007
【図４Ｂ】
2018206090.tif 000008
【図５Ａ】
2018206090.tif 000009
【図５Ｂ】
2018206090.tif 000010
【図５Ｃ】
2018206090.tif 000011
【図６Ａ】
2018206090.tif 000012
【図６Ｂ】
2018206090.tif 000013