特許 7530748 ＣＡＥ解析方法およびＣＡＥ解析装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-07-31

(45)【発行日】2024-08-08

(54)【発明の名称】ＣＡＥ解析方法およびＣＡＥ解析装置

(51)【国際特許分類】

   G06F 30/13 20200101AFI20240801BHJP

   G06F 30/20 20200101ALI20240801BHJP

【ＦＩ】

G06F30/13

G06F30/20

【請求項の数】 6

(21)【出願番号】P 2020100303

(22)【出願日】2020-06-09

(65)【公開番号】P2021196650

(43)【公開日】2021-12-27

【審査請求日】2023-04-04

(73)【特許権者】

【識別番号】504163612

【氏名又は名称】株式会社ＬＩＸＩＬ

(74)【代理人】

【識別番号】100106002

【弁理士】

【氏名又は名称】正林 真之

(74)【代理人】

【識別番号】100165157

【弁理士】

【氏名又は名称】芝 哲央

(74)【代理人】

【識別番号】100126000

【弁理士】

【氏名又は名称】岩池 満

(74)【代理人】

【識別番号】100160794

【弁理士】

【氏名又は名称】星野 寛明

(72)【発明者】

【氏名】安部 則弘

(72)【発明者】

【氏名】高橋 弥生

【審査官】合田 幸裕

(56)【参考文献】

【文献】特開２０１９－１６８８３８（ＪＰ，Ａ）

【文献】国際公開第２０１７／０４３５７１（ＷＯ，Ａ１）

【文献】中国特許出願公開第１０７０４５５６７（ＣＮ，Ａ）

【文献】中国特許出願公開第１０６０９６０７８（ＣＮ，Ａ）

【文献】吉川弘道 ほか，数値シミュレーションで考える構造解析 －ソフトで学ぶ非線形解析と応答解析－，第1版，日本，株式会社建通新聞社，2009年11月18日，pages 172-173

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０６Ｆ ３０／１０ － ３０／２８

ＩＥＥＥ Ｘｐｌｏｒｅ

ＪＳＴＰｌｕｓ（ＪＤｒｅａｍＩＩＩ）

(57)【特許請求の範囲】

【請求項1】

建具のフレーム部材に、前記フレーム部材と等価な曲げ剛性およびねじり剛性を有する板要素を適用してモデル化するモデル化工程と、
前記モデル化された建具のフレーム部材の強度解析を行う解析工程と、を有し、
前記モデル化工程は、
前記板要素を、前記板要素の短手方向の見付け幅が、前記フレーム部材の短手方向の見付け幅と等しい直方体状の板とし、さらに前記板の板厚を定義する形状定義過程を有する、ＣＡＥ解析方法。

【請求項2】

前記モデル化工程はさらに、
前記見付け幅と前記板厚から、前記板要素の断面２次モーメントおよびねじり定数を算出する第１算出過程と、
曲げ剛性をヤング率と断面２次モーメントの積で、ねじり剛性を横弾性係数とねじり定数の積で表し、前記フレーム部材と前記板要素について、曲げ剛性およびねじり剛性をそれぞれ等置して方程式を立式する等置立式過程と、
前記第１算出過程で算出された前記板要素の断面２次モーメントおよびねじり定数と前記方程式から、前記板要素のヤング率と横弾性係数を算出する第２算出過程と、
前記第２算出過程で算出されたヤング率と横弾性係数を板要素の材料物性として適用することで、板要素の曲げ剛性およびねじり剛性を定義する剛性定義過程と、を有する、請求項１に記載のＣＡＥ解析方法。

【請求項3】

前記解析工程は、
モデル化された前記板要素の曲げ応力をモデル化前の前記フレーム部材の断面係数で除することで、前記フレーム部材の曲げ応力を算出する曲げ応力換算過程を有する、請求項１または２に記載のＣＡＥ解析方法。

【請求項4】

建具のフレーム部材に、前記フレーム部材と等価な曲げ剛性およびねじり剛性を有する板要素を適用してモデル化するモデル適用部と、
前記モデル化された建具のフレーム部材の強度解析を行う解析部と、を有し、
前記モデル適用部は、
前記板要素を、前記板要素の短手方向の見付け幅が、前記フレーム部材の短手方向の見付け幅と等しい直方体状の板とし、さらに前記板の板厚を定義する形状定義部を有する、ＣＡＥ解析装置。

