4)外力の設定
Node 3,4 に垂直方向にかかる力は、
OpenSees Basic Examples (曲げモーメント-曲率計算) 部材の定義3
において、
# Set axial load
set P -180
と、定義しました、180です。これは、
「10% of axial capacity of columns」ということで、重力加速度による力として全体の 10%の力として180という値がセットされています。
Opensees のコマンドは、以下の通りです。
各コマンドの説明は、「OpenSeesコマンド集」を参照してください
# Define gravity loads
# --------------------
# Set a parameter for the axial load
set P 180.0; # 10% of axial capacity of columns
# Create a Plain load pattern with a Linear TimeSeries
pattern Plain 1 "Linear" {
# Create nodal loads at nodes 3 & 4
# nd FX FY MZ
load 3 0.0 [expr -$P] 0.0
load 4 0.0 [expr -$P] 0.0
}
5)解析の設定
解析する時の、連立方程式の解き方の手法などを設定します。
各コマンドの説明は、「OpenSeesコマンド集」を参照してください
# ------------------------------
# Start of analysis generation
# ------------------------------
# Create the system of equation, a sparse solver with partial pivoting
system BandGeneral
# Create the constraint handler, the transformation method
constraints Transformation
# Create the DOF numberer, the reverse Cuthill-McKee algorithm
numberer RCM
# Create the convergence test, the norm of the residual with a tolerance of
# 1e-12 and a max number of iterations of 10
test NormDispIncr 1.0e-12 10 3
# Create the solution algorithm, a Newton-Raphson algorithm
algorithm Newton
# Create the integration scheme, the LoadControl scheme using steps of 0.1
integrator LoadControl 0.1
# Create the analysis object
analysis Static
# ------------------------------
# End of analysis generation
# ------------------------------
10回分解析し、10回後のNode3 4 および 要素 1 の結果を表示します。
# ------------------------------
# Finally perform the analysis
# ------------------------------
# perform the gravity load analysis, requires 10 steps to reach the load level
analyze 10
# Print out the state of nodes 3 and 4
print node 3 4
# Print out the state of element 1
print ele 1
次のページ →
OpenSees Basic Examples (鉄筋コンクリート門型フレームの柱軸力のみを作用させた解析) コマンドのサマリ
Opensees コマンド 日本語解説 : 目次 |
0 件のコメント:
コメントを投稿