TOP > 国内特許検索 > 荷重伝達経路法に基づく数値構造解析装置

荷重伝達経路法に基づく数値構造解析装置

国内特許コード P09S000301
掲載日 2010年1月22日
出願番号 特願2007-542832
登録番号 特許第4310501号
出願日 平成18年11月6日(2006.11.6)
登録日 平成21年5月22日(2009.5.22)
国際出願番号 PCT/JP2006/322069
国際公開番号 WO 2007/052784
国際出願日 平成18年11月6日(2006.11.6)
国際公開日 平成19年5月10日(2007.5.10)
優先権データ
  • 特願2005-321695 (2005.11.7) JP
発明者
  • 高橋 邦弘
  • 櫻井 俊彰
出願人
  • 慶應義塾大学
発明の名称 荷重伝達経路法に基づく数値構造解析装置
発明の概要 【要約】荷重伝達経路法に基づく数値構造解析装置において、計算時間を大幅に短縮することを目的とする。解析対象構造物の支持点Bを固定し、特定負荷点Aに荷重をかけるようにパラメータを設定する。剛性行列保持手段1の全体剛性行列に基づいて、有限要素法計算手段2で、解析対象構造物の変形を計算して各点の変位量などの基本データを求める。特定負荷点Aと支持点Bを固定して、変化負荷点Cに3通りの検査荷重を与え、それぞれの変形を有限要素法計算手段2で計算して変位量を求める。部分剛性行列計算手段3で、解析対象構造物の内部剛性行列と荷重値と変位量に基づく多元連立一次方程式を解き、部分剛性行列KACを求める。剛性指標計算手段8で、部分剛性行列KACと基本データの変位量などから剛性指標U*の値を計算する。解析対象構造物の必要なすべての点を順次たどるように、変化負荷点Cを変更して各点のU*の値を計算する。
従来技術、競合技術の概要 【背景技術】構造設計においては、高剛性化や高強度化や軽量化などを目標にする。これらの目標を達成するためには、構造全体に関する検討が必要である。構造の強度剛性を論ずる際には、応力あるいはひずみを用いるのが普通である。しかし、構造設計においては、むしろ力の伝わる経路あるいは荷重の伝達に着目して構想をたてる必要がある。したがって、構造内部における力の伝達経路を知ることは重要な課題である。力の伝わり方に関して、従来は、応力による推定がなされてきた。応力とは、単位面積あたりに加わる力を示しており、応力の高い箇所に力が強く伝わると考えることもできる。応力は、実験的な測定あるいはコンピュータシミュレーションによって解析されてきた。主応力の値を等応力面に直交する矢印で示し、その分布図を構造物の図中に描き、矢印の分布によって力の伝わり方を表現することがしばしば行われる。しかし、この考え方には、いくつかの問題がある。応力集中により、力が有効に伝わっていなくても、大きな応力の生ずることがあり、応力を用いて力の伝達を考えると、誤った結論を導くことがある。例えば、構造内部に小さな円孔があると、その部分に応力集中が発生し、大きな主応力が現れる。しかし、円孔を設けることで、力を強く伝えることができると考えることは奇妙であるから、その箇所で大きな荷重の伝達が行われるようになったと考えるべきではない。もしそうであるなら、力を強く伝達させるためには円孔をあける、という誤った結論を導いてしまう。そうではなく、円孔をあけると、円孔からやや離れた部位で力を伝えるようになる、と結論されなくてはならない。そのためには、応力と異なる指標が必要である。そこで、本発明者らは、力の伝達とその経路に関するU*なる指標を提案した。U*は、応力と全く異なる発想に基づく指標である。その指標は直感的な概念に従って導いたものであるが、U*と力の伝達とが密接に関連していることを、U*の具体的計算例を通して明らかにし、構造設計におけるU*の有効性を示すことができる。さらに、U*ポテンシャルと剛性線ならびに剛性低下ベクトルなる概念を導入し、荷重伝達経路をU*で厳密に表現する方法を説明する。以下に、U*の概要を説明するが、具体的計算例などの詳細については、非特許文献1~9を参照されたい。U*とは、構造に力が加わったときに、構造内部の任意の1点を固定したときとしないときの反力の比から求められる。U*は、負荷点と任意点の結合の強さを表す指標であり、構造内部の荷重伝達を表す指標である。U*はポテンシャル量であり、U*分布を求めることで、構造において力を伝える向きを示すことが可能である。U*分布の勾配を求めることで、荷重伝達の方向である剛性低下ベクトルを求めることができ、剛性低下ベクトルを結ぶことで、力が伝達される道筋である剛性線を導出することができる。負荷点から支持点に向かう剛性線の中で、最も強く荷重を伝える経路を、荷重伝達経路と定義する。この荷重伝達経路を用いることで、構造内部の荷重伝達を定量的に評価することが可能となる。すなわち、指標U*を用いて、構造内部における荷重伝達の様子を定量的に表現し、構造内部のU*分布から、内部の主要な荷重伝達経路を見出す。得られた荷重伝達経路に対して評価を行うことにより、構造の改善点を明らかにすることが可能である。U*を用いて構造最適化を行ったところ、荷重伝達のみに着目したにもかかわらず、一般に行われているひずみエネルギー一様を目的とした最適解と同等の解を得ており、有効性を確認している。また、有限要素法を用いたコンピュータシミュレーション以外にも、実車車体を用いたU*測定を行っている。この構造解析の方法を、荷重伝達経路法とよぶことにする。U*の概念と定義を説明する。図8は、任意の弾性構造をモデル化したものである。図8(a)では、負荷点をAとし、支持点をBとし、構造物内部の任意の点をCとする。図8(b)は、任意点Cを固定した場合である。図8(a)と図8(b)では、共に負荷点Aに同一の強制変位dAを与えている。図8(a)の場合に、内部に貯えられたひずみエネルギーをUとする。図8(b)の場合の内部に貯えられたひずみエネルギーをU'とする。両者の比U'/Uの値は必ず1以上となり、内部の点Cと負荷点Aとの結合が強い程、この比は高い値を示すと考えられる。例えば、一本の線形ばねを考えたときには、点Cを変数としたU'/Uの分布は、ばねの軸に沿った曲線となる。U'/Uの代わりに、式(1)のように、(U'-U)/U'を使うと、この値は、負荷点Aで1をとり、支持点Bで0となる、一様に減少する直線で表示することができる。 U*≡(U'-U)/U' =1-(U'/U)-1・・・・・(1)式(1)がU*の定義であり、負荷点Aと点Cとの結合の強さを示していると考えられる。U*は、点Cの座標の関数である。荷重伝達経路について説明する。U*の分布は、数値計算あるいは実験によって求めることが可能である。構造物内のU*分布を等高線で表現したとき、その稜線に相当する曲線は、負荷点との結合の強い部分を順次結んだ線に相当している。力は剛性の高い部分をたどりつつ伝達されると考えることができるので、荷重は、この稜線を経路として伝達されるということができる。この稜線すなわち荷重伝達経路Sを、模式的に図8(c)に示す。荷重伝達経路に関して、好ましい構造が満たすべき3つの条件を以下にあげる。(1)荷重伝達経路におけるU*分布の一様性:荷重伝達経路におけるU*分布が、図9(a)に示される理想直線に近づくこと。(2)荷重伝達経路における内部剛性の連続性:荷重伝達経路におけるU*の曲率分布が、図9(b)のように理想直線に近づくこと。(3)荷重伝達経路の合致性:負荷点から支持点に向かう荷重伝達経路S1と負荷点と支持点を入れ替えたときに得られる荷重伝達経路S2が、図9(c)に示されるように合致すること。構造と荷重と変位の関係を説明する。構造を、3点A、B、Cで代表する。各点の荷重と変位の関係を示すと、式(2)となる。
【数1】・・・・・(2)
添字付きの各Kは、3×3の内部剛性行列である。添字付きのpは、荷重を表す列ベクトルである。添字付きのdは、変位を表す列ベクトルである。一見、有限要素法の初等的表記のように見えるが、そうではない。この式(2)で、3点に関する構造全体の効果を表現している。なお、荷重がゼロの時でも、剛体的な同一変位を許容する条件から、式(3)等の関係も成立する。 KAA=-(KAB+KAC)・・・・・(3)なお、K、p、dなどは行列やベクトルであるが、太字で表記できないので、通常表記とする。式(1)によるU*の定義を、式(2)によって表現する。Uは、図8(a)の支持点Bのみを固定した場合の変形エネルギーである。図8(a)において、支持点Bは固定されているので、dBは0となる。これより、式(2)を用いて、pAは、(KAAA+KACC)と表すことができる。この関係を用いて、図8(a)における変形エネルギーUは、式(4)で表される。 U=(1/2)pA・dA =(1/2)(KAAA+KACC)・dA・・・・・(4)積の演算には、ベクトルおよびテンソルの表記を用い、記号(・)は、ベクトルの内積を示す。一方、U'は、図9(b)の構造内の任意点Cを拘束した場合に貯えられる変形エネルギーである。このときも、負荷点Aの強制変位は変えず、dAのままとするから、そこに加わる力p'Aは、式(2)により、KAAAのみとなるので、U'は、式(5)で表される。 U'=(1/2)p'A・dA =(1/2)(KAAA)・dA・・・・・(5)式(4)と式(5)より、式(6)を得る。 (U'/U)={(KAAA)・dA}/{(KAAA+KACC)・dA}・・・・・(6)式(6)における分母はUであるため一定値である。U'/Uは、分子のKAAに依存していることが分かる。このKAAは、式(3)により、-(KAB+KAC)であるから、負荷点Aと任意点C間の結合の強さKACだけを示しているとはいえない。そのため、次のように式を変形する。式(4)と式(5)により、式(7)と書くことができる。 U'=U-(1/2)(KACC)・dA・・・・・(7)したがって、(U'/U)は、式(8)のように表現することができる。 (U'/U)=1-{(KACC)・dA}/(2U)・・・・・(8)さらに、式(1)のU*の定義に対して、式(8)を用いれば、最終的に、U*は、式(9)のようになる。 U*=1-{1-{(KACC)・dA}/(2U)}-1 ={1-(2U)/{(KACC)・dA}}-1 ={1-(pA・dA)/{(KACC)・dA}}-1・・・・・(9)式(9)のように表現することによって、U*は、負荷点Aと任意点Cの間の剛性を表す行列KACに依存した量であるということが明らかとなった。すなわち、任意点CのU*は、負荷点Aとの結合の強さに関する量であるといえる。U*の値の示す意味は、式(9)によって明らかにされたといえるであろう。なお、以上は、図8のモデルに基づく結果であるが、支持点が2個または複数個、さらには、ある領域が固定された場合についても、同様な結果が得られる。いずれの場合も、任意点Cを固定した場合と、固定しない場合とのエネルギー比によって定義されたU*は、式(9)と同一の形で表現され、負荷点Aと任意点C間の内部剛性行列KACによって表すことができる。特に、このKACを他の内部剛性行列と区別して呼ぶために、部分剛性行列ということにする。U*ポテンシャルと経路について説明する。このU*は、構造物の内部のあらゆる点において与えることができる。各点において定まるU*をポテンシャルと考え、U*ポテンシャルと呼ぶ。等ポテンシャル線に直交する曲線を、剛性線と名付ける。剛性線に沿う勾配gradU*は、剛性の低下量を表現している。この低下量を、式(10)で表す。 k(-)=-gradU*・・・・・(10)k(-)を剛性低下ベクトルと名付ける。k(-)は保存量である。U*の等高線と荷重経路について説明する。一般に、スカラー量分布の等高線に関する稜線とは、等高線に直交する線群の中で、最も勾配の小さい曲線として定義される。剛性線群の中で最も勾配の小さい曲線、すなわち、最も小さい剛性低下ベクトルをつないだ曲線が稜線である。力は剛性低下量の少ない部位をたどって伝わると考えるならば、この稜線を、荷重伝達経路と呼ぶことができる。荷重伝達経路をU*分布の稜線によって考えることにより、円孔がある場合には、円孔からやや離れた部位で力を伝えているということが結論される。U*を定義通りに計算する具体的計算法を説明する。定義どおりのままでU*を計算する場合には、次の手順が要る。(1)任意点Cを自由にしたまま、負荷点Aに一定の変位を与えるのに必要なエネルギーUを求める。(2)任意点Cを固定して、負荷点Aに一定の変位を与え、その時に必要なエネルギーを、有限要素法等の構造計算により求める。得られた値を式(1)に代入して、その点に関するU*を求める。(3)任意点Cの位置を順次変えて、この計算を繰り返す。これにより、全体にわたるU*の分布を求めることができる。U*を計算で求めるには、有限要素法(FEM)によるのが普通である。有限要素法の解法には、大きく分けて変位法と応力法がある。汎用の有限要素法計算ソフトでは変位法が用いられている。変位法においては、計算毎に別の荷重を与えて繰り返し計算を行う場合、すなわち荷重条件(力学的な境界条件)を変化させて繰り返し計算を行う場合には、全体剛性行列に関する連立方程式を1回解けばよい。しかし、計算毎に別の変位を与えながら繰り返し計算を行う場合、すなわち、変位条件(幾何学的な境界条件)を変化させて繰り返し計算を行う場合には、全体剛性行列に関する連立方程式を毎回繰り返して解く必要がある。そのため、繰り返し数の多い計算では、計算量が非現実的となる。有限要素法では、1解析中、全体剛性行列に関する連立方程式の計算にほとんどの時間を費やすため、繰り返し数の多い計算では、計算量が非現実的となる。U*の繰り返し計算では、任意点Cを拘束する。つまり、毎回異なる位置に選択した任意点Cにゼロの変位を与える。このように、幾何学的な条件を毎回変化させる必要があるので、そのままでは非現実的な計算量になる。全体剛性行列に関する連立1次方程式は、変数の数が何百万にもなる。何百万元の方程式を解くことになるから、1回解くだけでもかなりの時間を要する。しかも、U*分布を求めるには、それをまた何百万回も繰り返すので、その繰り返す回数を減らすことが重要である。U*を求める代わりに、KACなる行列を求めればよいことを説明する。KACは、負荷点Aと任意点Cの間のばね定数と考えることができる行列である。U*は、式(1)で与えられる。 U*=(1-(U'/U))・・・・・(1)式(1)は、次の式(9)のようにも表すことができる。 U*={1-(2U)/{(KACC)・dA}}-1・・・・・(9)ここで、Uは、任意点Cを固定しない場合の変形エネルギーである。変位dCと変位dAは、共に、任意点Cを固定しない場合の任意点Cと負荷点Aの変位である。そのため、UとdCとdAは、最初に一回計算しておけばよい。式(9)からU*を求めるには、KACなる行列を各点で求めればよいことになる。つまり、U*を求めるために必要な計算は、行列KACを各点で繰り返して求める計算でよいことになる。行列KACとは何かを説明する。そのために、構造の荷重と変位を表現する2種類の方法を述べる。まず、全体剛性行列による構造物の表現方法を説明する。有限要素法における有限要素に相当するメッシュ要素に分けて構造を表現する。要素数は100万に及ぶ例もある。これらの要素の各節点に順番に番号を付け、そこに加わる力と変位を、添え字つきのpとdで表す。いずれもベクトル量である。このとき、これらの力と変位が正比例関係にある場合に限って扱うとする。つまり、線形な力、変位関係を有している場合を扱うとする。これは、式(11)で表現できる。
【数2】・・・・・(11)
行列の成分は省略して記述してある。節点の数nは極めて大きくなることがあるので、この行列の成分はn×nの膨大な数となる。しかも、pもdもベクトルであるから3個のスカラー量から成っており、式(11)の行列の1個の成分は、3×3個のスカラー量から成っているから、実際にはさらに多くの成分から成っている。この行列は、有限要素法における剛性行列としてよく知られている。有限要素法においては、剛体変位を考慮した上で、この連立方程式を解くことが主要な計算となる。成分の数が多いので、この計算は時間を要する。通常は、汎用計算ソフトを用いて行われる。この剛性行列によって、構造物の荷重、変位の関係を表す方法が、第1の表現方法である。次に、負荷点Aと支持点Bと任意点Cを与えた場合において、3点のみにより構造物を表現する方法を説明する。ある構造を、任意点Cに関して、3点A、B、Cのみで表現した場合、有限要素のメッシュは無関係となり、単に3点に関する荷重、変位の関係になる。これらの関係式は簡単で、上述の式(2)で表される。この表現における添え字つきのKは、内部剛性行列と呼ばれる。この行列の成分は、式(2)に見るように、3×3である。ただし、その1個の成分は、また3×3個のスカラー量から成っている。すなわち、式(2)の行列は、9×9のスカラー量から成っている。この内部剛性行列によって、構造物の荷重、変位の関係を表す方法が、第2の表現方法である。この式(2)におけるKACが、式(9)におけるKACである。すなわち、このKACを求めれば、U*が式(9)により求まることになる。以上の説明に関する詳細は、以下に概要を示す本発明者らによるU*に関する論文を参照されたい。非特許文献1には、荷重伝達経路の概念における望ましい構造の条件が報告されている。新規な指標U*を導入して、全体構造の荷重伝達経路を表現する。構造中の負荷点と内部の任意の点との間の結合の度合いは、この指標を用いて定量的に表現できる。この概念に基づいて、望ましい構造の3つの条件、一様性と連続性と経路の合致性が導かれた。目的関数として、これらの3つの条件を導入して、数値計算で構造的最適化が計算された。応力や歪の概念を全く導入していないにもかかわらず、得られた構造は、合理的な形をしている。最終的に、単純構造の荷重伝達経路が、実験的に測定された。これらの値は、この指標が効果的に利用できることを示している。非特許文献2には、荷重伝達経路の概念を用いた大型トラック車室フロア部の振動低減に関する検討が報告されている。大型トラックの車室構造について、静的負荷のもとでの荷重伝達経路を検討した。その結果を、車室の振動低減に役立てる。単純モデルを用いた予備的なシミュレーションでは、床板振動は、床構造の前面交差部材に密接に関連している。有限要素法による荷重経路解析によると、低剛性の前面交差部材のために、荷重経路にはいくらかの不連続性と非一様性がある。非特許文献3には、大型トラックの車室床についてのADAMSによる振動解析とNASTRANによる構造評価が報告されている。大型トラックの車室構造の荷重伝達経路を、静的荷重の下で調べた。その結果を応用して、車室の振動を低減する。単純なモデルを使ったADAMS/Vibrationによる予備的シミュレーションでは、床板の振動は、床構造の前部交差部材の堅さと密接に関係している。NASTRANによる有限要素法での荷重伝達経路解析により、荷重伝達経路には、前部交差部材において不連続性とばらつきがあり、その部材の剛性を減殺している。非特許文献4には、大型トラックキャブのフロア振動に対するADAMSの適用とNASTRANによる負荷経路のU*解析が報告されている。大型トラックを軽量化し安価にすることは、構造設計において重要である。数値解析は、この点で重要な役割を果たしている。ADAMS/Vibrationによる予備的シミュレーションでは、床板の振動は、前部交差部材の剛性に密接に関係している。MSC/NASTRANによる荷重伝達経路解析により、荷重伝達経路には、前部交差部材において不連続性とばらつきがある。非特許文献5には、構造解析シミュレーションにおける荷重伝達の新たな表現法とその適用例が報告されている。新規な指標U*を導入して、荷重伝達経路を表現する。円孔のある板構造の荷重経路を、U*を適用して解析する。応力集中により孔の縁に強い応力がかかっていても、荷重伝達経路は、孔の縁を避けている。この結果は、直感的な予測と合致している。U*解析の計算法は、より複雑な境界条件を持つ構造に拡張できる。この方法の効果は、実際の大型トラック車室の有限要素法モデルを用いて証明された。非特許文献6には、前面衝突における乗用車室の荷重経路最適化とU*構造解析が報告されている。新規な指標U*を導入して、構造における荷重伝達を表現する。車室床構造のU*解析で、2つの場合を調べた。第1の場合として、U*の3つの条件が、目的関数として導入され、遺伝的アルゴリズムによる構造最適化に適用された。遺伝的アルゴリズム計算で得られた床構造は、独特の形をしている。部材は、床下部材の前面部材と側面部材の後部とに結合している。第2の場合として、衝突の際の床構造のU*の値と荷重経路が、PAM-CRASHで計算された。衝突が進行するとき、床下部材は第1の荷重伝達経路となる。衝突の最終段階では、床下部材が主な役割を担う。非特許文献7には、荷重伝達経路理論を適用した大型トラックの車室内の振動低減が報告されている。この目的は、静的荷重がかかった大型トラックの車室構造における荷重伝達経路を調べて、その結果を利用して車室の振動を低減することである。単純モデルを使った予備シミュレーションでは、床板振動は、床構造の前面交差部材の剛性に密接に関係している。有限要素法を使った荷重経路U*解析によると、前面交差部材の剛性の低さは、荷重経路における不連続性とばらつきに起因している。非特許文献8には、構造物内部における荷重伝達経路の新たな表現方法が報告されている。構造物における荷重伝達経路を表現するための指標U*がある。U*の値は、負荷点と内部の任意の点との間の結合度を表すものであることが、内部剛性行列からわかる。剛性場と、剛性線と、剛性低下ベクトルは、U*ポテンシャル線を使って定義される。荷重伝達経路は、剛性低下ベクトルが最小となる剛性線として表現できる。U*解析の適用例として、単純モデル構造について、有限要素法で計算した。U*の値の分布からみると、負荷点と支持点の間の対角部材は、荷重伝達において重要な役割を果たしている。非特許文献9には、実車大型トラックキャブ構造および縮尺模型を用いた荷重伝達U*解析の実験的検討が報告されている。U*の値の分布は、構造中における荷重伝達経路を示すことが知られている。実車大型トラックキャブ構造と、プラスチック縮尺模型について、2つの実験的なU*測定方法を開発した。これらの方法では、従来の方法とは違って、各部材の剛性データは必要ない。有限要素法では、U*解析において重要な役割を果たす実際の板の曲率の効果を表現できない。プラスチック縮尺模型を使うことで、U*の値の分布から強化の効果を直ちに測定できる。
【非特許文献1】Kunihiro Takahashi: Conditions for desirable structures based on a concept of load transfer courses, pp.699-702, 2001.
【非特許文献2】櫻井俊彰、星野裕昭、高橋邦弘:「荷重伝達経路の概念を用いた大型トラック客室フロア部の振動低減に関する検討」、社団法人自動車技術会学術講演会前刷集、No.36-02, pp.5-8, 2002.
【非特許文献3】Hiroaki Hoshino, Toshiaki Sakurai, Kunihiro Takahashi: Application of ADAMS for Vibration Analysis and Structure Evaluation By NASTRAN for Cab Floor of Heavy-Duty Truck, The 1st MSC.ADAMS European User Conference, London, November, 2002.
【非特許文献4】櫻井俊彰、星野裕昭、阿部正俊、高橋邦弘:「大型トラックキャブのフロア振動に対するADAMSの適用とNASTRANによる負荷経路U*の解析」(MSC.ADAMSユーザーコンファレンス2002).
【非特許文献5】櫻井俊彰、阿部正俊、翁祖栄、高橋邦弘:「構造解析シミュレーションにおける荷重伝達の新たな表現法とその適用例」、計算工学会計算工学講演会論文集、Vol.8, pp.401-404, 2003.
【非特許文献6】Toshiaki Sakurai, Junichi Tanaka, Akinori Otani, Changjun Zhang, Kunihiro Takahashi: Load Path Optimization and U* Structural Analysis for Passenger Car Compartments under Frontal Collision, International Body Engineering Conference 2003, pp.181-186, JSAE 20037007, SAE 2003-01-2734, 2003.
【非特許文献7】Hiroaki Hoshino, Toshiaki Sakurai, Kunihiro Takahashi: Vibration reduction in the cabins of heavy-duty trucks using the theory of load transfer paths, pp.165-171, JSAE Review Vol.24, 2003.
【非特許文献8】高橋邦弘、櫻井俊彰:「構造物内部における荷重伝達経路の新たな表現方法」、日本機械学会論文集(A編)Vol.71 No.708, pp.1097-1102, 2005.
【非特許文献9】井上健吾、一木優一郎、松田育磨、櫻井俊彰、石井英章、野原哲男、星野裕昭、高橋邦弘:「実車大型トラックキャブ構造および縮尺模型を用いた荷重伝達U*解析の実験的検討」、社団法人自動車技術会学術講演会前刷集、No.90-04, pp.27-30, 2004.
産業上の利用分野 本発明は、数値構造解析装置に関し、特に、有限要素法を利用して荷重伝達経路法により解析対象構造物の構造解析を行う数値構造解析装置に関する。
特許請求の範囲 【請求項1】解析対象構造物の全体剛性行列に基づいて有限要素法により解析対象構造物の変形を計算する有限要素法計算手段を備え、荷重伝達経路法により解析対象構造物の構造解析を行う数値構造解析装置において、解析対象構造物の必要なすべての点を順次たどるように変化負荷点を変更する位置変更手段と、特定負荷点と支持点を固定し、前記位置変更手段で選択した変化負荷点に3つの独立の検査荷重を順次かけるように計算パラメータを設定する検査荷重設定手段と、前記検査荷重をかけた状態について前記有限要素法計算手段によりそれぞれ変形を計算して前記変化負荷点の変位量と前記特定負荷点の荷重値を求める手段と、解析対象構造物の内部剛性行列と前記変化負荷点の変位量と前記特定負荷点の荷重値とに基づいて未知数が9個以下の多元連立一次方程式を解いて部分剛性行列KACを求める連立方程式計算手段とを備えたことを特徴とする数値構造解析装置。
【請求項2】解析対象構造物の全体剛性行列に基づいて有限要素法により解析対象構造物の変形を計算する有限要素法計算手段を備え、荷重伝達経路法により解析対象構造物の構造解析を行う数値構造解析装置において、支持点に関する要素を除いた全体剛性行列の逆行列を計算する第1逆行列計算手段と、解析対象構造物の必要なすべての点を順次たどるように変化負荷点を更新する位置更新手段と、前記逆行列から特定負荷点と前記変化負荷点に関する要素を抽出して小行列とする小行列抽出手段と、前記小行列の逆行列を計算して部分剛性行列KACを得る第2逆行列計算手段とを備えたことを特徴とする数値構造解析装置。
【請求項3】前記部分剛性行列KACを利用してスカラー化した剛性指標の値を計算する剛性指標計算手段を備えたことを特徴とする請求項1または2記載の数値構造解析装置。
【請求項4】解析対象構造物の全体剛性行列に基づいて有限要素法により有限要素法計算手段で解析対象構造物の変形を計算して、荷重伝達経路法により解析対象構造物の構造解析を行う数値構造解析方法において、特定負荷点と支持点を固定し、解析対象構造物の必要なすべての点を順次たどるように変化負荷点を位置変更手段で変更し、選択した変化負荷点に3つの独立の検査荷重をかけるように検査荷重設定手段で計算パラメータを設定し、前記検査荷重をかけた状態について前記有限要素法計算手段で有限要素法によりそれぞれ変形を計算して前記変化負荷点の変位量と前記特定負荷点の荷重値を求めて、解析対象構造物の内部剛性行列と前記変化負荷点の変位量と前記特定負荷点の荷重値とに基づいて連立方程式計算手段で未知数が9個以下の多元連立一次方程式を解いて部分剛性行列KACを求めることを特徴とする数値構造解析方法。
【請求項5】解析対象構造物の全体剛性行列に基づいて有限要素法により有限要素法計算手段で解析対象構造物の変形を計算して、荷重伝達経路法により解析対象構造物の構造解析を行う数値構造解析方法において、前記全体剛性行列から支持点に関する要素を除いた行列の逆行列を第1逆行列計算手段で計算し、解析対象構造物の必要なすべての点を順次たどるように変化負荷点を位置更新手段で更新し、前記逆行列から特定負荷点と前記変化負荷点に関する要素を小行列抽出手段で抽出した小行列の逆行列を第2逆行列計算手段で計算して部分剛性行列KACを求めることを特徴とする数値構造解析方法。
【請求項6】前記部分剛性行列KACを利用してスカラー化した剛性指標の値を剛性指標計算手段で計算し、解析対象構造物の必要なすべての点を順次たどるように位置更新手段で前記変化負荷点を更新して、全ての変化負荷点について剛性指標の値を計算することを特徴とする請求項4または5記載の数値構造解析方法。
【請求項7】請求項4~6のいずれかに記載の数値構造解析方法をコンピュータで実行するための処理手順を記述したコンピュータプログラム。
【請求項8】請求項7記載のコンピュータプログラムを記録した記録媒体。
産業区分
  • 計算機応用
国際特許分類(IPC)
Fターム
画像

※ 画像をクリックすると拡大します。

25924_01SUM.gif
出願権利状態 権利存続中
ライセンスをご希望の方、特許の内容に興味を持たれた方は、下記問合せ先にご相談下さい。


PAGE TOP

close
close
close
close
close
close
close