%總剛度矩陣
function?ASTIF=ASSEMBLE（NPOINNELEMNJIEZHICOORDLNODS）

ASTIF（1:NPOIN1:NPOIN）=0;%生成特定大小的總剛度矩陣，并賦值為0.總剛度矩陣，為節(jié)點數(shù)*節(jié)點數(shù)，代表了每一個節(jié)點與另一個節(jié)點間的關(guān)系矩陣
%?組裝總剛度矩陣
for?i=1:NELEM
????%?計算單剛
????[ESTIFSMATX]=FORMESTIF（iNJIEZHICOORDLNODS）;
????nt=LNODS（i:）;?%?臨時向量?nt用來記錄當(dāng)前單元的結(jié)點號，三個節(jié)點
????for?j=1:3
????????for?k=1:3???%兩個節(jié)點互相交錯行列所對應(yīng)的節(jié)點間的關(guān)系。
????????????ASTIF（nt（j）nt（k））=...
????????????ASTIF（nt（j）nt（k））+...%疊加剛度矩陣
????????????ESTIF（jk）;
????????????%?利用結(jié)點位移編號對應(yīng)關(guān)系將單元剛度分塊疊加到總剛矩陣
????????end
????end
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????660??2018-05-09?15:14??fem-eit\ASSEMBLE.m

?????文件???????8016??2018-05-09?15:09??fem-eit\data.txt

?????文件???????1149??2018-05-09?15:13??fem-eit\FORMESTIF.m

?????文件????????924??2018-05-09?15:09??fem-eit\gaus.m

?????文件???????1494??2018-05-09?15:12??fem-eit\main1.m

?????目錄??????????0??2018-05-09?15:17??fem-eit

-----------?---------??----------?-----??----

????????????????12243????????????????????6