Crowdsourcing: 2D axisymmetric rotor eigenbuckling analysis

BD
Bohlen, Dan (GE Aerospace, US)
Mon, Apr 15, 2024 11:56 AM

Hi All,

I have a puzzler.  Performing my first analysis of this type, but having run modal analysis I have some background.  I was able to run this analysis on 2 other parts successfully.  I had to fight through some bugs/issues in the model  (loose pieces, odd results in purely orthotropic regions), but I got results:  a bunch of eigenvalues and eigenvectors for a bunch of different nodal diameter.

So moving on to the third structure I'm all confident.  I solve for my first nodal diameter and the output file shows 20 reasonable eigenvalues (see below)  Same prestress at the previous structure - really the same .db - just a different portion of the rotor. The issue is Ansys only output 5 sets of results to the .rst file.  Those 5 results look fine.  I tried again running interactively (versus a script) and got the same thing.  I am stumped.

-----------  output from bucking solution....

USE SYM. BLOCK LANCZOS METHOD TO EXTRACT      20 MODES
EXTRACT EIGENVALUES IN THE RANGE  0.0000    TO  0.10000E+09

INCLUDE PRESTRESS EFFECTS FROM A PREVIOUS ANALYSIS

                   S O L U T I O N   O P T I O N S

PROBLEM DIMENSIONALITY. . . . . . . . . . . . .3-D
DEGREES OF FREEDOM. . . . . . UX  UY  UZ
ANALYSIS TYPE . . . . . . . . . . . . . . . . .BUCKLING
EXTRACTION METHOD. . . . . . . . . . . . . .BLOCK LANCZOS
PRESTRESS EFFECTS INCLUDED IF AVAILABLE . . . .YES
BUCKLING LOAD EXTRACTION RANGE. . . . . . . . .  0.0000    TO  0.10000E+09
GLOBALLY ASSEMBLED MATRIX . . . . . . . . . . .SYMMETRIC

                  L O A D   S T E P   O P T I O N S

LOAD STEP NUMBER. . . . . . . . . . . . . . . .    1
INERTIA LOADS                        X          Y          Z
OMEGA. . . . . . . . . . . .  0.0000    0.10472      0.0000
PRINT OUTPUT CONTROLS . . . . . . . . . . . . .NO PRINTOUT
DATABASE OUTPUT CONTROLS. . . . . . . . . . . .ALL DATA WRITTEN

BLOCK LANCZOS CALCULATION OF UP TO  20 EIGENVECTORS.
NUMBER OF EQUATIONS              =      141057
MAXIMUM WAVEFRONT                =          27
MAXIMUM MODES STORED            =          20
MINIMUM EIGENVALUE              =  0.00000E+00
MAXIMUM EIGENVALUE              =  0.10000E+09

***** EIGENVALUES (LOAD MULTIPLIERS FOR BUCKLING) *****
*** FROM BLOCK LANCZOS ITERATION ***

  SHAPE NUMBER   LOAD MULTIPLIER

        1           0.37446739
        2           0.40178640
        3           0.43843860
        4           0.44066720
        5           0.56409493
        6           0.57120969
        7           0.64626349
        8           0.68360993
        9           0.73609141
       10           0.74477059
       11           0.86760640
       12           0.90873043
       13           0.97509609
       14            1.1254436
       15            1.2205556
       16            1.3094186
       17            1.3385111
       18            1.4237132
       19            1.5651896
       20            1.6147387

---------------------  when the script goes to plot the last cases - no results case 6  -------------------------

CUMULATIVE DISPLAY NUMBER  1 WRITTEN TO FILE thbuc1004.png    - RASTER MODE.
DISPLAY TITLE=
prestress buckling ND=0 Mode=5

ALL ANNOTATION IS DELETED.

USE LOAD STEP    1  SUBSTEP    6  FOR LOAD CASE 0

*** ERROR ***                          CP =      11.812  TIME= 10:45:33
Load set not found on result file thbuc1.rst.

Dan Bohlen
Senior Engineer, Stress Analysis
STAR review chairman, military structures
GE Aerospace
1 Neumann Way
Evendale, OH  45215  USA

Building B90 Column  L 3.75
M/D H358  Cell  513-917-3402

