立体骨組や平面骨組の解析プログラム
【課題】従来の，柱梁接合部の応力や変形をも含めて解析する骨組解析法では，柱梁接合部と柱部材または梁部材が特定の位置関係にないと解析できず，かつ未知数の混入という解法上の問題があり，かつ立体骨組には適用できない，という課題があった。
【解決手段】以下の手段を併用する。（ア）節点モーメントに対応する外力項をパネルモーメントに対応する外力項で置き換える。（イ）パネルモーメントを応力でなく，自由度間の一対の外力として扱う事で，骨組全体の方程式中の対応外力項を既知数に変化させる。（ウ）柱梁接合部と柱部材または梁部材の位置関係が，特定の場合に制約されず，かつ立体骨組にも適用可能な汎用性のある変換マトリクスを作成する。これら手段の併用により，柱梁接合部の応力や変形をも含めて任意形状の立体骨組または平面骨組の解析が可能になる。
【解決手段】以下の手段を併用する。（ア）節点モーメントに対応する外力項をパネルモーメントに対応する外力項で置き換える。（イ）パネルモーメントを応力でなく，自由度間の一対の外力として扱う事で，骨組全体の方程式中の対応外力項を既知数に変化させる。（ウ）柱梁接合部と柱部材または梁部材の位置関係が，特定の場合に制約されず，かつ立体骨組にも適用可能な汎用性のある変換マトリクスを作成する。これら手段の併用により，柱梁接合部の応力や変形をも含めて任意形状の立体骨組または平面骨組の解析が可能になる。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は，柱や梁などの骨組部材を線材に置換して，立体骨組の応力や変形を計算するプログラムに関し，特に，柱梁接合部の変形をも計算するプログラムに関する。
【背景技術】
【０００２】
骨組解析では，節点を設けて，節点という点への作用力として，並進力やモーメントを考えている。一方，骨組中の柱梁接合部の応力は，例えば鋼構造骨組では非特許文献１に示されているパネルモーメントで評価をされる。パネルモーメントは，節点モーメントとは異なっていて，例えば，梁フェイスモーメント・柱せん断力・接合部寸法から計算されるか，または柱フェイスモーメント・梁せん断力・接合部寸法から計算されるが，いずれも複数部材の応力の関数になっている。このため，骨組解析中で柱梁接合部の応力や変形をも正しく計算するためには工夫が必要で，例えば非特許文献１や非特許文献２では，平面骨組の柱部材や梁部材の部材剛性マトリクス中において，外力項が節点モーメントに対応する列の他に，外力項がパネルモーメントに対応する列を新たに加える変換マトリクスを使用している。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】 日本建築学会，「鋼構造接合部設計指針」の５章：柱梁接合部パネル，丸善，２００６年３月発行，ｐ．２０６〜２３２
【非特許文献２】 小川厚治・多田元英，「柱・梁接合部パネルの変形を考慮した静的・動的応答解析プログラムの開発」，第１７回情報・システム・利用技術シンポジウム論文集，日本建築学会，１９９４年１２月発行，ｐ．７９〜８４
【非特許文献３】 百野泰樹ほか，「方杖ダンパーを用いた柱梁高力ボルト接合骨組の地震応答解析」，日本建築学会構造系論文集 第５８６号，日本建築学会，２００４年１２月発行，ｐ．２１９〜２２６
【発明の概要】
【発明が解決しようとする課題】
【０００４】
以上に述べた従来の骨組解析方法では，以下の点で課題があった。（１）柱梁接合部と柱部材または梁部材が特定の位置関係にある場合の変換マトリクスのみを扱っているので，任意形状の骨組に対する適応性に欠けてしまう。（２）パネルモーメントを与えるために柱部材と梁部材で全く異なる変換マトリクスを用いているので，本来は統一的に扱えるはずの柱部材と梁部材を分化させてしまい，任意形状の骨組に対する汎用性を低下させてしまう。（３）外力項に加えたパネルモーメントを部材応力のまま，骨組全体の剛性方程式中に紛れ込ませていて，本来は既知であるべき全体外力ベクトルの中に未知数が混入していて，部材応力を仮定しないと骨組全体の連立１次方程式をそのまま解けなくしてしまう。（４）平面骨組を対象としていて，立体骨組にはそのまま適用できない，という課題があった。
【課題を解決するための手段】
【０００５】
これら課題を解決するため，以下の（ア）〜（ウ）の手段を併用する。（ア）節点モーメントに対応する外力項をパネルモーメントに対応する外力項で置き換える。（イ）パネルモーメントを応力でなく，自由度間の一対の外力として扱う事で，骨組全体の方程式中の対応外力項を既知数に変化させる。（ウ）柱梁接合部と柱部材または梁部材の位置関係が，特定の場合に制約されず，かつ立体骨組にも適用可能な汎用性のある変換マトリクスを見出す。以上の併用する各手段（ア）〜（ウ）について，以下，順番に記す。
【０００６】
まず，前記手段（ア）の成立性を調べるため，節点モーメントの釣合条件とパネルモーメントの釣合条件の関係性を検討する。例えば平面骨組の柱梁接合部のパネルモーメントＭＰは，通常の形状の場合，左右の梁の略フェイスのモーメントをＭＬおよびＭＲ，上下の柱のせん断力をＱＵおよびＱＤ，梁のフランジ中心間距離をＷＢとして，符号の向きを図２の右手系方向にとると，次式のように表せる。
【数１】
【０００７】
なおＷＢは，鉄筋コンクリート構造では応力中心間距離などにできる。一方，上下の柱の略フェイスのモーメントをＭＵおよびＭＤ，左右の梁のせん断力をＱＬおよびＱＲ，柱のフランジ中心間距離をＷＣとして表現するなら，パネルモーメントは，次式のようにも表せる。
【数２】
【０００８】
上記の数式１と数式２を連立させると，パネルモーメントの釣合式が次式のように得られる。
【数３】
【０００９】
ところが数式３の右辺において，各括弧でまとめられた第１項〜第４項は，それぞれ，左右の梁および上下の柱の節点モーメントと同じになっている。つまり数式３は節点モーメントの釣合式そのものに等しい。したがって，節点モーメントの釣合式がパネルモーメントの釣合式で代替可能になって，手段（ア）に従い，骨組全体の方程式中の略任意の節点モーメントに対応する外力項を，パネルモーメントに対応する外力項で置き換えても骨組全体の方程式は成り立つ事になる。この手段は，物理的には，外力モーメントの作用領域を，節点という点に代えて，接合部やパネル等の大きさを持つ領域としてとらえ直す意味を持っている。
【００１０】
次に，前記手段（イ）については，（１）柱と梁が接合される節点の回転自由度について，柱側と梁側で異なる回転自由度を設定して，（２）両回転自由度間に，柱梁接合部の剛性や耐力を模擬する回転ばねまたは回転ばねを含む要素を設けて，更に（３）柱側の該回転自由度に対応する外力項には，単純な平面解析の場合には数式２に類する，柱側の回転自由度に係わるパネルモーメントが加わるように，柱および梁の部材剛性マトリクスを工夫して，同様に（４）梁側の該回転自由度に対応する外力項には，単純な平面解析の場合には数式１に類する，梁側の回転自由度に係わるパネルモーメントが加わるように，柱および梁の部材剛性マトリクスを工夫する。これらにより，パネルモーメントが自由度間の一対の外力として扱えて，骨組全体の方程式中で対応する外力項を未知数ではなく既知数にできる。
【００１１】
次に，前記手段（ウ）については，もし平面骨組の場合であれば図１で描けるような，柱梁接合部のせん断変形を追加する変換マトリクスを考える。まず，３次元中の節点変位ベクトルＤを，３方向の並進変位ＤＸ・ＤＹ・ＤＺ，自部材側（柱側または梁側で順不同）の３方向の回転変位ＲＣＸ・ＲＣＹ・ＲＣＺ，相手部材側（梁側または柱側で順不同）の３方向の回転変位ＲＢＸ・ＲＢＹ・ＲＢＺ，を使って次式でおく。
【数４】
【００１２】
また，該接合部のせん断変形を除外して考えた，自部材用の節点変位ベクトルＤＣは，図１の二重線の部分がせん断変形しなかった場合を考えれば良いから，３方向の自部材用の並進変位ＤＣＸ・ＤＣＹ・ＤＣＺ，を使って次式でおける。
【数５】
【００１３】
また，該せん断変形を除外して考えた，節点から自部材の略フェイス端への距離ベクトルＥＣを，３方向の成分ＥＣＸ・ＥＣＹ・ＥＣＺ，を使って次式でおく。
【数６】
【００１４】
このとき，ＤＣとＤの並進変位の差は，図１の二重線の部分のせん断変形を考えれば良いから，次式で表せる。
【数７】
【００１５】
上式を，せん断変形を付加するための変換マトリクスＴＣを用いて書き直すと，以下の式が得られる。
【数８】
【数９】
【００１６】
以上により，柱梁接合部のせん断変形を除外して考えた部材剛性マトリクスをＫＣとすると，柱梁接合部のせん断変形を考慮した部材剛性マトリクスＫは，次式で与えられる。
【数１０】
【００１７】
数式９の変換マトリクスＴＣは，距離ＥＣの３成分ＥＣＸ・ＥＣＹ・ＥＣＺを使って表されているので，柱梁接合部と柱部材または梁部材の位置関係には特に制約がなく，汎用性のある変換マトリクスが得られている。
【００１８】
なお立体骨組中の柱部材や梁部材において，材端の剛域のための変換マトリクスＴＥは，距離ＥＣに対して次式で書ける事が知られている。
【数１１】
【００１９】
このため，せん断変形を付加するための変換マトリクスを，数式１１の距離ＥＣに関する変換とまとめて，一つの変換マトリクスＴＥＴＣとすれば，次式で表せる。
【数１２】
【００２０】
上式のＴＥＴＣは，数式９のＴＣよりも簡単な形になって，実施時に都合が良い。なお以上の変換マトリクスは，例えば通常の骨組解析で使用する全体座標系−局所座標系間の変換マトリクスと積にしても表現可能であるから，マトリクスの具体的内容そのものは座標系の設定に伴って変化する性質がある。もちろん以上の変換マトリクスは，単純に２次元用に低減する事で，平面骨組解析にもそのまま適用できる。
【００２１】
得られた変換マトリクスとパネルモーメントの関係を調べておく。図２のような平面骨組の接合部を考えた場合，略フェイス端での部材外力の並進成分ＦＣＸ・ＦＣＹおよび回転成分ＭＣＺを，節点での部材外力の並進成分ＦＸ・ＦＹ，自部材側の回転成分ＭＰＣＺおよび相手部材側の回転成分ＭＰＢＺに変換するマトリクスは，数式１２のＴＥＴＣを２次元に低減して転置すれば良いから，次式で表せる。
【数１３】
【００２２】
上式の第３行と第４行について，図２のＭＬ・ＭＲ・ＭＵ・ＭＤ，ＱＬ・ＱＲ・ＱＵ・ＱＤ，ＷＢ・ＷＣを用いて，外力の向きを考慮して左右の梁および上下の柱の和をとると，柱側の回転角に関するモーメントＭＰＣ，および梁側の回転角に関するモーメントＭＰＢは以下の式のようになる。
【数１４】
【数１５】
【数１６】
【００２３】
つまり，ＭＰＣやＭＰＢが変換マトリクスによって節点モーメントではなく数式１や数式２で記載したパネルモーメントＭＰを与えている事がわかる。
【発明の効果】
【００２４】
本発明によれば，柱梁接合部の応力や変形をも含めて計算する骨組解析において，（１）柱梁接合部と柱部材または梁部材が特定の位置関係になくても計算ができ，（２）任意形状の骨組を計算でき，（３）立体骨組でも計算でき，（４）柱部材と梁部材で統一的な変換マトリクスを使うのでプログラムが簡明になり，（５）パネルモーメントを未知数のまま骨組全体の剛性方程式中に紛れ込ませる事がないので解法上の支障がなく，（６）簡単な形の変換マトリクスなので計算コストを増やさずに，（７）節点モーメントをパネルモーメントに替える事で節点を点から大きさを持つ領域へと拡張可能な物にして，（８）立体３方向のパネルモーメントが得られるため３次元的に剛性や耐力が相関した柱梁接合部の組み込みが容易になる，という効果がある。
【発明を実施するための形態】
【００２５】
以下，本発明の実施の最良の形態を図３に基づいて説明する。図３は最良の形態のフローチャートで，プログラム全体の流れを示す。以下，流れに沿って説明する。最初に入力データの取得を行う（Ｓ１）。該取得は，入力装置から行っても良く，既存のファイルから行っても良い。該入力データは，例えば，各節点の位置や自由度，各部材の種類や位置や剛性や耐力，与える外力や地動，等の骨組解析で用いる情報で良い。続いて，自由度の設定を行うが，柱側の回転自由度と梁側の回転自由度とは別個の自由度に設定する（Ｓ２）。
【００２６】
次に，柱や梁の部材剛性マトリクスを作成するが，最初に材端の剛域部分などを含まない剛性マトリクスを作成する（Ｓ３）。この作成では，基本的な線材理論や一般化塑性ヒンジ理論や弾塑性大変形理論や１次元有限要素法などを必要に応じて利用して，弾性または弾塑性の剛性マトリクスを作成する。続いて，数式１２の変換マトリクスを用いて，数式１０と同様の座標変換を行う（Ｓ４）。変換された剛性マトリクスは，対応自由度として柱側の回転自由度と梁側の回転自由度の両方を基本的に含むようになる。続いて，該変換後の剛性マトリクスを骨組全体の剛性マトリクスの対応位置へ足し込む（Ｓ５）。以上のステップＳ３〜Ｓ５を柱部材や梁部材について繰り返す。
【００２７】
次に，柱梁接合部について，例えば１〜３次元の回転ばね等の要素でモデル化して，該要素の剛性マトリクスを計算して，全体剛性マトリクスの対応位置へ足し込む（Ｓ６）。その足し込み位置は基本的には，柱側の回転自由度と梁側の回転自由度の間などとすれば良い。該回転ばね等の要素は，立体３方向で剛性や耐力の相関を有するばねとしても良く，軸剛性や軸耐力との相関を有するばねとしても良い。該ステップＳ６は柱梁接合部について繰り返す。次に，必要があれば，その他の部材等の剛性マトリクスを作成して，全体剛性マトリクスへ足し込む（Ｓ７）。該ステップＳ７はその他の部材等について繰り返す。
【００２８】
次に，荷重増分や地動などに応じた骨組全体の外力ベクトルを計算し（Ｓ８），連立１次方程式を解いて，変位増分を計算し，節点変位を更新する（Ｓ９）。該節点変位をもとに，柱部材や梁部材や柱梁接合部やその他部材等について，応力や歪を計算すると共に状態の更新を行う（Ｓ１０）。該ステップＳ１０は，柱部材や梁部材や柱梁接合部やその他の部材等について繰り返す。次に，必要に応じて，全部または特定の節点変位や部材応力や歪などの数値を記録装置のファイル等に出力すると共に，以上のステップＳ３〜Ｓ１０の主要な計算結果を内部メモリー等の記憶装置に保存する（Ｓ１１）。なお数値のファイル等への出力や，主要な計算結果の保存は，ステップＳ９やＳ１０における更新時に行っても良く，ステップＳ３〜Ｓ１０の中で随時に行っても良い。
【００２９】
以上の部分的繰り返しを含む一連のステップＳ３〜Ｓ１１を，増分解析のステップ毎に繰り返して，計算実行する。全ての増分解析ステップを終了するか，または増分解析の終了条件を満足したら，ステップＳ１１等で保存した計算結果を該記憶装置から読み出して，必要に応じた再構成を行って，変位や応力や歪などの計算結果のデータを記録装置のファイル等に出力する（Ｓ１２）。なお出力は別の時点で行っても良い。該出力ファイルは，別に用意されたプログラムで必要な数値や図を画面に表示したり，紙に印刷したりする時に利用される。
【産業上の利用可能性】
【００３０】
本発明は，骨組解析を利用する産業分野において利用可能で，建築分野や土木分野において建物や構築物設計のためのプログラムとして利用可能である。また本発明を利用したプログラムそのものを市販可能である。
【図面の簡単な説明】
【００３１】
【図１】 柱梁接合部のせん断変形の追加
【図２】 平面骨組での右手系の符号の向き
【図３】 最良の形態のフローチャート
【符号の説明】
【００３２】
１ 柱
２ 梁
３ 柱梁接合部
【技術分野】
【０００１】
本発明は，柱や梁などの骨組部材を線材に置換して，立体骨組の応力や変形を計算するプログラムに関し，特に，柱梁接合部の変形をも計算するプログラムに関する。
【背景技術】
【０００２】
骨組解析では，節点を設けて，節点という点への作用力として，並進力やモーメントを考えている。一方，骨組中の柱梁接合部の応力は，例えば鋼構造骨組では非特許文献１に示されているパネルモーメントで評価をされる。パネルモーメントは，節点モーメントとは異なっていて，例えば，梁フェイスモーメント・柱せん断力・接合部寸法から計算されるか，または柱フェイスモーメント・梁せん断力・接合部寸法から計算されるが，いずれも複数部材の応力の関数になっている。このため，骨組解析中で柱梁接合部の応力や変形をも正しく計算するためには工夫が必要で，例えば非特許文献１や非特許文献２では，平面骨組の柱部材や梁部材の部材剛性マトリクス中において，外力項が節点モーメントに対応する列の他に，外力項がパネルモーメントに対応する列を新たに加える変換マトリクスを使用している。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】 日本建築学会，「鋼構造接合部設計指針」の５章：柱梁接合部パネル，丸善，２００６年３月発行，ｐ．２０６〜２３２
【非特許文献２】 小川厚治・多田元英，「柱・梁接合部パネルの変形を考慮した静的・動的応答解析プログラムの開発」，第１７回情報・システム・利用技術シンポジウム論文集，日本建築学会，１９９４年１２月発行，ｐ．７９〜８４
【非特許文献３】 百野泰樹ほか，「方杖ダンパーを用いた柱梁高力ボルト接合骨組の地震応答解析」，日本建築学会構造系論文集 第５８６号，日本建築学会，２００４年１２月発行，ｐ．２１９〜２２６
【発明の概要】
【発明が解決しようとする課題】
【０００４】
以上に述べた従来の骨組解析方法では，以下の点で課題があった。（１）柱梁接合部と柱部材または梁部材が特定の位置関係にある場合の変換マトリクスのみを扱っているので，任意形状の骨組に対する適応性に欠けてしまう。（２）パネルモーメントを与えるために柱部材と梁部材で全く異なる変換マトリクスを用いているので，本来は統一的に扱えるはずの柱部材と梁部材を分化させてしまい，任意形状の骨組に対する汎用性を低下させてしまう。（３）外力項に加えたパネルモーメントを部材応力のまま，骨組全体の剛性方程式中に紛れ込ませていて，本来は既知であるべき全体外力ベクトルの中に未知数が混入していて，部材応力を仮定しないと骨組全体の連立１次方程式をそのまま解けなくしてしまう。（４）平面骨組を対象としていて，立体骨組にはそのまま適用できない，という課題があった。
【課題を解決するための手段】
【０００５】
これら課題を解決するため，以下の（ア）〜（ウ）の手段を併用する。（ア）節点モーメントに対応する外力項をパネルモーメントに対応する外力項で置き換える。（イ）パネルモーメントを応力でなく，自由度間の一対の外力として扱う事で，骨組全体の方程式中の対応外力項を既知数に変化させる。（ウ）柱梁接合部と柱部材または梁部材の位置関係が，特定の場合に制約されず，かつ立体骨組にも適用可能な汎用性のある変換マトリクスを見出す。以上の併用する各手段（ア）〜（ウ）について，以下，順番に記す。
【０００６】
まず，前記手段（ア）の成立性を調べるため，節点モーメントの釣合条件とパネルモーメントの釣合条件の関係性を検討する。例えば平面骨組の柱梁接合部のパネルモーメントＭＰは，通常の形状の場合，左右の梁の略フェイスのモーメントをＭＬおよびＭＲ，上下の柱のせん断力をＱＵおよびＱＤ，梁のフランジ中心間距離をＷＢとして，符号の向きを図２の右手系方向にとると，次式のように表せる。
【数１】
【０００７】
なおＷＢは，鉄筋コンクリート構造では応力中心間距離などにできる。一方，上下の柱の略フェイスのモーメントをＭＵおよびＭＤ，左右の梁のせん断力をＱＬおよびＱＲ，柱のフランジ中心間距離をＷＣとして表現するなら，パネルモーメントは，次式のようにも表せる。
【数２】
【０００８】
上記の数式１と数式２を連立させると，パネルモーメントの釣合式が次式のように得られる。
【数３】
【０００９】
ところが数式３の右辺において，各括弧でまとめられた第１項〜第４項は，それぞれ，左右の梁および上下の柱の節点モーメントと同じになっている。つまり数式３は節点モーメントの釣合式そのものに等しい。したがって，節点モーメントの釣合式がパネルモーメントの釣合式で代替可能になって，手段（ア）に従い，骨組全体の方程式中の略任意の節点モーメントに対応する外力項を，パネルモーメントに対応する外力項で置き換えても骨組全体の方程式は成り立つ事になる。この手段は，物理的には，外力モーメントの作用領域を，節点という点に代えて，接合部やパネル等の大きさを持つ領域としてとらえ直す意味を持っている。
【００１０】
次に，前記手段（イ）については，（１）柱と梁が接合される節点の回転自由度について，柱側と梁側で異なる回転自由度を設定して，（２）両回転自由度間に，柱梁接合部の剛性や耐力を模擬する回転ばねまたは回転ばねを含む要素を設けて，更に（３）柱側の該回転自由度に対応する外力項には，単純な平面解析の場合には数式２に類する，柱側の回転自由度に係わるパネルモーメントが加わるように，柱および梁の部材剛性マトリクスを工夫して，同様に（４）梁側の該回転自由度に対応する外力項には，単純な平面解析の場合には数式１に類する，梁側の回転自由度に係わるパネルモーメントが加わるように，柱および梁の部材剛性マトリクスを工夫する。これらにより，パネルモーメントが自由度間の一対の外力として扱えて，骨組全体の方程式中で対応する外力項を未知数ではなく既知数にできる。
【００１１】
次に，前記手段（ウ）については，もし平面骨組の場合であれば図１で描けるような，柱梁接合部のせん断変形を追加する変換マトリクスを考える。まず，３次元中の節点変位ベクトルＤを，３方向の並進変位ＤＸ・ＤＹ・ＤＺ，自部材側（柱側または梁側で順不同）の３方向の回転変位ＲＣＸ・ＲＣＹ・ＲＣＺ，相手部材側（梁側または柱側で順不同）の３方向の回転変位ＲＢＸ・ＲＢＹ・ＲＢＺ，を使って次式でおく。
【数４】
【００１２】
また，該接合部のせん断変形を除外して考えた，自部材用の節点変位ベクトルＤＣは，図１の二重線の部分がせん断変形しなかった場合を考えれば良いから，３方向の自部材用の並進変位ＤＣＸ・ＤＣＹ・ＤＣＺ，を使って次式でおける。
【数５】
【００１３】
また，該せん断変形を除外して考えた，節点から自部材の略フェイス端への距離ベクトルＥＣを，３方向の成分ＥＣＸ・ＥＣＹ・ＥＣＺ，を使って次式でおく。
【数６】
【００１４】
このとき，ＤＣとＤの並進変位の差は，図１の二重線の部分のせん断変形を考えれば良いから，次式で表せる。
【数７】
【００１５】
上式を，せん断変形を付加するための変換マトリクスＴＣを用いて書き直すと，以下の式が得られる。
【数８】
【数９】
【００１６】
以上により，柱梁接合部のせん断変形を除外して考えた部材剛性マトリクスをＫＣとすると，柱梁接合部のせん断変形を考慮した部材剛性マトリクスＫは，次式で与えられる。
【数１０】
【００１７】
数式９の変換マトリクスＴＣは，距離ＥＣの３成分ＥＣＸ・ＥＣＹ・ＥＣＺを使って表されているので，柱梁接合部と柱部材または梁部材の位置関係には特に制約がなく，汎用性のある変換マトリクスが得られている。
【００１８】
なお立体骨組中の柱部材や梁部材において，材端の剛域のための変換マトリクスＴＥは，距離ＥＣに対して次式で書ける事が知られている。
【数１１】
【００１９】
このため，せん断変形を付加するための変換マトリクスを，数式１１の距離ＥＣに関する変換とまとめて，一つの変換マトリクスＴＥＴＣとすれば，次式で表せる。
【数１２】
【００２０】
上式のＴＥＴＣは，数式９のＴＣよりも簡単な形になって，実施時に都合が良い。なお以上の変換マトリクスは，例えば通常の骨組解析で使用する全体座標系−局所座標系間の変換マトリクスと積にしても表現可能であるから，マトリクスの具体的内容そのものは座標系の設定に伴って変化する性質がある。もちろん以上の変換マトリクスは，単純に２次元用に低減する事で，平面骨組解析にもそのまま適用できる。
【００２１】
得られた変換マトリクスとパネルモーメントの関係を調べておく。図２のような平面骨組の接合部を考えた場合，略フェイス端での部材外力の並進成分ＦＣＸ・ＦＣＹおよび回転成分ＭＣＺを，節点での部材外力の並進成分ＦＸ・ＦＹ，自部材側の回転成分ＭＰＣＺおよび相手部材側の回転成分ＭＰＢＺに変換するマトリクスは，数式１２のＴＥＴＣを２次元に低減して転置すれば良いから，次式で表せる。
【数１３】
【００２２】
上式の第３行と第４行について，図２のＭＬ・ＭＲ・ＭＵ・ＭＤ，ＱＬ・ＱＲ・ＱＵ・ＱＤ，ＷＢ・ＷＣを用いて，外力の向きを考慮して左右の梁および上下の柱の和をとると，柱側の回転角に関するモーメントＭＰＣ，および梁側の回転角に関するモーメントＭＰＢは以下の式のようになる。
【数１４】
【数１５】
【数１６】
【００２３】
つまり，ＭＰＣやＭＰＢが変換マトリクスによって節点モーメントではなく数式１や数式２で記載したパネルモーメントＭＰを与えている事がわかる。
【発明の効果】
【００２４】
本発明によれば，柱梁接合部の応力や変形をも含めて計算する骨組解析において，（１）柱梁接合部と柱部材または梁部材が特定の位置関係になくても計算ができ，（２）任意形状の骨組を計算でき，（３）立体骨組でも計算でき，（４）柱部材と梁部材で統一的な変換マトリクスを使うのでプログラムが簡明になり，（５）パネルモーメントを未知数のまま骨組全体の剛性方程式中に紛れ込ませる事がないので解法上の支障がなく，（６）簡単な形の変換マトリクスなので計算コストを増やさずに，（７）節点モーメントをパネルモーメントに替える事で節点を点から大きさを持つ領域へと拡張可能な物にして，（８）立体３方向のパネルモーメントが得られるため３次元的に剛性や耐力が相関した柱梁接合部の組み込みが容易になる，という効果がある。
【発明を実施するための形態】
【００２５】
以下，本発明の実施の最良の形態を図３に基づいて説明する。図３は最良の形態のフローチャートで，プログラム全体の流れを示す。以下，流れに沿って説明する。最初に入力データの取得を行う（Ｓ１）。該取得は，入力装置から行っても良く，既存のファイルから行っても良い。該入力データは，例えば，各節点の位置や自由度，各部材の種類や位置や剛性や耐力，与える外力や地動，等の骨組解析で用いる情報で良い。続いて，自由度の設定を行うが，柱側の回転自由度と梁側の回転自由度とは別個の自由度に設定する（Ｓ２）。
【００２６】
次に，柱や梁の部材剛性マトリクスを作成するが，最初に材端の剛域部分などを含まない剛性マトリクスを作成する（Ｓ３）。この作成では，基本的な線材理論や一般化塑性ヒンジ理論や弾塑性大変形理論や１次元有限要素法などを必要に応じて利用して，弾性または弾塑性の剛性マトリクスを作成する。続いて，数式１２の変換マトリクスを用いて，数式１０と同様の座標変換を行う（Ｓ４）。変換された剛性マトリクスは，対応自由度として柱側の回転自由度と梁側の回転自由度の両方を基本的に含むようになる。続いて，該変換後の剛性マトリクスを骨組全体の剛性マトリクスの対応位置へ足し込む（Ｓ５）。以上のステップＳ３〜Ｓ５を柱部材や梁部材について繰り返す。
【００２７】
次に，柱梁接合部について，例えば１〜３次元の回転ばね等の要素でモデル化して，該要素の剛性マトリクスを計算して，全体剛性マトリクスの対応位置へ足し込む（Ｓ６）。その足し込み位置は基本的には，柱側の回転自由度と梁側の回転自由度の間などとすれば良い。該回転ばね等の要素は，立体３方向で剛性や耐力の相関を有するばねとしても良く，軸剛性や軸耐力との相関を有するばねとしても良い。該ステップＳ６は柱梁接合部について繰り返す。次に，必要があれば，その他の部材等の剛性マトリクスを作成して，全体剛性マトリクスへ足し込む（Ｓ７）。該ステップＳ７はその他の部材等について繰り返す。
【００２８】
次に，荷重増分や地動などに応じた骨組全体の外力ベクトルを計算し（Ｓ８），連立１次方程式を解いて，変位増分を計算し，節点変位を更新する（Ｓ９）。該節点変位をもとに，柱部材や梁部材や柱梁接合部やその他部材等について，応力や歪を計算すると共に状態の更新を行う（Ｓ１０）。該ステップＳ１０は，柱部材や梁部材や柱梁接合部やその他の部材等について繰り返す。次に，必要に応じて，全部または特定の節点変位や部材応力や歪などの数値を記録装置のファイル等に出力すると共に，以上のステップＳ３〜Ｓ１０の主要な計算結果を内部メモリー等の記憶装置に保存する（Ｓ１１）。なお数値のファイル等への出力や，主要な計算結果の保存は，ステップＳ９やＳ１０における更新時に行っても良く，ステップＳ３〜Ｓ１０の中で随時に行っても良い。
【００２９】
以上の部分的繰り返しを含む一連のステップＳ３〜Ｓ１１を，増分解析のステップ毎に繰り返して，計算実行する。全ての増分解析ステップを終了するか，または増分解析の終了条件を満足したら，ステップＳ１１等で保存した計算結果を該記憶装置から読み出して，必要に応じた再構成を行って，変位や応力や歪などの計算結果のデータを記録装置のファイル等に出力する（Ｓ１２）。なお出力は別の時点で行っても良い。該出力ファイルは，別に用意されたプログラムで必要な数値や図を画面に表示したり，紙に印刷したりする時に利用される。
【産業上の利用可能性】
【００３０】
本発明は，骨組解析を利用する産業分野において利用可能で，建築分野や土木分野において建物や構築物設計のためのプログラムとして利用可能である。また本発明を利用したプログラムそのものを市販可能である。
【図面の簡単な説明】
【００３１】
【図１】 柱梁接合部のせん断変形の追加
【図２】 平面骨組での右手系の符号の向き
【図３】 最良の形態のフローチャート
【符号の説明】
【００３２】
１ 柱
２ 梁
３ 柱梁接合部
【特許請求の範囲】
【請求項１】
柱や梁などの骨組部材を線材に置換して，立体骨組の応力や変形を計算するプログラムであって，以下の（ａ）〜（ｅ）の手順すなわち，（ａ）柱と梁が接合される節点の回転自由度について，柱側と梁側で異なる回転自由度を設定可能にした手順と，（ｂ）柱側の該回転自由度に対応する外力項を節点モーメントではなくパネルモーメントにして，該パネルモーメントは柱側の回転自由度に係わるパネルモーメントとなるように，柱および梁の部材剛性を座標変換する手順と，（ｃ）梁側の該回転自由度に対応する外力項を節点モーメントではなくパネルモーメントにして，該パネルモーメントは梁側の回転自由度に係わるパネルモーメントとなるように，柱および梁の部材剛性を座標変換する手順と，（ｄ）該柱側と該梁側の該回転自由度間に，柱梁接合部の剛性や耐力を模擬する回転ばねまたは回転ばねを含む要素を設定可能にした手順と，（ｅ）柱梁接合部を含めて骨組の応力や変形を計算する手順と，をコンピュータに実行させるためのプログラム。
【請求項２】
柱や梁などの骨組部材を線材に置換して，平面骨組の応力や変形を計算するプログラムであって，前記手順（ａ）〜（ｅ）をコンピュータに実行させるためのプログラム。
【請求項１】
柱や梁などの骨組部材を線材に置換して，立体骨組の応力や変形を計算するプログラムであって，以下の（ａ）〜（ｅ）の手順すなわち，（ａ）柱と梁が接合される節点の回転自由度について，柱側と梁側で異なる回転自由度を設定可能にした手順と，（ｂ）柱側の該回転自由度に対応する外力項を節点モーメントではなくパネルモーメントにして，該パネルモーメントは柱側の回転自由度に係わるパネルモーメントとなるように，柱および梁の部材剛性を座標変換する手順と，（ｃ）梁側の該回転自由度に対応する外力項を節点モーメントではなくパネルモーメントにして，該パネルモーメントは梁側の回転自由度に係わるパネルモーメントとなるように，柱および梁の部材剛性を座標変換する手順と，（ｄ）該柱側と該梁側の該回転自由度間に，柱梁接合部の剛性や耐力を模擬する回転ばねまたは回転ばねを含む要素を設定可能にした手順と，（ｅ）柱梁接合部を含めて骨組の応力や変形を計算する手順と，をコンピュータに実行させるためのプログラム。
【請求項２】
柱や梁などの骨組部材を線材に置換して，平面骨組の応力や変形を計算するプログラムであって，前記手順（ａ）〜（ｅ）をコンピュータに実行させるためのプログラム。
【図１】
【図２】
【図３】
【図２】
【図３】
【公開番号】特開２０１２−１１３６７９（Ｐ２０１２−１１３６７９Ａ）
【公開日】平成２４年６月１４日（２０１２．６．１４）
【国際特許分類】
【出願番号】特願２０１０−２７５０８４（Ｐ２０１０−２７５０８４）
【出願日】平成２２年１１月２２日（２０１０．１１．２２）
【特許番号】特許第４７６６３５５号（Ｐ４７６６３５５）
【特許公報発行日】平成２３年９月７日（２０１１．９．７）
【出願人】（５０７２７４３６４）株式会社テクニカルスルー (2)
【Ｆターム（参考）】
【公開日】平成２４年６月１４日（２０１２．６．１４）
【国際特許分類】
【出願日】平成２２年１１月２２日（２０１０．１１．２２）
【特許番号】特許第４７６６３５５号（Ｐ４７６６３５５）
【特許公報発行日】平成２３年９月７日（２０１１．９．７）
【出願人】（５０７２７４３６４）株式会社テクニカルスルー (2)
【Ｆターム（参考）】
[ Back to top ]