TAGGED: apdl, cyclic-symmetry, prestress
February 21, 2023 at 4:56 pmbiao.zhouSubscriber
I've got a problem with exporting the prestressed stiffness matrix of a sector model with cyclic symmetry. I'm looking for the way of recovering the condensed cyclic interface DOFs in the global stiffness matrix.
For a non-prestress sector model, I managed to export its full stiffness/mass matrix by running a modal analysis for the FREE model without constrained DOFs and cyclic interfaces and then by reading the matrix from the '.full' file.
For the sector model in the spinning condition, I have to firstly run a cyclic static analysis with the cyclic status on and constrained DOFs in the hub.
CYCLIC/SOLUTIONANTYPE,STATIC !STATIC ANALYSISOMEGA,0, 0, Vel_Rot_RADRESCONTRAL,LINEAR,ALL,1! BOUNDARY CONSTRAINTS: selected nodes at the disk hubALLSEL,ALLCMSEL,S,NSEL_BC,NODED, ALL, ALL, 0ALLSEL,ALLSOLVEFINISH
In the subsequent modal analysis, I tried to set the cyclic status off but failed. The constraints in the hub are deleted. Then with the cyclic status on, the exported stiffness/mass matrix only covers a part of the nodes in the sector model. It seems that the DOFs attached on one of the cyclic interfaces are condensed out in the global matrix. Moreover, it seems the node mapping with the exported matrice do not match the original sector model./SOLUANTYPE,STATIC, RESTART, , , PERTURBPERTURB,MODAL,,,NOKEEPSOLVE, ELFORMMODOPT,LANB,1,WRFULL,1! DELETE ALL NODE CONSTRAINTSALLSEL,ALLDDELE, ALL, ALLALLSEL,ALLSOLVEFINISHSo I'm looking for the way of exporting the full prestressed stiffness matrix with all the DOFs in the original sector model (free body) included.Any suggestions and comments will be highly appreciated.Thanks,
February 25, 2023 at 12:34 amBill BulatAnsys Employee
Are you using MAPDL directly of Mechanical with command objects?
I'm not perfectly clear on what you want to do or if it can be accomplished. When you SOLVE after the CYCLIC command has been issued, the first things that happen (before the solution is calculated) is a duplicate sector of the original sector is created and I think too that constraint equations are internally and transparently created that relate the nodal DOFs on the periodic sectors. Looking at the documentation of the CE command,
"The first unique degree of freedom in the equation is eliminated in terms of all other degrees of freedom in the equation."
So I would imagine that DOFs of at least some of the nodes on the periodic sectors will not be included in the system stiffness matrix.
Do you want a stiffness matrix of only the original sector or do you instead want one including both the original and duplicate sectors? I think the latter is possible, though my guess is that terms for some of the nodal DOFs on the periodic sectors will be missing.
February 25, 2023 at 2:42 pmbiao.zhouSubscriber
Thanks for your response. Yes, I am running the MAPDL codes directly. I want to obtain a stiffness matrix with all the dofs in the original base sector included.
Let me explain clearly how I attemped to do that.
% The reference sector model in Ansys firstly undergoes a cyclic static
% analysis in order to include the rotation-induced prestress effect.
% Then a linear perturbation modal analysis for the sector model with
% cyclic symmetry and without any constrained dof is performed in order to
% export the sector-level mass/stiffness matrix (of only the original base sector) for % the following model reduction research.
After the second step, Ansys gives the exported matrices of the original base sector. However, a part of the DOFs attached onto the cyclic symmetry interfaces are condensed in the exported K and M matrix. More specifically, both the high and low cyclic interface dofs are partially retained in the global matrix.
Now I'm trying to write some codes in order to read the imcomplete sector matrix and further to restore the terms related to the condensed cyclic interface dofs.
But this is not a easy task. Hope I can make it.
- You must be logged in to reply to this topic.
Boost Ansys Fluent Simulations with AWS
Computational Fluid Dynamics (CFD) helps engineers design products in which the flow of fluid components is a significant challenge. These different use cases often require large complex models to solve on a traditional workstation. Click here to join this event to learn how to leverage Ansys Fluids on the cloud, thanks to Ansys Gateway powered by AWS.
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.
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
- Solver Pivot Warning in Beam Element Model
- Understanding Force Convergence Solution Output
- An Unknown error occurred during solution. Check the Solver Output…..
- What is the difference between bonded contact region and fixed joint
- whether have the difference between using contact and target bodies
- The solver engine was unable to converge on a solution for the nonlinear problem as constrained.
- User manual
- Colors and Mesh Display
- material damping and modal analysis
© 2023 Copyright ANSYS, Inc. All rights reserved.