%% Rajni Chahal (rajni.chahal@mavs.uta.edu), PhD Student, UT Arlington 
% Dr Ashfaq Adnan, Dr Ajit Roy
%% November 20, 2016
% Please contact for any help in code

% Developing Filled Fiber & their Interfaces
n=50; % number of fibers
L=400; % x-dimension of RVE
M=300; % y-dimension of RVE
l=L/10; % length of fibers
r1=2; %radius of CNTs
d1=2*r1; 
r2=4; %radius of CNTs+Interface
d2=2*r2; 
axis([0 L 0 M]);
ax=gca;
th=[];
XX=[]; YY=[];
for i=1:n
    x=rand*L;
    y=rand*M;
    t=rand*2*pi;
    th=[th;t];
    X=[x-(l/2)*cos(t) x+(l/2)*cos(t)];
    Y=[y-(l/2)*sin(t) y+(l/2)*sin(t)];
    figure(1);
    line(X,Y,'Color',[0 1 1],'LineWidth',d1);
    axis off
%     box on
    set(gca,'units','pixels','position',[0 0 L M]);
    %store random numbers for fiber regeneration
    XX=[XX; X]; YY=[YY; Y];
    figure(1); % identities in figure(1) will have radius = fiber + interfaces radius
    box on
    line(X,Y,'Color',[1 0 0],'LineWidth',d2);
    axis off
    set(gca,'units','pixels','position',[0 0 L M]);
    
    %generating fiber portion to be removed from interface figure
    figure(2); % identities in figure(2) will have radius = fiber radius
    line(X,Y,'Color',[1 0 0],'LineWidth',d1);
    axis off
    box on
    set(gca,'units','pixels','position',[0 0 L M]);
end

%% subtract figure(2)from figure(1) to generate interfaces
% B=imread('b.jpg');
% A=imread('a.jpg');
% Z=imfuse(B,A);
% F=imshow(Z)