quadraticbarelementstiffness.m

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

function y = QuadraticBarElementStiffness(E,A,L)
%QuadraticBarElementStiffness   This function returns the element 
%                               stiffness matrix for a quadratic bar  
%                               with modulus of elasticity E,  
%                               cross-sectional area A, and length L.   
%                               The size of the element stiffness 
%                               matrix is 3 x 3.
y = E*A/(3*L)*[7 1 -8 ; 1 7 -8 ; -8 -8 16];