单元刚度矩阵源代码每一行含义
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];
你这没头没尾的鬼能知道啊 回复 2 # 16443 的帖子
这个函数是计算弹性模量E，横截面积A，惯性矩I，长度为L的平面钢架元的单元刚度矩阵。它返回6X6的单元刚度矩阵k。可以直接用，但不知道它怎么写出来的，每一步的意思不懂 回复 2 # 16443 的帖子
y = PlaneFrameElementStiffness(E,A,I,L,theta) matlab代码，先去学习一下吧。
页:
[1]