@@ -14,19 +14,19 @@ function plotmodel(modelvect::Vector{RheoModel}; ymaxG=nothing, ymaxJ=nothing)
1414 dσ = stressfunction (tstress, hstep ()) - stressfunction (tstress, hstep (offset= 5.0 ))
1515
1616 # Initialize Plots layout
17- plt = plot (layout= (3 ,1 ), size= (700 ,2000 ), left_margin = 15 Plots. mm)
17+ plt = plot (layout= (3 ,1 ), size= (700 ,1500 ), left_margin = 15 Plots. mm)
1818
1919 # --- Relaxation & Creep ---
2020 for i in 1 : length (modelvect)
2121 # Relaxation modulus
2222 dG = modelpredict (dϵ, modelvect[i])
2323 plot! (plt[1 ], dG. t, dG. σ, color= colplot[i], label= " model $i " , grid = true , linewidth = 3 , framestyle = :box )
24- plot! (plt[1 ], [- 5 * dt,0 ,dG. t[1 ]], [0.0 ,0.0 ,dG. σ[1 ]], color= colplot[i], label= " " )
24+ plot! (plt[1 ], [- 5 * dt,0 ,dG. t[1 ]], [0.0 ,0.0 ,dG. σ[1 ]], color= colplot[i], linewidth = 3 , label= " " )
2525
2626 # Creep modulus
2727 dJ = modelpredict (dσ, modelvect[i,1 ])
2828 plot! (plt[2 ], dJ. t, dJ. ϵ, color= colplot[i], label= " model $i " , grid = true , linewidth = 3 , framestyle = :box )
29- plot! (plt[2 ], [- 5 * dt,0 ,dJ. t[1 ]], [0.0 ,0.0 ,dJ. ϵ[1 ]], color= colplot[i], label= " " )
29+ plot! (plt[2 ], [- 5 * dt,0 ,dJ. t[1 ]], [0.0 ,0.0 ,dJ. ϵ[1 ]], color= colplot[i], linewidth = 3 , label= " " )
3030 end
3131
3232 # Set y-limits safely
@@ -60,7 +60,7 @@ function plotmodel(modelvect::Vector{RheoModel}; ymaxG=nothing, ymaxJ=nothing)
6060 xlabel! (plt[3 ], " Frequency" )
6161 ylabel! (plt[3 ], " Storage (—) and Loss (- -) moduli" )
6262 title! (plt[3 ], " Frequency response" )
63- xlims! (plt[3 ], 1e -2 , 1e2 )
63+ xlims! (plt[3 ], 0.9e -2 , 1.2e2 )
6464
6565 return plt
6666end
0 commit comments