%總剛度矩陣
function?ASTIF=ASSEMBLE（NPOINNELEMYOUNGPOISSTHICKCOORDLNODS）
%?張成特定大小總剛矩陣，并賦初值?0
%?ASTIF（1:2*NPOIN1:2*NPOIN）=0;
%?%?組裝總剛度矩陣
%?for?i=1:NELEM
%?????%?計算單剛
%?????[ESTIFSMATX]=FORMESTIF（iCOORDLNODSYOUNGPOISSTHICK）;
%?????nt=LNODS（i:）;?%?臨時向量?nt用來記錄當前單元的結點號
%?????for?j=1:3
%?????????for?k=1:3
%?????????????ASTIF（（nt（j）*2-1）:nt（j）*2（nt（k）*2-1）:nt（k）*2）=...
%?????????????ASTIF（（nt（j）*2-1）:nt（j）*2（nt（k）*2-1）:nt（k）*2）+...
%?????????????ESTIF（j*2-1:j*2k*2-1:k*2）;
%?????????????%?利用結點位移編號對應關系將單元剛度分塊疊加到總剛矩陣
%?????????end
%?????end
%?end
ASTIF（1:NPOIN1:NPOIN）=0;
%?組裝總剛度矩陣
for?i=1:NELEM
????%?計算單剛
????[ESTIFSMATX]=FORMESTIF（iCOORDLNODS）;
????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

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

?????文件???????1051??2017-06-27?17:39??FEM\ASSEMBLE.m

?????文件????????492??2017-06-26?19:36??FEM\BOUNDARY.m

?????文件????????697??2017-06-27?17:51??FEM\FORMESTIF.m

?????文件????????264??2017-06-26?16:11??FEM\FORMLOAD.m

?????文件????????785??2017-06-27?10:08??FEM\gaus.m

?????文件????????690??2017-06-26?20:17??FEM\Ke.txt

?????文件???????1727??2017-07-16?16:10??FEM\main1.m

?????文件????????135??2017-06-27?17:36??FEM\tri_in.txt

?????文件??????21924??2017-06-29?19:16??FEM\tri_in1.txt

?????文件????????426??2017-06-26?16:13??FEM\WRITESTRESS.m

?????目錄??????????0??2017-07-16?15:57??FEM

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

????????????????28191????????????????????11