*HEADING EXAMPLE PROBLEM: MODE 0 SPHERE USING BGT BOUNDARY CONDITION AND AC3D8R *************************** *PARAMETER f=1.0 T=4.0/f T4=T/4.0 omega=2*3.14159*f t0=-0.25/f *PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES *************************** *NODE,NSET=POLE 1, 0., 1. 31, 0., -1. 32, 0., 1.07958 62, 0., -1.07958 63, 0., 1.15915 93, 0., -1.15915 *NODE,NSET=INNER 2, 0.104528, 0.994522 3, 0.207912, 0.978148 4, 0.309017, 0.951057 5, 0.406737, 0.913545 6, 0.5, 0.866025 7, 0.587785, 0.809017 8, 0.669131, 0.743145 9, 0.743145, 0.669131 10, 0.809017, 0.587785 11, 0.866025, 0.5 12, 0.913545, 0.406737 13, 0.951057, 0.309017 14, 0.978148, 0.207912 15, 0.994522, 0.104528 16, 1., 0.0 17, 0.994522, -0.104528 18, 0.978148, -0.207912 19, 0.951057, -0.309017 20, 0.913545, -0.406737 21, 0.866025, -0.5 22, 0.809017, -0.587785 23, 0.743145, -0.669131 24, 0.669131, -0.743145 25, 0.587785, -0.809017 26, 0.5, -0.866025 27, 0.406737, -0.913545 28, 0.309017, -0.951057 29, 0.207912, -0.978148 30, 0.104528, -0.994522 *NODE,NSET=MID 33, 0.112847, 1.07367 34, 0.224457, 1.05599 35, 0.333609, 1.02674 36, 0.439105, 0.986245 37, 0.53979, 0.934944 38, 0.634561, 0.873399 39, 0.72238, 0.802284 40, 0.802284, 0.72238 41, 0.873399, 0.634561 42, 0.934944, 0.53979 43, 0.986245, 0.439105 44, 1.02674, 0.333609 45, 1.05599, 0.224457 46, 1.07367, 0.112847 47, 1.07958, 0.0 48, 1.07367, -0.112847 49, 1.05599, -0.224457 50, 1.02674, -0.333609 51, 0.986245, -0.439105 52, 0.934944, -0.53979 53, 0.873399, -0.634561 54, 0.802284, -0.72238 55, 0.72238, -0.802284 56, 0.634561, -0.873399 57, 0.53979, -0.934944 58, 0.439105, -0.986245 59, 0.333609, -1.02674 60, 0.224457, -1.05599 61, 0.112847, -1.07367 *NODE,NSET=OUTER 64, 0.121164, 1.1528 65, 0.241001, 1.13382 66, 0.358197, 1.10242 67, 0.471469, 1.05894 68, 0.579575, 1.00385 69, 0.681331, 0.937772 70, 0.775623, 0.861416 71, 0.861416, 0.775623 72, 0.937772, 0.681331 73, 1.00385, 0.579575 74, 1.05894, 0.471469 75, 1.10242, 0.358197 76, 1.13382, 0.241001 77, 1.1528, 0.121164 78, 1.15915, 0.0 79, 1.1528, -0.121164 80, 1.13382, -0.241001 81, 1.10242, -0.358197 82, 1.05894, -0.471469 83, 1.00385, -0.579575 84, 0.937772, -0.681331 85, 0.861416, -0.775623 86, 0.775623, -0.861416 87, 0.681331, -0.937772 88, 0.579575, -1.00385 89, 0.471469, -1.05894 90, 0.358197, -1.10242 91, 0.241001, -1.13382 92, 0.121164, -1.1528 *NSET, NSET=OUT 2, 16, 29, 34, 47, 60, 65, 78, 91, 79 *NCOPY,OLD SET=INNER,CHANGE NUMBER=100,NEW SET=INNER,SHIFT 0.,0.,0. 0.,0.,0., 0.,1.,0., 11.0 *NCOPY,OLD SET=MID,CHANGE NUMBER=100,SHIFT 0.,0.,0. 0.,0.,0., 0.,1.,0., 11.0 *NCOPY,OLD SET=OUTER,CHANGE NUMBER=100,SHIFT 0.,0.,0. 0.,0.,0., 0.,1.,0., 11.0 *ELEMENT, TYPE=AC3D8R,ELSET=LAYER_1 3, 3,5,105,103, 34, 36, 136, 134 *ELGEN, ELSET=LAYER_1 3,13,2,1 *ELEMENT, TYPE=AC3D8R,ELSET=LAYER_2 103, 34, 36, 136, 134, 65, 67, 167, 165 *ELGEN, ELSET=LAYER_2 103,13,2,1 *ELSET, ELSET=EALL LAYER_1, LAYER_2 *NCOPY, CHANGE NUMBER=1000, OLD SET=INNER, SHIFT, NEW SET=SHELL 0.,0.,0., 0., 0.0, 0., 0., 1.0, 0., 0.0 *ELEMENT, TYPE=S4R, ELSET=SHELL 2001, 1003, 1005, 1105, 1103 *ELGEN,ELSET=SHELL 2001, 13, 2, 1 *SURFACE,TYPE=ELEMENT,NAME=SECONDARY LAYER_1, S1 *SURFACE,TYPE=ELEMENT,NAME=MAIN SHELL,SPOS *TIE,NAME=COUPLING SECONDARY, MAIN ** ** FLUID MATERIAL PROPERTIES; ARTIFICIAL MEDIUM FOR CONVENIENCE ** USING RHO = 1, BULKMOD = (2*PI)^2, WE HAVE C = 2*PI AND K = F. ** *MATERIAL,NAME=AETHER *DENSITY 1.000000, *ACOUSTIC MEDIUM,BULK MODULUS 39.4784176, 2 *SOLID SECTION,MATERIAL=AETHER,ELSET=EALL ** *SHELL SECTION,MATERIAL=MAT1, ELSET=SHELL 0.01 *MATERIAL,NAME=MAT1 *ELASTIC 207E0,0.3 *DENSITY 7800, *NSET,NSET=RECRD 3, *TRANSFORM,NSET=SHELL,TYPE=S 0.,0.,0., 0.,0.,1. *AMPLITUDE,NAME=COSN,DEF=PERIODIC 1,,,0.0 1,0,0,0 *STEP *DYNAMIC,EXPLICIT , *IMPEDANCE,nonreflecting=spherical LAYER_2, I2,1.15915 ** *BOUNDARY,TYPE=VELOCITY,AMP=COSN SHELL, 1,1,0.159154943 *BOUNDARY SHELL, 2,3,0.0 SHELL, 4,6,0.0 *OUTPUT,FIELD,NUMBER=5 *NODE OUTPUT,NSET=SHELL V, *NODE OUTPUT POR *OUTPUT,HISTORY,FREQ=1 *NODE OUTPUT, NSET=OUT POR ** ** Output requests for qa testing ** *NSET,NSET=QA_TEST1 SHELL, *NSET,NSET=QA_TEST2,ELSET=EALL *OUTPUT,FIELD,NUM=1 *NODE OUTPUT,NSET=QA_TEST1 U, *NODE OUTPUT,NSET=QA_TEST2 POR, *OUTPUT,HISTORY,VAR=PRESELECT,FREQ=1 *END STEP