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

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