[SOLVED] error with TetrahedralTensorMassForceField

[shang]

dear all, did some have ever used TetrahedralTensorMassForceField before? I simulated my model by this force field , there does not show any tips but
“################ SIG11 – SIGSEGV: segfault ##########”
in the console, the xml is list bellow, does any know about this component and help me?

<Node name="root" dt="0.005" gravity="0 -9.81 0" >
    <VisualStyle displayFlags="hideBehaviorModels hideCollisionModels hideMappings hideForceFields" />
    <RequiredPlugin name="Geomagic plugin" pluginName="Geomagic" />
    
    <CollisionPipeline name="pipeline" depth="6" verbose="0"/>
    <BruteForceDetection name="detection" />
	<NewProximityIntersection name="Proximity" alarmDistance="0.3" contactDistance="0.2" />
    <CollisionResponse name="response" response="FrictionContact" />
    <CollisionGroup name="Group" />

    <FreeMotionAnimationLoop/>
    <GenericConstraintSolver tolerance="0.001" maxIt="1000"/>

    <Node name="Liver" gravity="0 -9.81 0">
		<EulerImplicitSolver tolerance="0.001" maxIt="1000"/>
        <CGLinearSolver name="linear solver" iterations="25" tolerance="1e-09" threshold="1e-09" />
        <MeshGmshLoader name="meshLoader" filename="mesh/artery-coarse-big2-shorten3.msh" />

        <Mesh src="@meshLoader"/>
	    <MechanicalObject name="dofs" src="@meshLoader" />
	    <TetrahedronSetTopologyContainer  name="Container" src="@meshLoader"/>
  	    <TetrahedronSetTopologyModifier   name="Modifier" />
        <TetrahedronSetTopologyAlgorithms name="TopoAlgo"   template="Vec3d" />
        <TetrahedronSetGeometryAlgorithms name="GeomAlgo"   template="Vec3d" />
        <DiagonalMass  name="computed using mass density" massDensity="0.001" />
        <TetrahedralTensorMassForceField name="TM" youngModulus="1000" poissonRatio="0.2" />
		
        <PrecomputedConstraintCorrection recompute="true"/>

        <FixedConstraint  name="FixedConstraint" indices="3 39 64" />

        <Node name="Visu" >
            <OglModel  name="VisualModel" fileMesh="mesh/artery-smooth-big2-shorten3.stl" />
		
            <BarycentricMapping name="visual mapping" input="@../dofs" output="@VisualModel" />
        </Node>
        <Node name="Surf2">
				<MeshSTLLoader name="loader" filename="mesh/artery-coarse-big2-shorten3.stl" />
                <Mesh src="@loader" />
                <MechanicalObject src="@loader" dx="0.5" />
                <Triangle />
                <BarycentricMapping />
        </Node>
    </Node>
</Node>

[shang]

I find there exists some bug in component ‘TetrahedralTensorMassForceField’ I simulated successfully with the tetrahedral mesh “mesh/liver.msh” but the exe crushed with only changing the gmsh file with “mesh/liver2.msh” ,so there exists with some problems with this component, does any know about this?

[shang]

exactly, this episode is right below

<Node name="root" dt="0.005" gravity="0 -9.81 0" >
    <VisualStyle displayFlags="hideBehaviorModels hideCollisionModels hideMappings hideForceFields" />
    <RequiredPlugin name="Geomagic plugin" pluginName="Geomagic" />
    
    <CollisionPipeline name="pipeline" depth="6" verbose="0"/>
    <BruteForceDetection name="detection" />
	<NewProximityIntersection name="Proximity" alarmDistance="0.3" contactDistance="0.2" />
    <CollisionResponse name="response" response="FrictionContact" />
    <CollisionGroup name="Group" />
    <FreeMotionAnimationLoop/>
    <GenericConstraintSolver tolerance="0.001" maxIt="1000"/>

    <Node name="Liver" gravity="0 -9.81 0">
        <EulerImplicitSolver name="cg_odesolver"  />
        <CGLinearSolver name="linear solver" iterations="25" tolerance="1e-09" threshold="1e-09" />
        <MeshGmshLoader name="meshLoader" filename="mesh/liver.msh" />
        <TetrahedronSetTopologyContainer name="topo" src="@meshLoader" />
        <MechanicalObject name="dofs" src="@meshLoader" />
        <TetrahedronSetGeometryAlgorithms template="Vec3d" name="GeomAlgo" />
        <DiagonalMass  name="computed using mass density" massDensity="1" />
        <TetrahedralTensorMassForceField name="TM" youngModulus="6000" poissonRatio="0.3" />
        <PrecomputedConstraintCorrection recompute="true"/>

        <FixedConstraint  name="FixedConstraint" indices="3 39 64" />
        <Node name="Visu" tags="Visual" gravity="0 -9.81 0">
            <OglModel  name="VisualModel" fileMesh="mesh/liver-smooth.obj" />
            <BarycentricMapping name="visual mapping" input="@../dofs" output="@VisualModel" />
        </Node>
        <Node name="Surf" gravity="0 -9.81 0">
            <SphereLoader filename="mesh/liver.sph" />
            <MechanicalObject name="spheres" position="@[-1].position" />
            <TSphereModel name="CollisionModel" listRadius="@[-2].listRadius"/>
            <BarycentricMapping name="sphere mapping" input="@../dofs" output="@spheres" />
        </Node>
    </Node>
