Hi,
I’m using AMMRV1.3.1 and the FreePostureMove model.
I’m trying to apply a time varying force to a dummy segment, which also moves in time but is rigidly connected to the thigh segment. I have included the code for the dummy segment and the force.
My problem is that the force doesn’t seem to be applied to the thigh segment. In the ChartFX I can see the force as it should be. But… I’ve only applied the force to the right thigh, but the muscle activities etc. are identical between left and right. I even magnified the force by a factor of 100 to try to highlight the difference.
Any help would be greatly appreciated!
Regards,
Chris
AnyFolder RThighCoP = {
AnySeg Seg= //Adds a dummy segment of zero mass etc
{
//AnyRefFrame &BodySeg = Main.HumanModel.BodyModel.Right.Leg.Seg.Thigh;
Mass = 0;
Jii = {0, 0, 0};
AnyDrawNode DRW=
{
Visible = On;
Opacity = 1;
RGB = {1, 0, 1};
ScaleXYZ = {0.025, 0.025, 0.025};
};
}; // *DummyThighSeg
AnyKinEq RotFix=
{
AnyKinRotational ROT=
{
Type = RotAxesAngles;
AnyRefFrame &Seg= ..Seg;
};
};// *RotFix
AnyKinEqInterPolDriver CoPDriver= {
Type = Bspline;
BsplineOrder = 4;
FileName = Main.MatlabInputParameters.SegForceFilePath +"\CoP_Thigh.txt";
AnyRefFrame &SegRef = .Seg;
AnyKinLinear LIN={
Ref = 0;//-1;
AnyRefFrame &SegRef = ..Seg;
AnyRefFrame &BodySeg = Main.HumanModel.BodyModel.Right.Leg.Seg.Thigh;
};
Reaction.Type={Off,Off,Off}; //Implies it has no influence on/from the environment
}; // *CoPDriver
AnyReacForce Reac= {
AnyKinLinear LIN= //Creates "infinitely" strong link between the segment and the dummy segment
{
Ref = 0;
AnyRefFrame &BodySeg = Main.HumanModel.BodyModel.Right.Leg.Seg.Thigh;
AnyRefFrame &Seg = ..Seg;
};
}; // *Reac
AnyFunInterpol ForceFun={
Type = Bspline;
BsplineOrder = 4;
FileName = Main.MatlabInputParameters.SegForceFilePath +"\Force_Thigh.txt"; //Read in Force on Seg
}; // *ForceFun
AnyForce3D Force= { //Apply the Force to the Seg
F = .ForceFun(t)*100;
//F={1000,0,0};
AnyRefFrame &Seg = .Seg;
}; // *Force
AnyDrawVector drwVector1 = {
AnyRefFrame &Seg= .Seg;
Visible = On;
Opacity = 1;
Vec = .Force.F *0.0001;
GlobalCoord = Off;
//Vec = {1, 1, 1};
Line = {
Thickness = 0.01;
RGB = {0.5, 1, 1};
End =
{
Style = Line3DCapStyleArrow;
RGB = {1, 1, 1};
Thickness = 0.025;
Length = 0.025;
};// *End
};// *Line
}; // *drwVector1
}; /* RThighCoP */