四面体（MATLAB）

四面体网格划分的有限元计算源代码，在将变形体划分为四面体后，输入节点坐标，程序利用弹塑性本构，将载荷和位移关系转变为线性方程组，是一个很好的求解简单三维变形体的程序，还可输出单元变形量，应力等。

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