%% Function name
%       VfVmrhoc

%% Revised
%       3 February 2014

%% Authors
%       Joshua Bigham, Trey Moore, & Autar Kaw
%       Section: All
%       Semester: Fall 2013

%% Purpose
%       Given densities of the fiber and matrix, as well as the weight volume
%       fraction, output the volume fractions for the fiber and matrix of the lamina

%% Usage
%       function [Vf,Vm,rhoc] = VfVmrhoc(rhof,rhom,Wf)
% Input Variables
%       rhof=fiber density
%       rhom=matrix density
%       Wf=weight fraction
% Output Variables
%       Vf=fiber volume fraction
%       Vm=matrix volume fraction
%       rhoc=density of lamina
% Keyword
%       volume fraction
%       density

%% License Agreement
%       http://www.eng.usf.edu/~kaw/OCW/composites/license/limiteduse.pdf

%% Code
function [Vf,Vm,rhoc] = VfVmrhoc(rhof,rhom,Wf)
% Density of lamina
rhoc=1/((Wf/rhof)+((1-Wf)/rhom));
% Volume fractions
Vf=Wf*(rhoc/rhof);
Vm=(1-Wf)*(rhoc/rhom);
end