My_coding – CellularBeam_Mesh ©
function OutputFile2(Nd,E)
% flange nodes to be set as z=0
nF = [find(Nd(:,2)<1) ; find(Nd(:,2)>max(Nd(:,2))-1); ...
find(Nd(:,3)<1) ; find(Nd(:,3)>max(Nd(:,3))-1)];
NdF = Nd(nF,1);
% flange element connectivity (grp2) and Loading
nL = [find(Nd(:,3)<1) ; find(Nd(:,3)>max(Nd(:,3))-1)];
NdL = Nd(nL,1);
Flange = find(sum(ismember(E,Nd(nL,1)),2));
EF = E(Flange,:);
E(Flange,:) = [];
% NEW FILE..
fid = fopen('fullCB_v1.dat', 'wt');
fprintf(fid, '');
fclose(fid);
% START WRITING IN FILE..
pr('#',0)
pr('analysis 3d.full statics',0)
pr('#',0)
pr('materials',0)