The Ansys Innovation Space website recently experienced a database corruption issue. While service has been restored there appears to have been some data loss from November 13. We are still investigating and apologize for any issues our users may have as a result.
General Mechanical

General Mechanical

Topics related to Mechanical Enterprise, Motion, Additive Print and more.

Is node order in the stiffness matrix different from the node order in solution?

    • C Deho
      Subscriber

      I performed a modal analysis and used APDL SMAT to extract the stiffness matrix. Then, in MATLAB, I applied a force to  just one node, but the results are significantly different from what I obtained from a static analysis. In ANSYS, excluding the fixed points, deformation occurred at all nodes. However, when performed in MATLAB, more than half of the nodes showed zero deformation. I used the governing equation F=Kx. 2Why might this be happening?

      Size of stiffness matrix is number of nodes*3 (X, Y, and Z directions). So, does this mean the order is node1x, 1y, 1z, ... nx, ny, nz? The analysis was conducted with default settings. 

       

       

       

    • mrife
      Ansys Employee

      Hi C Deho

      Please check the Ansys MAPDL help on the *SMAT command.  MAPDL does reorder the nodes/elements from what we as users see in the FEM.  And there is a mapping vector available that relates the FEM to Solving order.   

      • C Deho
        Subscriber

        Hi mike

        Do you mean that the node numbers seen in Workbench (Node ID) differ from those in the exported stiffness matrix? 

         

         

    • mrife
      Ansys Employee

      Hi C Deho

      I've replied a few times to this and received an odd error each time...so sorry if there ends up being multiple copies of basically the same answer.

      Anyway, the stiffness matrix does not use node IDs.  Nor does the force vector - it uses the degree of freedom IDs.  In any case yes, the IDs used to solve differ from those that WB Mechanical shows us for the mesh.  Hence the need for the vector of ID mapping.

      • C Deho
        Subscriber

        Hi mike

        I extracted a mapping vector using NOD2BCS and the resulting figure appeared as shown. It seems to map excluding fixed points starting from [2220,1]. What I'm curious about is what the values on both sides mean. Just like in the figure, there are fixed nodes at both the 200th and 800th positions, but in the mapping vector, all the values seem to be gathered towards the end. Isn't NOD2BCS supposed to provide the relationship between node IDs and the degrees of freedom (DOF) used by the solver?

        I'm using a student license, so I can't access education sources provided by Ansys. I've also searched the internet, but I can't find the information I'm looking for. It would be greatly appreciated if you could provide a more detailed explanation on how to match the DOF of the stiffness matrix with the node IDs for F=Kx calculation, though I know it might be a hassle.

         

        Regards,

        Deho

         

         

Viewing 2 reply threads
  • You must be logged in to reply to this topic.