Building 200  Desk Phone 3-8816

GE FOCUS:  Safety, Quality, Delivery, Cost

"In God we trust, all others bring data." W Edwards Deming

[IMG_0492]

Hi All, I have a puzzler. Performing my first analysis of this type, but having run modal analysis I have some background. I was able to run this analysis on 2 other parts successfully. I had to fight through some bugs/issues in the model (loose pieces, odd results in purely orthotropic regions), but I got results: a bunch of eigenvalues and eigenvectors for a bunch of different nodal diameter. So moving on to the third structure I'm all confident. I solve for my first nodal diameter and the output file shows 20 reasonable eigenvalues (see below) Same prestress at the previous structure - really the same .db - just a different portion of the rotor. The issue is Ansys only output 5 sets of results to the .rst file. Those 5 results look fine. I tried again running interactively (versus a script) and got the same thing. I am stumped. ----------- output from bucking solution.... USE SYM. BLOCK LANCZOS METHOD TO EXTRACT 20 MODES EXTRACT EIGENVALUES IN THE RANGE 0.0000 TO 0.10000E+09 INCLUDE PRESTRESS EFFECTS FROM A PREVIOUS ANALYSIS S O L U T I O N O P T I O N S PROBLEM DIMENSIONALITY. . . . . . . . . . . . .3-D DEGREES OF FREEDOM. . . . . . UX UY UZ ANALYSIS TYPE . . . . . . . . . . . . . . . . .BUCKLING EXTRACTION METHOD. . . . . . . . . . . . . .BLOCK LANCZOS PRESTRESS EFFECTS INCLUDED IF AVAILABLE . . . .YES BUCKLING LOAD EXTRACTION RANGE. . . . . . . . . 0.0000 TO 0.10000E+09 GLOBALLY ASSEMBLED MATRIX . . . . . . . . . . .SYMMETRIC L O A D S T E P O P T I O N S LOAD STEP NUMBER. . . . . . . . . . . . . . . . 1 INERTIA LOADS X Y Z OMEGA. . . . . . . . . . . . 0.0000 0.10472 0.0000 PRINT OUTPUT CONTROLS . . . . . . . . . . . . .NO PRINTOUT DATABASE OUTPUT CONTROLS. . . . . . . . . . . .ALL DATA WRITTEN BLOCK LANCZOS CALCULATION OF UP TO 20 EIGENVECTORS. NUMBER OF EQUATIONS = 141057 MAXIMUM WAVEFRONT = 27 MAXIMUM MODES STORED = 20 MINIMUM EIGENVALUE = 0.00000E+00 MAXIMUM EIGENVALUE = 0.10000E+09 ***** EIGENVALUES (LOAD MULTIPLIERS FOR BUCKLING) ***** *** FROM BLOCK LANCZOS ITERATION *** SHAPE NUMBER LOAD MULTIPLIER 1 0.37446739 2 0.40178640 3 0.43843860 4 0.44066720 5 0.56409493 6 0.57120969 7 0.64626349 8 0.68360993 9 0.73609141 10 0.74477059 11 0.86760640 12 0.90873043 13 0.97509609 14 1.1254436 15 1.2205556 16 1.3094186 17 1.3385111 18 1.4237132 19 1.5651896 20 1.6147387 --------------------- when the script goes to plot the last cases - no results case 6 ------------------------- CUMULATIVE DISPLAY NUMBER 1 WRITTEN TO FILE thbuc1004.png - RASTER MODE. DISPLAY TITLE= prestress buckling ND=0 Mode=5 ALL ANNOTATION IS DELETED. USE LOAD STEP 1 SUBSTEP 6 FOR LOAD CASE 0 *** ERROR *** CP = 11.812 TIME= 10:45:33 Load set not found on result file thbuc1.rst. Dan Bohlen Senior Engineer, Stress Analysis STAR review chairman, military structures GE Aerospace 1 Neumann Way Evendale, OH 45215 USA Building B90 Column L 3.75 M/D H358 Cell 513-917-3402 Building 200 Desk Phone 3-8816 GE FOCUS: Safety, Quality, Delivery, Cost "In God we trust, all others bring data." W Edwards Deming [IMG_0492]