(57)【要約】  管状体設計方法は、複数の繊維強化樹脂層を積層してなる管状体の設計変数の値を順次変化させながら、積層モデルを用いた離散的手法によって繰り返し計算を行って、前記管状体の複数の目的関数を算出する解析手順と、算出された前記複数の目的関数の中から最大値又は最小値となる前記設計変数の値を探索する探索手順と、を含む。
2017195898.tif 000002
【特許請求の範囲】
【請求項１】
複数の繊維強化樹脂層を積層してなる管状体の設計変数の値を順次変化させながら、積層モデルを用いた離散的手法によって繰り返し計算を行って、前記管状体の複数の目的関数を算出する解析手順と、
算出された前記複数の目的関数の中から最大値又は最小値となる前記設計変数の値を探索する探索手順と、
を含む管状体設計方法。
【請求項２】
前記目的関数は、複数の変数からなる
請求項１に記載の管状体設計方法。
【請求項３】
前記離散的手法は、有限要素法、境界要素法、個別要素法、Ｅｌｅｍｅｎｔ  Ｆｒｅｅ  Ｇａｌｅｒｋｉｎ（ＥＦＧ）法、ｅＸｔｅｎｄｅｄ  ＦＥＭ（ＸＦＥＭ）法、Ｓｍｏｏｔｈｅｄ  Ｐａｒｔｉｃｌｅ  Ｈｙｄｒｏｄｙｎａｍｉｃｓ（ＳＰＨ）法、Ｃａｐｓｕｌａｒ  Ｐａｒｔｉｃｌｅ  Ｍｅｔｈｏｄのうち少なくとも１つであり、
前記目的関数は、前記管状体の強度、重量、曲げ剛性、捩り剛性のうちいずれか１つ又は複数の線形結合で表されたものであり、
前記探索手順は、前記管状体の、積層された各層の積層角度、材料定数、パイプ形状のうち少なくとも１つを前記設計変数として選択する、
請求項１または請求項２に記載の管状体設計方法。
【請求項４】
前記管状体の積層される複数の層のうち少なくとも１つの層は、２つ以上の領域を備え、
前記解析手順は、前記層の領域毎に前記複数の目的関数を算出し、
前記探索手順は、前記層の領域毎に前記設計変数の値を探索する、
請求項１から請求項３のいずれか１項に記載の管状体設計方法。
【請求項５】
前記解析手順が、更に、前記設計変数である前記層の積層角度を離散的に選択するサンプリング手順、を含み、
前記解析手順は、前記サンプリング手順によって選択した、前記層の積層角度を、前記順次変化させた前記設計変数として、前記複数の目的関数を算出し、
前記探索手順は、前記解析手順によって算出した前記複数の目的関数に近似するよう算出した応答曲面を作成する応答曲面作成手順、を含み、
前記探索手順は、前記応答曲面作成手順によって作成した前記応答曲面中で前記目的関数が最大値又は最小値をとるような前記設計変数の値を探索する、
請求項１から請求項４のいずれか１項に記載の管状体設計方法。
【請求項６】
前記サンプリング手順は、サンプリング空間を繰り返しのある無限空間に拡張し、各サンプリング点同士の距離の最小値を最大化するサンプリングを行う、
請求項５に記載の管状体設計方法。
【請求項７】
前記応答曲面作成手順は、フーリエ級数近似にて応答曲面の作成を行う、
請求項５または請求項６に記載の管状体設計方法。
【請求項８】
前記フーリエ級数近似は、一次近似又は二次近似である、
請求項７に記載の管状体設計方法。
【請求項９】
前記管状体は、ゴルフシャフトであり、
請求項１から請求項８のいずれか１項に記載の管状体設計方法によってゴルフシャフトを設計する、ゴルフシャフト設計方法。
【請求項１０】
最内層のアングル層２層を除いた各層の積層角度の和である垂直積層成分が、バット側からチップ側に向かうにしたがって小さくなるゴルフシャフト。
【請求項１１】
少なくとも１層が積層角度を変えた複数のプリプレグからなる
請求項１０に記載のゴルフシャフト。
【請求項１２】
バット端部を基準とし、０ｍｍ以上かつ３２５ｍｍ未満を位置Ｃ近傍、３２５ｍｍ以上かつ６２５ｍｍ未満を位置Ｂ近傍、６２５ｍｍ以上〜９７５ｍｍを位置Ａ近傍、９７５ｍｍ以上かつチップ端までを位置Ｔ近傍と定義したとき、１）〜４）を満たす、
請求項１０または請求項１１に記載のゴルフシャフト。
１）位置Ｃ近傍での垂直積層成分が、１９５以上、２２５未満である
２）位置Ｂ近傍での垂直積層成分が、１８５以上であり、Ｃ点の垂直積層成分よりも小さい
３）位置Ａ近傍での垂直積層成分が、９０以上、１８０未満である
４）位置Ｔ近傍での垂直積層成分が、０である
【請求項１３】
位置Ｃ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、９０°、１５°〜４５°、９０°であり、位置Ｂ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、９０°、５°〜１５°、９０°であり、位置Ａ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、０°、０°、９０°であり、位置Ｔ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、０°、０°、０°である、
請求項１０から１２のいずれか１項に記載のゴルフシャフト。
【請求項１４】
位置Ｃ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、９０°、＋１５°〜４５°、−１５°〜−４５°、９０°であり、位置Ｂ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、９０°、−５°〜−１５°、＋５°〜＋１５°、９０°であり、位置Ａ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、０°、０°、９０°であり、位置Ｔ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、０°、０°、０°である、
請求項１０または請求項１１に記載のゴルフシャフト。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、管状体設計方法、ゴルフシャフト設計方法およびゴルフシャフトに関する。
本願は、２０１６年５月１３日に出願された日本国特願２０１６−０９７２６８号に基づき優先権を主張し、その内容をここに援用する。
【背景技術】
【０００２】
２００８年１月１日より、公益財団法人日本ゴルフ協会（ＪＧＡ）の規制により高反発クラブがゴルフ規則に不適合となった。この規則では、クラブフェースの反発係数（ＣＯＲ）が０．８３０以内でなくてはならない。これにより、ヘッド単体で飛距離を追及することが難しくなったため、ゴルフクラブの製造会社は、シャフトに着目するようになった。シャフトで飛距離を伸ばす手段として、軽量化が挙げられる。しかしながら現状を鑑みると、ゴルフシャフトの設計は、試行錯誤（トライアル・アンド・エラー）の繰り返しによって設計されており、真に最軽量の積層構成がなされているか不明である。また、試行錯誤のため、試作時間もコストもかかる。仮に試作時間も予算も十分にあったとしても、全組み合わせを試作することは現実的に不可能である。これは、ゴルフシャフトだけではなく、釣竿、自転車のフレーム、テニスラケット等のスポーツ用途部材、さらにはロボットアームやカーボンロール、ドライブシャフト等の産業用途部材にも共通の課題である。【０００３】
従来、試作時間やコスト削減のためにシミュレーションと最適化計算を用いることが提案されている。
例えば、特許文献１に記載の技術分野はゴルフシャフトの設計方法及び製造方法に関するもので、シャフトの設計用パラメータを関数表現して離散的手法を用いた解析と最適化計算を行うものである。また、特許文献２に記載の技術もゴルフシャフトに関するものであり、用具によるスイング動作の変化に着目し、シャフトの性能を最適化したものである。【０００４】
また、ゴルフシャフト以外の分野に関しては、特許文献３、特許文献４に記載の技術が挙げられる。特許文献３に記載の技術は、樹脂の物性を大域的最適化によって決定しようとするものである。特許文献４に記載の技術は、繊維強化樹脂と金属のハイブリッド材に関するものである。特許文献４に記載の技術は、繊維強化樹脂層と金属の間にある接着層が強度発現に大きく起因するため、各層を積層シェル要素でモデリングし、界面となる接着層のパラメータを最適化する手法である。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】日本国特開２００８−３０２０１７号
【特許文献２】日本国特開２０１１−０００４２５号
【特許文献３】日本国国際公開第２０１３／０４２６００号
【特許文献４】日本国特開２００３−９４５５９号
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、特許文献１に記載の技術では、シャフトの最適な曲げ剛性値を知ることはできるが、それを達成するための具体的な積層構成を知ることができない。特許文献２に記載の技術においても同様にシャフトの最適な性能を知ることはできるが、具体的な積層構成を知ることができない。特許文献３に記載の技術では、特許文献１、２と同じく最適な樹脂の性能を知ることはできても、具体的な構成を知ることはできない。特許文献４に記載の技術では、界面のパラメータを知ることはできても繊維強化樹脂層の最適構成を知ることはできない。このように、特許文献１〜４に記載の技術においては、複数の繊維強化樹脂層からなる管状体の積層構成まで知ることができない。【０００７】
本発明に係る態様は、上記の課題を鑑みてなされたものであって、目的に応じて最適な積層構成を設計することができる管状体設計方法およびゴルフシャフト設計方法を提供することを目的とする。併せて、管状体の一例として、ゴルフシャフトを提供することを目的とする。
【課題を解決するための手段】
【０００８】
［１］上記目的を達成するため、本発明の一態様に係る管状体設計方法は、複数の繊維強化樹脂層を積層してなる管状体の設計変数の値を順次変化させながら、積層モデルを用いた離散的手法によって繰り返し計算を行って、前記管状体の複数の目的関数を算出する解析手順と、算出された前記複数の目的関数の中から最大値又は最小値となる前記設計変数の値を探索する探索手順と、を含む。
［２］上記［１］の態様において、前記目的関数は、複数の変数からなってもよい。
［３］上記［１］または［２］の態様において、前記離散的手法は、有限要素法、境界要素法、個別要素法、Ｅｌｅｍｅｎｔ  Ｆｒｅｅ  Ｇａｌｅｒｋｉｎ（ＥＦＧ）法、ｅＸｔｅｎｄｅｄ  ＦＥＭ（ＸＦＥＭ）法、Ｓｍｏｏｔｈｅｄ  Ｐａｒｔｉｃｌｅ  Ｈｙｄｒｏｄｙｎａｍｉｃｓ（ＳＰＨ）法、Ｃａｐｓｕｌａｒ  Ｐａｒｔｉｃｌｅ  Ｍｅｔｈｏｄのうち少なくとも１つであり、前記目的関数は、前記管状体の強度、重量、曲げ剛性、捩り剛性のうちいずれか１つ又は複数の線形結合で表されたものであり、前記探索手順は、前記管状体の、積層された各層の積層角度、材料定数、パイプ形状のうち少なくとも１つを前記設計変数として選択してもよい。
［４］上記［１］から［３］のいずれか１つの態様において、前記管状体の積層される複数の層のうち少なくとも１つの層は、２つ以上の領域を備え、前記解析手順は、前記層の領域毎に前記複数の目的関数を算出し、前記探索手順は、前記層の領域毎に前記設計変数の値を探索してもよい。
［５］上記［１］から［４］のいずれか１つの態様において、前記解析手順が、更に、前記設計変数である前記層の積層角度を離散的に選択するサンプリング手順、を含み、前記解析手順は、前記サンプリング手順によって選択した、前記層の積層角度を、前記順次変化させた前記設計変数として、前記複数の目的関数を算出し、前記探索手順は、前記解析手順によって算出した前記複数の目的関数に近似するよう算出した応答曲面を作成する応答曲面作成手順、を含み、前記探索手順は、前記応答曲面作成手順によって作成した前記応答曲面中で前記目的関数が最大値又は最小値をとるような前記設計変数の値を探索してもよい。
［６］上記［５］の態様において、前記サンプリング手順は、サンプリング空間を繰り返しのある無限空間に拡張し、各サンプリング点同士の距離の最小値を最大化するサンプリングを行ってもよい。
［７］上記［５］または［６］の態様において、前記応答曲面作成手順は、フーリエ級数近似にて応答曲面の作成を行ってもよい。
［８］上記［７］の態様において、前記フーリエ級数近似は、一次近似又は二次近似であってもよい。【０００９】
［９］上記目的を達成するため、本発明の一態様に係るゴルフシャフト設計方法は、前記管状体が、ゴルフシャフトであり、上記［１］から［８］のいずれか１つの態様によってゴルフシャフトを設計する。【００１０】
［１０］上記目的を達成するため、本発明の一態様に係るゴルフシャフトは、最内層のアングル層２層を除いた各層の積層角度の和である垂直積層成分が、バット側からチップ側に向かうにしたがって小さくなる。
［１１］上記［１０］の態様において、少なくとも１層が積層角度を変えた複数のプリプレグからなってもよい。
［１２］上記［１０］または［１１］の態様において、バット端部を基準とし、０ｍｍ以上かつ３２５ｍｍ未満を位置Ｃ近傍、３２５ｍｍ以上かつ６２５ｍｍ未満を位置Ｂ近傍、６２５ｍｍ以上〜９７５ｍｍを位置Ａ近傍、９７５ｍｍ以上かつチップ端までを位置Ｔ近傍と定義したとき、１）〜４）を満たしてもよい。
１）位置Ｃ近傍での垂直積層成分が、１９５以上、２２５未満である
２）位置Ｂ近傍での垂直積層成分が、１８５以上であり、Ｃ点の垂直積層成分よりも小さい
３）位置Ａ近傍での垂直積層成分が、９０以上、１８０未満である
４）位置Ｔ近傍での垂直積層成分が、０である
［１３］上記［１０］から［１２］のいずれか１つの態様において、位置Ｃ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、９０°、１５°〜４５°、９０°であり、位置Ｂ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、９０°、５°〜１５°、９０°であり、位置Ａ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、０°、０°、９０°であり、位置Ｔ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、０°、０°、０°であってもよい。
［１４］上記［１０］または［１１］の態様において、位置Ｃ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、９０°、＋１５°〜４５°、−１５°〜−４５°、９０°であり、位置Ｂ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、９０°、−５°〜−１５°、＋５°〜＋１５°、９０°であり、位置Ａ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、０°、０°、９０°であり、位置Ｔ近傍での積層角度が前記管状体の内側から順に、＋４５°、−４５°、０°、０°、０°、０°、０°であってもよい。
【発明の効果】
【００１１】
本発明によれば、所望の物性を目的関数として選択し、前記目的関数に材料物性と構成を設計変数として入力した後、前記設計変数を順次変化させながら繰り返し計算を実施することで、最適解として所望の物性を満たす材料構成を得ることができる。具体的には、複数の繊維強化樹脂層からなる管状体において、目的に応じて各層の最適な積層構成を設計することができる。特に、軽量化や高強度化に応用することができ、複数の繊維強化樹脂層からなる管状体に特有の課題を解決できる。また、前記応用事例においては、従来の管状体にはない、軽量かつ高剛性のものが得られる。
【図面の簡単な説明】
【００１２】
【図１】実施形態に係る設計装置のブロック図である。
【図２】実施例に係る設計処理のフローチャートである。
【図３】実施例に係る離散的手法を用いた解析のための積層シェル要素モデルを示す図である。
【図４】ＳＧ基準  ＣＰＳＡ  ００９８における３点曲げ試験を示す図である。
【図５】実施例に係る３点曲げ試験に基づくモデルを示す図である。
【図６】実施例に係る設計処理のフローチャートである。
【図７】実施例に係るサンプリング手順を示す図である。
【図８】実施例に係る６層目と７層目のサンプリング点を用いて算出した３点曲げ強度に基づく応答曲面の一例を示す図である。
【図９】実施例に係る５点の場合のサンプリング模式図である。
【図１０】実施例に係る１次のフーリエ級数近似を行った応答曲面を示す図である。
【図１１】実施例に係る２次のフーリエ級数近似を行った応答曲面を示す図である。
【図１２】実施例に係る制約条件を課したときの位置Ｃにおける最適構成例を示す図である。
【図１３】実施例に係る制約条件を課したときの位置Ｂ、Ａ、Ｔにおける最適構成例を示す図である。
【図１４】本実施形態の手法を用いて求めた強度最大化のゴルフシャフトの最適積層の構成例を示す図である。
【発明を実施するための形態】
【００１３】
以下、図面を参照しながら本発明の実施形態について説明する。
なお、本実施形態では、複数の層（例えば繊維強化樹脂層）からなる管状体の例として、ゴルフシャフトを例に説明するが、これに限られない。本実施形態は、釣竿、自転車のフレーム、テニスラケット等のスポーツ用途部材、さらにはロボットアームやカーボンロール、ドライブシャフト等の産業用途部材の複数の繊維強化樹脂層からなる管状体に適用することができる。
【００１４】
＜設計装置１の構成＞
図１は、実施形態に係る設計装置１のブロック図である。図１に示すように、設計装置１は、目的関数選択部１１、設計変数選択部１２、設計変数変更部１３、解析部１４、目的関数最大化・最小化部１５、および記憶部１６を備える。
【００１５】
設計装置１は、利用者（設計者、製造者等）によって設計装置１に入力される入力情報を用いて、複数の繊維強化樹脂層を積層してなる管状体であるゴルフシャフトを設計し、設計した情報を表示部や印刷装置等の外部装置に出力する。ここで、入力情報とは、目的関数と設計変数の情報である。また、出力される情報は、例えば、積層数、材料の情報、各層の積層角度（配向角ともいう）、重量、強度等である。なお、設計装置１は、目的関数と設計変数の情報のうちいくつかの情報を、予め記憶部１６に記憶させておいてもよい。目的関数と設計変数については、後述する。
なお、設計装置１の一部または全ては、パーソナルコンピュータおよびパーソナルコンピュータにインストールされているソフトウェアによって実現するようにしてもよい。
【００１６】
目的関数選択部１１は、記憶部１６が記憶する情報、または入力された入力情報から少なくとも１つの目的関数を選択する。ここで、目的関数とは、例えば、ゴルフシャフトの強度、ゴルフシャフトの重量、ゴルフシャフトの曲げ剛性、ゴルフシャフトの捩り剛性、３点曲げ強度等である。なお、目的関数選択部１１は、複数の目的関数を選択してもよい。また、目的関数は、選択した複数の目的関数を線形結合して表した情報、または複数の目的関数が線形結合されて表された情報であってもよい。目的関数選択部１１は、選択した目的関数を解析部１４に出力する。
【００１７】
設計変数選択部１２は、記憶部１６が記憶する情報、または入力された入力情報から少なくとも１つの設計変数を選択する。ここで、設計変数とは、管状体であるゴルフシャフトを設計するために必要な情報であり、ゴルフシャフトの積層構成、ゴルフシャフトの材料定数、ゴルフシャフトの各層の積層角度、ゴルフシャフトの各部の径、ゴルフシャフトの各部の厚み等である。なお、設計変数選択部１２は、複数の設計変数を選択してもよい。また、設計変数選択部１２は、選択しなかった変数を、任意の固定値として解析部１４に出力するようにしてもよい。任意の固定値は、例えば記憶部１６が記憶するようにしてもよい。設計変数選択部１２は、例えば、設計変数に積層角度のみを選択し、材料定数、径、厚みを選択しなかったとする。この場合、設計変数選択部１２は、ゴルフシャフトに一般的に使用される材料定数、径、厚みを固定値とする。そして、解析部１４は、これらの固定値を用いて、例えば全てのシャフトの平均的な物性等の解析を実行するようにしてもよい。また、軽量シャフトにおいて強度を最大化することが目的であった場合、設計変数選択部１２は、母集団を軽量シャフトに用いられる材料の物性値とするようにしてもよい。
【００１８】
設計変数変更部１３は、設計変数選択部１２が選択した設計変数を順次変更する。設計変数変更部１３は、例えば、設計変数として積層角度を選んだときは、積層角度を所定の角度（例えば４５°）刻みで変更する。設計変数変更部１３は、変更した設計変数を解析部１４に出力する。設計変数変更部１３は、管状体であるゴルフシャフトを作成するために、例えば繊維強化樹脂層を複数層積層するため各層ごとに積層角度を変更する。なお、材料の積層方向を変更すると、ヤング率が積層角度に応じて変化する。設計変数変更部１３は、積層角度毎のヤング率を周知の手法を用いて算出する。なお、変更時に、設計変数変更部１３は、計算時間短縮のため、実験計画法、数理計画法等を用いてなるべく効率的に設計変数を変更することが好ましい。
【００１９】
解析部１４は、目的関数選択部１１が選択した目的関数と、設計変数変更部１３が順次変更した設計変数と、積層シェル要素モデルを用いて、離散的手法を用いて解析を行う。
離散的手法は、有限要素法、境界要素法、個別要素法、Ｅｌｅｍｅｎｔ  Ｆｒｅｅ  Ｇａｌｅｒｋｉｎ（ＥＦＧ）法、ｅＸｔｅｎｄｅｄ  ＦＥＭ（ＸＦＥＭ）法、Ｓｍｏｏｔｈｅｄ  Ｐａｒｔｉｃｌｅ  Ｈｙｄｒｏｄｙｎａｍｉｃｓ（ＳＰＨ）法、Ｃａｐｓｕｌａｒ  Ｐａｒｔｉｃｌｅ  Ｍｅｔｈｏｄ（ＣＰＭ）のうち少なくとも１つである。一般的に普及しており、信頼性が高いため、好ましくは有限要素法が良い。解析部１４は、例えば構造解析ソフトウェアを用いて解析を行う。構造解析ソフトウェアは、例えばＬｉｖｅｒｍｏｒｅ  Ｓｏｆｔｗａｒｅ  Ｔｅｃｈｎｏｌｏｇｙ  Ｃｏｒｐｏｒａｔｉｏｎ社のＬＳ−ＤＹＮＡ（登録商標）である。例えば、解析部１４は、ＬＳ−ＤＹＮＡ（登録商標）の解法設定において、時間発展を動的陽解法に設定し、解析ソルバーを非線形行動解析に設定し、離散的手法を有限要素法に設定する。また、解析部１４は、ＬＳ−ＤＹＮＡ（登録商標）の要素設定において、要素形態を肉厚シェル要素に設定し、要素定式化を２次元平面応力要素に設定し、アワグラスコントロールを粘性型に設定する。なお、上述した例は一例であり、解析部１４は、他の設定を用いてもよい。なお、積層シェル要素モデルについては、後述する。
【００２０】
目的関数最大化・最小化部１５は、解析部１４が解析した結果において、目的関数の最大値又は最小値を探索することで最適点を求める。目的関数最大化・最小化部１５は、目的関数が、例えばゴルフシャフトの強度、ゴルフシャフトの曲げ剛性、ゴルフシャフトの捩り剛性の場合、最大値を探索し、目的関数が、例えばゴルフシャフトの重量の場合、最小値を探索する。実サンプルによる強度試験は、結果のばらつきが大きいことから、シミュレーションによる理論解を知ることが有用であることから、目的関数として強度を選択することは極めて有意義である。なお、目的関数最大化・最小化部１５は、極大値又は極小値を探索するようにしてもよい。また、目的関数最大化・最小化部１５は、最大化または最小化に、例えばＳｉｍｕｌａｔｅｄ  Ａｎｎｅａｌｉｎｇ法（焼きなまし法）等を用いる。ここで、Ｓｉｍｕｌａｔｅｄ  Ａｎｎｅａｌｉｎｇ法とは、大域的最適化問題への汎用の乱択アルゴリズムであり、広大な探索空間内の与えられた関数の大域的最適解に対して、よい近似を与える。なお、目的関数最大化・最小化部１５は、Ｓｉｍｕｌａｔｅｄ  Ａｎｎｅａｌｉｎｇ法以外の手法を用いて最大化または最小化を行うようにしてもよい。目的関数最大化・最小化部１５は、例えば、解析部１４が解析した結果を用いて例えば放射基底関数（ＲＢＦ；Ｒａｄｉａｌ  Ｂａｓｉｓ  Ｆｕｎｃｔｉｏｎ）等を用いて、応答曲面を作成するようにしてもよい。そして、目的関数最大化・最小化部１５は、作成した応答曲面を用いて、解空間の近似関数を算出してもよく、最急勾配法、動的計画法等の各種最適化手法を用いて最適点のみを導出してもよい。
【００２１】
記憶部１６は、目的関数と設計変数の情報のうちいくつかの情報を、予め記憶する。
【実施例１】
【００２２】
実施例１における設計処理手順の一例を説明する。
図２は、実施例に係る設計処理のフローチャートである。
（ステップＳ１）目的関数選択部１１は、記憶部１６が記憶する情報、または入力された入力情報から少なくとも１つの目的関数を選択する（目的関数選択手順、解析手順）。【００２３】
（ステップＳ２）設計変数選択部１２は、記憶部１６が記憶する情報、または入力された入力情報から少なくとも１つの設計変数を選択する（設計変数選択手順、解析手順）。
（ステップＳ３）設計変数変更部１３は、設計変数選択部１２が選択した設計変数を順次変更する（設計変数変更手順、解析手順）。【００２４】
（ステップＳ４）解析部１４は、積層シェル要素モデルを用いて、例えば有限要素法で解析を行う（解析手順）。
（ステップＳ５）目的関数最大化・最小化部１５は、ステップＳ１で目的関数選択部１１が選択した目的関数の最大値又は最小値を探索することで最適点を求める（目的関数最大化・最小化手順、探索手順）。【００２５】
次に、離散的手法を用いた解析のための積層シェル要素モデルの一例を説明する。
図３は、実施例に係る離散的手法を用いた解析のための積層シェル要素モデルを示す図である。なお、図３に示す図では、半断面図で示している。
図３に示すように、積層シェル要素モデル１００は、内径が１４．１［ｍｍ］、長さが３００［ｍｍ］、繊維強化樹脂層を第１層Ｌｙ１〜第７層Ｌｙ７による７層積層し、各層の積層角度が０°（度）（長手方向に繊維が沿う方向）である。
また、積層シェル要素モデルでは、有限個の要素を節点で接続することにより、略棒状のモデルとして定義する。また、各要素それぞれには、例えば、要素の長さ、曲げ弾性率、ヤング率Ｅ、断面二次モーメントＩ等のパラメータが含まれる。【００２６】
次に、記憶部１６が記憶する材料の情報、または設計装置１に入力される材料の情報の一例を説明する。
表１は、記憶部１６が記憶する材料の情報、または設計装置１に入力される材料の情報の一例を示す。【００２７】
【表１】
2017195898.tif 000003
【００２８】
表１に示す例では、繊維に樹脂等を含浸させたシート状のものであるプリプレグ（ｐｒｅｐｒｅｇ）Ａ〜Ｉの情報を、記憶部１６が記憶する。記憶部１６は、プリプレグ（χ）（χはＡ〜Ｉのいずれか）に、品番、繊維の引張弾性率［ＧＰａ］、目付［ｇ／ｍ２］、樹脂含有率［質量％］、および厚さ［ｍｍ］を対応付けて記憶する。例えば、記憶部１６は、プリプレグＡに、品番がＴＲ３５０Ｃ０５０Ｓ、繊維の引張弾性率が２３５［ＧＰａ］、目付が７５［ｇ／ｍ２］、樹脂含有率が２５［質量％］、および厚さが０．０６２［ｍｍ］を対応付けて記憶する。ここで、目付とは、単位面積当たりの重さである。なお、表１では、プリプレグＣを除き、後の解析を考慮して糸の弾性率以外を同一となるように調整してある。
なお、表１に示した例は一例であり、記憶部１６は、他の品番の材料、０．０６２［ｍｍ］以外の厚さ等の材料の情報も記憶していてもよい。【００２９】
ここで、設計装置１が行う計算例を、図２を参照しつつ説明する。
目的関数選択部１１は、記憶部１６が記憶する情報から、ゴルフシャフトの曲げ剛性を選択する（ステップＳ１）。続けて、設計変数選択部１２は、記憶部１６が記憶する情報から、プリプレグの種類を選択する（ステップＳ２）。続けて、設計変数変更部１３は繊維強化樹脂層を複数層積層するため各層ごとに積層角度プリプレグの種類を変更する（ステップＳ３）。解析部１４は、積層シェル要素モデルを用いて、例えば有限要素法によって解析を行う（ステップＳ４）。目的関数最大化・最小化部１５は、ステップＳ４で解析した結果において、目的関数である剛性の最大値を探索することで最適点を求める（ステップＳ５）。
この結果、最適点に近い材料としては、プリプレグＩ（表１）であった。
【実施例２】
【００３０】
次に、実施例２について説明する。
一般に、繊維強化樹脂層は、剛性が高いほど脆い。このため、弾性率だけでなく強度も考慮する必要がある。特に、０°層のみで構成した場合、周方向に潰れが発生するため強度が大きく低下する場合がある。この場合は、積層角度を変更することで強度を高めることができる。このため、実施例２では、このような強度問題を考慮するため目的関数として３点曲げ強度を選択し、設計変数として積層角度を選択する例を説明する。【００３１】
ここで、ＳＧ基準  ＣＰＳＡ  ００９８、ゴルフクラブ用シャフトの認定基準及び基準確認方法（一般財団法人製品安全協会、通商産業大臣承認５産第２０８７号、平成５年１０月４日）に記載の『Ｃ型シャフトの強度（１）３点曲げ試験』の概要について説明する。ここで、Ｃ型とは、シャフトの材質が繊維強化プラスティック製等の非金属製であるものである。【００３２】
図４は、ＳＧ基準  ＣＰＳＡ  ００９８における３点曲げ試験を示す図である。ここで、位置Ｔとは、先端から９０［ｍｍ］の位置である。位置Ａとは、先端から１７５［ｍｍ］の位置である。位置Ｂとは、先端から５２５［ｍｍ］の位置である。位置Ｃとは、グリップエンドから１７５［ｍｍ］の位置である。なお、図４において、符号Ｓは、シャフトである。
ＳＧ基準  ＣＰＳＡ  ００９８によれば、「Ｃ型シャフトにあたっては、３点曲げ試験を行ったときに破壊しないこと」が認定基準となっている。【００３３】
この試験に基づいて、実施例では、以下のモデルを用いる。
図５は、実施例に係る３点曲げ試験に基づくモデルを示す図である。図５に示すようにシャフトの積層シェル要素モデル１００は、内径が１４．１［ｍｍ］、長さが３００［ｍｍ］、繊維強化樹脂層を第１層Ｌｙ１〜第７層Ｌｙ７による７層積層したものである。
３点の位置は、図５に示すように、一端から１５ｍｍの位置、両端（一端と他端）の中心の１５０［ｍｍ］に位置、他端から１５［ｍｍ］の位置である。ＳＧ基準  ＣＰＳＡ  ００９８に従って、一端および他端それぞれから１５［ｍｍ］の位置のＲが１２．５［ｍｍ］であり、両端（一端と他端）の中心の１５０［ｍｍ］に位置における圧子２０１のＲが７５［ｍｍ］である。このモデルでは、中央の圧子２０１が、上から一定の速度で降りてくるものである。【００３４】
設計装置１の構成は、実施例１と同じである。
設計変数変更部１３は、設計変数選択部１２が選択した設計変数を順次変更する。設計変数変更部１３は、設計変数を変更する処理において、設計変数を離散的に決めるサンプリングを行う。なお、設計変数変更部１３が行うサンプリングについては、後述する。
目的関数最大化・最小化部１５は、目的関数の最大値又は最小値を探索することで最適点を求める。目的関数最大化・最小化部１５は、目的関数の最大値又は最小値を探索する処理において、目的関数を近似的に算出して応答曲面を作成する。なお、目的関数最大化・最小化部１５が行う応答曲面の作成については、後述する。【００３５】
実施例２における設計処理手順の一例を説明する。
図６は、実施例に係る設計処理のフローチャートである。なお、実施例１と同じ処理については、同じ符号を用いて説明を省略する。
（ステップＳ３１）設計変数変更部１３は、設計変数選択部１２が選択した設計変数を順次変更する（設計変数変更手順、解析手順）。なお、設計変数変更部１３は、ステップＳ３１において、設計変数を離散的に決めるサンプリングを行う（ステップＳ３１１、サンプリング手順）。【００３６】
（ステップＳ５１）目的関数最大化・最小化部１５は、目的関数の最大値又は最小値を探索することで最適点を求める（目的関数最大化・最小化手順、探索手順）。目的関数最大化・最小化部１５は、ステップＳ５１において、目的関数を近似的に算出して応答曲面を作成する（ステップＳ５１１、応答曲面作成手順）。【００３７】
次に、ステップＳ３１１のサンプリング手順について説明する。
図７は、実施例に係るサンプリング手順を示す図である。なお、図７では、説明を単純化するために、外側から２層（６層目Ｌｙ６、７層目Ｌｙ７）の積層角度を最適化する場合を考える。このように、本実施例では、積層のうち、隣接する２つの層（例えば、外側から２層）を選択する。【００３８】
図７において、横軸は６層目の積層角度を示し、縦軸は７層目（最外層）の積層角度を示している。図７において、６層目の積層角度Ａ［°（度）］と７層目の積層角度Ｂ［°（度）］を（Ａ，Ｂ）で表すと、丸印に示す（−９０，−９０）、（−９０，９０）、（０，０）、（９０，−９０）、（９０，９０）がサンプリング点である。なお、積層シェルによる強度解析は１回の解析時間が長いため、なるべく解析回数を減らすことが実用上重要になる。このため、サンプリングは、なるべく少ないサンプリング点で全体を網羅できるようにするのが好ましい。【００３９】
本実施形態では、サンプリング手順で用いる手法の一例として、Ｓｐａｃｅ−ｆｉｌｌｉｎｇ法で５点をサンプリングする。
ここで、Ｓｐａｃｅ−Ｆｉｌｌｉｎｇ法とは、空間をバランスよく埋める手法であり、「各点の距離の最小値」を最大化するアルゴリズムで構成される。設計変数変更部１３は、このように、サンプリングしたサンプリング点を変更した設計変数とする。なお、設計変数変更部１３は、サンプリングの手法としてＳｐａｃｅ−Ｆｉｌｌｉｎｇ法以外の手法を用いてもよい。【００４０】
次に、ステップＳ５１１の応答曲面作成手順について説明する。
図８は、実施例に係る６層目と７層目のサンプリング点を用いて算出した３点曲げ強度に基づく応答曲面の一例を示す図である。
図８において、水平面の第１の軸が第６層の積層角度（ａｎｇｌｅ６）を示し、水平面の第２の軸が第７層の積層角度（ａｎｇｌｅ７）を示し、縦軸が強度［Ｎ］の大きさを示す。なお、積層角度は、−９０°〜＋９０°の範囲とする。
図８に示す例では、６層目の積層角度を−９０°、７層目の積層角度を−９０°にしたときが最大値であり、強度（剛性）が最も大きくなることを示している。【００４１】
ここで、設計装置１が行う計算例を、図６を参照しつつ説明する。
目的関数選択部１１は、記憶部１６が記憶する情報から、３点曲げ強度を選択する（ステップＳ１）。続けて、設計変数選択部１２は、記憶部１６が記憶する情報から、積層角度を選択する（ステップＳ２）。続けて、設計変数変更部１３は、設計変数に対してサンプリングを行う（ステップＳ３１１）。続けて、設計変数変更部１３は、サンプリングした点を、例えば外側から２層（６層目、７層目）それぞれの変更した積層角度とする（ステップＳ３１）。続けて、解析部１４は、積層シェル要素を用いて、例えば有限要素法によって解析を行う（ステップＳ４）。続けて、目的関数最大化・最小化部１５は、ステップＳ４１で解析した結果を用いて応答曲面を作成する（ステップＳ５１１）。続けて、目的関数最大化・最小化部１５は、作成した応答曲面を用いて、目的関数の最大値を探索することで最適点を求める（ステップＳ５１）。
【実施例３】
【００４２】
実施例２の図８に示した応答曲面において、極大値となる６層目の積層角度Ａ［°（度）］と６層目の積層角度Ｂ［°（度）］は、（−９０，−９０）、（−９０，９０）、（９０，−９０）、（９０，９０）である。管状体の積層角度は、−９０°と＋９０°とが同一の積層を示す。このため、４つの極大点の強度（縦軸）は、一致しなければならない。しかしながら、図８に示した例では、４つの極大点の強度（縦軸）が一致していない。
この理由は、サンプリング手順の時点で、単純に空間を埋めるような手法を用いた場合に−９０°と＋９０°が同一であるという周期性が考慮されていないためである。実施例３では、このような周期性を考慮した例を説明する。【００４３】
周期性を考慮するには、サンプリング手順の時点で工夫する必要がある。周期性を考慮する一つの方法として、サンプリング空間を周期的に無限に拡張する方法が考えられる。
すなわち、実施例２では、−９０°と＋９０°付近が２重にサンプリングされていることになり、５点サンプリングしているにもかかわらず、２点分と同等である。そこで、本実施例では、サンプリング空間を拡張することで、２重サンプリングを防ぐ。【００４４】
図９は、実施例に係る５点の場合のサンプリング模式図である。図９において、丸印は実施例２におけるサンプリング手法によって選択されたサンプリング点を示し、Ｘ印が実施例３におけるサンプリング手法によって選択されたサンプリング点を示す。実施例３では、図９に示すように、サンプリング空間を拡張することで、２重サンプリングを防ぐ。【００４５】
次に、サンプリング空間を拡張して周期性を考慮したＳｐａｃｅ−Ｆｉｌｌｉｎｇ法について説明する。定義域の設計変数空間をｖ１として、前記設計変数空間ｖ１を、設計変数空間ｖ２，ｖ３，……に拡張すると、無限大の設計変数空間Ｖが生成する。ｖ１に一つ点Ｎ１を生成すると、拡張された前記設計変数空間ｖ２，ｖ３，……にも同様に点Ｎ２，Ｎ３，……が生成される。設計変数空間ｖ１を拡張して無限大の設計変数空間Ｖとするため、ｖ１に生成された点Ｎ１ｉ、Ｎ１ｊの２点間の距離ｄｉｊは、同一の設計変数空間の中の距離を示すとは限らない。他の設計変数空間にも生成されたＮ２ｊ，Ｎ３ｊ…．．とＮ１ｉとの距離の中で最小の距離が、２点間の距離になる。このため、解析部１４は、前記設計変数空間をＶ中のｎ個の点を一様に分布させるために、次式（１）を解くことでサンプリング点を求める。【００４６】
【数１】
2017195898.tif 000004
【００４７】
なお、式（１）において、ｉは、１層毎のサンプル数を表し、例えば２以上５０以下の値である。２点以上であれば、信頼性が高い結果が得られる。５０点以下であれば、計算時間を短くできる。また、ｊは、積層の総数を表し、例えば１以上１５以下であり、実施形態では７の例を説明している。また、Ｌは、複数のサンプル点のうち、任意の２点間の距離を表す。また、ｆは、関数を表す。
解析部１４は、式（１）において、ｍａｘ（ｍｉｎ（Ｌ））によって新たなサンプル点を求める。ｍａｘ（ｍｉｎ（Ｌ））の意味合いは、最小のＬが最大になるように配置することで、空間全体にバランスよくサンプリング点をばらまくことを意味している。
このように算出されたサンプリング点は、図７または図９である。なお、図７および図９では、サンプリング点の数が５つの例を示したが、数はこれに限られない。サンプリング点の数は、所定の値であってもよく、設計装置１の利用者が設定または変更するようにしてもよい。【００４８】
目的関数最大化・最小化部１５は、算出したサンプリング点を用いて、例えば放射基底関数（ＲＢＦ）によって応答曲面を作成する。または、目的関数最大化・最小化部１５は、応答曲面作成にも周期性を考慮するためにフーリエ級数近似を用いて周期性を担保するとさらによい。目的関数最大化・最小化部１５は、フーリエ級数近似を、次式（２）を用いて行う。【００４９】
【数２】
2017195898.tif 000005
【００５０】
なお、式（２）において、ｎはサンプル数を表す。図７または図９に示した例では、ｎ＝５となる。
図１０は、実施例に係る１次のフーリエ級数近似を行った応答曲面を示す図である。なお、図１０に示す例は、実施例２と同様に、外側から２層（６層目、７層目）を選択した例である。図１０において、水平面の第１の軸が第６層の積層角度（ａｎｇｌｅ６）を示し、水平面の第２の軸が第７層の積層角度（ａｎｇｌｅ７）を示し、縦軸が強度［Ｎ］の大きさを示す。
図１０に示す応答曲面により、解空間全体を知ることができ、強度を最大化するための最適な積層構成（ここでは設計変数を積層角度としたため各層の積層角度）を知ることができる。
【実施例４】
【００５１】
実施例３では、応答曲面作成に周期性を考慮するために１次のフーリエ級数近似を行う例を説明したが、実施例４では２次のフーリエ級数近似を行う例を説明する。
１次近似の場合、全体が平滑化されてしまい局所的な最適点が存在する場合に見つけ出すことができない可能性がある。これを２次近似に変更することで局所的な最適点も見つけることができる。【００５２】
図１１は、実施例に係る２次のフーリエ級数近似を行った応答曲面を示す図である。なお、図１１に示す例は、実施例２と同様に、外側から２層（６層目、７層目）を選択した例である。図１１において、水平面の第１の軸が第６層の積層角度（ａｎｇｌｅ６）を示し、水平面の第２の軸が第７層の積層角度（ａｎｇｌｅ７）を示し、縦軸が強度［Ｎ］の大きさを示す。【００５３】
次に、総積層数が７の場合について説明する。
実施例４では、前述した位置Ｃに対して、次式（３）の制約条件を用いる。【００５４】
【数３】
2017195898.tif 000006
【００５５】
なお、式（３）において、ＥＩは、シャフトの曲げ剛性、Ｅは、ヤング率、Ｉは、シャフトの断面二次モーメントである。ＧＪは、シャフトの捩り剛性、Ｇは、せん断弾性係数、Ｊは、ねじり定数（断面ねじりモーメント）である。また、ｉは、層の順番であり、１層目が１、７層目が７である。また、Ｅ（ａｉ）ｘｘは、Ｅ（ａ１）ｘｘが最内側の層Ｌｙ１のヤング率、Ｅ（ａ２）ｘｘが最内側から２層目の層Ｌｙ２のヤング率、・・・を表す。Ｉａｉは、Ｉａ１が最内側の層Ｌｙ１の断面二次モーメント、Ｉａ２が最内側から２層目の層Ｌｙ２の断面二次モーメント、・・・を表す。Ｅ（ａｉ）ｓｓは、Ｅ（ａ１）ｓｓが最内側の層Ｌｙ１のせん断弾性係数、Ｅ（ａ２）ｓｓが最内側から２層目の層Ｌｙ２のせん断弾性係数、・・・を表す。Ｊａｉは、Ｊａ１が最内側の層Ｌｙ１のねじり定数、Ｊａ２が最内側から２層目の層Ｌｙ２のねじり定数、・・・を表す。また、ｆ、ｇそれぞれは、関数を表す。また、θ１は１層目の積層角度、θ２は２層目の積層角度、・・・、θ７は７層目の積層角度を表す。
なお、ＥＩはシャフトとして機能する程度の硬さを有することを制約条件とし、３５［Ｎ・ｍ２］以上とした。ＧＪはシャフトとして機能する程度の捩れにくさを有することを制約条件とし、１２［Ｎ・ｍ２］以上とした。【００５６】
式（３）に示した制約条件を課したときの位置Ｃにおける最適構成例を説明する。
図１２は、実施例に係る制約条件を課したときの位置Ｃにおける最適構成例を示す図である。
図１２に示すように、最外層から最内層の順に、７層目Ｌｙ７の積層角度が９０°、６層目Ｌｙ６の積層角度が１５°（計算上は１４．１８°だが１５°とした）、５層目Ｌｙ５の積層角度が９０°、４層目Ｌｙ４の積層角度が０°、３層目Ｌｙ３の積層角度が０°、２層目Ｌｙ１の積層角度が−４５°、１層目Ｌｙ１の積層角度が４５°である。
なお、図１２に示した構成は、強度が最大であり、重量が最も軽い構成の例である。【００５７】
同様に、位置Ｂに次式（４）の制約条件を課し、位置Ａに次式（５）の制約条件を課し、位置Ｔに次式（６）の制約条件を課し積層構造を算出した。なお、位置Ｂ、ＡおよびＴそれぞれの位置では、必要な性能が異なるため各位置の内径を１３．０［ｍｍ］、９．８［ｍｍ］、８．４［ｍｍ］とした。なお、上述した位置Ｃの内径は１４．１［ｍｍ］とした。【００５８】
【数４】
2017195898.tif 000007
【００５９】
【数５】
2017195898.tif 000008
【００６０】
【数６】
2017195898.tif 000009
【００６１】
なお、式（４）〜（６）それぞれの値は、ＥＩがシャフトとして機能する程度の硬さを有する内径の大きさに応じた制約条件、ＧＪがシャフトとして機能する程度の捩れにくさを有する内径の大きさに応じた制約条件として設定した値である。【００６２】
式（３）〜（６）に示した制約条件を課したときの位置Ｃ、Ｂ、Ａ、Ｔにおける最適構成例を説明する。
図１３は、実施例に係る制約条件を課したときの位置Ｂ、Ａ、Ｔにおける最適構成例を示す図である。
図１３の符号ｇ１０１が示す領域に示すように、位置Ｂでは、最外層から最内層の順に、７層目Ｌｙ７の積層角度が９０°、６層目Ｌｙ６の積層角度が５°、５層目Ｌｙ５の積層角度が９０°、４層目Ｌｙ４の積層角度が０°、３層目Ｌｙ３の積層角度が０°、２層目Ｌｙ１の積層角度が−４５°、１層目Ｌｙ１の積層角度が４５°である。【００６３】
図１３の符号ｇ１０２が示す領域に示すように、位置Ａでは、最外層から最内層の順に、７層目Ｌｙ７の積層角度が９０°、６層目Ｌｙ６の積層角度が０°、５層目Ｌｙ５の積層角度が０°、４層目Ｌｙ４の積層角度が０°、３層目Ｌｙ３の積層角度が０°、２層目Ｌｙ１の積層角度が−４５°、１層目Ｌｙ１の積層角度が４５°である。【００６４】
図１３の符号ｇ１０３が示す領域に示すように、位置Ｔでは、最外層から最内層の順に、７層目Ｌｙ７の積層角度が０°、６層目Ｌｙ６の積層角度が０°、５層目Ｌｙ５の積層角度が０°、４層目Ｌｙ４の積層角度が０°、３層目Ｌｙ３の積層角度が０°、２層目Ｌｙ１の積層角度が−４５°、１層目Ｌｙ１の積層角度が４５°である。【００６５】
表２に、各位置における比較例と実施例の最適積層の強度評価の結果を示す。表２は、計算から求めたものを実際に作成して、各位置の強度を測定した結果である。このように、実施例によれば、重量が同じであっても、比較例に対して強度が高まることが確認された。【００６６】
【表２】
2017195898.tif 000010
【００６７】
なお、図１２および図１３に示した計算結果において、位置Ｃ、位置Ｂ点における外側から２層目の層Ｌｙ６がそれぞれ１５°、５°である。この設計値をゴルフシャフトに適用すると、ルール違反となる可能性がある。「シャフトが非対称の特性を持つようにデザインされることを効果的に制限」すなわち「故意に左右非対称にしてはならない」というルールがあるためである（例えば、クラブと球についての規則ガイド２０１６年（第５版）、日本ゴルフ協会）。そこで、材料を約半分の厚みのプリプレグＣ（表１）を用いて±１５°、±５°とすることが望ましい。また、比較例と最適積層のシャフト重量は同等であり、剛性は比較例の方が高かった。この場合では、剛性が小さい方が強度に有利であることが分かる。
剛性評価は１［ｋｇ］の錘を用いた片持ち曲げ試験にて行った。ゴルフシャフトの片持ち曲げ試験は、例えば前述したＳＧ基準  ＣＰＳＡ  ００９８に従うものであり、詳細は省略する。【００６８】
＜ゴルフシャフトの一例＞
ここで、強度を最大化したゴルフシャフトの最適積層の例について説明する。設計条件は、位置Ｃの内径が１４．１［ｍｍ］、位置Ｂの内径が１３．０［ｍｍ］、位置Ａの内径が９．８［ｍｍ］、位置Ｔの内径が８．４［ｍｍ］である。
図１４は、強度を最大化したゴルフシャフトの最適積層の構成例を示す図である。【００６９】
本実施形態におけるゴルフシャフトは、マンドレルと呼ばれる芯金に、各プリプレグを巻き付け、加熱硬化後にマンドレルを引き抜き、外層を研磨することで得られる。一般的なゴルフシャフト製造方法と同一である。なお、図１４では、研磨しろを図示していない。また、図１４では、チップ部（位置Ｔ側）に外径調整層が存在するが、図示していない。また、位置Ｃ側はバット部である。計算の結果、各層に用いた材料は、全てプリプレグＥ（表１）であった。【００７０】
最内層であるアングル層３１の位置Ｃ、Ｂ、ＡおよびＴそれぞれの積層角度が＋４５°である。
アングル層３２は、位置Ｃ、Ｂ、ＡおよびＴそれぞれの積層角度が−４５°であり、アングル層３１に積層して構成される。【００７１】
ストレート層３３は、位置Ｃ、Ｂ、ＡおよびＴそれぞれの積層角度が０°であり、アングル層３２に積層して構成される。
ストレート層３４は、位置Ｃ、Ｂ、ＡおよびＴそれぞれの積層角度が０°であり、ストレート層３３に積層して構成される。【００７２】
変化層３５ｃは、積層角度が＋９０°であり、ストレート層３４の位置Ｃに積層して構成される。
変化層３６ｃは、積層角度が＋１５°であり、変化層３５ｃに積層して構成される。
変化層３７ｃは、積層角度が＋９０°であり、変化層３６ｃに積層して構成される。【００７３】
変化層３５ｂは、積層角度が＋９０°であり、ストレート層３４の位置Ｂに積層して構成される。
変化層３６ｂは、積層角度が５°であり、変化層３５ｂに積層して構成される。
変化層３７ｂは、積層角度が＋９０°であり、変化層３６ｂに積層して構成される。【００７４】
変化層３５ａは、積層角度が０°であり、ストレート層３４の位置Ａに積層して構成される。
変化層３６ａは、積層角度が０°であり、変化層３５ａに積層して構成される。
変化層３７ａは、積層角度が＋９０°であり、変化層３６ａに積層して構成される。【００７５】
変化層３５ｔは、積層角度が０°であり、ストレート層３４の位置Ａに積層して構成される。
変化層３６ｔは、積層角度が０°であり、変化層３５ｔに積層して構成される。
変化層３７ｔは、積層角度が０°であり、変化層３６ｔに積層して構成される。【００７６】
なお、各層は、他の隣接する層に巻き付けることで構成するようにしてもよい。
なお、変化層３５ｃ，変化層３５ｂ，変化層３５ａ，変化層３５ｔ、変化層３６ｃ，変化層３６ｂ，変化層３６ａ，変化層３６ｔ、変化層３７ｃ，変化層３７ｂ，変化層３７ａ，変化層３７ｔそれぞれは、各端部が突き合わされるように、すなわち重ね合わせることなく且つ隙間なく配置されるように形成されることが望ましいが、それぞれが単体で巻き付けられてもよい。ここで、図１４の変化層３５ｃ、３５ｂに代表される斜めに突き合わされている部分は応力集中を避けるための逃がし部である。なお、突き合わせた部位に関しては、重ならない限り、例えば、製造上不可避な程度の間隔を有して配置されていてもよい。【００７７】
変化層３５ｃ，変化層３６ｃ，変化層３７ｃそれぞれは、バット端部を起点として０−３２５ｍｍ（以下、位置Ｃ近傍と呼ぶ）配置される。ここで、長さの基準は、逃がし部の中央までとする。変化層３５ｂ、変化層３６ｂ、変化層３７ｂそれぞれは、３２５−６２５ｍｍ（以下、位置Ｂ近傍と呼ぶ）配置される。また、変化層３５ａ，変化層３６ａ，変化層３７ａそれぞれは、６２５−９７５ｍｍ（以下、位置Ａ近傍と呼ぶ）配置される。変化層３５ｔ、変化層３６ｔ、変化層３７ｔそれぞれは、９７５−チップ端まで（以下、位置Ｔ近傍と呼ぶ）に配置される。【００７８】
以上の構成は、以下のＩ〜ＶＩの特徴を備えている。
Ｉ．バット側からチップ側に向かうに従って、順次垂直積層成分が小さくなる。
ここで垂直積層成分とは、最内層のアングル層の２層を除いた各層の積層角度の和と定義した。なお、アングル層とは、シャフト長手方向に対する繊維の配向が＋２０°〜＋７５°、−２０°〜−７５°である層を指す。また、垂直積層成分は絶対値であり、無次元量とした。
位置Ｃ近傍、位置Ｂ近傍、位置Ａ近傍、位置Ｔ近傍の垂直積層成分を算出する場合は、それぞれの領域において、各層の垂直積層成分を加算した値を各層のプリプレグの枚数で除算してから、求めた各層の垂直積層成分を全層に亘って加算して求める。例えば、０°層と１５°層を比較した場合、垂直積層成分は、０°層では０であり、１５°層では１５になるため、垂直積層成分は１５°層の方が大きい。また、位置Ｃ近傍の積層数が３層であり、第一の層が９０°層の１枚であり、第二の層が４５°と１５°の２枚を突き合わせており、第三の層が９０°を２枚突き合わせている場合、垂直積層成分は９０／１＋（４５＋１５）／２＋（９０＋９０）／２＝２１０となる。
ＩＩ．少なくとも１層が積層角度を変えた複数のプリプレグからなる。複数のプリプレグからなる層は、全長に亘る層であってもよく、部分的に配置された層でもよい。【００７９】
ＩＩＩ．位置Ｃ近傍では積層構成が内側から順に、＋４５°、−４５°、０°、０°、９０°、１５°〜４５°、９０°である。垂直積層成分は、１９５〜２２５となる。位置Ｃ近傍は剛性が高いため、曲げによるたわみが少ないうちに、潰れて破壊する。すなわち、潰し剛性が支配的になる。潰し剛性を、疑似的かつ相対的に数値化したものである垂直積層成分が１９５以上であれば、潰し剛性が高くなるため好ましい。垂直積層成分が２２５以下であれば、曲げ剛性が十分であるため、好ましい。【００８０】
ＩＶ．位置Ｂ近傍では積層構成が内側から順に、＋４５°、−４５°、０°、０°、９０°、５°〜１５°、９０°である。垂直積層成分は１８５〜１９５であり、位置Ｃの垂直積層成分よりも小さい。位置Ｂ近傍は、位置Ｃ近傍よりも細く柔らかいため、破壊時に位置Ｃ近傍よりも曲げ成分が大きくなる。垂直積層成分を、位置Ｃ近傍の垂直積層成分よりも小さい前記範囲とすることで、曲げ剛性を高めることができる。【００８１】
Ｖ．位置Ａ近傍では、積層構成が内側から順に、＋４５°、−４５°、０°、０°、０°、０°、９０°である。垂直積層成分は９０である。位置Ｂの垂直積層成分より小さい。仮に位置Ａでの垂直積層成分が大きすぎる場合、剛性不足に陥る可能性高いため、垂直積層成分が１８０未満であることが好ましい。位置Ａ近傍は、位置Ｂ近傍よりも細く柔らかいため、破壊時に位置Ｂ近傍よりも曲げ成分が大きくなる。そのため、曲げ剛性を高めるために、垂直積層成分を位置Ｂ近傍の垂直積層成分よりも小さい上記範囲とすることが好ましい。【００８２】
ＶＩ．位置Ｔ近傍では、積層構成が内側から順に、＋４５°、−４５°、０°、０°、０°、０°、０°である。垂直積層成分は０である。位置Ｔ近傍では、破壊はほぼ曲げによって起こるため、垂直積層成分は０が好ましい。【００８３】
上記特徴は、バット側ほど曲げ剛性が高いため、３点曲げ実施時には曲げ変形はあまり生じず、潰しで壊れやすいことに依存する。すなわち、円筒形のシャフトが楕円に変形して壊れる。そのため、潰し剛性を高めるための９０°層を多く配置すべきという結論に達する。
一方、チップ側は曲げ剛性が低く、３点曲げ実施時には曲げ変形が大きく生じ、潰し変形はあまり生じない。よって、曲げ剛性を高めるための０°層を多く配置すべきという結論に達する。【００８４】
以上のように、本実施形態の設計手法により、上述した構成（図１４）を有するシャフトが得られた。なお、ここでは層数は７層としプリプレグＥで構成したが、上記特徴を満たせば同じように強度を高めることが可能となる。各変化層も同様に調節可能である。設計装置１が設計したシャフトは概ね３０［ｇ］以下であった（実測は２９．５［ｇ］であった）。剛性は１［ｋｇ］の錘を用いた片持ち曲げ試験にて９０［ｍｍ］以上になるよう設計した。片持ち曲げ試験の実測値は９１［ｍｍ］であった。【００８５】
以上のように、実施形態によれば、設計したい目的関数（例えば剛性）が最大又は最小となる設計変数（例えば材料）を求めることができる。
また、実施形態によれば、設計したい目的関数（例えば３点曲げ強度）が最大又は最小となる設計変数（例えば積層角度）を層毎に求めることができる。【００８６】
また、実施形態によれば、Ｓｐａｃｅ−Ｆｉｌｌｉｎｇ法等を用いたサンプリング手順を行うことで、少ないサンプリング点で全体を網羅でき、なるべく少ないサンプリング点で全体を網羅でき、解析時間を短縮することができる。
また、実施形態によれば、フーリエ級数近似によるサンプリング手順において、サンプリング空間を拡張したので、２重サンプリングを防ぐことができる。
また、実施形態によれば、サンプリング手順において、フーリエ級数近似が一次の場合、解空間全体を知ることができ、強度を最大化するための最適な積層構成（ここでは設計変数を積層角度としたため各層の積層角度）を知ることができる。また、フーリエ級数近似が二次の場合、局所的な最適点も見つけることができる。【００８７】
また、実施形態によれば、積層のうち少なくとも１つの層が複数の領域（例えば、位置Ｃ近傍、位置Ｂ近傍、位置Ａ近傍、位置Ｔ近傍）を備えるため、層毎かつ領域毎に設計したい目的関数が最大又は最小となる設計変数を求めることができる。これにより、実施形態によれば、例えば、重量が同じであっても、比較例に対して強度が高めることができる。
【産業上の利用可能性】
【００８８】
軽量化や高強度化に応用することができ、複数の繊維強化樹脂層からなる管状体に特有の課題を解決できる。
【符号の説明】
【００８９】
１…設計装置、１１…目的関数選択部、１２…設計変数選択部、１３…設計変数変更部、１４…解析部、１５…目的関数最大化・最小化部、１６…記憶部、Ｓ１…目的関数選択手順（解析手順）、Ｓ２…設計変数選択手順（解析手順）、Ｓ３、Ｓ３１…設計変数選択手順（解析手順）、Ｓ４…解析手順、Ｓ５、Ｓ５１…目的関数最大化・最小化手順（探索手順）、Ｓ３１１…サンプリング手順、Ｓ５１１…応答曲面作成手順、３１、３２・・・アングル層、３３、３４・・・ストレート層、３５、３６、３７、３５ａ、３６ａ、３７ａ、３５ｂ、３６ｂ、３７ｂ、３５ｃ、３６ｃ、３７ｃ、３５ｔ、３６ｔ、３７ｔ・・・変化層
【図１】
2017195898.tif 000011
【図２】
2017195898.tif 000012
【図３】
2017195898.tif 000013
【図４】
2017195898.tif 000014
【図５】
2017195898.tif 000015
【図６】
2017195898.tif 000016
【図７】
2017195898.tif 000017
【図８】
2017195898.tif 000018
【図９】
2017195898.tif 000019
【図１０】
2017195898.tif 000020
【図１１】
2017195898.tif 000021
【図１２】
2017195898.tif 000022
【図１３】
2017195898.tif 000023
【図１４】
2017195898.tif 000024
【国際調査報告】
2017195898.tif 000025
2017195898.tif 000026
2017195898.tif 000027