On Sun, Apr 24, 2016 at 5:10 PM, Bhalla, Amneet Pal S <amne...@live.unc.edu> wrote:
> > > On Apr 22, 2016, at 6:24 PM, David Knezevic <david.kneze...@akselos.com> > wrote: > > That term is omitted in the code, but it should be easy to add in by > adding a surface integral term to the residual. Also, note that the > traction term \int_\Gamma g_i v_i ds will not contribute to the Jacobian > (similar to how the volumetric load doesn't contribute to the Jacobian). > > > Yes, I have added the term \int_\Gamma g_i v_i ds in the residual > evaluation. > OK, if that's all you want then it should be straightforward. You can follow the approach for assembling the traction term in systems_of_equations_ex6, for example. > > But I'm not sure if the formulation above is what you're asking about? > > I guess it won't contribute to Jacobian, as it should be treated > independent of the unknown solution field. > > P.S: I was able to convert that example into a 2d case for a known > displacement field (from which I constructed F, PK1 etc) and > verified that I am getting correct solution. The field is (u,v) = (XY, > Y^2), which works with a zero Dirichlet BC at Y = 0. To provide > a nonzero Dirichlet BC at some other location, I did something like: > > http://paste.ofcode.org/i2348da2BnaDcEZhkXFZa3 > > Does this make sense? > Yes, that makes sense, that is how you attach a non-zero DirichletBoundary. David ------------------------------------------------------------------------------ Find and fix application performance issues faster with Applications Manager Applications Manager provides deep performance insights into multiple tiers of your business applications. It resolves application problems quickly and reduces your MTTR. Get your free trial! https://ad.doubleclick.net/ddm/clk/302982198;130105516;z _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
