*HEADING
Quarter model: FAN=HUB/4+BLADE fan_cyclicsymmodel
Cyclic symmetry model option and steady state dynamics
**
*RESTART,WRITE,FREQUENCY=1
**
** NODE DEFINITION: HUB/4+BLADE
** Nodes for the hub: 1--90
**
*NODE
1, 3.50000, 3.00000, 4.00000
3, 3.46638, 3.18024, 4.00000
7, 3.18024, 3.46638, 4.00000
9, 3.00000, 3.50000, 4.00000
**
28, 6.00000, 3.00000, 4.00000
30, 5.79830, 4.08144, 4.00000
34, 4.08144, 5.79830, 4.00000
36, 3.00000, 6.00000, 4.00000
**
82, 6.00000, 3.00000, 0.00000
84, 5.79830, 4.08144, 0.00000
88, 4.08144, 5.79830, 0.00000
90, 3.00000, 6.00000, 0.00000
**
9998, 3.00000, 3.00000, 4.00000
9999, 3.00000, 3.00000, 0.00000
*NGEN,LINE=C,NSET=NHUB
1,3,1,9998
3,7,1,9998
7,9,1,9998
*NGEN,LINE=C,NSET=RINGF
28,30,1,9998
30,34,1,9998
34,36,1,9998
*NGEN,LINE=C,NSET=RINGB
82,84,1,9999
84,88,1,9999
88,90,1,9999
*NFILL,NSET=HUB
NHUB,RINGF,3,9
RINGF,RINGB,6,9
**
** Total nodes for blade: 91--120,
** plus common nodes 57--61
**
*NODE,NSET=BLADE
91, 6.34825, 4.61208, 2.00000
92, 6.04542, 5.17363, 2.04275
93, 5.64799, 5.65929, 2.11107
94, 5.16740, 6.05420, 2.19011
95, 4.61648, 6.34954, 2.25973
96, 6.97098, 5.04709, 2.00000
97, 6.60952, 5.64419, 2.08975
98, 6.14537, 6.18965, 2.22829
99, 5.61424, 6.64337, 2.39622
100, 5.04534, 6.97776, 2.55559
101, 7.71735, 5.22982, 2.00000
102, 7.23607, 5.92235, 2.13475
103, 6.57818, 6.65456, 2.33376
104, 5.87275, 7.29336, 2.58456
105, 5.23136, 7.72660, 2.83877
106, 8.44058, 5.48119, 2.00000
107, 7.86154, 6.25188, 2.17913
108, 7.04773, 7.13818, 2.43316
109, 6.20063, 7.93303, 2.76013
110, 5.49725, 8.45477, 3.10513
111, 8.94199, 6.04637, 2.00000
112, 8.35346, 6.80121, 2.24584
113, 7.55742, 7.64823, 2.59979
114, 6.73147, 8.41155, 3.04013
115, 6.03050, 8.93642, 3.47975
116, 8.88910, 6.78398, 2.00000
117, 8.41737, 7.43225, 2.26100
118, 7.89643, 8.00249, 2.59810
119, 7.35361, 8.48209, 3.02957
120, 6.78399, 8.88910, 3.50000
**
** NODE SETS
**
*NSET,NSET=XAXIS,GENERATE
10,82,9
*NSET,NSET=YAXIS,GENERATE
18,90,9
*surface,type=node,name=secondary1
YAXIS,1.5
*surface,type=node,name=main1
XAXIS,1.5
*tie,cyclicsym,position=0.5,name=tie
secondary1,main1
*cyclic symmetry model,n=4
3., 3., 0., 3., 3., 10.
*ELEMENT,TYPE=S4R
1,1,10,11,2
*ELGEN,ELSET=HUB
1,8,1,1,9,9,8
*ELEMENT,TYPE=S4R
73,57,91,92,58
77,91,96,97,92
*ELGEN,ELSET=CYLBLADE
73,4,1,1
77,4,1,1,5,5,4
**
*ELSET,ELSET=HUBF,GENERATE
1, 24, 1
*ELSET,ELSET=HUBT,GENERATE
25, 72, 1
**
*SHELL SECTION,ELSET=HUBF,MATERIAL=MAT
0.250, 3
*SHELL SECTION,ELSET=HUBT,MATERIAL=MAT
0.125, 3
*SHELL SECTION,ELSET=CYLBLADE,MATERIAL=MAT
0.125, 3
*MATERIAL,NAME=MAT
*DENSITY
7850.,
*ELASTIC,TYPE=ISOTROPIC
2.E+11, 0.29
*NSET,NSET=NPRINT,GENERATE
28,36,1
55,63,1
**
*NSET,NSET=PRINT,GEN
116,120,1
*BOUNDARY
NHUB,ENCASTRE
**
*STEP,nlgeom,INC=1000
Step 1: centrifugal loading
*STATIC
.1,1.
**
*DLOAD
CYLBLADE, CENT, 6.3e6, 3., 3., 0., 0., 0., -1.
HUBF, CENT, 6.3e6, 3., 3., 0., 0., 0., -1.
HUBT, CENT, 6.3e6, 3., 3., 0., 0., 0., -1.
**
*EL PRINT,FREQ=0
*NODE FILE,FREQ=2,NSET=PRINT
U
*NODE PRINT,FREQ=0
*NODE FILE,FREQ=0
**
*END STEP
**
*STEP
** *FREQUENCY,eigensolver=lanczos
*FREQUENCY,eigensolver=lanczos,normalization=mass
26,,,,8
*SELECT CYCLIC SYMMETRY MODES, NMIN=0, NMAX=0
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*MODAL FILE
*END STEP
*STEP
*STEADY STATE DYNAMICS,FREQUENCY SCALE=LINEAR,INTERVAL=RANGE
6.0,137.0,10,4.
*DLOAD, LOAD CASE=1, OP=NEW, CYCLIC MODE=0
CYLBLADE, P, 6.3e6
HUBF, P, 6.3e6
HUBT, P, 6.3e6
*DLOAD, LOAD CASE=2, OP=NEW, CYCLIC MODE=0
CYLBLADE, P, 6.3e6
HUBF, P, 6.3e6
HUBT, P, 6.3e6
*MODAL DAMPING,MODAL=DIRECT
1,26,0.03
*SELECT EIGENMODES, GENERATE
1,26,1
*EL PRINT,ELSET=CYLBLADE,FREQUENCY=10
S11,E11
*NODE PRINT,NSET=PRINT,FREQUENCY=10
U,
*NODE FILE,NSET=PRINT,FREQUENCY=10
U,V,A,RF
*MODAL PRINT,FREQUENCY=10
GU,
GA,
GPU,
*EL FILE,ELSET=CYLBLADE,FREQUENCY=10
S,
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S,
*OUTPUT,HISTORY,FREQUENCY=10
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S,
*END STEP
*STEP
** *FREQUENCY,eigensolver=lanczos
*FREQUENCY,eigensolver=lanczos,normalization=mass
56,,,,8
*SELECT CYCLIC SYMMETRY MODES, NMIN=1, NMAX=1
*OUTPUT, FIELD, FREQUENCY=0
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*NODE FILE,FREQ=0
*EL FILE,FREQ=0
*MODAL FILE
**
*END STEP
*STEP
*STEADY STATE DYNAMICS,FREQUENCY SCALE=LINEAR,INTERVAL=RANGE
6.0,137.0,10,4.
*DLOAD, LOAD CASE=1, OP=NEW, CYCLIC MODE=1
CYLBLADE, P, 6.3e6
HUBF, P, 6.3e6
HUBT, P, 6.3e6
*DLOAD, LOAD CASE=2, OP=NEW, CYCLIC MODE=1
CYLBLADE, P, 6.3e6
HUBF, P, 6.3e6
HUBT, P, 6.3e6
*MODAL DAMPING,MODAL=DIRECT
1,26,0.03
*SELECT EIGENMODES, GENERATE
1,26,1
*EL PRINT,ELSET=CYLBLADE,FREQUENCY=10
S11,E11
*NODE PRINT,NSET=PRINT,FREQUENCY=10
U,
*NODE FILE,NSET=PRINT,FREQUENCY=10
U,V,A,RF
*MODAL PRINT,FREQUENCY=10
GU,
GA,
GPU,
*EL FILE,ELSET=CYLBLADE,FREQUENCY=10
S,
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S,
*OUTPUT,HISTORY,FREQUENCY=10
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S,
*END STEP
*STEP
** *FREQUENCY,eigensolver=lanczos
*FREQUENCY,eigensolver=lanczos,normalization=mass
56,,,,8
*SELECT CYCLIC SYMMETRY MODES, NMIN=2, NMAX=2
*OUTPUT, FIELD, FREQUENCY=0
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*EL FILE,FREQ=0
*NODE FILE,FREQ=0
*MODAL FILE
**
*END STEP
*STEP
*STEADY STATE DYNAMICS,FREQUENCY SCALE=LINEAR,INTERVAL=RANGE
6.0,137.0,10,4.
*DLOAD, LOAD CASE=1, OP=NEW, CYCLIC MODE=2
CYLBLADE, P, 6.3e6
HUBF, P, 6.3e6
HUBT, P, 6.3e6
*DLOAD, LOAD CASE=2, OP=NEW, CYCLIC MODE=2
CYLBLADE, P, 6.3e6
HUBF, P, 6.3e6
HUBT, P, 6.3e6
*MODAL DAMPING,MODAL=DIRECT
1,26,0.03
*SELECT EIGENMODES, GENERATE
1,26,1
*EL PRINT,ELSET=CYLBLADE,FREQUENCY=10
S11,E11
*NODE PRINT,NSET=PRINT,FREQUENCY=10
U,
*NODE FILE,NSET=PRINT,FREQUENCY=10
U,V,A,RF
*MODAL PRINT,FREQUENCY=10
GU,
GA,
GPU,
*EL FILE,ELSET=CYLBLADE,FREQUENCY=10
S,
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S,
*OUTPUT,HISTORY,FREQUENCY=10
*NODE OUTPUT,NSET=PRINT
U,V,A,RF
*ELEMENT OUTPUT,ELSET=CYLBLADE
S,
*END STEP