【請求項5】

前記モデル適用部はさらに、
前記見付け幅と前記板厚に基づいて前記板要素のヤング率と横弾性係数を算出する算出部と、
前記算出部で算出されたヤング率と横弾性係数を前記板要素の材料物性として適用して、前記板要素の曲げ剛性およびねじり剛性を定義する剛性定義部と、を有する、請求項４に記載のＣＡＥ解析装置。

【請求項6】

前記解析部は、
モデル化された前記板要素の曲げ応力をモデル化前の前記フレーム部材の断面係数で除することで、前記フレーム部材の曲げ応力を算出する曲げ応力換算部を有する、請求項４または５に記載のＣＡＥ解析装置。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、ＣＡＥ解析方法およびＣＡＥ解析装置に関する。

【背景技術】

【0002】

近年、サッシ業界においてはサッシフレームの断熱化・スリム化が進められており、サッシ設計にコンピュータ支援による解析（以下、ＣＡＥ解析）が欠かせない技術となっている。強風・突風に対するサッシの耐風圧強度性能をＣＡＥ解析で評価するには、フレーム部材の曲げやねじれ挙動を再現し、サッシフレーム部材やガラスの変形と応力、固定部品の反力などを計算する技術が求められる。

【0003】

サッシフレーム部材や金物部品を精細にモデル化した立体要素による３次元モデルを作成してもよいが、サッシ全体の精細な３次元モデルを作成した場合にはモデリング作業や計算時間に多大な工数を要する。

【0004】

特許文献１には、複数の部材間の溶接部についてシェル要素、ソリッド要素およびビーム要素を用いてモデル化するＣＡＥ解析手法が開示されている。これにより、部材間の溶接部についてモーメントの伝達を表現することを可能としながら簡易に適切なモデル化を行うことができるとされている。

【先行技術文献】

【特許文献】

【0005】

【文献】特開２０１５－１０９０２１号公報

【発明の概要】

【発明が解決しようとする課題】

【0006】

特許文献１の手法は、部材間の溶接部についてモーメントの伝達を表現しつつ簡易にモデル化を行うものであり、部材の曲げ剛性及びねじれ剛性に基づいて、部材自体に生じる曲げやねじれ挙動を再現するものではない。

【課題を解決するための手段】

【0007】

本開示は、建具のフレーム部材に、前記フレーム部材と等価な曲げ剛性およびねじり剛性を有する板要素を適用してモデル化するモデル化工程と、前記モデル化された建具のフレーム部材の強度解析を行う解析工程と、を有するＣＡＥ解析方法を提供する。

【図面の簡単な説明】

【0008】

【図1】本実施形態の建具を示す斜視図である。

【図2】本実施形態の建具の下框部材を示す斜視図である。

【図3】下框部材を板要素でモデル化したものを示す図である。

【図4】本実施形態の建具を板要素でモデル化したものを示す図である。

【図5】耐風圧試験のシミュレーションを示す図である。

【図6】本開示の解析方法を示すフローチャートである。

【図7】本開示の解析装置の構成を示すブロック図である。

【0009】

（解析対象）
本開示のＣＡＥ解析方法は、例えば図１に示すような建具１００に適用することができる。本実施形態の建具１００は引き違い窓であり、例えば面材６および面材６を保持するサッシ部分の耐風圧強度性能の評価に用いることができる。本開示のＣＡＥ解析方法の適用先は、引き違い窓に限定されない。框などのフレーム部材のような細長形状の部材を有する建具であれば、本開示のＣＡＥ解析方法を適用できる。

【0010】

建具１００は、枠体１と、上框２と、下框３と、縦框４と、召合せ框５と、面材６と、を有して構成される。上框２と、下框３と、縦框４と、召合せ框５は、面材６の四方を囲むように保持して引き違い窓の障子部分を構成し、枠体１は、引き違い窓の障子部分を枠内に保持する。枠体１や各框は、例えばアルミなどの金属や樹脂など、種々の材質で構成されるものを使用できる。面材６は、例えばガラスで構成される。

【0011】