</Node>

and this episode caused crush with “################ SIG11 – SIGSEGV: segfault ##########”

<Node name="root" dt="0.005" gravity="0 -9.81 0" >
    <VisualStyle displayFlags="hideBehaviorModels hideCollisionModels hideMappings hideForceFields" />
    <RequiredPlugin name="Geomagic plugin" pluginName="Geomagic" />
    
    <CollisionPipeline name="pipeline" depth="6" verbose="0"/>
    <BruteForceDetection name="detection" />
	<NewProximityIntersection name="Proximity" alarmDistance="0.3" contactDistance="0.2" />
    <CollisionResponse name="response" response="FrictionContact" />
    <CollisionGroup name="Group" />
    <FreeMotionAnimationLoop/>
    <GenericConstraintSolver tolerance="0.001" maxIt="1000"/>

    <Node name="Liver" gravity="0 -9.81 0">
        <EulerImplicitSolver name="cg_odesolver"  />
        <CGLinearSolver name="linear solver" iterations="25" tolerance="1e-09" threshold="1e-09" />
        <MeshGmshLoader name="meshLoader" filename="mesh/liver2.msh" />
        <TetrahedronSetTopologyContainer name="topo" src="@meshLoader" />
        <MechanicalObject name="dofs" src="@meshLoader" />
        <TetrahedronSetGeometryAlgorithms template="Vec3d" name="GeomAlgo" />
        <DiagonalMass  name="computed using mass density" massDensity="1" />
        <TetrahedralTensorMassForceField name="TM" youngModulus="6000" poissonRatio="0.3" />
        <PrecomputedConstraintCorrection recompute="true"/>

        <FixedConstraint  name="FixedConstraint" indices="3 39 64" />
        <Node name="Visu" tags="Visual" gravity="0 -9.81 0">
            <OglModel  name="VisualModel" fileMesh="mesh/liver-smooth.obj" />
            <BarycentricMapping name="visual mapping" input="@../dofs" output="@VisualModel" />
        </Node>
        <Node name="Surf" gravity="0 -9.81 0">
            <SphereLoader filename="mesh/liver.sph" />
            <MechanicalObject name="spheres" position="@[-1].position" />
            <TSphereModel name="CollisionModel" listRadius="@[-2].listRadius"/>
            <BarycentricMapping name="sphere mapping" input="@../dofs" output="@spheres" />
        </Node>
    </Node>
</Node>

[Hugo]

Hey Xiaojuan,

I checked the file “mesh/liver2.msh” it looks pretty strange (look ids are not starting from 1). I wouldn’t use it.

In your last post (#10128), the first scene is working fine for me.

Best,

Hugo

[shang]

yes, Hugo, you’re right, the first scene works fine, I noticed that the bug in the .msh file too, I guess that I have to reformat my msh file like the format of file ” mesh/liver”. I will have a try,thank you , I will update the issue if it is solved!

[shang]

hi Hugo, I have reformated my .msh file with like the format, namely the text style as the “mesh/liver”, and it worked fine, I successfully load the msh file and simulated the scene with TetrahedralTensorMassForceField using the .msh file.
I will make this issue as solved, thank you!

[Author]

Posts