JP7374689B2

JP7374689B2 - 設計システム及び設計方法

Info

Publication number: JP7374689B2
Application number: JP2019177127A
Authority: JP
Inventors: 竜介池田; 康男新田; 悠歩河本; 敬章津嘉田; 義輝田中; 孝誠谷内; 聡子後藤; 将武松本
Original assignee: Shimizu Corp
Current assignee: Shimizu Corp
Priority date: 2019-09-27
Filing date: 2019-09-27
Publication date: 2023-11-07
Anticipated expiration: 2039-09-27
Also published as: JP2021056619A

Description

本発明は、構造物の設計のための応力計算を行う設計システム及び設計方法に関する。

海外の原子力発電所等の構造物の設計では３次元の有限要素法（３Ｄ－ＦＥＭ：Finite Element Method）を用いた地震応答解析（動解）から直接得られる要素応力や複数要素で構成される部材応力を用いた設計が行われている（例えば、特許文献１参照）。

３Ｄ－ＦＥＭを用いた解析により応力時刻歴データを計算すると、出力されるデータ量が膨大なので全時刻で構造物の断面算定をすると、計算時間が膨大となる。そこで、設計上クリティカル（重大）となる応力データを抽出し、抽出した応力データに基づいて設計が行われている。

特開２０１１－１０７０４０号公報

しかしながら、３Ｄ－ＦＥＭを用いた解析結果に基づいて、時刻に関係なく応力の最大値を用いた設計を行うと、例えば軸力の最大値と曲げモーメントの最大値が同時に発生する保守的な設計となってしまうという課題がある。

本発明は、上述する問題点に鑑みてなされたものであり、構造物に作用する応力の膨大なデータから簡易に必要なデータを抽出しつつも合理的な設計を行える設計システム及び設計方法を提供することを目的とする。

上記の目的を達するために、本発明は、３次元の有限要素法を用いて設計対象の構造物の設計を支援する設計システムであって、前記構造物の構造及び荷重に関する設計データに基づいて、地震入力に対して前記構造物に作用する応力の成分を３次元解析し、前記構造物に作用する軸力と曲げモーメントとの関係を示す軌跡を算出し、前記軌跡が含まれる領域の各象限毎に最大値を抽出して、前記構造物の断面に作用する第１設計用応力を算出する演算部を備えることを特徴とする、設計システムである。

本発明によれば、地震入力により構造物の断面に作用する動的な応力の３次元の有限要素法を用いた解析において、各象限毎に解析結果の最大値を設定することにより、設計上クリティカルとなる応力データを合理的に抽出すると共に、設計に用いるデータの数を大幅に削減することができる。

また、本発明は、前記演算部は、前記設計データに基づいて、固定荷重を含む静的荷重に対して前記構造物に作用する応力の成分を３次元解析し、前記構造物の断面に作用する第２設計用応力を算出し、前記第１設計用応力と前記第２設計用応力とに基づいて、前記構造物の断面に作用する応力を算出するように構成されていてもよい。

本発明によれば、固定荷重を含む静的荷重により構造物の断面に作用する応力の３次元の有限要素法を用いた解析をすることにより、第１設計用応力と第２設計用応力とを組み合わせた応力に基づいて構造物の断面算定することができる。

３次元の有限要素法を用いて設計対象の構造物の設計を支援する設計方法であって、前記構造物の構造及び荷重に関する設計データに基づいて、地震入力に対して前記構造物に作用する応力の成分を３次元解析し、前記構造物に作用する軸力と曲げモーメントとの関係を示す軌跡を算出し、各象限の最大値を設計用応力空間を設定し、前記構造物の断面に作用する第１設計用応力を算出することを特徴とする、設計方法である。

本発明によれば、地震入力により構造物の断面に作用する動的な応力の３次元の有限要素法を用いた解析において、解析結果の各象限毎の最大値を設定することにより、設計上クリティカルとなる応力データを合理的に抽出すると共に、設計に用いるデータの数を大幅に削減することができる。

本発明によれば、構造物に作用する応力の膨大なデータから簡易に必要なデータを抽出しつつも合理的な設計を行うことができる。

