資源簡介
四面體網格劃分的有限元計算源代碼,在將變形體劃分為四面體后,輸入節點坐標,程序利用彈塑性本構,將載荷和位移關系轉變為線性方程組,是一個很好的求解簡單三維變形體的程序,還可輸出單元變形量,應力等。

代碼片段和文件信息
function?y?=?TetrahedronAssemble(Kkijmn)
%TetrahedronAssemble???This?function?assembles?the?element?stiffness
%??????????????????????matrix?k?of?the?linear?tetrahedral?(solid)
%??????????????????????element?with?nodes?i?j?m?and?n?into?the
%??????????????????????global?stiffness?matrix?K.
%??????????????????????This?function?returns?the?global?stiffness??
%??????????????????????matrix?K?after?the?element?stiffness?matrix??
%??????????????????????k?is?assembled.
K(3*i-23*i-2)?=?K(3*i-23*i-2)?+?k(11);
K(3*i-23*i-1)?=?K(3*i-23*i-1)?+?k(12);
K(3*i-23*i)?=?K(3*i-23*i)?+?k(13);
K(3*i-23*j-2)?=?K(3*i-23*j-2)?+?k(14);
K(3*i-23*j-1)?=?K(3*i-23*j-1)?+?k(15);
K(3*i-23*j)?=?K(3*i-23*j)?+?k(16);
K(3*i-23*m-2)?=?K(3*i-23*m-2)?+?k(17);
K(3*i-23*m-1)?=?K(3*i-23*m-1)?+?k(18);
K(3*i-23*m)?=?K(3*i-23*m)?+?k(19);
K(3*i-23*n-2)?=?K(3*i-23*n-2)?+?k(110);
K(3*i-23*n-1)?=?K(3*i-23*n-1)?+?k(111);
K(3*i-23*n)?=?K(3*i-23*n)?+?k(112);
K(3*i-13*i-2)?=?K(3*i-13*i-2)?+?k(21);
K(3*i-13*i-1)?=?K(3*i-13*i-1)?+?k(22);
K(3*i-13*i)?=?K(3*i-13*i)?+?k(23);
K(3*i-13*j-2)?=?K(3*i-13*j-2)?+?k(24);
K(3*i-13*j-1)?=?K(3*i-13*j-1)?+?k(25);
K(3*i-13*j)?=?K(3*i-13*j)?+?k(26);
K(3*i-13*m-2)?=?K(3*i-13*m-2)?+?k(27);
K(3*i-13*m-1)?=?K(3*i-13*m-1)?+?k(28);
K(3*i-13*m)?=?K(3*i-13*m)?+?k(29);
K(3*i-13*n-2)?=?K(3*i-13*n-2)?+?k(210);
K(3*i-13*n-1)?=?K(3*i-13*n-1)?+?k(211);
K(3*i-13*n)?=?K(3*i-13*n)?+?k(212);
K(3*i3*i-2)?=?K(3*i3*i-2)?+?k(31);
K(3*i3*i-1)?=?K(3*i3*i-1)?+?k(32);
K(3*i3*i)?=?K(3*i3*i)?+?k(33);
K(3*i3*j-2)?=?K(3*i3*j-2)?+?k(34);
K(3*i3*j-1)?=?K(3*i3*j-1)?+?k(35);
K(3*i3*j)?=?K(3*i3*j)?+?k(36);
K(3*i3*m-2)?=?K(3*i3*m-2)?+?k(37);
K(3*i3*m-1)?=?K(3*i3*m-1)?+?k(38);
K(3*i3*m)?=?K(3*i3*m)?+?k(39);
K(3*i3*n-2)?=?K(3*i3*n-2)?+?k(310);
K(3*i3*n-1)?=?K(3*i3*n-1)?+?k(311);
K(3*i3*n)?=?K(3*i3*n)?+?k(312);
K(3*j-23*i-2)?=?K(3*j-23*i-2)?+?k(41);
K(3*j-23*i-1)?=?K(3*j-23*i-1)?+?k(42);
K(3*j-23*i)?=?K(3*j-23*i)?+?k(43);
K(3*j-23*j-2)?=?K(3*j-23*j-2)?+?k(44);
K(3*j-23*j-1)?=?K(3*j-23*j-1)?+?k(45);
K(3*j-23*j)?=?K(3*j-23*j)?+?k(46);
K(3*j-23*m-2)?=?K(3*j-23*m-2)?+?k(47);
K(3*j-23*m-1)?=?K(3*j-23*m-1)?+?k(48);
K(3*j-23*m)?=?K(3*j-23*m)?+?k(49);
K(3*j-23*n-2)?=?K(3*j-23*n-2)?+?k(410);
K(3*j-23*n-1)?=?K(3*j-23*n-1)?+?k(411);
K(3*j-23*n)?=?K(3*j-23*n)?+?k(412);
K(3*j-13*i-2)?=?K(3*j-13*i-2)?+?k(51);
K(3*j-13*i-1)?=?K(3*j-13*i-1)?+?k(52);
K(3*j-13*i)?=?K(3*j-13*i)?+?k(53);
K(3*j-13*j-2)?=?K(3*j-13*j-2)?+?k(54);
K(3*j-13*j-1)?=?K(3*j-13*j-1)?+?k(55);
K(3*j-13*j)?=?K(3*j-13*j)?+?k(56);
K(3*j-13*m-2)?=?K(3*j-13*m-2)?+?k(57);
K(3*j-13*m-1)?=?K(3*j-13*m-1)?+?k(58);
K(3*j-13*m)?=?K(3*j-13*m)?+?k(59);
K(3*j-13*n-2)?=?K(3*j-13*n-2)?+?k(510);
K(3*j-13*n-1)?=?K(3*j-13*n-1)?+?k(511);
K(3*j-13*n)?=?K(3*j-13*n)?+?k(512);
K(3*j3*i-2)?=?K(3*j3*i-2)?+?k(61);
K(3*j3*i-1)?
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????631??2002-01-25?04:49??cj\TetrahedronElementVolume.m
?????文件???????6375??2002-01-25?06:43??cj\TetrahedronAssemble.m
?????文件????????615??2002-01-25?05:24??cj\TetrahedronElementPStresses.m
?????文件???????2364??2002-01-25?05:11??cj\TetrahedronElementStiffness.m
?????文件???????2404??2002-01-25?05:10??cj\TetrahedronElementStresses.m
?????目錄??????????0??2008-03-31?20:41??cj
-----------?---------??----------?-----??----
????????????????12607????????????????????7
- 上一篇:matlab讀取tdms文件
- 下一篇:風力發電機matlab仿真模型
評論
共有 條評論