beamk.m

来自「matlab编程实现有限元计算,对有限元学习很有帮助」· M 代码 · 共 8 行

function y=beamK(E,I,L)
syms x
A=[1 0 0 0;0 1 0 0;1 L L*L L*L*L;0 1 2*L 3*L*L];
H=[0 0 2 6*x];
B=-H*inv(A);
Btr=transpose(B);
K=E*I*int(Btr*B,x,0,L);
y=double(K);