张拉整体梭型结构及其MATLAB构型程序

for i=1:nc
plot3(N(1,i),N(2,i),N(3,i),'o');
hold on
end
xlabel('x[ ]'); ylabel('y[ ]'); zlabel('z[ ]'); %---杆矩阵 B---
CB=[eye(n) -eye(n) zeros(n,2)];
CB=[CB;zeros(1,2*n) -1 1];
elseif i==n-1
ຫໍສະໝຸດ Baidu
e1i=[1 zeros(1,n-2) -1];
end
e1=[e1;e1i];
end
for i=0:n-1
if i==0
e2i=[zeros(1,n-1) 1];
elseif i>0
e2i=[zeros(1,i-1) 1 zeros(1,n-i)]; end e2=[e2;e2i]; end CS1=[-eye(n) zeros(n) ones(n,1) zeros(n,1)]; CS2=[zeros(n) -eye(n) zeros(n,1) ones(n,1)]; CS3=[-eye(n) e2 zeros(n,2)]; %---结构索矩阵--CS=[CS1;CS2;CS3]; CST=CS'; S=N*CST; [nsr nsc]=size(CST); for i=1:nsc; CSTi=CST(:,i); n2n=find(~(CSTi+ones(nsr,1))); n1=N(:,n2n); n2=N(:,n2n)+S(:,i); s_b=line([n1(1) n2(1)],[n1(2) n2(2)],[n1(3) n2(3)]); set([s_b],'LineWidth',1,'color','r'); hold on; end axis equal
end
for i=1:n
theta1=2*pi*i/n;
n2=[ru*cos(theta1+phi);ru*sin(theta1+phi);h0+h];
N=[N n2];
end
nd=[0;0;0];
%底部节点
nu=[0;0;2*h0+h]; %底部节点
N=[N nd nu];
[nr nc]=size(N);
张拉整体梭型结构及其 MATLAB 构型程序
30
25
20
z[ ]
y[ ]
15
10
5
0 4 2 0-2-4 y[ ]
-4 -2 0 2 4 x[ ]
4 3 2 1 0 -1 -2 -3 -4
-4 -3 -2 -1 0 1 2 3 4 5 x[ ]
n=5; h=10;
%---棱柱杆数--%---棱柱高度---
h0=10;
rd=5;
%---下底面半径---
ru=5;
%---上底面半径---
phi=30*pi/180; %---扭转角---
%---节点矩阵---
N=[];
for i=1:n
theta=2*pi*(i-1)/n;
n1=[rd*cos(theta);rd*sin(theta);h0];
N=[N n1];
set([g_b],'LineWidth',4,'color','b');
hold on;
end %---索矩阵 S---
e1=[];e2=[];
for i=0:n-1
if i==0
e1i=[-1 1 zeros(1,n-2)];
elseif i>0&i<n-1
e1i=[zeros(1,i) -1 1 zeros(1,n-i-2)];
CBT=CB';
B=N*CBT;
[nbr nbc]=size(CBT);
for i=1:nbc;
CBTi=CBT(:,i);
n1n=find(~(CBTi+ones(nbr,1)));
n1=N(:,n1n); n2=N(:,n1n)+B(:,i);
g_b=line([n1(1) n2(1)],[n1(2) n2(2)],[n1(3) n2(3)]);