有限元求解柏松方程。本文采用ForTran語言編制程序。程序中大部分變量采用有名公共區(qū)存儲方式存儲,這樣可以減少內(nèi)存占用量。
IFG:生成有限元網(wǎng)格信息,即元素節(jié)點局部編碼與總體編碼對照表,節(jié)點實際坐標,邊界節(jié)點編碼與邊界點上的已知值
GKD:生成總剛一維存儲對角元的地址,計算總剛一維存儲長度
FIXP:設(shè)置已知節(jié)點函數(shù)值
GK(NI,NJ,ADJ,AIJ):單元剛度矩陣計算
GF(NI,N,M,LE,YI,FE):單元列陣的計算
AK(I,J,AIJ):總剛度矩陣元素迭加
QEB:總剛度矩陣和總列陣合成
BDE:邊界條件處理
SOLGS:Gauss-Seidel迭代法求解方程組
UDIFF(NI,NFLAG,UDIF,LE,ADJ):標準元素內(nèi)形狀函數(shù)導數(shù)計算
DIFF:節(jié)點上 , 加權(quán)平均
標簽:
ForTran
程序
有限元
方程
上傳時間:
2017-09-12
上傳用戶:erkuizhang
有限元求解柏松方程。本文采用ForTran語言編制程序。程序中大部分變量采用有名公共區(qū)存儲方式存儲,這樣可以減少內(nèi)存占用量。
IFG:生成有限元網(wǎng)格信息,即元素節(jié)點局部編碼與總體編碼對照表,節(jié)點實際坐標,邊界節(jié)點編碼與邊界點上的已知值
GKD:生成總剛一維存儲對角元的地址,計算總剛一維存儲長度
FIXP:設(shè)置已知節(jié)點函數(shù)值
GK(NI,NJ,ADJ,AIJ):單元剛度矩陣計算
GF(NI,N,M,LE,YI,FE):單元列陣的計算
AK(I,J,AIJ):總剛度矩陣元素迭加
QEB:總剛度矩陣和總列陣合成
BDE:邊界條件處理
SOLGS:Gauss-Seidel迭代法求解方程組
UDIFF(NI,NFLAG,UDIF,LE,ADJ):標準元素內(nèi)形狀函數(shù)導數(shù)計算
DIFF:節(jié)點上 , 加權(quán)平均
標簽:
ForTran
程序
有限元
方程
上傳時間:
2017-09-12
上傳用戶:問題問題