Add new file; main dgibi code to calculate the principal stresses in transient method

Written main dgibi code to calculate the principal stresses in transient method. First version, with stresses approaching zero at the end and peaking in the middle.