*HEADING
Simple loudspeaker analysys
**
** Node defintion for Air
**
*PREPRINT,MODEL=YES,HISTORY=YES
*NODE,NSET=FRONT
  1,-0.2500,-.2000,-.3000
  2,-0.1250,-.2000,-.3000
  3,-0.0000,-.2000,-.3000
  4, 0.1250,-.2000,-.3000
  5, 0.2500,-.2000,-.3000
 41,-0.2500,-.2000,-.1500
 42,-0.1218,-.2000,-.1218
 43, 0.0000,-.2000,-.1725
 44, 0.1218,-.2000,-.1218
 45, 0.2500,-.2000,-.1500
 81,-0.2500,-.2000, .0000
 82,-0.1725,-.2000, .0000
 83,-0.0000,-.2000, .0000
 84, 0.1725,-.2000, .0000
 85, 0.2500,-.2000, .0000
121,-0.2500,-.2000, .1500
122,-0.1218,-.2000, .1218
123, 0.0000,-.2000, .1725
124, 0.1218,-.2000, .1218
125, 0.2500,-.2000, .1500
161,-0.2500,-.2000, .3000
162,-0.1250,-.2000, .3000
163,-0.0000,-.2000, .3000
164, 0.1250,-.2000, .3000
165, 0.2500,-.2000, .3000
*NODE,NSET=REAR
 31,-0.2500, .2000,-.3000
 32,-0.1250, .2000,-.3000
 33,-0.0000, .2000,-.3000
 34, 0.1250, .2000,-.3000
 35, 0.2500, .2000,-.3000
 71,-0.2500, .2000,-.1500
 72,-0.1218, .2000,-.1218
 73, 0.0000, .2000,-.1725
 74, 0.1218, .2000,-.1218
 75, 0.2500, .2000,-.1500
