compute_P

compute_P(X_at_nodes, U_at_nodes, S_at_int_points, element_type)[source]

Compute the first Piola-Kirchhoff stress tensors for num_elem elements of the given element type. This element type has n nodes and ips integration points. 2d elements have d=2 dimensions. 3d elements have d=3 dimensions.

See also

Module conforce_gen.element_definitions for a list of supported element types

Parameters:

  • X_at_nodes – Array of shape (num_elem, n, d) containing the coordinates at n nodes of num_elem elements.

  • U_at_nodes – Array of shape (num_elem, n, d) containing the displacements at n nodes of num_elem elements.

  • S_at_int_points – Array of shape (num_elem, ips, d, d) containing the symmetric stress tensors at ips integration points for num_elem elements.

  • element_type – str, name of an element type

Returns:

P_at_int_points: Array of shape (num_elem, ips, d, d) containing the 1. Piola-Kirchhoff stress tensors evaluated on ips integration points for num_elem element.