syms EI l a EA = 1/2*EI; k = [EA/EI 0 0 -EA/EI 0 0;0 12/l^2 -6 0 -12/l^2 -6/l;0 -6/l 4 0 6/l 2;-EA/EI 0 0 EA/EI 0 0;0 -12/l^2 6/l 0 12/l^2 6/l;0 -6/l 2 0 6/l 4]; k = EI/l*k; Q = [cos(a) sin(a) 0 0 0 0;-sin(a) cos(a) 0 0 0 0;0 0 1 0 0 0;0 0 0 cos(a) sin(a) 0;0 0 0 -sin(a) cos(a) 0;0 0 0 0 0 1]; l_2 = 5; format short alpha_2 = atan(-3/4); alpha_2 = alpha_2*180/pi; k_2_lokal = subs(k,l,l_2); Q_2 = subs(Q,a,alpha_2); k_2 = transpose(Q_2).* k_2_lokal.*Q_2 此段代码的原理为直接刚度法。出现的k_2结果非常复杂,cos这串数字的结果用计算器计算为4/5,不知道该如何简化。 使用过simplfy(),double()和format都无法简化。 Matlab版本为64位R2020b |
最佳答案
1 条回复