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 (out-of-balance 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 --only-names=res

to see the nodal forces (in the last step).

Does it help?

r.