IndexOutOfRangeException in model.Solve()

[AmitR72]

Hi there,
Thanks for this lovely library!

I'm trying to apply the examples in the documentation to a model I draw in Rhino & GrassHopper with a plugin I write in c# using the library.
I get the following error on model.Solve(): IndexOutOfRangeException: 'Index was outside the bounds of the array.'
when I tried Solve_MPC, it didn't show an error, but all forces on the bars (elements) were 0. I tried to use the frame and different loads examples too.
At first, I got an error that the matrix isn't symmetric and not positive definite. And after changing the constraints following this #11, I got over this error, but now getting this one.

this is my model, where the long vector is ~10 and the shorter is ~5 in length. the other vectors are just marks of the constraints, but as I mentioned, I applied Fixed for all the nodes.

[epsi1on]

Hi
Thanks. Note that Model.Solve() is deprecated and will be removed. Model.Solve_MPC() should be used instead.
It seems you are trying to model a 2d truss, if yes, then maybe all nodes are totally fixed to ground ...

[epsi1on]

so it seems you've fixed DoFs more than enough so all loads are directly directed to ground and no element have internal force. Can you please tell me which nodes in your visualized model are connected to the ground?
Is it a 2D truss?
Here is more info on how to track the notposdef exception:
https://github.com/BriefFiniteElementNet/BriefFiniteElement.Net/wiki/How-to-fix-NotPosDef-error
I think it is a little more complicated than is seems.
Thanks

[AmitR72]

Thank you! It worked!
now, another question: if I use GetInternalForceAt(0) on each bar after Solve_MPC(), I get:
curve at (0, 5, 0) --- (5, 5, 0) with internal force: {F: -12.5, 0, 0, M: 0, 0, 0}
curve at (0, 5, 0) --- (10, 5, 0) with internal force: {F: 17.5, 0, 0, M: 0, 0, 0}
curve at (5, 5, 0) --- (7.5, 0, 0) with internal force: {F: -16.7705098312485, 0, 0, M: 0, 0, 0}
curve at (0, 5, 0) --- (2.5, 0, 0) with internal force: {F: -11.180339887499, 0, 0, M: 0, 0, 0}
curve at (5, 5, 0) --- (2.5, 0, 0) with internal force: {F: 11.180339887499, 0, 0, M: 0, 0, 0}
curve at (10, 5, 0) --- (7.5, 0, 0) with internal force: {F: -39.1311896062463, 0, 0, M: 0, 0, 0}
curve at (2.5, 0, 0) --- (7.5, 0, 0) with internal force: {F: -10, 0, 0, M: 0, 0, 0}

Why do I have only forces at x? Not what I would expcet. And, doest negative force means the bar is in strech, while positive force means it's in stress? I want to get this info for each bar, stress/strech.
Thanks again!

[epsi1on]

BarElement.GetInternalForceAt(0) wil give you the internal force at the middle of element.
For truss elements we have only axial force (in local coordination system, it is along element's X axis)
About the sign, as long as i remember, positive should meamn stretch and negative mean compress.

[epsi1on]

Also if the comment is solution to your problem, please use the button (mark as answer) to mark it as answer.
Thanks

[AmitR72]

Sorry for the many questions...
If I want to apply kind of a G force on a structure.. Should I just use the Dead loadCase with CoordinationSystem.Global for each bar? including the vertical bars? I want to make sure that the self-weight of the vertical bars are also being taking in account in the forces caculations, and I'm not sure that the uniformLoad suits this. (pretty similar to the 4th example here https://www.codeproject.com/Articles/794983/Finite-Element-Method-Programming-in-Csharp only without xy forces and z forces with self-mass of the verticals bars)
Thanks!

[epsi1on]

I think it would work, UniformLoad, with Vector.K direction and DeadLoad load case and CoordinationSystem.Global direction on every bar (even verticals). Anyways need to be tested.

There are also 7 other well documented examples on documentation. They could help too. link

Example list is in the left panel: