本資料のうち，枠囲みの内容は，営業秘密又は防護上の観点から公開できません。

東海第二発電所 工事計画審査資料
資料番号

工認-387 改3

提出年月日

平成30年7月19日

NT2 補③ Ｖ-5-3 R0

Ⅴ-5-3 計算機プログラム（解析コード）の概要・ＳＡＰ-Ⅳ

目次

NT2 補③ Ｖ-5-3 R1

1.

はじめに

1

1.1 使用状況一覧

22.

3解析コードの概要

2.1 SAP-Ⅳ CNDYN Ver.4.1

3

2.2 SAP-Ⅳ Ver.1.00

43.

5SAP-Ⅳ Ver.1.00の解析手法について

3.1 一般事項

5

3.2 解析コードの特徴

5

3.3 解析手法

53.4 解析フローチャート

7

3.5 検証（Verification）及び妥当性確認（Validation）

8

1.

はじめに
本資料は，添付書類において使用した計算機プログラム（解析コード）SAP-Ⅳについて説明するものである。

NT2 補③ Ｖ-5-3 R0

本解析コードを使用した添付書類を示す使用状況一覧，解析コードの概要を以降に記載する。

1

1.1 使用状況一覧
使用添付書類バージョン

Ⅴ-2-4-2-2

使用済燃料貯蔵ラックの耐震性についての計算書

CNDYN Ver.4.1

Ⅴ-2-5-4-1-1

残留熱除去系熱交換器の耐震性についての計算書

CNDYN Ver.4.1

Ⅴ-2-6-3-2-1

水圧制御ユニットの耐震性についての計算書

CNDYN Ver.4.1

Ⅴ-2-6-5-1

起動領域計装の耐震性についての計算書

CNDYN Ver.4.1

Ⅴ-2-11-2-1

燃料取替機の耐震性についての計算書

CNDYN Ver.4.1

Ⅴ-2-11-2-4

チャンネル着脱機の耐震性についての計算書

CNDYN Ver.4.1

Ⅴ-2-11-2-7

制御棒貯蔵ハンガの耐震性についての計算書

CNDYN Ver.4.1

Ⅴ-2-8-3-3-1

緊急時対策所換気系ダクトの耐震性についての計算書

Ver.1.00

Ⅴ-2-別添1-6

供給配管の耐震計算書(統合版)Ver.8.0
rev3

Ⅴ-2-別添2-2

溢水源としない耐震B，Cクラス機器の耐震性についての (統合版)Ver.8.0
計算書

NT2 補③ Ｖ-5-3 R1

rev32

2.

解析コードの概要
2.1 SAP-Ⅳ CNDYN Ver.4.1
コード名
項目
使用目的

SAP-Ⅳ
3次元有限要素法（シェルモデル）による固有値解析及び応力解析
3次元有限要素法（はりモデル）による固有値解析及び地震応答解析

開発機関

米国カリフォルニア大学

開発時期

1973年

使用したバージョン

CNDYN Ver.4.1

コードの概要SAP-Ⅳ（以下，「本解析コード」という。）は，任意形状の三次元モデル（主にはり要素及びシェル要素）に対して，有限要素法を用いて静的解析及び動的解析を行うもので，主として，機器の固有値計算並びに自重，運転時荷重及び地震力による応力計算等に用いる。

NT2 補③ Ｖ-5-3 R2

本解析コードは，機械工学，土木工学，航空工学等の分野において，多くの実績を有している。
【検証（Verification）】
・両持ちはりの単純支持円筒モデルについて，SAP-Ⅳによる解析結果と理論解とを比較して検討し，SAP-Ⅳによる解析結果が妥当であることを確認している。
・平板のモデルについて，シェルモデルによる固有値解析及び応力解析を行い,SAP-Ⅳによる解析結果と理論解とを比較して検討し，SAPⅣによる解析結果が妥当であることを確認している。
・本解析コードの運用環境について、開発機関から提示された要件を検証満足していることを確認している。

（Verification）【妥当性確認（Validation）】

及び妥当性確認（Validation）