本発明の実施形態に係る設計システムの構成を示すブロック図である。３次元の有限要素法によりモデル化された建物を示す斜視図である。有限要素法の要素に加わる応力を示す図である。解析結果の軌跡を包絡する方法を示す図である。設計システムにおいて実行される処理の流れを示すフローチャートである。

以下、図面を参照しつつ、本発明に係る設計システム１の実施形態について説明する。設計システム１は、地震力により建物の断面に作用する応力を３次元の有限要素法（３Ｄ－ＦＥＭ）を用いて解析する設計支援装置である。

図１に示されるように、設計システム１は、設計データが入力される入力部２と、入力されたデータに基づいて設計値を算出する演算部４と、演算部４の算出結果を表示する表示部６と、演算部４の演算に必要なデータを記憶する記憶部８と、を備える。

設計システム１は、例えば、パーソナルコンピュータ、タブレット型端末、スマートフォン等の端末装置により実現される。設計システム１は、ネットワークを通じて演算結果を出力するサーバ装置であってもよい。

入力部２は、キーボード、タッチパネル等により実現されるデータ入力のためのユーザインタフェースである。入力部２は、タブレット型端末やスマートフォンにより無線又は有線等により接続される別体の端末装置であってもよい。入力部２からは、設計対象の建物の構造及び荷重等の設計に関する設計データが入力される。入力された設計データは、記憶部８に記憶される。設計データは、例えば、設計対象物の寸法、間取り、部材の重量、材料、地震波の波形、風荷重等の固定荷重等の各種データが含まれる。

記憶部８は、フラッシュメモリやＨＤＤ（Hard Disk Drive）等の記憶媒体により構成された記憶装置である。記憶部８は、入力部２により入力された設計データの他、３Ｄ－ＦＥＭの解析に必要な数式を実行するプログラム等のデータを記憶する。記憶部８は、設計システム１に内蔵されている。記憶部８は、設計システム１に着脱自在な記憶装置であってもよいし、ネットワークを通じて接続されるサーバ装置に内蔵されていてもよい。

演算部４は、メモリや記憶部８に記憶されたデータに基づいて、建物の設計に必要な３Ｄ－ＦＥＭ等の演算を実行する。演算部４は、ＣＰＵ（Central Processing Unit）、ＧＰＵ（Graphics Processing Unit）等のプロセッサがプログラム（ソフトウェア）を実行することで実現される。これらの各機能部のうち一部または全部は、ＬＳＩ（Large Scale Integration）やＡＳＩＣ（Application Specific Integrated Circuit）、ＦＰＧＡ（Field-Programmable Gate Array）等のハードウェアによって実現されてもよいし、ソフトウェアとハードウェアの協働によって実現されてもよい。

プログラムは、予めＨＤＤ（Hard Disk Drive）やフラッシュメモリなどの記憶装置に格納されていてもよいし、ＤＶＤやＣＤ－ＲＯＭなどの着脱可能な記憶媒体に格納されており、記憶媒体がドライブ装置に装着されることで記憶装置にインストールされてもよい。プログラムは、ネットワークを通じて通じた外部サーバから実行されるものであってもよい。

表示部６は、例えば、ＬＣＤ（Liquid Crystal Display）、有機ＥＬ（Electro Luminescence）ディスプレイ、ＬＥＤ（Light Emitting Diode）ディスプレイ等の表示装置である。表示部６は、必ずしも設計システム１に設けられていなくてもよく、設計システム１と無線又は有線で接続されるパーソナルコンピュータ、タブレット型端末、スマートフォン等の他の端末装置により実現されてもよい。

次に、演算部４の具体的な処理の内容について説明する。ユーザは、入力部２を介して設計対象物である建物等の構造物の設計データを入力する。設計データは、記憶部８に記憶される。

図２に示されるように、演算部４は、記憶部８から設計データを読み出し、建物の３次元モデルを生成する。演算部４は、例えば、設計データに基づいて原子炉建屋等の構造物の３次元モデルを生成する。

演算部４は、設計データに基づいて有限要素法（Finite Element Method：ＦＥＭ）モデルを用いて、構造物を無数の要素に分割し、各要素に作用する応力成分を算出する。演算部４は、弾塑性地震応答解析を行う。演算部４は、例えば、地震入力により構造物に作用する動的なｎ個（ｎは自然数）の応力成分を算出する。演算部４は、地震時以外の他に、固定荷重を含む静的荷重により構造物に作用する応力成分を算出する。静的荷重は、例えば、Ｄ：固定荷重、Ｌ：積載荷重、Ｔ：温度荷重、Ｓ：積雪荷重、Ｗ：風圧力、Ｈ：土圧および水圧等の荷重が含まれる。