本実施形態において解析対象としての建具１００は、細長形状の各部材で構成されるフレーム部材と、フレーム部材に囲まれる面材６と、に分けられる。本開示の解析方法は、細長形状のフレーム部材を、前記フレーム部材と等価な曲げ剛性およびねじり剛性を有する板要素でモデル化するＣＡＥ解析方法である。

【0012】

本実施形態において建具１００の細長形状の各部材とは、枠体１と、上框２と、下框３と、縦框４と、召合せ框５である。例えば下框３は図２に示すような細長形状の部材であり、その延びる方向に沿って、ガラス溝３１、レール溝３２、中間部３３を有する。ガラス溝３１は面材６の下端と対向する部分である。ガラス溝３１は面材６の端部を挿入されて、面材６を保持する。レール溝３２は枠体１の下部と対向する部分である。レール溝３２は枠体１から上方に向けて延びるレールが挿入され、枠体１にスライド可能に取り付けられる。中間部３３はガラス溝３１およびレール溝３２の中間に位置する部分であり、中空構造を有する。本実施形態においては、上框２と、下框３と、縦框４と、召合せ框５と、面材６についての解析を考える。

【0013】

（部材のモデル化）
建具１００の各框は、いずれも板要素でモデル化することができる。各框は、曲げ剛性およびねじり剛性が等価な板要素でモデル化され、各部材の構造特性を再現する。例えば図３に示すように、下框３は、板要素３ａでモデル化される。板要素３ａの見付け幅ｂは元の下框３の見付け幅に等しい。板厚ｔは見付け幅ｂに比例する値で定義される。例えば板厚ｔは、ｔ＝ｂ／１０で定義される。

【0014】

ガラス溝３１やレール溝３２、および中空構造を有する中間部３３を有し、複雑な断面形状を呈する下框３は、曲げ剛性（Ｅ・Ｉ）とねじり剛性（Ｇ・Ｊ）が等価な１枚の板要素３ａに単純化してモデル化することができる。以下にその手順を示す。

【0015】

実際の下框３の曲げ剛性とねじり剛性を、それぞれＥ・Ｉ、Ｇ・Ｊとする。ここで、ＥおよびＧは、それぞれ下框３のヤング率および横弾性係数を表し、Ｉは当該部材の断面二次モーメント、Ｊはねじり定数を表す。ＥおよびＧは下框３を構成する材料による固有の値である。ＩおよびＪは下框３の形状等から算出することができる。

【0016】

一方、モデル化した板要素３ａの曲げ剛性とねじり剛性をそれぞれＥ'・Ｉ'、Ｇ'・Ｊ'とする。ここで、Ｅ'およびＧ'は、モデル化した板要素３ａのヤング率および横弾性係数を表し、Ｉ'はモデル化した板要素３ａの断面二次モーメント、Ｊ'はねじり定数を表す。

【0017】

下框３の見付け幅ｂを用いて、板の肉厚ｔを例えばｔ＝ｂ／１０と定義する。板における断面二次モーメントおよびねじり定数の定義と、板要素３ａの形状変数ｂおよびｔから、Ｉ'とＧ'はそれぞれ
Ｉ’＝（ｂ・ｔ^３）／１２＝ｂ^４／１２０００ …（１）
Ｊ’＝（ｂ・ｔ^３）／３＝ｂ^４／３０００ …（２）
で算定できる。

【0018】

モデル化前後で曲げ剛性およびねじり剛性が等しいとすれば、Ｅ・Ｉ＝Ｅ'・Ｉ'、Ｇ・Ｊ＝Ｇ’・Ｊ'である。したがって残る未知数Ｅ'およびＧ'はそれぞれ、（１）式、（２）式および既知の変数Ｅ、Ｇ、Ｉ、Ｊ、ｂ、ｔ（＝ｂ／１０）を用いて、
Ｅ’＝Ｅ・Ｉ／Ｉ'＝１２（Ｅ・Ｉ）／（ｂ・ｔ^３）＝１２０００（Ｅ・Ｉ）／ｂ^４
…（３）
Ｇ’＝Ｇ・Ｊ／Ｊ'＝３（Ｇ・Ｊ）／（ｂ・ｔ^３）＝３０００（Ｇ・Ｊ）／ｂ^４
…（４）
で表される。

【0019】

