%-constraints--------------------------------------------%
%-----Node, Dof, Value
% axial displacement and torsion are restricted
for k=1:size([Link],1)
if k==1
[Link]=[[Link](k,1), 1, 0
[Link](k,1), 4, 0
];
else
[Link]=[[Link]
[Link](k,1), 1, 0
[Link](k,1), 4, 0
];
end