Comparing strain rates between initial and final states with plotmodel

swellsmo

Hi,

I'm trying to plot principal strain rates from the initial model velocities and compare them to the solved stress balance principal strain rates. Ideally I would like the plot command to look something like this: plotmodel(md,'data',strainrate_initial,'data',strainrate_final)

I think this means I need to have both the initial and final strain rate calculations stored in md. I am using the mechanicalproperties function, but running it on the initial velocities overwrites the results from the stress balance results velocities. My next thought was to modify md.stressbalance.requested_outputs to get the final strain rates in md.results.StressbalanceSolution, but I'm struggling to find a list of possible requested outputs, and I'm not sure if requested_outputs will compute principal strains specifically, which mechanicalproperties does. Is there a requested_output I'm not finding for principal strain rates, or is there a way to have two mechanicalproperties runs stored in md.results?

Thanks for any help!

mathieumorlighem

Hello @swellsmo

md.results is a struct so you can rename anything you want. For example, if you want to stick with mechanicalproperties, you can do this:

md=mechanicalproperties(md, vx1, vy1);
md.results.strainrate1 = md.results.strainrate;
 md=mechanicalproperties(md, vx2, vy2);
md.results.strainrate2 = md.results.strainrate;
plotmodel(md, 'data', md.results.strainrate1.effectivevalue, 'data', md.results.strainrate2.effectivevalue);

You can do the same with md.results.SB1 = md.results.StressbalanceSolution, etc. You can request strain rates like this md.stressbalance.requested_outputs = {'default', 'StrainRatexx', 'StrainRateyy', 'StrainRatexy'}; but you will need to diagonalize the matrix yourself (although we could also add this as a requested output at some point if you would prefer).
I hope this answers your questions!
Best
Mathieu

swellsmo

Thank you!