上記算定式にて板要素３ａのヤング率Ｅ’と横弾性係数Ｇ’を求め、これを板に置換された下框３の材料物性として定義すれば、曲げ剛性とねじり剛性が等価な板要素３ａに置換された下框３のＣＡＥ解析モデルを作成することができる。

【0020】

図４に示すように、上框２、縦框４および召合せ框５についても同様にして、それぞれ板要素２ａ、４ａおよび５ａとしてモデル化することができる。板要素２ａ、４ａおよび５ａの見付け幅をｂ_２、ｂ_４およびｂ_５とし、板要素３ａの見付け幅をｂ_３とする。各板要素について、板厚ｔ_ｉ＝ｂ_ｉ／１０の関係が成り立つものとする（ｉは２から５までの整数）。

【0021】

図５に示すように、建具１００に板要素を適用したモデルには、板要素上にグリッドで示す複数の節面が定義される。この節面ごとに載荷される荷重の大きさや変位をコンピュータ解析することで、建具全体の強度性能を得ることができる。

【0022】

モデリングを板要素のような単純な要素で行い、当該板要素に元の部材の曲げ剛性及びねじり剛性を適用することで、部材の曲げやねじれなどの挙動を再現可能となっている。これにより、各框の変形による影響を加味しつつ、精細な３次元モデルの作成にかかる時間と工数を短縮し建具１００の耐風圧強度性能を短時間で計算・確認することができる。

【0023】

一度モデル化した建具の部材を一部変更する場合には、変更した部材の断面性能値を適用して修正することで、新たな物体について簡単にモデル化することができる。すなわち、引き違い窓において框の断面形状や肉厚などを設計変更した場合には、作成済みの引き違い窓モデルについて変更した框の断面性能値を適用すれば、瞬時に新たなモデルに反映されるので、モデリング時間や工数が削減され好適である。

【0024】

本実施形態のＣＡＥ解析方法では、フレーム部材のたわみやねじれなどの変形量はＣＡＥ解析の計算結果をそのまま採用することができるが、曲げ応力の算定には以下の換算が必要となる。

【0025】

実際のフレーム部材の曲げ応力を求めるには、先ず、モデル化した板の曲げ応力σ'から、曲げモーメントＭ＝ｂ・ｔ^２／６・σ'を求める。これを実際のフレーム部材の断面係数Ｚで除すれば、曲げ応力σ＝Ｍ／Ｚが算定できる。これにより、フレーム部材の実機において、塑性変形リスクまたは許容応力に対する余裕度を考慮できるとともに、過剰品質を避けて経済設計を行うことができる。

【0026】

（解析装置）
本実施形態に係るＣＡＥ解析を実行する装置を、ＣＡＥ解析装置２０と称する。図７に示すようにＣＡＥ解析装置２０は、ＣＡＤデータ取得部２１と、モデル適用部２２と、解析部２３と、を有して構成される。

【0027】

ＣＡＤデータ取得部２１は、予め用意された建具１００の各部材の解析対象領域を表すＣＡＤデータを取得し、モデル化処理用のデータに変換する。モデル適用部２２は、取得した建具１００の解析領域を示すデータに基づき、図５に示すように当該解析領域に存在する各部材のモデルとして板要素を適用する。解析部２３は、モデル化された建具１００全体について解析を行う。

【0028】

モデル適用部２２は、板要素の形状および寸法を定義して設定する形状定義部２２１と、板要素の見付け幅と板厚に基づいて、フレーム部材と曲げ剛性及びねじり剛性が等しくなる板要素のヤング率と横弾性係数を算出する算出部２２２と、算出部２２２で算出されたヤング率と横弾性係数を板要素の材料物性として適用することで板要素の曲げ剛性およびねじり剛性を定義する剛性定義部２２６と、を有する。

【0029】

算出部２２２は、板要素の見付け幅と板厚から板要素の断面２次モーメントおよびねじり定数を算出する第１演算部２２３と、フレーム部材と板要素について、曲げ剛性およびねじり剛性をそれぞれ等置して方程式を立式する等置立式部２２４と、前記方程式に基づいて板要素のヤング率と横弾性係数を算出する第２演算部２２５と、を有する。

【0030】

解析部２３は、実際のフレーム部材の曲げ応力を換算する曲げ応力換算部２４を有する。曲げ応力換算部２４は、モデル化された板要素の曲げ応力σ'から求めた曲げモーメントＭを実際のフレーム部材の断面係数Ｚで除することによって、実際のフレーム部材の曲げ応力σ＝Ｍ／Ｚを換算する。

【0031】

（解析方法）
本開示のＣＡＥ解析方法は、細長形状のフレーム部材を有する建具の強度性能を解析する方法であって、モデル化工程と、解析工程と、を有する。

【0032】

モデル化工程は、板要素の形状および寸法を定義して設定する形状定義過程と、板要素の断面２次モーメントおよびねじり定数を算出する第１算出過程と、フレーム部材と板要素と、の曲げ剛性およびねじり剛性を等置して方程式を立式する等置立式過程と、板要素のヤング率と横弾性係数を算出する第２算出過程と、第２算出過程で算出されたヤング率と横弾性係数を板要素の材料物性として適用することで、板要素の曲げ剛性およびねじり剛性を定義する剛性定義過程と、を有する。解析工程は、フレーム部材の曲げ応力を算出する曲げ応力換算過程を有する。各工程および過程について、図６に示すフローチャートに沿って説明する。

【0033】

モデル化工程では、ＣＡＤデータ取得部２１が解析対象物体のＣＡＤデータを受け取ってモデル化処理用のデータに変換する（ステップＳ１）。次いで、取得したＣＡＤデータに基づいて、図７に示すように当該解析領域に存在する各部材のモデルとして板要素を適用する（ステップＳ２）。このとき、次の過程に従って板要素の寸法を定義するとともに、もとのフレーム部材と同じ曲げ剛性およびねじり剛性を適用する。

【0034】

形状定義過程では、形状定義部２２１が、フレーム部材を構成する框を、同じ見付け幅および長さを有する直方体状の板で置換し、さらに板厚を見付け幅に比例する値として仮定して定義する。

【0035】

第１算出過程では、第１演算部２２３が、形状定義過程で設定した板要素の寸法と、上記（１）式および（２）式に基づいて、前記板要素の見付け幅と板厚から、板要素の断面２次モーメントおよびねじり定数を算出する。この時点で、板要素のヤング率と横弾性係数は未知数である。

【0036】

等置立式過程では、等置立式部２２４が、フレーム部材および対応する板要素の、曲げ剛性をヤング率と断面２次モーメントの積で、ねじり剛性を横弾性係数とねじり定数の積で表す。これを用いて、フレーム部材と板要素について、曲げ剛性およびねじり剛性をそれぞれ等置して方程式（上記（３）および（４）式）を立式する。

【0037】

第２算出過程では、第２演算部２２５が、第１算出過程で算出された板要素の断面２次モーメントおよびねじり定数と前記方程式から、未知数であった板要素のヤング率と横弾性係数を算出する。

【0038】

剛性定義過程では、剛性定義部２２６が、算出された板要素のヤング率と横弾性係数をモデル化された板要素の材料物性として適用する。これにより、板要素にフレーム部材と等しい曲げ剛性およびねじり剛性が定義される。

【0039】

解析工程では、解析部２３が、モデル化された建具１００について風圧載荷時の荷重伝達を再現して解析を行い、結果を出力する（ステップＳ５）。図５に示す解析シミュレーション結果では、グリッドの色が明るくなるほど風圧載荷時のモデルの変位が大きいことを示している。

【0040】

フレーム部材の曲げ剛性およびねじり剛性をモデルに反映したことにより、曲げやねじれなどの挙動を再現した解析が可能になるとともに、解析モデルを簡易化して短時間で強度解析できる。また、部材変更の際には変更後の部材の断面性能値を反映させるだけで、瞬時に新たなモデルに反映される。これにより、実機試験の試験数を削減し、開発期間の短縮が期待できる。

【0041】

さらに解析工程は、実際のフレーム部材の曲げ応力を換算する曲げ応力換算過程を有する。曲げ応力換算過程では、モデル化された板要素の曲げ応力σ'から求めた曲げモーメントＭを実際のフレーム部材の断面係数Ｚで除することによって、実際のフレーム部材の曲げ応力σ＝Ｍ／Ｚを換算する(ステップＳ６)。

【0042】

本実施形態のＣＡＥ解析方法について説明した。本実施形態の解析方法によれば、以下の効果を奏する。

【0043】