演算部４は、地震時及び固定荷重などにより作用する応力を３次元ＦＥＭ応答解析モデルで組み合わせて算出する。演算部４は、算出した組合せ応力に基づいて設計用応力を算出する。演算部４は、算出した設計用応力を用いて構造物の断面算定を実施する。

図３に示されるように、演算部４は、各要素に作用する応力を成分毎に時刻歴に基づいて算出する。各要素は、構造物を構成する部位の構造に応じてシェル要素と梁要素とに分けられる。シェル要素は、板やシェルの様な形状の連続体からなる薄板形状の部材のモデル化に用いられる要素である。シェル要素は、見かけ上において厚みがゼロの面で構成され、計算上は板厚分の剛性を持つ。梁要素は、梁などの様な形状の連続体からなる棒状のような形状の部材のモデル化に使用される要素である。梁要素は、見かけ上において線だけの要素で、構成され、計算上は指定した断面の剛性を持つ。

演算部４は、部材をシェル要素に基づいて解析する。演算部４は、シェル要素に作用する８成分の応力時刻歴データを算出する。演算部４は、例えば、部材を梁要素に基づいて解析する。演算部４は、梁要素に作用する６成分の応力時刻歴データを算出する。

シェル要素の断面設計では、例えば、膜力と曲げの応力６成分（Ｎｘ，Ｎｙ，Ｎｘｙ，Ｍｘ，Ｍｙ，Ｍｘｙ）のつり合いを計算する断面設計が行われる。以下、地震応答解析について説明する。

演算部４は、設計データに基づいて地震入力に対して構造物に作用する動的な応力を３次元解析（地震応答解析）する。演算部４は、地震が入力される所定時間内において作用する応力の全データを応力時刻歴データとして出力する。

演算部４は、設計用データに基づいて、固定荷重などに対して構造物に作用する静的な応力を３次元解析する。演算部４は、地震入力により構造物に作用する応力と、固定荷重などにより構造物に作用する応力とを組み合わせた組合せ応力の応力時刻歴データを算出する。

全ての応力時刻歴データは、応力成分を算出するための壁、床などに作用する軸力、せん断力、曲げモーメント等のｎ個の応力の関係を示すｎ次元空間に描かれる計算結果が含まれており、膨大なデータとなる。そのため、演算部４は、膨大な応力時刻歴データから断面設計に用いる設計用応力を抽出する。演算部４は、膨大な応力時刻歴データの中から断面設計上、クリティカル（重大）となる応力データを抽出する。

演算部４は、例えば、軸力と曲げモーメントとの関係を示す計算結果の軌跡が含まれる領域を座標系の各象限毎の最大値又は最小値を設計値として算出する。演算部４は、膜力と曲げの応力６成分について設計用応力空間を設定する。以下２次元の領域を例に説明する。

図４に示されるように、演算部４は、所定時間における軸力と曲げモーメントとの関係を示す解析結果の全ての応力時刻歴データＤを示す軌跡の中から、座標系における第１象限から第４象限毎に対応する符号に応じて応力の最大値又は最小値を抽出する。演算部４は、抽出した応力の最大値又は最小値を設定する。

図５に示されるように、演算部４は、座標系において全ての符号組合せ毎に最大応力を抽出する。

演算部４は、例えば、座標系において応力時刻歴データＤの第１象限Ｑ１における最大値又は最小値を抽出する。第１象限Ｑ１の符号は、（＋，＋）であるので演算部４は、ｘ軸の最大値Ｐ１を通りｙ軸に平行な直線Ｌ１と、ｙ軸の最大値Ｐ２を通りｘ軸に平行な直線Ｌ２との交点となる点Ｐ３とにより包絡された領域Ｒ１を設定する。

