Triangle2D3Node有限元基础教程曾攀算例4.3.2

主命令函数：

%%-----主程序
clc
clear all;
 %---物理参数------------------
E=1;              %弹性模量,单位Pa
t=1;              %单元厚度,单位m
% a=1;                 %尺寸a，单位m
% q=10e6;              %均布力，单位N/m2
NU=0.25;             %泊松比
%-------------------------------
x1=0;
y1=0;
x2=1;
y2=0;
x3=1;
y3=1;
x4=0;
y4=1;
Triangle2D3Node=Triangle2D3Node_Func;
k1=Triangle2D3Node.Stiffness(E,NU,t,x1,y1,x2,y2,x4,y4,1)
k2=Triangle2D3Node.Stiffness(E,NU,t,x2,y2,x3,y3,x4,y4,1)
% 两单元刚度矩阵相同
KK=zeros(8);%（8*8）初始KK==0提供位置
KK=Triangle2D3Node.Assembly(KK,k1,1,2,4); %杆1
KK=Triangle2D3Node.Assembly(KK,k2,2,3,4)  %加杆2
% KK=Triangle2D3Node.Assembly(KK,k1,3,4,2)  %或者用这个节点顺序与k1一样
k=KK([3 4 5 6 8],[3 4 5 6 8])
p=[-1; 0; 1; 0 ;0]  %5*1  力边界条件
u=k\p
q=[0,0,u(1),u(2),u(3),u(4),0,u(5)]'
P=KK*q  %节点力
fp=[0 0 p(1) p(2) p(3) p(4) 0 p(5)]'
F=P-f