MATLAB中计算等腰直角三角形单元且μ=14记杨氏模量为E厚度为t球形函数矩阵N应变矩阵B应力矩阵S单元刚度矩阵K

以下是MATLAB代码：

% 定义材料参数 E = 1e5; % 杨氏模量 t = 0.1; % 厚度 mu = 1/4; % 泊松比

% 定义节点坐标 x = [0,0;1,0;0,1];

% 计算单元面积 A = 0.5*det([1,1,1; x(1,:),1; x(2,:),1]);

% 计算球形函数矩阵[N] N = [1-x(1,1)-x(1,2), x(1,1), x(1,2); ... 1-x(2,1)-x(2,2), x(2,1), x(2,2); ... 1-x(3,1)-x(3,2), x(3,1), x(3,2)];

% 计算应变矩阵[B] B = 1/(2At) * [-x(2,2)+x(3,2), x(1,2)-x(3,2), -x(1,2)+x(2,2); ... -x(3,1)+x(2,1), -x(1,1)+x(3,1), x(1,1)-x(2,1)];

% 计算应力矩阵[S] S = E/(1-mu^2) * [1, mu, 0; mu, 1, 0; 0, 0, (1-mu)/2];

% 计算单元刚度矩阵[K] K = tAB'SB;