演算部４は、例えば、座標系において応力時刻歴データＤの第２象限Ｑ２における最大値又は最小値を抽出する。第２象限Ｑ２の符号は、（－，＋）であるので演算部４は、ｘ軸の最小値Ｐ４を通りｙ軸に平行な直線Ｌ３と、ｙ軸の最大値Ｐ５を通りｘ軸に平行な直線Ｌ４との交点となる点Ｐ６とにより包絡された領域Ｒ２を設定する。

演算部４は、例えば、座標系において応力時刻歴データＤの第３象限Ｑ３における最大値又は最小値を抽出する。第３象限Ｑ３の符号は、（－，－）であるので演算部４は、ｘ軸の最小値Ｐ７を通りｙ軸に平行な直線Ｌ５と、ｙ軸の最小値Ｐ８を通りｘ軸に平行な直線Ｌ６との交点となる点Ｐ９とにより包絡された領域Ｒ３を設定する。

演算部４は、例えば、座標系において応力時刻歴データＤの第４象限Ｑ４における最大値又は最小値を抽出する。第４象限Ｑ４の符号は、（＋，－）であるので演算部４は、ｘ軸の最大値Ｐ１０を通りｙ軸に平行な直線Ｌ７と、ｙ軸の最小値Ｐ１１を通りｘ軸に平行な直線Ｌ８との交点となる点Ｐ１２とにより包絡された領域Ｒ４を設定する。

演算部４は、設定した設計用応力空間である領域Ｒ１からＲ４に基づいて建物の断面算定を行う。即ち、演算部４は、Ｐ３，Ｐ６，Ｐ９，Ｐ１２における応力値を用いて建物の断面算定を行う。

図６は、設計システム１において実行される設計方法の処理の流れを示すフローチャートである。演算部４は、入力部２に入力された設計用データに基づいて、３Ｄ－ＦＥＭを用いて建物の３Ｄモデルを構築し、地震入力により３Ｄモデルの各要素に作用する応力をそれぞれ算出し、建物に作用する動的な応力を解析する（ステップＳ１０）。演算部４は、設計用データに基づいて、３Ｄ－ＦＥＭを用いて固定荷重を含む静的荷重により建物に作用する応力を算出する（ステップＳ１２）。

演算部４は、地震時及び固定荷重を含む静的荷重により建物に作用する応力を組み合わせた応力時刻歴データを算出する（ステップＳ１４）。演算部４は、応力時刻歴データを示す軌跡の中から、全データを象限毎の座標の絶対値の最大値を設計用応力として抽出する（ステップＳ１６）。演算部４は、算出した設計用応力用いて建物の断面算定を行う（ステップＳ１８）。

上述したように設計システム１によれば、上記処理により全時刻歴データから象限毎に最大値又は最小値を抽出することにより、データ量を大幅に低減することができる。

以上、本発明の一実施形態について説明したが、本発明は上記の一実施形態に限定されるものではなく、その趣旨を逸脱しない範囲で適宜変更可能である。

１設計システム
２入力部
４演算部
６表示部
８記憶部

Claims

３次元の有限要素法を用いて設計対象の構造物の設計を支援する設計システムであって、
前記構造物の構造及び荷重に関する設計データに基づいて、地震時及び固定荷重を含む静的荷重に対して前記構造物に作用する応力が組み合わされた成分を３次元解析し、前記構造物に作用する軸力と曲げモーメントとの関係を示す軌跡を算出し、前記軌跡が含まれる領域について座標系の各象限毎の最大値又は最小値を抽出し、前記最大値又は前記最小値に基づいて前記構造物の断面に作用する応力を算出する演算部を備えることを特徴とする、
設計システム。
記憶部に記憶されたプログラムを演算部が実行することにより３次元の有限要素法を用いて設計対象の構造物の設計を支援する設計方法であって、
前記記憶部に記憶された前記構造物の構造及び荷重に関する設計データに基づいて、地震時及び固定荷重を含む静的荷重に対して前記構造物に作用する応力が組み合わされた成分を前記演算部で３次元解析し、
前記構造物に作用する軸力と曲げモーメントとの関係を示す軌跡を前記演算部で算出し、
前記軌跡が含まれる領域について座標系の各象限毎の最大値又は最小値を前記演算部で抽出し、
前記最大値又は前記最小値に基づいて前記構造物の断面に作用する応力を前記演算部で算出することを特徴とする、
設計方法。