*HEADING EXAMPLE PROBLEM: MODE 0 SPHERE USING BGT BOUNDARY CONDITION AND ACAX4R. *************************** *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=INNER 1, 0., 1. 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 31, 0., -1. *NODE,NSET=MID 32, 0., 1.07958 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 62, 0., -1.07958 *NODE,NSET=OUTER 63, 0., 1.15915 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 93, 0., -1.15915 *NSET, NSET=OUT 1, 8, 15, 16, 17, 30, 31 32, 39, 47, 54, 62 63, 70, 78, 85, 93 *NCOPY, CHANGE NUMBER=100, OLD SET=INNER, SHIFT, NEW SET=BEAMS 0.,0.,0., 0., 0.0, 0., 0., 1.0, 0., 0.0 *ELEMENT, TYPE=ACAX4R, ELSET=LAYER_1 1, 1, 2, 33, 32 *ELEMENT, TYPE=ACAX4R, ELSET=LAYER_2 100, 32, 33, 64, 63 *ELGEN, ELSET=LAYER_1 1,30,1,1 *ELGEN, ELSET=LAYER_2 100,30,1,1 *ELSET, ELSET=EALL LAYER_1, LAYER_2 *ELEMENT, TYPE=B21,ELSET=BEAMS 2001, 101, 102 *ELGEN,ELSET=BEAMS 2001, 30, 1, 1 *SURFACE,TYPE=ELEMENT,NAME=SLAVE LAYER_1, S1 *SURFACE,TYPE=ELEMENT,NAME=MASTER BEAMS,SPOS *TIE,NAME=COUPLING SLAVE, MASTER *********************************** *ELEMENT, TYPE=ACINAX2 10000,63,64 *ELGEN, ELSET=INF 10000,30,1,1 *SOLID SECTION, ELSET=INF, MATERIAL=AETHER 1.,0.,0.,0. *********************************** ** ** 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 ** *BEAM SECTION,MATERIAL=MAT1,ELSET=BEAMS,SECTION=RECT 0.1,0.1 *MATERIAL,NAME=MAT1 *ELASTIC 207E0,0.3 *DENSITY 7800, *NSET,NSET=RECRD 1, *TRANSFORM,NSET=BEAMS,TYPE=C 0.,0.,0., 0.,0.,1. *AMPLITUDE,NAME=COSN,DEF=PERIODIC 1,,,0.0 1,0,0,0 *STEP *DYNAMIC,EXPLICIT , ** *BOUNDARY,TYPE=VELOCITY,AMP=COSN BEAMS, 1,1,0.159154943 *BOUNDARY BEAMS, 2,2,0.0 BEAMS, 6,6,0.0 *OUTPUT,FIELD,NUMBER=5 *NODE OUTPUT,NSET=BEAMS V, *NODE OUTPUT POR *OUTPUT,HISTORY,FREQ=1 *NODE OUTPUT, NSET=OUT POR ** ** Output requests for qa testing ** *NSET,NSET=QA_TEST1 BEAMS, *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