36 lines (30 with data), 1.6 kB

AnyFolder FractureForces ={
  
  AnyForceMomentMeasure RigidJointReactionMeasure =  {
    AnyForceBase &ref1 = Main.RatModel.Jnts.Schnittlager.Constraints.Reaction; //this is a reference to the constraints of the joint
    AnyRefFrame &ref2 = Main.RatModel.Seg.Femur_p.Cut_Node;
  };
  
  AnyVector Fcut_dist = - RigidJointReactionMeasure.F;
  AnyVector Fcut_dist_local = - RigidJointReactionMeasure.F * Main.RatModel.Seg.Femur_p.Axes * (-1); 
  AnyVector Fcut_prox = RigidJointReactionMeasure.F;
  AnyVector Fcut_prox_local = RigidJointReactionMeasure.F * Main.RatModel.Seg.Femur_p.Axes * (-1); 
  AnyVar Fcut_total = vnorm(Fcut_dist_local);
  AnyVector Mcut_dist = - RigidJointReactionMeasure.M;
  AnyVector Mcut_dist_local = - RigidJointReactionMeasure.M * Main.RatModel.Seg.Femur_p.Axes * (-1); 
  AnyVector Mcut_prox = RigidJointReactionMeasure.M;
  AnyVector Mcut_prox_local = RigidJointReactionMeasure.M * Main.RatModel.Seg.Femur_p.Axes * (-1); 
  AnyVar Mcut_total = vnorm(Mcut_prox);
  
  AnyForceMomentMeasure2 masseffect_d = {
    AnyRefFrame &r = Main.RatModel.Seg.Femur_d.Cut_Node;
    AnySeg &s = Main.RatModel.Seg.Femur_d;
    IncludeGravity = On;
    IncludeInertiaForces = On;
  };
  AnyForceMomentMeasure2 masseffect_p = {
    AnyRefFrame &r = Main.RatModel.Seg.Femur_d.Cut_Node;
    AnySeg &s = Main.RatModel.Seg.Femur_p;
    IncludeGravity = On;
    IncludeInertiaForces = On;
  };
  
  AnyFloat mass_effect_F_sum = masseffect_d.F + masseffect_p.F;
  AnyFloat mass_effect_M_sum = masseffect_d.M + masseffect_p.M;
  
};