本実施形態のＣＡＥ解析方法は、建具１００のフレーム部材（上框２、下框３、縦框４、召合わせ框５）に、前記フレーム部材と等価な曲げ剛性およびねじり剛性を有する板要素を適用してモデル化するモデル化工程と、前記モデル化された建具のフレーム部材の強度解析を行う解析工程と、を有する。これにより、曲げやねじれなどの挙動を再現した解析が可能になるとともに、解析モデルを簡易化して短時間で強度解析できる。実機試験の試験数の削減や、開発期間の短縮が期待できる。

【0044】

前記モデル化工程は、前記板要素を、その見付け幅が前記フレーム部材の見付け幅と等しい直方体状の板とし、さらに前記板の板厚を定義する形状定義過程と、前記見付け幅と前記板厚から、前記板要素の断面２次モーメントおよびねじり定数を算出する第１算出過程と、曲げ剛性をヤング率と断面２次モーメントの積で、ねじり剛性を横弾性係数とねじり定数の積で表し、前記フレーム部材と前記板要素について、曲げ剛性およびねじり剛性をそれぞれ等置して方程式を立式する等置立式過程と、前記第１算出過程で算出された前記板要素の断面２次モーメントおよびねじり定数と前記方程式から、前記板要素のヤング率と横弾性係数を算出する第２算出過程と、前記第２算出過程で算出されたヤング率と横弾性係数を板要素の材料物性として適用することで、板要素の曲げ剛性およびねじり剛性を定義する剛性定義過程と、を有する。これにより、算出した曲げ剛性およびねじり剛性を板要素に反映することができる。

【0045】

前記解析工程は、モデル化された前記板要素の曲げ応力をモデル化前の前記フレーム部材の断面係数で除することで、前記フレーム部材の曲げ応力を算出する曲げ応力換算過程を有する。これにより、実際のフレーム部材の曲げ応力を算出することができる。

【0046】

本実施形態のＣＡＥ解析装置２０は、建具１００のフレーム部材に、前記フレーム部材と等価な曲げ剛性およびねじり剛性を有する板要素を適用してモデル化するモデル適用部２２と、モデル化された建具１００のフレーム部材の強度解析を行う解析部２３と、を有する。これにより、曲げやねじれなどの挙動を再現した解析が可能になるとともに、解析モデルを簡易化して短時間で強度解析できる。実機試験の試験数の削減や、開発期間の短縮が期待できる。

【0047】

モデル適用部２２は、前記板要素を、その見付け幅が前記フレーム部材の見付け幅と等しい直方体状の板とし、さらに前記板の板厚を定義する形状定義部２２１と、前記見付け幅と前記板厚に基づいて前記板要素のヤング率と横弾性係数を算出する算出部２２２と、算出部２２２で算出されたヤング率と横弾性係数を前記板要素の材料物性として適用して、前記板要素の曲げ剛性およびねじり剛性を定義する剛性定義部２２６と、を有する。これにより、算出した曲げ剛性およびねじり剛性を板要素に反映することができる。

【0048】

解析部２３は、モデル化された前記板要素の曲げ応力をモデル化前の前記フレーム部材の断面係数で除することで、前記フレーム部材の曲げ応力を算出する曲げ応力換算部２４を有する。これにより、実際のフレーム部材の曲げ応力を算出することができる。

【0049】

本開示のＣＡＥ解析方法およびＣＡＥ解析装置は、上記の実施形態に限定されない。例えば建具１００以外にも、細長形状の部材を有する物体であれば、本開示のＣＡＥ解析方法を適用することができる。また、本開示のＣＡＥ解析方法を耐風圧強度以外の強度解析に適用してもよい。

【符号の説明】

【0050】

１００…建具、１…枠体、２…上框、３…下框、３１…ガラス溝、３２…レール溝、３３…中間部、４…縦框、５…召合せ框、６…面材、２ａ…モデル化後の上框、３ａ…モデル化後の下框、４ａ…モデル化後の縦框、５ａ…モデル化後の召合せ框、２０…ＣＡＥ解析装置、２１…ＣＡＤデータ取得部、２２…モデル適用部、２２１…形状定義部、２２２…算出部、２２３…第１演算部、２２４…等置立式部、２２５…第２演算部、２２６…剛性定義部、２３…解析部、２４…曲げ応力換算部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】