<div dir="ltr"><div dir="ltr">On Mon, Dec 19, 2022 at 6:12 AM TARDIEU Nicolas via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg2247495761693931989">




<div dir="ltr">
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">Dear PETSc users,</span></div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)"><br>
</span></div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">I plan to solve nonlinear quasi-static problems with PETSc. More
 precisely, these are solid mechanics problems with elasto-plasticity. </span></div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">So they do not involve  "physical time", rather  "pseudo time",
 which is mandatory to describe the stepping of the loading application.</span></div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">In general, the loading vector F(x, t) is expressed
 as the following product  F(x, t)=F0(x)*g(t), where g is a scalar function of the pseudo-time.</span></div>
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I see how to use a SNES in order to solve a certain step of the loading history but I wonder if a TS can be used to deal with the loading history through the definition of this g(t) function ? </div></div></div></div></blockquote><div><br></div><div>I believe so. We would want you to formulate it as a differential equation,</div><div><br></div><div>  F(x, \dot x, t) = G(x)</div><div><br></div><div>which is your case is easy</div><div><br></div><div>  F(x, t) = 0</div><div><br></div><div>so you would just put your function completely into the IFunction. Since there is no \dot x term, this is a DAE,</div><div>and you need to use one of the solvers for that.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg2247495761693931989"><div dir="ltr"><div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Furthermore, since too large load steps can lead to non-convergence, a stepping strategy is almost always required to restart a load step that failed. Does TS offer such a feature ?</div></div></div></div></blockquote><div><br></div><div>I think you can use PETSc adaptation. There is a PID option, which might be able to capture your behavior.</div><div><br></div><div>  Thanks,</div><div><br></div><div>      Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg2247495761693931989"><div dir="ltr"><div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Thank you for your answers.</div>
<div id="m_3468023608163333585Signature">
<div name="divtagdefaultwrapper">
<font size="3"></font>
<div style="margin:0px"><font size="3" face="Calibri,sans-serif">Regards, </font>
</div>
<font size="3"></font>
<div style="margin:0px"><font face="Calibri,sans-serif"><font size="3">Nicolas<br>
<span style="font-size:11pt">--</span></font><br>
</font><font size="2" face="Arial,sans-serif" color="#ff8100"><span style="font-size:9pt"><b>Nicolas Tardieu</b></span></font><br>
<div style="margin:0px"><font size="2" face="Calibri,sans-serif"><span style="font-size:11pt"><font size="1" face="Arial,sans-serif" color="#ff8100"><span style="font-size:8pt"><b>Ing PhD Computational Mechanics</b></span></font><font size="1" face="Arial,sans-serif" color="#ff8100"><span style="font-size:8pt"><b><br>
</b></span></font></span></font></div>
<div style="margin:0px"><font size="2" face="Calibri,sans-serif"><span style="font-size:11pt"><font size="1" face="Arial,sans-serif" color="#0062e1"><span lang="en-US" style="font-size:8pt">EDF - R&D Dpt ERMES</span></font><font size="1" face="Helv,sans-serif" color="#0062e1"><span lang="en-US" style="font-size:2pt"><br>
</span></font></span></font></div>
<div style="margin:0px"><span style="font-family:Arial,sans-serif;font-size:8pt;color:rgb(0,98,225);background-color:rgba(0,0,0,0)">PARIS-SACLAY, FRANCE</span></div>
</div>
</div>
</div>
</div>
<p></p>

<p><br>
Ce message et toutes les pièces jointes (ci-après le 'Message') sont établis à l'intention exclusive des destinataires et les informations qui y figurent sont strictement confidentielles. Toute utilisation de ce Message non conforme à sa destination, toute diffusion ou toute publication totale ou partielle, est interdite sauf autorisation expresse.</p>

<p>Si vous n'êtes pas le destinataire de ce Message, il vous est interdit de le copier, de le faire suivre, de le divulguer ou d'en utiliser tout ou partie. Si vous avez reçu ce Message par erreur, merci de le supprimer de votre système, ainsi que toutes ses copies, et de n'en garder aucune trace sur quelque support que ce soit. Nous vous remercions également d'en avertir immédiatement l'expéditeur par retour du message.</p>

<p>Il est impossible de garantir que les communications par messagerie électronique arrivent en temps utile, sont sécurisées ou dénuées de toute erreur ou virus.<br>
____________________________________________________</p>

<p>This message and any attachments (the 'Message') are intended solely for the addressees. The information contained in this Message is confidential. Any use of information contained in this Message not in accord with its purpose, any dissemination or disclosure, either whole or partial, is prohibited except formal approval.</p>

<p>If you are not the addressee, you may not copy, forward, disclose or use any part of it. If you have received this message in error, please delete it and all copies from your system and notify the sender immediately by return message.</p>

<p>E-mail communication cannot be guaranteed to be timely secure, error or virus-free.</p></div>

</div></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>