function tm = taumZLP(lamm, lamf, PROPS) %This function calculates the tension taum for the given lamdas and the %material constants cnontained in PROPS. It differs from taum because it %does not include INPUTS, only PROPS. b = PROPS{5}; LAMMAX = PROPS{2}/75; CBETA = PROPS{1}*75/LAMMAX/b; NMOLEC = PROPS{3}; calf = PROPS{4}; calf=(6*LAMMAX^2 - 9*LAMMAX + 4)/(LAMMAX - 1)^2; x0 = 1/LAMMAX; %The ZLP approximate consitutive model is coded here. % term1 = (1/4/(1-lamm*x0)^2 - 0.25 + lamm*x0)/lamm/x0; % rtterm = sqrt(lamm^2+3*lamf^2)*x0/2; % rtterm2 = sqrt(lamf^2+3*lamm^2)*x0/2; % term2 = (1/4/(1-rtterm)^2 - 0.25 + rtterm)/2/rtterm; % term3 = 3*(1/4/(1-rtterm2)^2 - 0.25 + rtterm2)/2/rtterm2; % tm = CBETA*(4*lamm/3/lamf*(term1+term2+term3) - calf/(lamm*lamf)^2); % THE INTERPOLATED (SIMPLIFIED VERSION IS CODED HERE: lama = sqrt(lamm*lamf); %d0 = -LAMMAX^2/(4*lama*(lama-LAMMAX)^2); d0 = -1/(4*lama*x0^2*(lama-1/x0)^2); d1n = 48*(lama*x0)^4 - 153*(lama*x0)^3 + 171*(lama*x0)^2 - 71*lama*x0 + 1; d1 = d1n/4/lama/(lama*x0 - 1)^3; dinf = 1/4/lama; lams = sqrt(lamm/lamf); %e0 = -lams^3/4/(lams-LAMMAX)^2; %einf = lams^3*x0^2/4; %mu0ovrcbeta = calf + (3*LAMMAX^2 - LAMMAX)/4/(LAMMAX - 1)^2; %e1 = 6*mu0ovrcbeta*x0 - lams^3*(lams + LAMMAX)/4/(LAMMAX - lams)^3 - 6*x0*calf; taus = 2*CBETA/3/x0*(d0 + d1*(lams - 1) + dinf*lams/(1 - lams*lama*x0)^2); %tau0 = 4*CBETA/3/x0*(3/4/(1 - x0)^2 - 3/4 + 3*x0); %taua1 = 2*CBETA/3/x0*(e0+e1*(lama - 1) + einf*lama/(1 - lama*lams*x0)^2); %taua = taua1 + tau0 - CBETA*calf/lamm^2/lamf^2; %tm = taua + taus; tm = taus;