Bug fixes for ustruct Guccione model and follower pressure load

Code/Source/svFSI/mat_models.cpp

@@ -690,10 +690,10 @@ void get_pk2cc_dev(const ComMod& com_mod, const CepMod& cep_mod, const dmnType&
  2.0*g3*(Es(0,1)*RmRm.slice(3) + Es(0,2)*RmRm.slice(5));
 
  auto CCb = 2.0*ten_dyad_prod(Sb, Sb, nsd);
- Sb += Sb * r2;
+ Sb = Sb * r2;
 
  // Fiber reinforcement/active stress
- Sb += Sb + Tfa*mat_dyad_prod(fl.col(0), fl.col(0), nsd);
+ Sb += Tfa*mat_dyad_prod(fl.col(0), fl.col(0), nsd);
 
  double r1 = J2d*mat_ddot(C, Sb, nsd) / nd;
  S = J2d*Sb - r1*Ci;

Code/Source/svFSI/ustruct.cpp

@@ -142,7 +142,7 @@ void b_ustruct_3d(const ComMod& com_mod, const int eNoN, const double w, const V
  for (int a = 0; a < eNoN; a++) {
  NxFi(0,a) = Nx(0,a)*Fi(0,0) + Nx(1,a)*Fi(1,0) + Nx(2,a)*Fi(2,0);
  NxFi(1,a) = Nx(0,a)*Fi(0,1) + Nx(1,a)*Fi(1,1) + Nx(2,a)*Fi(2,1);
- NxFi(2,a) = Nx(1,a)*Fi(1,2) + Nx(1,a)*Fi(1,2) + Nx(2,a)*Fi(2,2);
+ NxFi(2,a) = Nx(0,a)*Fi(0,2) + Nx(1,a)*Fi(1,2) + Nx(2,a)*Fi(2,2);
  }
 
  nFi(0) = nV(0)*Fi(0,0) + nV(1)*Fi(1,0) + nV(2)*Fi(2,0);