-
-
August 27, 2018 at 4:42 pm
ericzhng
SubscriberI used a set of simple APDL commands (attached below) trying to get the local stiffness matrix (8 X
. The results I got is somehow different from what I compute by hand. This is for a 2D plane element (1 element), no constraints, applied with a zero force. This way we can get the full local stiffness matrix.
Also when I tried to derive mass matrix, which I have to specify a transient analysis, the local stiffness matrix changes to a 25 X 25 size matrix. And still I couldn't output the mass matrix, not sure why.
-
August 29, 2018 at 8:48 pm
John Doyle
Ansys EmployeeTry using the "EMATWRITE,yes" before the SOLVE.
and then after the solve, add the following:
/POST1
*DMAT, stiffness, D, Import, EMAT,plane2d.emat, STIFF, 1
*PRINT,stiffness,stiffness_182.txt
*DMAT, mass, D, Import, EMAT,plane2d.emat, MASS, 1
*PRINT,mass,mass_182.txt
-
August 31, 2018 at 3:39 pm
ericzhng
SubscriberThanks Doyle, That is great help! I was wondering why I can't extract mass matrix.
But same issue exists that the stiffness matrix I obtained from ANSYS and the one I computed are different. I double checked my formulation of 2d element, it should be correct.
The matK is what I imported from ANSYS, the other k is what I computed. Could you give me some hints on this, please?
Thank you!
-
August 31, 2018 at 3:40 pm
ericzhng
SubscriberSorry, I forgot to reply to you in last post.
-
November 11, 2019 at 2:43 pm
emix
SubscriberI know this post is old, but in case anybody is still interested, the difference in stiffness matrix between plane182 element and the standard Q4 elements is that the first has enhanced strain formulations to prevent shear locking.
To deactivate this function just use the following:
KEYOPT,1,1,0
By setting KEYOPT(1)=0, the stiffness matrix obtained from ANSYS would be the same as the one computed.
-
- You must be logged in to reply to this topic.

Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.

Ansys Blog
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.
- Saving & sharing of Working project files in .wbpz format
- An Unknown error occurred during solution. Check the Solver Output…..
- Understanding Force Convergence Solution Output
- Solver Pivot Warning in Beam Element Model
- Colors and Mesh Display
- How to calculate the residual stress on a coating by Vickers indentation?
- whether have the difference between using contact and target bodies
- What is the difference between bonded contact region and fixed joint
- The solver engine was unable to converge on a solution for the nonlinear problem as constrained.
- User manual
-
2600
-
2088
-
1319
-
1108
-
459
© 2023 Copyright ANSYS, Inc. All rights reserved.