$ TEST No IC8 - SHEAR DIFFUSION $ $ Short description: $ For case 1 the system is generated with an element mesh consisting of $ 24x2 quadrilateral elements and the respective beam elements (type TRUS). $ In the subsequent cases the number of elements are doubled, up to $ case 3 having an element mesh of 96x8 quadrilateral elements. $ The calculation is made with the program ASE. $ $ Results: $ - Comparative value Sxy in point E: Desired value 27.8 MPa analytical, 34.3 MPa with $ 8-node quad-element of second order and 3-node beam element $ obtained values: $ 24x2: 25.463 MPa $ 48x4: 28.804 MPa $ 96x8: 28.939 MPa $ $ - further results of the ASE-calculation: values in MPa $ $ Shear stress curve in the middle of the panel with 24x2 elements is depicted $ inadequately, especially in the region from point E up to the end of the panel. $ The outcome of the shear stress curve in the middle of the panel for 96x8 $ elements is very similar (minor deviations, 5%) to the results of the FE-calculation $ of the benchmark - larger deviations to the analytical result can only be found in $ the middle of the panel-end (analytical -0.001 MPa, Benchmark 8.88 MPa, ASE 9.577 MPa). $ $ Normal stress curve in Beam 1 $ - Maximum value at x=1.8m, difference of ASE to benchmark: for 24x2 elements $ deviation of 10%, for 48x4 elements deviation of 3%, for 96x8 elements $ deviation 0.9% $ (Benchmark: 46.8 MPa, ASE: 47.2 MPa = 9.4kN/0.2E-3m2, analytical: 44.1 MPa) $ - In the middle of the beam the results of the Benchmark, ASE (48x4 and 96x8 $ elements) and the analytical calculation correspond: 22 MPA or 4.4 kN/0.2E-3 $ - larger deviations for ASE at the beginning of the beam (x=0 or 0.0375 m): $ Benchmark 0.115 MPA, ASE 24x2: 3.5 MPa, 48x4: 1.5 MPa, 96x8: 1.0 MPA, $ analytical: -0.244 MPa $ +PROG GENF urs:1 HEAD TEST No IC8 - SHEAR DIFFUSION TXB element mesh:24x2 elements SYST FRAM GDIV 10000 MAT 1 E 210e+6 MUE 0.3 NODE 1 0.0 0.0 25 1.8 0.0 101 0.0 0.075 125 1.8 0.075 MESH 1 101 125 25 2 24 1 T1 0.712E-3 SVAL 1 1 A 0.2e-3 SVAL 2 1 A 0.4e-3 TRUS (1001 1024 1) (1 1) (2 1) NCS 1 TRUS (2001 2024 1) (101 1) (102 1) NCS 1 TRUS (3001 3002 1) (1 50) (51 50) NCS 2 NODE (1 25 1) FIX py 101 FIX px END +PROG ASE urs:2 HEAD TEST No IC8 - SHEAR DIFFUSION TXB Elementnetz: 24x2 Elemente CTRL QTYP 0 LC 1 POLO 1 1.8 0. 0. PX 10 END +PROG WING urs:3 HEAD CTRL MEAS 2D PAGE FIRS 11 UNIO 0 UNII 0 MOVE SIZE TYPE LP SC 0 W - H - MARG FULL FORM STAN HLEG 2.50000 $$ WLEG 0 SCAL - SCHH H1 .21765 H2 .21765 H3 .24500 H4 .10000 $$ H5 .17150 H6 .17150 H7 .17150 HS .56943 MS 0 $$ FACW 1.00000 REVI 100.000 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 NEGY 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 QUAD TYPE T UNIT 1.00000 SCHH .15000 ND 6 LOAD TYPE FULL UNIT -2.50000 LC NO 1 T 10.00000E+09 DESI 1 QUAD TYPE SXM UNIT 1.00000 SCHH .10000 STYP NODE ND 3 QUAD TYPE SYM UNIT 1.00000 SCHH .10000 STYP NODE ND 3 QUAD TYPE TXYM UNIT 1.00000 SCHH .10000 STYP NODE ND 3 BEAM N UNIT 10.00000 SCHH .10000 ND 3 STYP TRUS MOVE 0. 0.0375 DRAW 1.8 0.0375 DTYP ZZ SECT TXYM ND 3 END $---------------------------------------------------------------------------------- +PROG GENF urs:4 HEAD TEST No IC8 - SHEAR DIFFUSION TXB element mesh:48x4 elements SYST FRAM GDIV 10000 MAT 1 E 210e+6 MUE 0.3 NODE 1 0.0 0.0 49 1.8 0.0 201 0.0 0.075 249 1.8 0.075 MESH 1 201 249 49 4 48 1 T1 0.712E-3 SVAL 1 1 A 0.2e-3 SVAL 2 1 A 0.4e-3 TRUS (1001 1048 1) (1 1) (2 1) NCS 1 TRUS (2001 2048 1) (201 1) (202 1) NCS 1 TRUS (3001 3004 1) (1 50) (51 50) NCS 2 NODE (1 49 1) FIX py 201 FIX px END +PROG ASE urs:5 HEAD TEST No IC8 - SHEAR DIFFUSION TXB element mesh:48x4 elements CTRL QTYP 0 LC 2 POLO 1 1.8 0. 0. PX 10 END +PROG WING urs:6 HEAD CTRL MEAS 2D PAGE FIRS 11 UNIO 0 UNII 0 MOVE SIZE TYPE LP SC 0 W - H - MARG FULL FORM STAN HLEG 2.50000 $$ WLEG 0 SCAL - SCHH H1 .21765 H2 .21765 H3 .24500 H4 .05000 $$ H5 .17150 H6 .17150 H7 .17150 HS .56943 MS 0 $$ FACW 1.00000 REVI 100.000 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 NEGY 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 QUAD TYPE T UNIT 1.00000 SCHH .05000 ND 6 LOAD TYPE FULL UNIT -2.50000 LC NO 2 T 10.00000E+09 DESI 1 QUAD TYPE SXM UNIT 1.00000 SCHH .05000 STYP NODE ND 3 QUAD TYPE SYM UNIT 1.00000 SCHH .05000 STYP NODE ND 3 QUAD TYPE TXYM UNIT 1.00000 SCHH .05000 STYP NODE ND 3 BEAM N UNIT 7.50000 SCHH .10000 ND 3 STYP TRUS MOVE 0. 0.0375 DRAW 1.8 0.0375 DTYP ZZ SECT TXYM ND 3 END $-------------------------------------------------------------------------- +PROG GENF urs:7 HEAD TEST No IC8 - SHEAR DIFFUSION TXB element mesh:96x8 elements SYST FRAM GDIV 10000 MAT 1 E 210e+6 MUE 0.3 NODE 1 0.0 0.0 97 1.8 0.0 801 0.0 0.075 897 1.8 0.075 MESH 1 801 897 97 8 96 1 T1 0.712E-3 SVAL 1 1 A 0.2e-3 SVAL 2 1 A 0.4e-3 TRUS (1001 1096 1) (1 1) (2 1) NCS 1 TRUS (2001 2096 1) (801 1) (802 1) NCS 1 TRUS (3001 3008 1) (1 100) (101 100) NCS 2 NODE (1 97 1) FIX py 801 FIX px END +PROG ASE urs:8 HEAD TEST No IC8 - SHEAR DIFFUSION TXB element mesh:96x8 elements CTRL QTYP 0 LC 3 POLO 1 1.8 0. 0. PX 10 END +PROG WING urs:9 HEAD CTRL MEAS 2D PAGE FIRS 11 UNIO 0 UNII 0 MOVE SIZE TYPE LP SC 0 W - H - MARG FULL FORM STAN HLEG 2.50000 $$ WLEG 0 SCAL - SCHH H1 .21765 H2 .21765 H3 .24500 H4 .05000 $$ H5 .17150 H6 .17150 H7 .17150 HS .56943 MS 0 $$ FACW 1.00000 REVI 100.000 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 3 T 10.00000E+09 DESI 1 VIEW TYPE STAN X 0 Y 0 Z 1.00000 AXIS NEGY 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 QUAD TYPE T UNIT 1.00000 SCHH .05000 ND 6 LOAD TYPE FULL UNIT -2.50000 LC NO 3 T 10.00000E+09 DESI 1 QUAD TYPE SXM UNIT 1.00000 SCHH .05000 STYP NODE ND 3 QUAD TYPE SYM UNIT 1.00000 SCHH .05000 STYP NODE ND 3 QUAD TYPE TXYM UNIT 1.00000 SCHH .05000 STYP NODE ND 3 BEAM N UNIT 7.50000 SCHH .05000 ND 3 STYP TRUS MOVE 0. 0.0375 DRAW 1.8 0.0375 DTYP ZZ SECT TXYM ND 3 END