111,-0.2500, .2000, .0000
112,-0.1725, .2000, .0000
113,-0.0000, .2000, .0000
114, 0.1725, .2000, .0000
115, 0.2500, .2000, .0000
151,-0.2500, .2000, .1500
152,-0.1218, .2000, .1218
153, 0.0000, .2000, .1725
154, 0.1218, .2000, .1218
155, 0.2500, .2000, .1500
191,-0.2500, .2000, .3000
192,-0.1250, .2000, .3000
193,-0.0000, .2000, .3000
194, 0.1250, .2000, .3000
195, 0.2500, .2000, .3000
**
*NFILL,NSET=AIR
FRONT,REAR,3,10
**
*NODE 
83,-0.0000,-.1600, .0000
**
** ELEMENT DEFINTION FOR AIR
**
*ELEMENT,TYPE=AC3D8
1,1,2,12,11,41,42,52,51
31,41,42,52,51,81,82,92,91
36,44,45,55,54,84,85,95,94
91,121,122,132,131,161,162,172,171
*ELGEN,ELSET=AIR 
1,4,1,1,3,10,10,1,1,1
31,1,1,1,3,10,10,2,40,30,1,1,1
36,1,1,1,3,10,10,2,40,30,1,1,1
91,4,1,1,3,10,10,1,1,1
**
*ELEMENT,TYPE=AC3D6,ELSET=AIR
32,92,52,93,82,42,83
33,52,53,93,42,43,83
34,53,54,93,43,44,83
35,54,94,93,44,84,83
62,132,92,93,122,82,83
63,133,132,93,123,122,83
64,134,133,93,124,123,83
65,94,134,93,84,124,83
**
42,102,62,103,92,52,93
43,62,63,103,52,53,93
44,63,64,103,53,54,93
45,64,104,103,54,94,93
72,142,102,103,132,92,93
73,143,142,103,133,132,93
74,144,143,103,134,133,93
75,104,144,103,94,134,93
**
52,112,72,113,102,62,103
53,72,73,113,62,63,103
54,73,74,113,63,64,103
55,74,114,113,64,104,103
82,152,112,113,142,102,103
83,153,152,113,143,142,103
84,154,153,113,144,143,103
85,114,154,113,104,144,103
**
*SOLID SECTION,ELSET=AIR,MAT=ACOUSTIC
*MATERIAL,NAME=ACOUSTIC
*DENSITY
1.115
*ACOUSTIC MEDIUM,BULK MODULUS
1.3194E5
**
** SURFACE DEFINITIONS
**
*ELSET,ELSET=AIR_FRONT
1,2,3,4
31,36,61,66
91,92,93,94
*SURFACE,NAME=AIR_FRONT
AIR_FRONT,S3
*ELSET,ELSET=AIR_REAR_PLANE,GEN
21,24,1
51,56,5
81,86,5
111,114,1
*ELSET,ELSET=AIR_REAR_CONE,GEN
52,55,1,
82,85,1
*SURFACE, NAME=AIR_REAR
AIR_REAR_PLANE,S5
AIR_REAR_CONE,S1
*ELSET,ELSET=AIR_LEFT,GEN
1,111,10
*SURFACE, NAME=AIR_LEFT
AIR_LEFT,S6
*ELSET,ELSET=AIR_RIGHT,GEN
4,24,10
36,56,10
66,86,10
94,114,10
*SURFACE, NAME=AIR_RIGHT
AIR_RIGHT, S4
*ELSET,ELSET=AIR_TOP,GEN
91,94,1
101,104,1
111,114,1
*SURFACE, NAME=AIR_TOP
AIR_TOP,S2
*ELSET,ELSET=AIR_BOT,GEN
1,4,1
11,14,1
21,24,1
*SURFACE, NAME=AIR_BOT
AIR_BOT,S1
*ELSET,ELSET=AIR_CONE,GEN
32,35,1
62,65,1
*SURFACE, NAME=AIR_CONE
AIR_CONE,S2
**
**
** NODE DEFINITION FOR BOX AND CONE
**
*NCOPY,OLDSET=AIR,CHANGE NUMBER=1000,SHIFT,NEW SET=BOX
0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,1.0,0.0
**
** ELEMENT DEFINITION FOR CONE
**
*ELEMENT,TYPE=S3R,ELSET=CONE
1032,1082,1042,1083
1033,1042,1043,1083
1034,1043,1044,1083
1035,1044,1084,1083
1062,1122,1082,1083
1063,1123,1122,1083
1064,1124,1123,1083
1065,1084,1124,1083
**
*SHELL SECTION,ELSET=CONE,MAT=POLYETHYLENE 
0.1E-3
*MATERIAL,NAME=POLYETHYLENE
*ELASTIC
3.400E9,0.3
*DENSITY
450
**
** ELEMENT DEFINITION FOR BOX
** 
*ELEMENT,TYPE=S4R,ELSET=FRONT
1001,1001,1002,1042,1041
1002,1002,1003,1043,1042
1003,1003,1004,1044,1043
1004,1004,1005,1045,1044
1031,1041,1042,1082,1081
1036,1044,1045,1085,1084
1061,1081,1082,1122,1121
1066,1084,1085,1125,1124
1091,1121,1122,1162,1161
1092,1122,1123,1163,1162
1093,1123,1124,1164,1163
1094,1124,1125,1165,1164
**
*ELCOPY,ELEMENT SHIFT=20,OLDSET=FRONT,SHIFTNODES=30,NEWSET=REAR 
*ELCOPY,ELEMENT SHIFT=20,OLDSET=CONE,SHIFTNODES=30,NEWSET=REAR
**
*NODE
1113,0.0,0.2,0.0
**
*ELEMENT,TYPE=S4R
2001,1001,1041,1051,1011
*ELGEN,ELSET=LEFT
2001,4,40,10,3,10,1
*ELEMENT,TYPE=S4R
2101,1005,1015,1055,1045
*ELGEN,ELSET=RIGHT
2101,4,40,10,3,10,1
*ELEMENT,TYPE=S4R
3001,1161,1162,1172,1171
*ELGEN,ELSET=TOP
3001,4,1,1,3,10,10
*ELEMENT,TYPE=S4R
3101,1001,1011,1012,1002
*ELGEN,ELSET=BOT
3101,4,1,1,3,10,10
*ELSET,ELSET=BOX
FRONT,REAR,LEFT,RIGHT,TOP,BOT
*SHELL SECTION,ELSET=BOX,MAT=WOOD
0.015
*MATERIAL,NAME=WOOD
*ELASTIC
11.61E9,0.3
*DENSITY
562
**
**
** SURFACE DEFINITION FOR THE BOX AND CONE
**
*SURFACE, NAME=BOX_FRONT
FRONT,SPOS
*SURFACE, NAME=BOX_REAR
REAR,SNEG
*SURFACE, NAME=BOX_RIGHT
RIGHT,SPOS
*SURFACE, NAME=BOX_LEFT
LEFT,SPOS
*SURFACE, NAME=BOX_TOP
TOP,SPOS
*SURFACE, NAME=BOX_BOT
BOT,SPOS
*SURFACE, NAME=CONE
CONE,SPOS
**
** CONTACT DEFINITION
**
*TIE, NAME=COUPLED
AIR_FRONT,BOX_FRONT
AIR_REAR,BOX_REAR
AIR_LEFT,BOX_LEFT
AIR_RIGHT,BOX_RIGHT
AIR_TOP,BOX_TOP
AIR_BOT,BOX_BOT
*TIE, NAME=COUPLEDCONE, TYPE=SURFACE TO SURFACE
AIR_CONE,CONE
**
** BOUNDARY CONDITIONS
**
*NSET,NSET=CORNERS
1001,1005,1031,1035
*NSET,NSET=MID
1012,1013,1014,1022,1023,1024
*NSET,NSET=CONECENTER
1083
*NSET,NSET=AIR-CONECENTER
83
*BOUNDARY
MID,1,3,0
**
*STEP,NAME=COUPLED-FREQUENCY
Default coupled acoustic-structural frequency extraction
*FREQUENCY,EIGENSOLVER=LANCZOS,norm=mass
150,,,-1
*OUTPUT,FIELD,freq=0
*END STEP
**------------------------------------------------
*STEP,NAME=SSD-DIRECT
*STEADY STATE DYNAMICS, DIRECT,FREQUENCY SCALE=LINEAR,INTERVAL=RANGE
300,400,50,1
*CLOAD
1083,2,1.0
*OUTPUT,FIELD,FREQ=0
*OUTPUT, HISTORY
*NODE OUTPUT,NSET=MID
RF,
*END STEP
**------------------------------------------------
***STEP,NAME=SSD-SUBSPACE
**Steadystatedynamics based on coupled eigenmodes
***STEADY STATE DYNAMICS,SUBSPACE PROJECTION, FREQUENCY SCALE=LINEAR,INTERVAL=RANGE
**300,400,100,1
***CLOAD
**1083,2,1.
***CLOAD, LOAD CASE=2
**1083,2,0.06
***OUTPUT, FIELD
***NODE OUTPUT
**U,POR
***OUTPUT, HISTORY
***NODE OUTPUT,NSET=CONECENTER
**U,
***NODE OUTPUT,NSET=AIR-CONECENTER
**POR,
***END STEP
***STEP,NAME=SSD-DIRECT
**Steadystatedynamics direct
***STEADY STATE DYNAMICS,DIRECT, FREQUENCY SCALE=LINEAR,INTERVAL=RANGE
**300,400,100,1
***CLOAD
**1083,2,1.
***CLOAD, LOAD CASE=2
**1083,2,0.06
***END STEP
***STEP,NAME=UNCOUPLED-FREQUENCY
**Uncoupled acoustic-structural frequency extraction
***FREQUENCY,EIGENSOLVER=LANCZOS,acoustic coupling=off
**15,,,,7,15
***OUTPUT,FIELD
***NODE OUTPUT
**U,POR
***END STEP
***STEP,NAME=SSD-SUBSPACE-UNCOUPLED
**Steadystatedynamics based on uncoupled eigenmodes
***STEADY STATE DYNAMICS,SUBSPACE PROJECTION, FREQUENCY SCALE=LINEAR,INTERVAL=RANGE
**300,400,100,1
***CLOAD
**1083,2,1.
***CLOAD, LOAD CASE=2
**1083,2,0.06
***OUTPUT, FIELD
***NODE OUTPUT
**U,POR
***OUTPUT, HISTORY
***NODE OUTPUT,NSET=CONECENTER
**U,
***NODE OUTPUT,NSET=AIR-CONECENTER
**POR,
***END STEP