

August 27, 2018 at 4:42 pmericzhngSubscriber
I 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 pmJohn DoyleAnsys Employee
Try 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 pmericzhngSubscriber
Thanks 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 pmericzhngSubscriber
Sorry, I forgot to reply to you in last post.

November 11, 2019 at 2:43 pmemixSubscriber
I 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.