・本解析コードは，機械工学，土木工学，航空工学等の分野において，多くの実績を有しており，妥当性は十分に確認されている。
・SAP-Ⅳ CNDYN Ver.4.1のマニュアルにより，今回の工認申請で使用する3次元有限要素法（シェルモデル）による固有値解析及び応力解析並びに3次元有限要素法（はりモデル）による固有値解析及び地震応答解析に，本解析コードが適用できることを確認している。
・今回の工事計画における構造に対し使用する要素及び解析については，工事の計画において使用された実績がある。
・今回の工事計画認可申請において使用するバージョンは，工事の計画において使用されているものと同じであることを確認している。

32.2 SAP-Ⅳ Ver.1.00
コード名
項目

SAP-Ⅳ

使用目的

三次元有限要素法（はり要素）による固有値解析及び応力解析

開発機関

米国カリフォルニア大学

開発時期

1973 年使用したバージョン

Ver.1.00
SAP-Ⅳ（以下，「本解析コード」という。）は，任意形状の三次元モコードの概要

デル（主にはり要素及びシェル要素）に対して，有限要素法を用いて静的解析及び動的解析を行うもので，主として，機器の固有値計算並びに自重，運転時荷重及び地震力による応力計算等に用いる。
本解析コードは，機械工学，土木工学，航空工学等の分野において，多くの実績を有している。

NT2 補③ Ｖ-5-3 R0【検証（Verification）】
・入力されたデータはインプットデータとして解析結果とともに出力され，入力データと一致することを確認している。
・片持ちばりの自重による固定端モーメント及び自由端たわみ，固有振動数を本解析コードの静的解析結果及び固有値解析結果と理論解を比較して検討し，本解析コードによる解析結果が妥当であることを確認している。
・本解析コードの運用環境について，開発機関から提示された要件を検証満足していることを確認している。

（Verification）【妥当性確認（Validation）】

及び・本解析コードは，機械工学，土木工学，航空工学等の様々な分野に妥当性確認（Validation）おける使用実績を有しており，妥当性は十分に確認されている。
・原子力産業界において，工認申請範囲外の支持構造物に対する耐震性評価に本解析コードが使用された実績があることを確認している。
・今回の工認申請で行うはりモデルの固有値解析及び静的解析という解析の使用目的に照らして，用途及び適用範囲が上述の妥当性確認範囲であることを確認している。
・開発機関が提示するマニュアルにより今回の工認申請で使用する三次元有限要素法（はり要素）による固有値解析及び応力解析に本解析コードが適用できることを確認している。

4

3.

SAP-Ⅳの解析手法について
3.1 一般事項
本書は，米国カリフォルニア大学において開発された汎用解析コード SAP-Ⅳの説明書である。
本解析コードは，有限要素法を用いて任意形状の三次元系の静的解析及び動的解析を行うもので，主として，機器・配管系の自重，運転時荷重及び地震力による応力計算に用いる。
3.2 解析コードの特徴
SAP-Ⅳは，有限要素法を用いて静的解析及び動的解析を行う三次元構造解析用の計算機プ
ログラムである。
主な特徴を下記に示す。
・二次元及び三次元有限要素プログラムである。
・モデル要素として九種類の要素を扱うことができる。
・静的解析後に，同じモデルを使い，解析条件及び荷重条件を変更し動的解析を行う

NT2 補③ Ｖ-5-3 R0

ことができる。
3.3 解析手法
3.3.1 静的解析
要素の平衡方程式を組み立て，構造物全体に対して次の平衡方程式を作り，これを解く。
Ｋｕ＝Ｒ
ここで，
Ｋ：剛性マトリックス
ｕ：変位ベクトル
Ｒ：荷重ベクトル
計算された節点の変位から，要素の変形及び要素の応力を求める。

5

3.3.2 動的解析
(1) 固有値解析
次の固有値方程式を解いて，固有振動数と振動モードを求める。
Ｋφ＝ω2 Ｍφ
ここで，Ｋ：剛性マトリックス
φ：固有ベクトル
ω：固有振動数
Ｍ：質量マトリックス
上記方程式の解法として，剛性マトリックスと質量マトリックスが小さい場合にはデターミナント法を用い，剛性マトリックスと質量マトリックスが大きい場合はサブスペース反復法が用いられる。

