Skip to content

Instantly share code, notes, and snippets.

@swordlql
swordlql / odbPrinStressAngle.py
Created April 15, 2017 12:24 — forked from crmccreary/odbPrinStressAngle.py
Python script for use within Abaqus to create a fieldOutput of the angle between a specified vector and the normal to the maximum principal stress plane
import sys, getopt, os, string
import math
from odbAccess import *
from abaqusConstants import *
import numpy as np
from numpy import linalg as LA
def get_p1_vector(s):
'''
11, 22, 33, 12, 13, 23