On Apr 22, 2016, at 6:24 PM, David Knezevic <[email protected]<mailto:[email protected]>> 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. 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? ------------------------------------------------------------------------------ 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 [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
