solid axisymmetric

Three Dimensional Finite Element Method
Solid Axisymmetric-9

■プログラムSTATIC8QFXAXISYMM-NEW.FOR■
入力データを元に変位u(r,z)とw(r,z)を解析するプログラムです。本プログラムの構成は、2次元の静解析と殆ど同じです。ただし、積分される内容はかなり違います。とりあえずプログラムで使われているサブプログラムとその機能を紹介しておきます。

Table 1 out of 2

Subroutine名	与える変数	送り返す変数	役割/目的
GRULE	N	SAI,WT	Gauss-Legendre積分法のサンプリング座標値と重み値を返す
DERIV	ND,INTEPT,SAI	BPP	形状関数の微分を Gauss-Legendre積分法のサンプリング座標値で計算する 4節点1次、8節点2次、9節点2次、 12節点3次要素が使える
ISOPARA	ND,E1,E2	F	形状関数の値を座標(E1,E2)で計算する 4節点1次、8節点2次、9節点2次、 12節点3次要素が使える
SHAPEF	ND,INTEPT,SAI	SF	形状関数の値を Gauss-Legendre積分法のサンプリング座標値で計算する
GSM	VISCO,FLMDA,GAMMA MXE,MXN,MXNA,MXWA INTEPT,ND,BPP WT,NE,NNODE,NBWA XCOORD,NODEX	AM	剛性マトリクス[K]を作るための積分を行う
SOLVE	MXE,MXN,MXBA,NBWA ND,NE,NNODE,NODEX AM,RHS,MXNA,MXWA NFORCE,IBFORCE BVFORCE,NBF IBNDF,BVF	U	サブプログラムNUEMANN、FORM、 SYSTEMを callし変位を計算する
NUEMANN	MXNA,MXBA,RHS NFORCE,IBFORCE BVFORCE	RHS	ベクトル{RHS}にNuemann型境界条件をはめ込む
FROM	MXNA,MXB,MXWA NNODE2,NBWA,AM RHS,NBF,BVF,IBNDF	AM RHS	[K]マトリクスにDirichlet型境界条件をはめ込む
SYSTEM	MXNA,MXWA,NNODE2 NBWA,AM,RHS	RHS	[K]{U}={RHS}の連立方程式を解き結果を{RHS}に入れる

BACK	NEXT