*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=slave1 YAXIS,1.5 *surface,type=node,name=master1 XAXIS,1.5 *tie,cyclicsym,position=0.5,name=tie slave1,master1 *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