<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Aug 25, 2016 at 8:46 AM, Karin&NiKo <span dir="ltr"><<a href="mailto:niko.karin@gmail.com" target="_blank">niko.karin@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear PETSc gurus,</div><div><br></div><div>Thanks to the help of Matthew, I have been able to reproduce in PETSc some tests of the paper of Peter Brune et al. entitled "Composing scalable nonlinear algebraic solvers", with special attention to the elasticity test.</div><div>I have also tried to reproduce it in a widely used mechanics finite element solver and I cannot obtain the same results, mainly because of my lack of undestanding of the boundary conditions and of the loading. </div><blockquote style="margin-right:0px" dir="ltr"><div>According to the paper, I undestand that the edges in red in the attached image are fully clamped. If I do that, I do observe a rotation of the grey  face (see attached image). </div></blockquote><blockquote style="margin-right:0px" dir="ltr"><div>If I clamp the over-mentioned edges plus I forbid the normal displacement of the grey face, I get a quite similar warped shape but the details of the deformation near the clamped faces are very different</div></blockquote><blockquote style="margin-right:0px" dir="ltr"><div>In the paper, the loading is defined as a volume force applied to the whole structure whereas in Wriggers' book, it is defined as a nodal force. </div><div><br></div></blockquote><div dir="ltr">Could you please give me some details on these points?</div></div></blockquote><div><br></div><div>I am looking at the code:</div><div><br></div><div><a href="https://bitbucket.org/petsc/petsc/src/bd8b9ddea6dc39522a77cdc0585cbede38eb66bb/src/snes/examples/tutorials/ex16.c?at=master&fileviewer=file-view-default#ex16.c-204">https://bitbucket.org/petsc/petsc/src/bd8b9ddea6dc39522a77cdc0585cbede38eb66bb/src/snes/examples/tutorials/ex16.c?at=master&fileviewer=file-view-default#ex16.c-204</a><br></div><div><br></div><div>and to me it clearly looks like we are fixing just the top edge on each side. I do not see anything forcing the end faces to</div><div>stay aligned with the y-z plane, but maybe I am missing something. Can you just compare the residual/Jacobian for the</div><div>initial guess?</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 dir="ltr"><div>Best regards,</div><div>Nicolas</div><div class="gmail_extra"><br></div><div class="gmail_quote">2016-08-23 20:25 GMT+02:00 Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left:1px solid rgb(204,204,204)"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span>On Tue, Aug 23, 2016 at 10:54 AM, Karin&NiKo <span dir="ltr"><<a href="mailto:niko.karin@gmail.com" target="_blank">niko.karin@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left:1px solid rgb(204,204,204)"><div dir="ltr"><div><div><div><div><div>Dear PETSc team,<br><br></div>I have read with high interest the paper of Peter Brune et al. entitled "Composing scalable nonlinear algebraic solvers".<br></div>Nevertheless I would like to be able to reproduce the tests that are presented within (mainly the elasticity problem, ex16). <br><br></div>Could you please provide us with the command lines of these tests?<br></div></div></div></blockquote><div><br></div></span><div>I believe Peter used the attached script.</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;padding-left:1ex;border-left:1px solid rgb(204,204,204)"><div dir="ltr"><div><div></div>Best regards,<br></div>Nicolas<span><font color="#888888"><br></font></span></div><span><font color="#888888">
</font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><span class=""><font color="#888888"><div><br></div>-- <br><div data-smartmail="gmail_signature">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>
</font></span></font></span></div></div>
</blockquote></div><div class="gmail_extra"><br></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">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></div>