做桁架静力学类和我继续这个错误在运行。三重检查我的工作我没有任何负整数
指数在位置1是无效的。数组
必须是正整数或指数
逻辑值。
错误在CEE210_CP3_verification(线
70)
PosVec (j) = NodeP (EN,:) -NodeP (SN:);
nNodesP = 6;
nMembersP = 9;
NodeP = 0 (nNodesP, 3);%初始化内存存储
MemberP = 0 (nMembersP, 2);
% ....输入节点坐标
% 9节点3节点4成员
:NodeP (1) = (0, 0, 0);
:NodeP (2) = (4 3 0);
:NodeP (3) = (8 3 0);
:NodeP (4) = (12, 0, 0);
:NodeP (5) = (4 0 0);
:NodeP (6) = (8, 0, 0);
成员(1:)= [1,2];% 1节点1到节点2
成员(2:)= (1、5);% 2节点2节点3
成员(3:)= (2、5);% 3节点2节点4成员
成员(4:)= (2、3);% 4节点2节点5成员
成员(5)= (2,6);% 5节点1到节点3
:成员(6日)= (6,5);% 6节点1到节点4成员
:成员(7日)= (3、6);% 7节点3节点6成员
:成员(8日)= (3、4);% 8节点3节点5成员
:成员(9日)= (4、6);% 9节点3节点4成员
PosVecP = 0 (nMembersP, 3);%初始化内存存储
UnitVecP = 0 (nMembersP, 3);
OppPosVecP = 0 (nMembersP, 3);
OppUnitVecP = 0 (nMembersP, 3);
Memberlength = 0 (nMembersP, 1);
(j = 1: nMembersP)
SN = MemberP (j, 1);% SN =开始节点
EN = MemberP (j, 2);
PosVecP (j) = NodeP (EN,:) -NodeP (SN:);
OppPosVecP (j) = -PosVecP (j:); % EN =结束节点
Memberlength (j) =规范(PosVecP (j,:));
UnitVecP (j) = PosVecP (j:)。/ Memberlength (j, 1);
OppUnitVecP (j) = -UnitVecP (j:);
%完成本节计算位置向量,长度、单位向量,
%和消极的单位向量桁架的每个成员。
结束
% %共创输出命令窗口
流(' % s \ n ', ' - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -”)
流(' % s \ n ', ' - - - - - - - - - - - -平面桁架- - - - - - - - - - - - - - - - - -”)
流(' % s \ n ', ' - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -”)
(i = 1: nNodesP)
流(' % s % 8我% 8.3 f % 8.3 f % 8.3 f \ n’,‘节点:’,我,NodeP(我:))
结束
流(“\ n”)
(i = 1: nMembersP)
流(' % s % 8我% 8.3 f % 8.3 f % 8.3 f \ n”,“单位向量:”,我,UnitVecP(我:))
结束
流(“\ n \ n \ n”)
% %共创的阴谋
图一=图(1);clf;网格;轴相等;抓住;
包含(“X”);ylabel (Y);标题(“平面桁架”);
m = 1: nMembersP
SN = MemberP (m, 1);
EN = MemberP (m, 2);
X = [NodeP (SN, 1);NodeP(嗯,1)];
Y = [NodeP (SN, 2);NodeP (EN, 2)];
Z = [NodeP (SN, 3);NodeP (EN), 3)];
p =情节(X, Y);
集(p,“颜色”,“k”,“线宽”,2);
结束% m, nMembersP
散射(NodeP (: 1) NodeP(:, 2),“填满”,“红”);
