$ $ TEST No IC37 - AXISYMMETRIC STIFFENED CYLINDER - PRESSURE $ $ Note print error: false OUTPUT Direct stress Syy on inner surface point C $ correct(?) OUTPUT Hoop stress Stt on inner surface point C $ Reason: Compare the analytical result with $ the benchmark results $ $ Short decription: $ For case 1 the system is generated with an element mesh of 12x1 quadrilateral elements, $ whereby the element width is equivalent to the shell thickness of 0.01 m, and two $ additional quad-elements on the outer side of the cylinder, for the thicker cross section, $ are generated. $ In the subsequent case 2 the number of elements are doubled for the shell axis, $ therefore 24x1 elements are used, with four additional quad-elements on the outer side $ of the cylinder (for the thicker cross section). $ The calculation was made with the program TALPA, $ because only this program can analyse an axisymmetrical state. $ Note: Rotational axis for benchmark test IC 37 is the Y-axis $ Rotational axis for program TALPA is the X-axis $ Stress comparison: Benchmark = TALPA $ ---------------------------- $ Syy,Saa = Sxx $ Stt = Szz $ $ Results: $ - Comparative value Stt in point C inner side of cylinder: Desired 49.11 MPa analytical, $ 1.62 MPa Benchmark fine element mesh $ with TALPA obtained values (from GRAF = Szz): $ 12x1: 45.43 MPa $ 24x1: 44.93 MPa $ $ - further results of the TALPA-calculation: values in MPa $ $ Point C inside Point C outside $ Saa Stt Saa Stt $ ------------------------------------------------------------------- $ 12x1 elem. 21.45 45.43 -2.57 36.6 $ 24x1 elem. 21.17 44.93 -3.65 35.25 $ $ Desired: $ Benchmark $ fine mesh 1.62 29.9 48.4 43.9 $ $ $ Point B inside Point B outside $ Saa Stt Saa Stt $ ------------------------------------------------------------------- $ 12x1 elem. 11.96 48.48 13.04 47.56 $ 24x1 elem. 9.73 48.59 15.25 48.97 $ $ Desired: $ Benchmark $ fine mesh 29.6 49.9 20.4 47.1 $ $ $ Displacement in point C inside: values in m $ $ horizontal vertical $ -------------------------------------------------------------------- $ 12x1 elem. 4.87E-6 9.256E-5 $ 24x1 elem. 4.87E-6 9.141E-5 $ $ Benchmark 0.2483E-4 0.7004E-4 PROG GENF urs:1 HEAD TEST No IC37 - AXISYMMETRIC STIFFENED CYLINDER - PRESSURE TXB Element mesh: coarse SYST AXIA MAT 1 E 210e+6 MUE 0.3 $ generate bottom thin part NODE 1 0.0 0.5 2 0.0 0.51 9 0.4 0.5 10 0.4 0.51 MESH 1 2 10 9 1 4 1 $ generate mid-part NODE 11 0.4888 0.5 12 0.4888 0.51 15 0.5112 0.5 16 0.5112 0.51 17 0.6 0.5 18 0.6 0.51 MESH 9 10 12 11 1 1 1 MESH 11 12 16 15 1 2 1 MESH 15 16 18 17 1 1 1 NODE 101 0.4888 0.5224 103 0.5112 0.5224 MESH 12 101 103 16 1 2 1 $ generate top thin part NODE 25 1.0 0.5 26 1.0 0.51 MESH 17 18 26 25 1 4 1 NODE 1,2 FIX pxmm NODE 25,26 FIX mm END +PROG TALPA urs:2 HEAD TEST No IC37 - AXISYMMETRIC STIFFENED CYLINDER - PRESSURE TXB Element mesh: coarse LC 1 BOLO 1 25 2 P 1000 BOLO 25 26 1 PX 12500 $ single load 125kN as line load over cylinder thickness 0.01m END +PROG WING urs:3 HEAD $ BATCH-AUSGABE DES INTERAKTIVEN GRAF; DATUM: 8.02.99; ZEIT: 11: 1:14 PAGE FIRS 20 UNIO 0 UNII 0 SIZE TYPE LP SC 0 W - H - MARG YES FORM STAN HLEG 2.50000 WLEG 0 $$ SCAL - SCHH H1 .29375 H2 .29375 H3 .24500 H4 .17150 $$ H5 .17150 H6 .17150 H7 .17150 HS .35925 MS 0 $$ FACW 1.00000 REVI 25.0000 MINH .15000 MINW .60000 COLO C1 1001 C2 3001 C3 6001 C4 8001 C5 5002 C6 7001 $$ C7 7502 C8 7501 C10 6004 C11 3101 C12 3301 C13 4001 $$ C14 2001 C15 5006 C16 1001 FILL TYPE COLO FROM - TO - INC - NUMB 10 FILL TYPE 1 FROM - TO 13 INC 1 FILL TYPE 2 FROM - TO 4 INC 1 FILL TYPE 3 FROM - TO 14 INC 1 FILL TYPE 4 FROM - TO 7 INC 1 FILL TYPE 5 FROM - TO 12 INC 1 FILL TYPE 6 FROM - TO 3 INC 1 FILL TYPE 7 FROM - TO 5 INC 1 FILL TYPE 8 FROM - TO 6 INC 1 FILL TYPE 9 FROM - TO 2 INC 1 FILL TYPE 10 FROM - TO 11 INC 1 GRP NUMB (0 1023 1) OPTI YES ; MOVE GRP NUMB ADDI OPTI YES GRP NUMB OPTI ETYP (0 1023 1) FULL - $GRP NUMB SBRI OPTI OFF GRP NUMB KINE OPTI FULL GDIV 1000 BOX LC NO 1 T 10.00000E+09 DESI 1 VIEW TYPE STAN X 0 Y 0 Z -1.00000 AXIS NEGX ROTA 0 PERS A INFI BX 0 BY 0 BZ 0 ANGL 99.0000 VISI HIDD NO DMAX -10.0000000E+11 DMIN 10.0000000E+11 DEFO TYPE NO STRU NUME 0 NUMN 1.00000 MARK 0 MFIX 5.00000 LOAD TYPE FULL UNIT -2.50000 QUAD TYPE SX UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE SY UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE SZ UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE TXY UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 BOX XMIN 0.29 YMIN 0.49 XMAX 0.61 YMAX 0.6 TYPE WIND STRU NUME 0 NUMN 1.00000 MARK 0 MFIX 5.00000 QUAD TYPE SX UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE SY UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE SZ UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE TXY UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 NODE UX UNIT 100 SCHH 0.1 ND 5 NODE UY UNIT 100 SCHH 0.1 ND 5 END $------------------------------------------------------------------------------------ PROG GENF urs:4 HEAD TEST No IC37 - AXISYMMETRIC STIFFENED CYLINDER - PRESSURE TXB Element mesh: fine SYST AXIA MAT 1 E 210e+6 MUE 0.3 $ generierung unterer dünner teil NODE 1 0.0 0.5 2 0.0 0.51 17 0.4 0.5 18 0.4 0.51 MESH 1 2 18 17 1 8 1 $ generierung mittelteil NODE 21 0.4888 0.5 22 0.4888 0.51 29 0.5112 0.5 30 0.5112 0.51 33 0.6 0.5 34 0.6 0.51 MESH 17 18 22 21 1 2 1 MESH 21 22 30 29 1 4 1 MESH 29 30 34 33 1 2 1 NODE 101 0.4888 0.5224 105 0.5112 0.5224 MESH 22 101 105 30 1 4 1 $ generierung oberer dünner teil NODE 49 1.0 0.5 50 1.0 0.51 MESH 33 34 50 49 1 8 1 NODE 1,2 FIX pxmm NODE 49,50 FIX mm END +PROG TALPA urs:5 HEAD TEST No IC37 - AXISYMMETRIC STIFFENED CYLINDER - PRESSURE TXB Element mesh: fine LC 2 BOLO 1 49 2 P 1000 BOLO 49 50 1 PX 12500 $ einzellast 125kN als linienlast über zylinderdicke 0.01m END +PROG WING urs:6 HEAD $ BATCH-AUSGABE DES INTERAKTIVEN GRAF; DATUM: 8.02.99; ZEIT: 11: 1:14 PAGE FIRS 20 UNIO 0 UNII 0 SIZE TYPE LP SC 0 W - H - MARG YES FORM STAN HLEG 2.50000 WLEG 0 $$ SCAL - SCHH H1 .29375 H2 .29375 H3 .24500 H4 .17150 $$ H5 .17150 H6 .17150 H7 .17150 HS .35925 MS 0 $$ FACW 1.00000 REVI 25.0000 MINH .15000 MINW .60000 COLO C1 1001 C2 3001 C3 6001 C4 8001 C5 5002 C6 7001 $$ C7 7502 C8 7501 C10 6004 C11 3101 C12 3301 C13 4001 $$ C14 2001 C15 5006 C16 1001 FILL TYPE COLO FROM - TO - INC - NUMB 10 FILL TYPE 1 FROM - TO 13 INC 1 FILL TYPE 2 FROM - TO 4 INC 1 FILL TYPE 3 FROM - TO 14 INC 1 FILL TYPE 4 FROM - TO 7 INC 1 FILL TYPE 5 FROM - TO 12 INC 1 FILL TYPE 6 FROM - TO 3 INC 1 FILL TYPE 7 FROM - TO 5 INC 1 FILL TYPE 8 FROM - TO 6 INC 1 FILL TYPE 9 FROM - TO 2 INC 1 FILL TYPE 10 FROM - TO 11 INC 1 GRP NUMB (0 1023 1) OPTI YES ; MOVE GRP NUMB ADDI OPTI YES GRP NUMB OPTI ETYP (0 1023 1) FULL - $GRP NUMB SBRI OPTI OFF GRP NUMB KINE OPTI FULL GDIV 1000 BOX LC NO 2 T 10.00000E+09 DESI 1 VIEW TYPE STAN X 0 Y 0 Z -1.00000 AXIS NEGX ROTA 0 PERS A INFI BX 0 BY 0 BZ 0 ANGL 99.0000 VISI HIDD NO DMAX -10.0000000E+11 DMIN 10.0000000E+11 DEFO TYPE NO STRU NUME 0 NUMN 1.00000 MARK 0 MFIX 5.00000 LOAD TYPE FULL UNIT -2.50000 QUAD TYPE SX UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE SY UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE SZ UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE TAU UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 BOX XMIN 0.29 YMIN 0.49 XMAX 0.61 YMAX 0.6 TYPE WIND STRU NUME 0 NUMN 1.00000 MARK 0 MFIX 5.00000 QUAD TYPE SX UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE SY UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE SZ UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 QUAD TYPE TXY UNIT 1.00000 SCHH .10000 STYP NODE ND 2.00000 NODE UX UNIT 100 SCHH 0.1 ND 5 NODE UY UNIT 100 SCHH 0.1 ND 5 END