OpenSees Blog 日本語 目次

最初に「このブログをみるためのガイド」をご覧ください。

Update中OpenSees コマンド 日本語解説 の 目次 OpenSeesコマンド はじめての方は「ここから
このblogで使用しているOpenSeesコマンド集は「ここ」 Update中
OpenSees のソースコード 解析に挑戦してみる 目次最初のページは「ソースコードのダウンロードとビルド」
Appendix:C言語での 「有限要素法における平面トラスの解析」目次最初のページは「Microsoft Visual Studioの導入方法」







目次の中で、更新したページにはNewがついています

このブログ内の単語を検索したい場合は、左上OpenSeesロゴの上に検索窓から検索できます。


2017年1月18日水曜日

OpenSees Basic Examples ( D2: 2階建て単スパンフレームでの固有値解析) 解説(Node, Element)



0)必要となる数字をシンボルに代入しておき、実行環境を生成します

#input
set m  [expr 100.0/386.0]
set numModes 2      # 解析する固有値モードの数

#material
set A 63.41    # Element 1,2,3,4,5,6の断面積
set I 320.0       #Element 2,4,6の慣性モーメント
set E 29000.0  # ヤング率

#geometry
set L 240.
set h  120.

# create data directory "mkdir"のディレクトリを生成します
file mkdir modes; 



1)Model, Node(節点), Element(要素)の設定
解析する時の、次元、節点(Node)、要素(Element)の
Opensees のコマンドは、以下の通りです。

1.1) 2次元、自由度3 で、解析するよう設定します
model BasicBuilder [次元数] [自由度数] 

model BasicBuilder -ndm 2 -ndf 3

(参照ページはここ

1.2) Node(節点) の 座標を設定します
Node [ nodeの番号 ] [ x座標] [y座標] -mass [節点質量(自由度の第1成分:u )]
[節点質量(自由度の第2成分:v )] [節点質量(自由度の第3成分:θ )]

$がついた変数は、0)の所で定義したシンボルの数字がはいります。


# nodal coordinates:
node 1   0.  0. ;
node 2   $L  0. ;
node 3   0.  $h ;
node 4   $L  $h ;
node 5   0.  [expr 2*$h];    #  = 2 x h 
node 6   $L  [expr 2*$h];       #  = 2 x h 

1.3) Nodeのどこを固定するかを設定します
ここでは、nodeの1と2を固定します。

fix [ nodeの番号 ] [ uの固定の定義]  [ vの固定の定義]  [ θの固定の定義] 
0: 固定されていない
1: 固定されている


# Single point constraints -- Boundary Conditions
fix 1 1 1 1;
fix 2 1 1 1;


(参照ページはここ

1.4) Nodeの質量の定義

Nodeの質量を指定します。

Node3,4  Node5,6 の水平方向のみの質量が定義されています

mass [ nodeの番号 ] [ uの質量]  [ vの質量]  θ回転慣性質量] 


# assign mass
mass 3 $m 0. 0. ;
mass 4 $m 0. 0. ;
mass 5 [expr $m/2.] 0. 0. ;
mass 6 [expr $m/2.] 0. 0. ;

(参照ページはここ

1.5) 座標変換の定義

geomTransf [座標変換の方法] [座標変換方法のID]


# define geometric transformation:
set TransfTag 1;
geomTransf Linear $TransfTag ;


(参照ページはここ



1.6) Element(要素)を設定します。

Element コマンド の中の elasticBeamColum を使用します。

element elasticBeamColumn [ elementの番号 ] [ 片方のnode番号] [片方のnode番号] [ 断面積] [ヤング率] [Z軸の慣性モーメント] [座標変換方法のID]


# define elements:
# columns
element elasticBeamColumn 1 1 3 $A $E [expr 2.*$I] $TransfTag;
element elasticBeamColumn 2 3 5 $A $E $I           $TransfTag;
element elasticBeamColumn 3 2 4 $A $E [expr 2.*$I] $TransfTag;
element elasticBeamColumn 4 4 6 $A $E $I           $TransfTag;
# beams
element elasticBeamColumn 5 3 4 $A $E [expr 2.*$I] $TransfTag;
element elasticBeamColumn 6 5 6 $A $E $I           $TransfTag;


(参照ページはここ



次のページ →
OpenSees Basic Examples ( D3: 2階建て単スパンフレームでの固有値解析)) 解説(固有値解析)


Opensees コマンド 日本語解説 : 目次

0 件のコメント:

コメントを投稿