入力された外力(荷重データ)から、固定されていない成分のみのデータを作成します。
これは、<有限要素法>固定条件の処理 で作成したconstraint_sequence[ ] を再利用します。
最初に、外部から入力されたデータから、固定されているNodeの外力x成分loading_point_x[ ];、及び y成分loading_point_x[ ];をシーケンシャルに並べますloading_force_TF[ ]。
for (point_index=0; point_index<point_num; point_index++) {
loading_force_TF[2*point_index] = loading_point_x[point_index];
loading_force_TF[2*point_index+1] = loading_point_y[point_index];
}
constraint_sequence[ ]を使用して、計算する必要のない成分(=固定されている成分)の外力成分を抜き去っていき、valid_loading_force[ ] を作成します。
for (int i = 0; i < valid_point_num; i++)
valid_loading_force[i]=loading_force_TF[constraint_sequence[i]];
以上、 <有限要素法>荷重データの処理 のプログラムをまとめたものは、以下の様になります。
(途中、いらなくなったメモリを解放しています)
次のページ →
<有限要素法>変位の計算
「有限要素法における平面トラスの解析」プログラムをC言語で書いてみる : 目次 |
0 件のコメント:
コメントを投稿