<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 25, 2016 at 9: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:0 0 0 .8ex;border-left:1px #ccc solid;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></div></blockquote><div>This does not sound stable.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><blockquote style="margin-right:0px" dir="ltr"><div> </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></div></blockquote><div>Your picture looks like this problem. </div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><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 dir="ltr"><br></div><div>Best regards,</div><div>Nicolas</div><div><div class="h5"><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-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid"><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-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid"><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-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid"><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"><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></div></div>
</blockquote></div><div class="gmail_extra"><br></div></div></div></div>
</blockquote></div><br></div></div>