<div dir="ltr"><div>Hello</div><div><br></div><div>Thank you for the response. To give you a better idea of what I am doing, allow me to describe the problem:</div><div><br></div><div>A cantilever plate measuring 4cm x 4cm with thickness 0.1cm in fixed into the wall along its left edge. We want to analyze the plane stress situation (i.e, in the 4cmx4cm plane). The plate is made of steel and is loaded at the top right corner by a force say P Newton. <br></div><div><br></div><div>Hence, we are interested in Txx, Tyy, Txy stresses and corresponding strains and displacements.<br><div><img src="cid:ii_jhycm1z30_163c3fbdc2352600" width="389" height="481"></div><div><br></div><div>I have looked at the function 'g' that you have mentioned; it seems that this function assigns to the local variable 'y', no matter the coordinates. However, I am interested in imposing Neumann conditions only at the end of the plate. (Apart from the "Tip loading")<br></div><div><br></div><div>(The fig. is from "Finite Element Procedures By K.J Bathe)</div><div><br></div><div>Thanks</div><div><br></div><div>Nageswararao E.<br></div>​<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jun 3, 2018 at 2:52 AM, Christian Engwer <span dir="ltr"><<a href="mailto:christian.engwer@uni-muenster.de" target="_blank">christian.engwer@uni-muenster.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Nagesh,<br>
<br>
the current implementation does not include inhomogeneous Neumann<br>
boundary conditions. They should be fairly easy to implement, but<br>
nobody needed them up to now and thus nobody implemented them.<br>
<br>
Please have a look at linearelasticity.hh:310. This is where you might<br>
query a function for the surface force. The you have to compute the<br>
resulting residual contribution into `y`. The function `g` is supposed<br>
to return the surface force. The first question you have to answer<br>
your self is, which type of surface forces do you want to model here,<br>
i.e. which diemensionality do they have?<br>
<br>
Best<br>
Christian Engwer<br>
<div><div class="h5"><br>
On Sat, Jun 02, 2018 at 07:59:48PM +0530, Nagesh Eranki wrote:<br>
> Hello<br>
> <br>
> I am a mechanical engineering undergraduate trying to learn the ropes of<br>
> Finite Element Analysis. To facilitate and practice the theory, I am<br>
> currently going through some test cases in:<br>
> <br>
> /dune-2.6.0/dune-pdelab/dune/<wbr>pdelab/test<br>
> <br>
> In particular, I am interested in the file "testelasticity.cc". This file<br>
> contains "#include <dune/pdelab/localoperator/<wbr>linearelasticity.hh>" which I<br>
> have tracked down and studied. That in turn requires "#include<br>
> "linearelasticityparameter.hh"<wbr>, also which   I have studied.<br>
> <br>
> I have noticed that there is a function "isDirichlet" that returns the type<br>
> bool.<br>
> <br>
> I am interested to impose non-homogenous Neumann conditions, in addition to<br>
> adding a concentrated "tip force" at the end of the cantilever.<br>
> <br>
> Please advise on how to implement this.<br>
> <br>
> Thanks<br>
> <br>
> Nageswararao E.<br>
<br>
</div></div>> ______________________________<wbr>_________________<br>
> dune-pdelab mailing list<br>
> <a href="mailto:dune-pdelab@lists.dune-project.org">dune-pdelab@lists.dune-<wbr>project.org</a><br>
> <a href="https://lists.dune-project.org/mailman/listinfo/dune-pdelab" rel="noreferrer" target="_blank">https://lists.dune-project.<wbr>org/mailman/listinfo/dune-<wbr>pdelab</a><br>
</blockquote></div><br></div>