compute_P_for_C3D15

compute_P_for_C3D15(X_at_nodes, U_at_nodes, S_at_int_points)[source]

Compute the first Piola-Kirchhoff stress tensors for num_elem elements of typ C3D15. Each element has n=15 nodes and ips=9 integration points.

Parameters:

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

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

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

Returns:

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