单元刚度矩阵源代码每一行含义

wm12345

function y = PlaneFrameElementStiffness(E,A,I,L,theta)
%PlaneFrameElementStiffness   This function returns the element
%                             stiffness matrix for a plane frame   
%                             element with modulus of elasticity E,  
%                             cross-sectional area A, moment of
%                             inertia I, length L, and angle
%                             theta (in degrees).
%                             The size of the element stiffness
%                             matrix is 6 x 6.
x = theta*pi/180;
C = cos(x);
S = sin(x);
w1 = A*C*C + 12*I*S*S/(L*L);
w2 = A*S*S + 12*I*C*C/(L*L);
w3 = (A-12*I/(L*L))*C*S;
w4 = 6*I*S/L;
w5 = 6*I*C/L;
y = E/L*[w1 w3 -w4 -w1 -w3 -w4 ; w3 w2 w5 -w3 -w2 w5 ;
   -w4 w5 4*I w4 -w5 2*I ; -w1 -w3 w4 w1 w3 w4 ;
   -w3 -w2 -w5 w3 w2 -w5 ; -w4 w5 2*I w4 -w5 4*I];

16443

你这没头没尾的鬼能知道啊

wm12345

回复 2 # 16443 的帖子

这个函数是计算弹性模量E，横截面积A，惯性矩I，长度为L的平面钢架元的单元刚度矩阵。它返回6X6的单元刚度矩阵k。可以直接用，但不知道它怎么写出来的，每一步的意思不懂

wm12345

回复 2 # 16443 的帖子

y = PlaneFrameElementStiffness(E,A,I,L,theta)

16443

matlab代码，先去学习一下吧。