fke_1.m

来自「在有限元分析中计算结构整体刚度矩阵、单元刚度矩阵、单元结点位移、应力应变等」· M 代码 · 共 14 行

function Ke=fKe_1(E,nu,t,ECC)
%三角单元刚度矩阵函数
%%%%%%%%%%%%%%%%参数说明
%弹性模量E单位MPa
%泊松比nu
%厚度t单位cm
%结点编号向量e单位cm
%%%%%%%%%%%%%%%%%%%%%%%控制程序
E=E*10^2;%弹性模量量纲转换为N/cm^2
%%%%%%%%%%%%%%%%%%%%%%%%产生结点坐标
Ae=fAe_1(ECC);
D=E/(1-nu^2)*[1 nu 0;nu 1 0;0 0 (1-nu)/2];
Be=fBe_1(ECC);
Ke=Be'*D*Be*t*Ae;