Code used to fit cyclic data

Critical times (tB-tDD (called tD' in paper)), maximum strain (emax) and rate of strain (ν)

tB=1/2; (*time of peak of first cycle*)
tD=1; (*time of end of first cycle*)
tBB=3/2; (*time of peak of second cycle - called tB' in paper*)
tDD=2; (*time of end of second cycle - called tD' in paper*) 
emax=6/100; (*peak fascicle strain*)
ν=emax/tB; (*strain rate*)

(*imposed fascicle strain*)
ee1[t_]:=ν t
ee2[t_]:=ν tB-ν(t-tB)
ee3[t_]:=ν(t-tD)
ee4[t_]:=ν(tBB-tD)-ν(t-tBB)
e[t_]:=If[t<0,0,If[t<tB,ee1[t],If[t<tD,ee2[t],If[t<tBB,ee3[t],If[t<tDD,ee4[t],0]]]]]

Critical time at which fibril becomes taut

tA[ec1_]:=ec1/ν

Relaxation function and its derivative

EE[t_]:=E0+E1 Exp[-t/τ1]+E2 Exp[-t/τ2]

DE[t_]:=Evaluate[D[EE[t],t]]

Creep function and its derivative

J[t_]:=J0+J1 Exp[-t/τc1]+J2 Exp[-t/τc2]

DJ[t_]:=Evaluate[D[J[t],t]]

Expressions for σ1 (these and all following expressions were determined by analytically evaluating relevant integrals in Mathematica)

σ11[t_,ec1_]:=((-E0 ec1+E0 t ν+E1 ν (τ1-(E^((ec1-t ν)/(ν τ1))) τ1)+E2 ν (τ2-(E^((ec1-t ν)/(ν τ2))) τ2))/(1+ec1))

σ12[t_,ec1_]:=((-E0 ec1+E0 tB ν+(E^(-(t/(τ1)))) ((E^(tB/τ1))-(E^(ec1/(ν τ1)))) E1 ν τ1+(E^(-(t/(τ2)))) ((E^(tB/τ2))-(E^(ec1/(ν τ2)))) E2 ν τ2)/(1+ec1))

σ13[t_,ec1_]:=((ν (E0 (-t+tB)+(-1+(E^((-t+tB)/(τ1)))) E1 τ1+(-1+(E^((-t+tB)/(τ2)))) E2 τ2))/(1+ec1))

Find tC

σ1[t_,ec1_]:=If[t<tA[ec1],0,If[t<tB,σ11[t,ec1],σ12[t,ec1]+σ13[t,ec1]]]

tC[ec1_]:=tC/.FindRoot[σ12[tC,ec1]+σ13[tC,ec1]==0,{tC,tB}]

Expressions for e1

e11[t_,ec1_]:=(1/(1+ec1))((E^((-t+tB)/(τc1))) J1 (ν (E1 τ1 (-1+(((E^((ec1-tB ν)/(ν τ1))) τ1)/(τ1-τc1)))+E2 τ2 (-1+(((E^((ec1-tB ν)/(ν τ2))) τ2)/(τ2-τc1))))+E0 (ec1+ν (-tB+τc1)))-(E^((-t+(ec1/ν))/(τc1))) J1 (ν (E1 τ1 (-1+((τ1)/(τ1-τc1)))+E2 τ2 (-1+((τ2)/(τ2-τc1))))+E0 (ec1+ν (-(ec1/ν)+τc1)))+(E^((-t+tB)/(τc2))) J2 (ν (E1 τ1 (-1+(((E^((ec1-tB ν)/(ν τ1))) τ1)/(τ1-τc2)))+E2 τ2 (-1+(((E^((ec1-tB ν)/(ν τ2))) τ2)/(τ2-τc2))))+E0 (ec1+ν (-tB+τc2)))-(E^((-t+(ec1/ν))/(τc2))) J2 (ν (E1 τ1 (-1+((τ1)/(τ1-τc2)))+E2 τ2 (-1+((τ2)/(τ2-τc2))))+E0 (ec1+ν (-(ec1/ν)+τc2))))

e12[t_,ec1_,tc_]:=-(1/(1+ec1))(-(E^((-t+tB)/(τc1))) J1 (-E0 ec1+E0 tB ν-E1 ν τ1-E2 ν τ2+(((E^(-(tB/(τ1)))) (2 (E^(tB/τ1))-(E^(ec1/(ν τ1)))) E1 ν ((τ1)^2))/(τ1-τc1))+(((E^(-(tB/(τ2)))) (2 (E^(tB/τ2))-(E^(ec1/(ν τ2)))) E2 ν ((τ2)^2))/(τ2-τc1))+E0 ν τc1)+(E^((-t+tc)/(τc1))) J1 (-E0 ec1+2 E0 tB ν-E0 tc ν-E1 ν τ1-E2 ν τ2+(((E^(-(tc/(τ1)))) (2 (E^(tB/τ1))-(E^(ec1/(ν τ1)))) E1 ν ((τ1)^2))/(τ1-τc1))+(((E^(-(tc/(τ2)))) (2 (E^(tB/τ2))-(E^(ec1/(ν τ2)))) E2 ν ((τ2)^2))/(τ2-τc1))+E0 ν τc1)-(E^((-t+tB)/(τc2))) J2 (-E0 ec1+E0 tB ν-E1 ν τ1-E2 ν τ2+(((E^(-(tB/(τ1)))) (2 (E^(tB/τ1))-(E^(ec1/(ν τ1)))) E1 ν ((τ1)^2))/(τ1-τc2))+(((E^(-(tB/(τ2)))) (2 (E^(tB/τ2))-(E^(ec1/(ν τ2)))) E2 ν ((τ2)^2))/(τ2-τc2))+E0 ν τc2)+(E^((-t+tc)/(τc2))) J2 (-E0 ec1+2 E0 tB ν-E0 tc ν-E1 ν τ1-E2 ν τ2+(((E^(-(tc/(τ1)))) (2 (E^(tB/τ1))-(E^(ec1/(ν τ1)))) E1 ν ((τ1)^2))/(τ1-τc2))+(((E^(-(tc/(τ2)))) (2 (E^(tB/τ2))-(E^(ec1/(ν τ2)))) E2 ν ((τ2)^2))/(τ2-τc2))+E0 ν τc2))

Find tAA (called tA' in paper)

e1[t_,ec1_]:=If[t<tA[ec1],0,If[t<tB,1/(1+ec1)(ee1[t]-ec1),If[t<tC[ec1],1/(1+ec1)(ee2[t]-ec1),e11[t,ec1]+e12[t,ec1,tC[ec1]]]]]

tAA[ec1_]:=tAA/.FindRoot[(1+ec1)(e11[tAA,ec1]+e12[tAA,ec1,tC[ec1]])+ec1-ee3[tAA]==0,{tAA,tD}]

Find tCC (called tC' in paper)

σ2[t_,ec1_]:=If[t<tA[ec1],0,If[t<tB,σ11[t,ec1],If[t<tC[ec1],σ12[t,ec1]+σ13[t,ec1],0]]]

e2[t_,ec1_]:=If[t<tA[ec1],0,If[t<tB,1/(1+ec1)(ee1[t]-ec1),If[t<tC[ec1],1/(1+ec1)(ee2[t]-ec1),If[t<tAA[ec1],e11[t,ec1]+e12[t,ec1,tC[ec1]],If[t<tBB,1/(1+ec1)(ee3[t]-ec1),1/(1+ec1)(ee4[t]-ec1)]]]]]

Expressions for σ3

σ31[t_,ec1_]:=σ12[t,ec1]

σ32[t_,ec1_,tc_]:=((ν (E0 tB-E0 tc+(E^(-(t/(τ1)))) ((E^(tB/τ1))-(E^(tc/τ1))) E1 τ1+(E^(-(t/(τ2)))) ((E^(tB/τ2))-(E^(tc/τ2))) E2 τ2))/(1+ec1))

σ33[t_,ec1_,tc_,taa_]:=(1/(1+ec1))(1/((((τ1-τc1))^2) (((τ2-τc1))^2))(E^(-(((t+tc) (τ1+τ2))/(τ1 τ2))-(taa/(τc1)))) J1 ((E^((t/(τ1))+(taa/(τ2)))) E2 τ2 (τ1-τc1)+(E^((taa/(τ1))+(t/(τ2)))) E1 τ1 (τ2-τc1)+(E^(t ((1/(τ1))+(1/(τ2))))) E0 (τ1-τc1) (τ2-τc1)) (-2 (E^((tB/(τ2))+tc ((1/(τ1))+(1/(τc1))))) E2 ν ((τ2)^2) (τ1-τc1)+(E^((ec1/(ν τ2))+tc ((1/(τ1))+(1/(τc1))))) E2 ν ((τ2)^2) (τ1-τc1)-2 (E^((tB/(τ1))+tc ((1/(τ2))+(1/(τc1))))) E1 ν ((τ1)^2) (τ2-τc1)+(E^((ec1/(ν τ1))+tc ((1/(τ2))+(1/(τc1))))) E1 ν ((τ1)^2) (τ2-τc1)+(E^(tc ((1/(τ1))+(1/(τ2))+(1/(τc1))))) (ν (E1 τ1+E2 τ2)+E0 (ec1+ν (-2 tB+tc-τc1))) (τ1-τc1) (τ2-τc1)+2 (E^(tc ((1/(τ1))+(1/(τ2)))+(tB/(τc1)))) ν (E2 τ2 (τ1-τc1)+E1 τ1 (τ2-τc1)+E0 (τ1-τc1) (τ2-τc1)) τc1-(E^(tc ((1/(τ1))+(1/(τ2)))+(ec1/(ν τc1)))) ν (E2 τ2 (τ1-τc1)+E1 τ1 (τ2-τc1)+E0 (τ1-τc1) (τ2-τc1)) τc1)-1/((((τ1-τc1))^2) (((τ2-τc1))^2))(E^(-t ((1/(τ1))+(1/(τ2)))-tc ((1/(τ1))+(1/(τ2))+(1/(τc1))))) J1 ((E^((t/(τ1))+(tc/(τ2)))) E2 τ2 (τ1-τc1)+(E^((tc/(τ1))+(t/(τ2)))) E1 τ1 (τ2-τc1)+(E^(t ((1/(τ1))+(1/(τ2))))) E0 (τ1-τc1) (τ2-τc1)) (-2 (E^((tB/(τ2))+tc ((1/(τ1))+(1/(τc1))))) E2 ν ((τ2)^2) (τ1-τc1)+(E^((ec1/(ν τ2))+tc ((1/(τ1))+(1/(τc1))))) E2 ν ((τ2)^2) (τ1-τc1)-2 (E^((tB/(τ1))+tc ((1/(τ2))+(1/(τc1))))) E1 ν ((τ1)^2) (τ2-τc1)+(E^((ec1/(ν τ1))+tc ((1/(τ2))+(1/(τc1))))) E1 ν ((τ1)^2) (τ2-τc1)+(E^(tc ((1/(τ1))+(1/(τ2))+(1/(τc1))))) (ν (E1 τ1+E2 τ2)+E0 (ec1+ν (-2 tB+tc-τc1))) (τ1-τc1) (τ2-τc1)+2 (E^(tc ((1/(τ1))+(1/(τ2)))+(tB/(τc1)))) ν (E2 τ2 (τ1-τc1)+E1 τ1 (τ2-τc1)+E0 (τ1-τc1) (τ2-τc1)) τc1-(E^(tc ((1/(τ1))+(1/(τ2)))+(ec1/(ν τc1)))) ν (E2 τ2 (τ1-τc1)+E1 τ1 (τ2-τc1)+E0 (τ1-τc1) (τ2-τc1)) τc1)+1/((((τ1-τc2))^2) (((τ2-τc2))^2))(E^(-(((t+tc) (τ1+τ2))/(τ1 τ2))-(taa/(τc2)))) J2 ((E^((t/(τ1))+(taa/(τ2)))) E2 τ2 (τ1-τc2)+(E^((taa/(τ1))+(t/(τ2)))) E1 τ1 (τ2-τc2)+(E^(t ((1/(τ1))+(1/(τ2))))) E0 (τ1-τc2) (τ2-τc2)) (-2 (E^((tB/(τ2))+tc ((1/(τ1))+(1/(τc2))))) E2 ν ((τ2)^2) (τ1-τc2)+(E^((ec1/(ν τ2))+tc ((1/(τ1))+(1/(τc2))))) E2 ν ((τ2)^2) (τ1-τc2)-2 (E^((tB/(τ1))+tc ((1/(τ2))+(1/(τc2))))) E1 ν ((τ1)^2) (τ2-τc2)+(E^((ec1/(ν τ1))+tc ((1/(τ2))+(1/(τc2))))) E1 ν ((τ1)^2) (τ2-τc2)+(E^(tc ((1/(τ1))+(1/(τ2))+(1/(τc2))))) (ν (E1 τ1+E2 τ2)+E0 (ec1+ν (-2 tB+tc-τc2))) (τ1-τc2) (τ2-τc2)+2 (E^(tc ((1/(τ1))+(1/(τ2)))+(tB/(τc2)))) ν (E2 τ2 (τ1-τc2)+E1 τ1 (τ2-τc2)+E0 (τ1-τc2) (τ2-τc2)) τc2-(E^(tc ((1/(τ1))+(1/(τ2)))+(ec1/(ν τc2)))) ν (E2 τ2 (τ1-τc2)+E1 τ1 (τ2-τc2)+E0 (τ1-τc2) (τ2-τc2)) τc2)-1/((((τ1-τc2))^2) (((τ2-τc2))^2))(E^(-t ((1/(τ1))+(1/(τ2)))-tc ((1/(τ1))+(1/(τ2))+(1/(τc2))))) J2 ((E^((t/(τ1))+(tc/(τ2)))) E2 τ2 (τ1-τc2)+(E^((tc/(τ1))+(t/(τ2)))) E1 τ1 (τ2-τc2)+(E^(t ((1/(τ1))+(1/(τ2))))) E0 (τ1-τc2) (τ2-τc2)) (-2 (E^((tB/(τ2))+tc ((1/(τ1))+(1/(τc2))))) E2 ν ((τ2)^2) (τ1-τc2)+(E^((ec1/(ν τ2))+tc ((1/(τ1))+(1/(τc2))))) E2 ν ((τ2)^2) (τ1-τc2)-2 (E^((tB/(τ1))+tc ((1/(τ2))+(1/(τc2))))) E1 ν ((τ1)^2) (τ2-τc2)+(E^((ec1/(ν τ1))+tc ((1/(τ2))+(1/(τc2))))) E1 ν ((τ1)^2) (τ2-τc2)+(E^(tc ((1/(τ1))+(1/(τ2))+(1/(τc2))))) (ν (E1 τ1+E2 τ2)+E0 (ec1+ν (-2 tB+tc-τc2))) (τ1-τc2) (τ2-τc2)+2 (E^(tc ((1/(τ1))+(1/(τ2)))+(tB/(τc2)))) ν (E2 τ2 (τ1-τc2)+E1 τ1 (τ2-τc2)+E0 (τ1-τc2) (τ2-τc2)) τc2-(E^(tc ((1/(τ1))+(1/(τ2)))+(ec1/(ν τc2)))) ν (E2 τ2 (τ1-τc2)+E1 τ1 (τ2-τc2)+E0 (τ1-τc2) (τ2-τc2)) τc2))

σ34[t_,ec1_,taa_]:=((ν (E0 (t-taa)+E1 (τ1-(E^((-t+taa)/(τ1))) τ1)+(1-(E^((-t+taa)/(τ2)))) E2 τ2))/(1+ec1))

σ35[t_,ec1_,taa_]:=((ν (-E0 taa+E0 tBB+(E^(-(t/(τ1)))) (-(E^(taa/τ1))+(E^(tBB/τ1))) E1 τ1+(E^(-(t/(τ2)))) (-(E^(taa/τ2))+(E^(tBB/τ2))) E2 τ2))/(1+ec1))

σ36[t_,ec1_]:=((ν (E0 (-t+tBB)+(-1+(E^((-t+tBB)/(τ1)))) E1 τ1+(-1+(E^((-t+tBB)/(τ2)))) E2 τ2))/(1+ec1))

σ3[t_,ec1_]:=If[t<tA[ec1],0,If[t<tB,σ11[t,ec1],If[t<tC[ec1],σ12[t,ec1]+σ13[t,ec1],If[t<tAA[ec1],0,If[t<tBB,σ31[t,ec1]+σ32[t,ec1,tC[ec1]]+σ33[t,ec1,tC[ec1],tAA[ec1]]+σ34[t,ec1,tAA[ec1]],σ31[t,ec1]+σ32[t,ec1,tC[ec1]]+σ33[t,ec1,tC[ec1],tAA[ec1]]+σ35[t,ec1,tAA[ec1]]+σ36[t,ec1]]]]]]

tCC[ec1_]:=tCC/.FindRoot[σ31[tCC,ec1]+σ32[tCC,ec1,tC[ec1]]+σ33[tCC,ec1,tC[ec1],tAA[ec1]]+σ35[tCC,ec1,tAA[ec1]]+σ36[tCC,ec1]==0,{tCC,tBB}]

Fibril stress

σf[t_,ec1_]:=If[t<tA[ec1],0,If[t<tB,σ11[t,ec1],If[t<tC[ec1],σ12[t,ec1]+σ13[t,ec1],If[t<tAA[ec1],0,If[t<tBB,σ31[t,ec1]+σ32[t,ec1,tC[ec1]]+σ33[t,ec1,tC[ec1],tAA[ec1]]+σ34[t,ec1,tAA[ec1]],If[t<tCC[ec1],σ31[t,ec1]+σ32[t,ec1,tC[ec1]]+σ33[t,ec1,tC[ec1],tAA[ec1]]+σ35[t,ec1,tAA[ec1]]+σ36[t,ec1],0]]]]]]

Total fascicle stress

NN[x_,μ_,σ_]:=PDF[TruncatedDistribution[{0,Infinity},NormalDistribution[μ,σ]],x](*Define critical strain distribution*)

σF[t_,μ_,σ_,ϕ_]:=ϕ NIntegrate[Hold[NN[ec1,μ,σ]σf[t,ec1]],{ec1,0,0.06},AccuracyGoal->8](*Note-the matrix term has been removed*)

Creep function parameters

J0=(1/E0);
J1=(E0 (E1-E2) (τ1-τ2)+(E1+E2) (E1 τ1+E2 τ2-(Sqrt[(E0^2) (((τ1-τ2))^2)+2 E0 (τ1-τ2) (E1 τ1-E2 τ2)+(((E1 τ1+E2 τ2))^2)])))/(2 E0 (E0+E1+E2) (Sqrt[-4 E0 (E0+E1+E2) τ1 τ2+(((E1 τ1+E2 τ2+E0 (τ1+τ2)))^2)]));
J2=(-E0 (E1-E2) (τ1-τ2)-(E1+E2) (E1 τ1+E2 τ2+(Sqrt[(E0^2) (((τ1-τ2))^2)+2 E0 (τ1-τ2) (E1 τ1-E2 τ2)+(((E1 τ1+E2 τ2))^2)])))/(2 E0 (E0+E1+E2) (Sqrt[-4 E0 (E0+E1+E2) τ1 τ2+(((E1 τ1+E2 τ2+E0 (τ1+τ2)))^2)]));
τc1=1/((1/(2 (E0+E1+E2) τ1 τ2)) (E1 τ1+E2 τ2+E0 (τ1+τ2)+(Sqrt[(E0^2) (((τ1-τ2))^2)+2 E0 (τ1-τ2) (E1 τ1-E2 τ2)+(((E1 τ1+E2 τ2))^2)])));
τc2=1/((1/(2 (E0+E1+E2) τ1 τ2))(E1 τ1+E2 τ2+E0 (τ1+τ2)-(Sqrt[(E0^2) (((τ1-τ2))^2)+2 E0 (τ1-τ2) (E1 τ1-E2 τ2)+(((E1 τ1+E2 τ2))^2)])));

Relaxation function parameters

E0=3000;(*MPa*)
E1=0.2 3000/0.63;(*MPa*)
E2=0.17 3000/0.63;(*MPa*)
τ1=1.9;(*s*)
τ2=52;(*s*)

Comparison with experimental data (First two load cycles only)

Import Cycle1

Cycle1=Import["INSERT FILE PATH HERE"];
Cycle1Time=Drop[Cycle1[[All,1]],4]-Cycle1[[5,1]];
Cycle1Strain=Drop[Cycle1[[All,10]],4];
Cycle1Stress=Drop[Cycle1[[All,11]],4];
Cycle1StressTime=Transpose[{Cycle1Time,Cycle1Stress}];
Cycle1StrainTime=Transpose[{Cycle1Time,Cycle1Strain}];

CyclesAveraged1=Join[Transpose[{MovingAverage[Join[Table[-0.07+0.01t,{t,0,6}],Cycle1Time[[1;;17]]],15],MovingAverage[Join[Table[0,{t,0,6}],Cycle1Stress[[1;;17]]],15]}],Transpose[{MovingAverage[Cycle1Time[[6;;25]],13],MovingAverage[Cycle1Stress[[6;;25]],13]}],Transpose[{MovingAverage[Cycle1Time[[17;;34]],11],MovingAverage[Cycle1Stress[[17;;34]],11]}],Transpose[{MovingAverage[Cycle1Time[[28;;43]],9],MovingAverage[Cycle1Stress[[28;;43]],9]}],Transpose[{MovingAverage[Cycle1Time[[40;;60]],1],MovingAverage[Cycle1Stress[[40;;60]],1]}],Transpose[{MovingAverage[Cycle1Time[[58;;73]],9],MovingAverage[Cycle1Stress[[58;;73]],9]}],Transpose[{MovingAverage[Cycle1Time[[67;;84]],11],MovingAverage[Cycle1Stress[[67;;84]],11]}],Transpose[{MovingAverage[Cycle1Time[[76;;96]],13],MovingAverage[Cycle1Stress[[76;;96]],13]}],Transpose[{MovingAverage[Cycle1Time[[84;;117]],15],MovingAverage[Cycle1Stress[[84;;117]],15]}],Transpose[{MovingAverage[Cycle1Time[[106;;125]],13],MovingAverage[Cycle1Stress[[106;;125]],13]}],Transpose[{MovingAverage[Cycle1Time[[117;;134]],11],MovingAverage[Cycle1Stress[[117;;134]],11]}],
Transpose[{MovingAverage[Cycle1Time[[128;;143]],9],MovingAverage[Cycle1Stress[[128;;143]],9]}],Transpose[{MovingAverage[Cycle1Time[[140;;150]],1],MovingAverage[Cycle1Stress[[140;;150]],1]}],Transpose[{MovingAverage[Cycle1Time[[158;;173]],9],MovingAverage[Cycle1Stress[[158;;173]],9]}],Transpose[{MovingAverage[Cycle1Time[[167;;184]],11],MovingAverage[Cycle1Stress[[167;;184]],11]}],Transpose[{MovingAverage[Cycle1Time[[176;;196]],13],MovingAverage[Cycle1Stress[[176;;196]],13]}],Transpose[{MovingAverage[Cycle1Time[[184;;217]],15],MovingAverage[Cycle1Stress[[184;;217]],15]}]];

Theoretical stress - critical strain distribution parameters were determined from fit to relaxation data

StressScaled=Quiet[Table[{t,σF[t,0,0.114,0.8]},{t,0,2,0.02}]];

CP1B=ListPlot[{CyclesAveraged1,StressScaled},Joined->True,PlotStyle->{{Red},{Black}},AxesLabel->{t,s},PlotRange->{{0,2},{-5,80}},PlotLabel->d]

Import Cycle2

Cycle2=Import["INSERT FILE PATH HERE"];
Cycle2Time=Drop[Cycle2[[All,1]],4]-Cycle2[[5,1]];
Cycle2Strain=Drop[Cycle2[[All,10]],4];
Cycle2Stress=Drop[Cycle2[[All,11]],4];
Cycle2StressTime=Transpose[{Cycle2Time,Cycle2Stress}];
Cycle2StrainTime=Transpose[{Cycle2Time,Cycle2Strain}];

CyclesAveraged2=Join[Transpose[{MovingAverage[Join[Table[-0.07+0.01t,{t,0,6}],Cycle2Time[[1;;17]]],15],MovingAverage[Join[Table[0,{t,0,6}],Cycle2Stress[[1;;17]]],15]}],Transpose[{MovingAverage[Cycle2Time[[6;;25]],13],MovingAverage[Cycle2Stress[[6;;25]],13]}],Transpose[{MovingAverage[Cycle2Time[[17;;34]],11],MovingAverage[Cycle2Stress[[17;;34]],11]}],Transpose[{MovingAverage[Cycle2Time[[28;;43]],9],MovingAverage[Cycle2Stress[[28;;43]],9]}],Transpose[{MovingAverage[Cycle2Time[[40;;60]],1],MovingAverage[Cycle2Stress[[40;;60]],1]}],Transpose[{MovingAverage[Cycle2Time[[58;;73]],9],MovingAverage[Cycle2Stress[[58;;73]],9]}],Transpose[{MovingAverage[Cycle2Time[[67;;84]],11],MovingAverage[Cycle2Stress[[67;;84]],11]}],Transpose[{MovingAverage[Cycle2Time[[76;;96]],13],MovingAverage[Cycle2Stress[[76;;96]],13]}],Transpose[{MovingAverage[Cycle2Time[[84;;117]],15],MovingAverage[Cycle2Stress[[84;;117]],15]}],Transpose[{MovingAverage[Cycle2Time[[106;;125]],13],MovingAverage[Cycle2Stress[[106;;125]],13]}],Transpose[{MovingAverage[Cycle2Time[[117;;134]],11],MovingAverage[Cycle2Stress[[117;;134]],11]}],
Transpose[{MovingAverage[Cycle2Time[[128;;143]],9],MovingAverage[Cycle2Stress[[128;;143]],9]}],Transpose[{MovingAverage[Cycle2Time[[140;;150]],1],MovingAverage[Cycle2Stress[[140;;150]],1]}],Transpose[{MovingAverage[Cycle2Time[[158;;173]],9],MovingAverage[Cycle2Stress[[158;;173]],9]}],Transpose[{MovingAverage[Cycle2Time[[167;;184]],11],MovingAverage[Cycle2Stress[[167;;184]],11]}],Transpose[{MovingAverage[Cycle2Time[[176;;196]],13],MovingAverage[Cycle2Stress[[176;;196]],13]}],Transpose[{MovingAverage[Cycle2Time[[184;;217]],15],MovingAverage[Cycle2Stress[[184;;217]],15]}]];

Theoretical stress

StressScaled2=Quiet[Table[{t,σF[t,0.039,0.016,0.8]},{t,0,2,0.02}]];

CP2B=ListPlot[{CyclesAveraged2,StressScaled2},Joined->True,PlotStyle->{{Green},{Black}},AxesLabel->{t,s},PlotRange->{{0,2},{-5,80}},PlotLabel->e]

Import Cycle 3

Cycle3=Import["INSERT FILE PATH HERE"];
Cycle3Time=Drop[Cycle3[[All,1]],4]-Cycle3[[5,1]];
Cycle3Strain=Drop[Cycle3[[All,10]],4];
Cycle3Stress=Drop[Cycle3[[All,11]],4];
Cycle3StressTime=Transpose[{Cycle3Time,Cycle3Stress}];
Cycle3StrainTime=Transpose[{Cycle3Time,Cycle3Strain}];

CyclesAveraged3=Join[Transpose[{MovingAverage[Join[Table[-0.07+0.01t,{t,0,6}],Cycle3Time[[1;;17]]],15],MovingAverage[Join[Table[0,{t,0,6}],Cycle3Stress[[1;;17]]],15]}],Transpose[{MovingAverage[Cycle3Time[[6;;25]],13],MovingAverage[Cycle3Stress[[6;;25]],13]}],Transpose[{MovingAverage[Cycle3Time[[17;;34]],11],MovingAverage[Cycle3Stress[[17;;34]],11]}],Transpose[{MovingAverage[Cycle3Time[[28;;43]],9],MovingAverage[Cycle3Stress[[28;;43]],9]}],Transpose[{MovingAverage[Cycle3Time[[40;;60]],1],MovingAverage[Cycle3Stress[[40;;60]],1]}],Transpose[{MovingAverage[Cycle3Time[[58;;73]],9],MovingAverage[Cycle3Stress[[58;;73]],9]}],Transpose[{MovingAverage[Cycle3Time[[67;;84]],11],MovingAverage[Cycle3Stress[[67;;84]],11]}],Transpose[{MovingAverage[Cycle3Time[[76;;96]],13],MovingAverage[Cycle3Stress[[76;;96]],13]}],Transpose[{MovingAverage[Cycle3Time[[84;;117]],15],MovingAverage[Cycle3Stress[[84;;117]],15]}],Transpose[{MovingAverage[Cycle3Time[[106;;125]],13],MovingAverage[Cycle3Stress[[106;;125]],13]}],Transpose[{MovingAverage[Cycle3Time[[117;;134]],11],MovingAverage[Cycle3Stress[[117;;134]],11]}],
Transpose[{MovingAverage[Cycle3Time[[128;;143]],9],MovingAverage[Cycle3Stress[[128;;143]],9]}],Transpose[{MovingAverage[Cycle3Time[[140;;150]],1],MovingAverage[Cycle3Stress[[140;;150]],1]}],Transpose[{MovingAverage[Cycle3Time[[158;;173]],9],MovingAverage[Cycle3Stress[[158;;173]],9]}],Transpose[{MovingAverage[Cycle3Time[[167;;184]],11],MovingAverage[Cycle3Stress[[167;;184]],11]}],Transpose[{MovingAverage[Cycle3Time[[176;;196]],13],MovingAverage[Cycle3Stress[[176;;196]],13]}],Transpose[{MovingAverage[Cycle3Time[[184;;217]],15],MovingAverage[Cycle3Stress[[184;;217]],15]}]];

Theoretical stress

StressScaled3=Quiet[Table[{t,σF[t,0.047,0.011,0.8]},{t,0,2,0.02}]];

CP3B=ListPlot[{CyclesAveraged3,StressScaled3},Joined->True,PlotStyle->{{Blue},{Black}},PlotRange->{{0,2},{-5,80}},AxesLabel->{t,s},PlotLabel->f]

GR2=GraphicsRow[{CP1B,CP2B,CP3B},ImageSize->Full]