constrained nodal force output
Hi,
Is there any way to output constrained nodal force of "examples/linear_elasticity/material_nonlinearity.py"? I would like to confirm material nonlinearity with nodal force.
Any help or clarification would be greatly appreciated. Best regards,
Takuo Fujita
Hi Takuo,
On 09/22/2015 02:13 AM, pyontaku14 wrote:
Hi,
Is there any way to output constrained nodal force of "examples/linear_elasticity/material_nonlinearity.py"? I would like to confirm material nonlinearity with nodal force.
By nodal force you mean the residuals (outofbalance forces) in FE nodes, right? If so, those can be obtained by evaluating the residual vector with removed boundary conditions, so that the forces in fixed nodes are computed. Add the following into post_process():
def post_process(out, pb, state, extend=False):
...
ebcs = pb.ebcs
pb.time_update(ebcs={})
vec = pb.get_variables().create_state_vector()
pb.equations.evaluate(mode='weak', dw_mode='vector', asm_obj=vec)
out['res'] = Struct(name='res', mode='vertex', data=vec.reshape((1, 3)),
dofs=None)
pb.time_update(ebcs=ebcs)
...
and run the example. Then use
./postproc.py cylinder.h5 b step 1 onlynames=res
to see the nodal forces (in the last step).
Does it help?
r.
participants (2)

pyontaku14

Robert Cimrman