5.7. Example Damped Modal Cyclic Symmetry Analysis

This example damped modal cyclic symmetry analysis presents a simplified ring-strut-ring structure used in many rotating-machinery applications.

5.7.1. Problem Description

This example is very similar to Example Modal Cyclic Symmetry Analysis, except that a damped modal analysis is performed. The component is a simplified fan inlet case for a military aircraft engine.

5.7.2. Problem Specifications

The problem specifications follow those of Example Modal Cyclic Symmetry Analysis, with the addition of alpha (mass proportional) and structural damping coefficients equal to 0.1.

5.7.3. Input File for the Analysis

Use the input file (input_cyclicExample07.dat, download: input_cyclicExample07.zip) given below to perform this example damped modal cyclic symmetry analysis. The file contains the complete geometry, material properties, and solution options for the finite element model.

! Damped modal cyclic symmetry analysis for a ring-strut-ring configuration.

! Plot settings.
/view,1,1,1,2
/plopts,minm,0
/plopts,date,0
/plopts,logo,on
/pnum,real,1
/number,1

/prep7

! Define geometry, mesh, and boundary conditions or read the CDB file.
generateGeom = 0                    ! Set value to 1 to generate model.
*if,generateGeom,eq,1,then
    ! Define geometry parameters.
    r1 = 5
    r2 = 10
    d1 = 2
    nsect = 24
    alpha_deg = 360 / nsect
    alpha_rad = 2*acos(-1) / nsect

    ! Set coordinate system to cylindrical.
    csys,1

    ! Create keypoints, lines, and areas.
    ! Define keypoints and lines.
    k,1,0,0,0
    k,2,0,0,d1
    k,3,r1,0,0
    k,4,r1,0,d1
    l,3,4

    ! Rotate line to create area.
    arotat,1,,,,,,1,2,alpha_deg/2

    k,7,r2,0,0
    k,8,r2,0,d1
    l,7,8
    arotat,5,,,,,,1,2,alpha_deg/2
    arotat,2,,,,,,1,2,alpha_deg/2
    arotat,6,,,,,,1,2,alpha_deg/2

    a,5,6,10,9

    ! Mesh with shell-181 elements.
    mshkey,1
    et,1,181
    r,1,0.20     ! Thickness of the inner and outer rings.
    r,2,0.1      ! Thickness of the strut.
    mp,ex,1,10e6
    mp,prxy,1,0.3
    mp,dens,1,1e-4
    esize,0.25

    ! Select areas and assign attributes.
    asel,,,,1,4
    aatt,,1
    asel,,,,5
    aatt,,2
    allsel
    amesh,all

    ! Apply boundary conditions.
    dk,5,uz,0

    ! Save the database.
    cdwrite,comb,cyclicExample07,cdb
*else
    cdread,comb,cyclicExample07,cdb,,,
*endif

! Configure cyclic symmetry.
cyclic

! Plot elements.
eplot

! List cyclic status.
cyclic,status
! List cyclic solution option settings.
cycopt,status

finish

! Perform a damped modal solve.
/solution
csys,1
antype,modal
modopt,damp,6       ! Select damped eigensolver (unity norm by default).
mxpand,6,,,yes
alphad,0.1          ! Apply global alpha damping.
dmpstr,0.1          ! Apply structural damping.
solve
finish
/post1
! Print eigenvalues.
set, list,,, ,,, ,order
! Set cyclic expansion
/cycexpand,,on

rsys,1

! Read results for load step 1, substep 3
! Plot the real and imaginary tangential displacement contours.
set,1,3,,real
plns,u,y

set,1,3,,imag
plns,u,y

! Read results for load step 13, substep 1
! Plot the real tangential displacement contour.
set,13,1,,real
plns,u,y

! Read results for load step 3, substep 5
! Plot the real and imaginary tangential displacement contours.
set,3,5,,real
plns,u,y

set,3,5,,imag
plns,u,y
fini

5.7.4. Analysis Steps

The following table describes the input listing and the steps involved in this example damped modal cyclic symmetry analysis in more detail.

Step Description Command
1.

Create the Ring-Strut-Ring model.

---

2.

Specify a cyclic symmetry analysis and configure the database accordingly.

CYCLIC
3.

Specify the damped eigensolver.

MODOPT,DAMP
4.

Apply global alpha damping.

ALPHAD,0.1
5.Apply structural damping. DMPSTR,0.1
6.

Solve the modal cyclic symmetry analysis.

7.List complex eigenvalues. SET,LIST
8.

Specify the global cylindrical coordinate system.

9.

Read results for "load step 1 - substep 3 - harmonic index 0." This can be done for both real and imaginary results.

SET,1,3,,REAL
SET,1,3,,IMAG
10.

Plot the real and imaginary tangential displacement contour.

Executing this step causes the struts of the assembly to bend "in phase."

PLNSOL,U,Y
11.Read results for "load step 13 - substep 1 - harmonic index 12." SET,13,1,REAL
12.

Plot the real tangential displacement contour.

Executing this step causes the struts of the assembly to bend "out of phase."

PLNSOL,U,Y
13.

Read results for "load step 3 - substep 5 - harmonic index 1."

SET,3,5,,REAL
SET,3,5,,IMAG
14.

Plot the imaginary tangential displacement contour.

----

This step completes the example damped modal cyclic symmetry analysis.

PLNSOL,U,Y

The results of your analysis should match those shown below.

Table 5.2: Sample Complex Eigenvalues

SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE HRM-INDEX
1-1.04E-041110
2-5.40E-051110
3-1.58E-021220
45.40E-051220
5-17.821330
6357.131330
717.8041440
8-357.131440
9-38.8351550
10778.481550
1138.8191660
12-778.481660
13-6.00E-122171
147.57E-152171
15-3.47E-052281
163.87E-062281
17-1.59E-022391
18-3.87E-062391
19-1.59E-0224101
201.63E-1724101
21-33.08925111
22663.2725111
2333.07326121
24-663.2726121
2510.16431132
26-203.9531132
2710.16432142
...... ...... ... ...

Figure 5.22: Damped Modal Cyclic Symmetry Analysis Results - Real

Damped Modal Cyclic Symmetry Analysis Results - Real

Figure 5.23: Damped Modal Cyclic Symmetry Analysis Results - Imaginary

Damped Modal Cyclic Symmetry Analysis Results - Imaginary


Note:  Mode shape values may vary slightly depending on your computer system.