NT2 補③ Ｖ-5-3 R0(2) 動的応答解析
動的応答解析においては，次の方程式を用いる。
・・・
Ｍｕ＋Ｃｕ＋Ｋｕ＝Ｒ(ｔ)ここで，Ｍ：質量マトリックス
・・
ｕ：加速度ベクトル・ｕ：速度ベクトル
ｕ：変位ベクトル
Ｃ：減衰マトリックス
Ｋ：剛性マトリックス
Ｒ(ｔ)：時刻歴の外荷重ベクトル
上記方程式を用いて，モーダル解析法による時刻歴応答解析，直接積分法による時刻歴応答解析ができる。

6

3.4 解析フローチャート
本解析コードを用いた解析フローチャートを図 3-1 に示す。

ＳＴＡＲＴ

入力データ整理NT2 補③ Ｖ-5-3 R0

解析種別

静的解析

固有値解析

節点変位，荷重の計算固有振動数の計算

出力

ＥＮＤ

図3-1 解析フローチャート

7

3.5 検証（Verification）及び妥当性確認（Validation）
3.5.1 検証（Verification）
今回の解析に用いた解析コード SAP-Ⅳの検証として，SAP-Ⅳによる解析結果と手計算による理論解との比較検証を行う。
(1) 解析ケース
表 3-1 に示す 2 ケースについて解析を行い，理論解と比較する。
表 3-1 検証ケース
ケース

解析種別

内容

要素条件

はり要素

角形鋼管，一様密度

片持ちばりの自重による1

静的解析

固定端モーメント及び自由端たわみ

2

固有値解析

片持ちばりの固有振動数

NT2 補③ Ｖ-5-3 R0(2) 解析条件
使用状況一覧における支持構造物の代表的な部材を使用した解析モデルを使用する。
図 3-2 に解析条件を，図 3-3 に解析モデルを示す。
長さＬ(mm)

断面形状

−(mm)

断面積

Ａ(mm2)

Ｉ(mm4)

縦弾性係数

Ｅ(N/mm2)

ｔ：板厚質量密度ρ(ton/mm3)

材料：

重力加速度g(m/s2)

片端部固定断面二次モーメントＬ

図3-2 解析条件

図3-3 解析モデル（使用要素：はり要素）

8

(3) 解析結果
a. 解析結果と理論解との比較（ケース1）
表 3-2 に解析結果と理論解との比較結果を示す。表 3-2 に示すとおり，解析結果と理論解はよく一致しており，解析コード SAP-Ⅳが検証されていることを確認した。
表 3-2 解析結果と理論解の比較
理論解＊

比較項目

解析結果

誤差

Ｍｍａｘ
（固定端モーメント）
δｍａｘ（自由端たわみ）σ（発生応力）NT2 補③ Ｖ-5-3 R0

注記＊：理論解は下記の式で計算する。
（出典：機械工学便覧）
ｗ Ｌ ，ｗＭｍａｘ
2ｗ Ｌdｍａｘ
8 Ｅ Ｉr Ａ ｇg 1000

（出典：設計・建設規格）
sｘ ＋3 ｘｙb. 解析結果と理論解との比較（ケース2）
表 3-3 に解析結果と理論解の固有振動数の比較結果を示す。表 3-3 に示すとおり，解析結果を理論解はよく一致しており，解析コード SAP-Ⅳが検証されていることを確認した。
表 3-3 解析結果と理論解の比較（単位：Hz）
次数

理論解＊

解析結果

誤差

注記＊：理論解は下記の式で計算する。
（出典：機械工学便覧）
ｆl

Ｅ Ｉ

2 p Ｌ

Ａ r

9

3.5.2 妥当性確認（Validation）
本解析コードの妥当性確認の内容は，以下のとおりである。
本解析コードは，機械工学，土木工学，航空工学等の様々な分野における使用実績を有しており，妥当性は十分に確認されている。
原子力産業界において，工認申請範囲外の支持構造物に対する耐震性評価に本解析コードが使用された実績があることを確認している。
今回の工認申請で行うはりモデルの固有値解析及び静的解析という解析の使用目的に照らして，用途，適用範囲が上述の妥当性確認範囲であることを確認している。
開発機関が提示するマニュアルにより今回の工認申請で使用する三次元有限要素法（はり要素）による固有値解析，応力解析に本解析コードが適用できることを確認している。
3.5.3 評価結果
3.5.1 及び 3.5.2 より，本解析コードを使用目的に示す固有値解析及び応力解析に用NT2 補③ Ｖ-5-3 R0Eいることは妥当である。

10