*HEADING
 TENNIS RACKET AND BALL IMPACT
RIGID BODY VERIFICATION : tennis_rig_gcont.inp
*PREPRINT,MODEL=YES,HISTORY=YES,CONTACT=YES
********************
**
**  Strings.
**
********************
*NODE
**  Bottom curve.
 104, -2.700,-6.625,0.
 109,  0.000,-8.500,0.
 114,  2.700,-6.625,0.
**  Left side curve.
 501, -4.020,-3.750,0.
1201, -4.520, 0.000,0.
1901, -3.780, 3.750,0.
**  Right side curve.
 517,  4.020,-3.750,0.
1217,  4.520, 0.000,0.
1917,  3.780, 3.750,0.
**  Top curve.
2103, -3.240, 4.500,0.
2109,  0.000, 5.500,0.
2115,  3.240, 4.500,0.
**  Rectangular portion.
 202, -3.780,-6.625,0.
 216,  3.780,-6.625,0.
 502, -3.780,-3.750,0.
 516,  3.780,-3.750,0.
1902, -3.780, 3.750,0.
1916,  3.780, 3.750,0.
2002, -3.780, 4.500,0.
2016,  3.780, 4.500,0.
*NGEN,NSET=N200
202,216,1
*NGEN,NSET=N500
502,516,1
*NGEN,NSET=N1900
1902,1916,1
*NGEN
2002,2016,1
*NFILL,NSET=STRINGS
N200,N500,3,100
N500,N1900,14,100
*NGEN,NSET=STRINGS,LINE=P
104,114,1,109
517,1917,100,1217
501,1901,100,1201,0.,0.,1.
2103,2115,1,2109,0.,0.,1.
*NSET,NSET=MIDPLANE,GEN
105,113,1
2104,2114,1
501,1801,100
517,1817,100
*NSET,NSET=MIDPLANE
204,303,402,214,315,416
1902,2003,1916,2015
**
*ELEMENT,TYPE=T3D2,ELSET=STRINGS
   1, 204, 205
  11, 303, 304
  23, 402, 403
 247,2003,2004
1001, 402, 502
1016, 303, 403
1033, 204, 304
1231, 315, 415
1248, 416, 516
2001, 501, 502
2015, 516, 517
2029, 105, 205
2038,2004,2104
*ELGEN,ELSET=STRINGS
   1, 10,1,1
  11, 12,1,1
  23, 14,1,1, 16,100,14
 247, 12,1,1
1001, 15,100,1
1016, 17,100,1
1033, 18,100,1, 11,1,18
1231, 17,100,1
1248, 15,100,1
2001, 14,100,1
2015, 14,100,1
2029,  9,1,1
2038, 11,1,1 
*SOLID SECTION,ELSET=STRINGS,MATERIAL=STRING
**  Diameter of strings is 0.1 in.
7.854E-3, 
**
*NSET, NSET=NCONTACT, GENERATE
205,213,1
304,314,1
403,415,1
502,516,1
602,616,1
702,716,1
802,816,1
902,916,1
1002,1016,1
1102,1116,1
1202,1216,1
1302,1316,1
1402,1416,1
1502,1516,1
1602,1616,1
1702,1716,1
1802,1816,1
1903,1915,1
2004,2014,1
********************
**
**  Frame.
**
********************
*NCOPY,CHANGE NUMBER=10000,OLD SET=MIDPLANE,SHIFT,NEW SET=FRONT
0.,0.,.25
0.,0.,0., 0.,0.,1., 0.0
*NCOPY,CHANGE NUMBER=20000,OLD SET=MIDPLANE,SHIFT,NEW SET=BACK
0.,0.,-.25
0.,0.,0., 0.,0.,1., 0.0
*NODE,NSET=REFNODE
10000, 0,0.,0.
*NODE,NSET=HANDLE,INPUT=tennis_rig1.inp
**
**
*ELEMENT, TYPE=B31, ELSET=BEAM
   50061,   50051,   50115
   50062,   50050,   50116
   50063,   50049,   50117
   50064,   50115,   50010
   50065,   50116,   50023
   50066,   50117,   50036
*ELEMENT, TYPE=C3D8R, ELSET=BRICKS
   50049,   50079,   50080,   50083,   50082,   50088,   50089,
   50092,   50091
   50050,   50080,   50081,   50084,   50083,   50089,   50090,
   50093,   50092
   50051,   50082,   50083,   50086,   50085,   50091,   50092,
   50095,   50094
   50052,   50083,   50084,   50087,   50086,   50092,   50093,
   50096,   50095
   50053,   50088,   50089,   50092,   50091,   50097,   50098,
   50101,   50100
   50054,   50089,   50090,   50093,   50092,   50098,   50099,
   50102,   50101
   50055,   50091,   50092,   50095,   50094,   50100,   50101,
   50104,   50103
   50056,   50092,   50093,   50096,   50095,   50101,   50102,
   50105,   50104
   50057,   50097,   50098,   50101,   50100,   50106,   50107,
   50110,   50109
   50058,   50098,   50099,   50102,   50101,   50107,   50108,
   50111,   50110
   50059,   50100,   50101,   50104,   50103,   50109,   50110,
   50113,   50112
   50060,   50101,   50102,   50105,   50104,   50110,   50111,
   50114,   50113
   50067,   50044,   50043,   50039,   50026,   50080,   50081,
   50084,   50083
   50068,   50045,   50044,   50026,   50013,   50079,   50080,
   50083,   50082
   50069,   50013,   50026,   50025,   50012,   50082,   50083,
   50086,   50085
   50070,   50026,   50039,   50038,   50025,   50083,   50084,
   50087,   50086
*ELEMENT,TYPE=S4R,ELSET=FRAME, INPUT=tennis_rig2.inp
**
*ELGEN,ELSET=FRAME
3001,  8,1,1
3013, 13,100,1
3039, 10,1,1
3052, 13,100,1
4001,  8,1,1
4013, 13,100,1
4039, 10,1,1
4052, 13,100,1
*SHELL SECTION, MATERIAL=MAT1, ELSET=FRAME
0.08, 
*SOLID SECTION, MATERIAL=MAT1, ELSET=BRICKS
*BEAM SECTION, MATERIAL=MAT1, ELSET=BEAM, SECTION=CIRC
0.08, 
*ELSET, ELSET=RACKET
FRAME, BRICKS, BEAM
*MATERIAL,NAME=MAT1
*DENSITY
1.07E-3, 
*ELASTIC
1.E6, 
*NODE,NSET=SPHERE,INPUT=tennis_rig3.inp
**
*ELEMENT,TYPE=S4R,ELSET=SPHERE,INPUT=tennis_rig4.inp
*SHELL SECTION,ELSET=SPHERE,MATERIAL=RUBBER
.2,3
**
*NSET,NSET=CAVREF
 30001,
*SURFACE,TYPE=ELEMENT,NAME=CAVSURF
 SPHERE, SNEG
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16, UNIVERSAL GAS CONSTANT=73.562
*FLUID CAVITY, NAME=BALL, REF NODE=CAVREF, SURFACE=CAVSURF, BEHAVIOR=GAS, AMBIENT PRESSURE=14.7

*FLUID BEHAVIOR, NAME=GAS
*MOLECULAR WEIGHT
 0.1367
*CAPACITY, TYPE=TABULAR
 112.847
********************
**
** Material definitions.
**
********************
*MATERIAL,NAME=RUBBER
*DENSITY
1.E-4, 
*HYPERELASTIC,N=1
100.,25.,1.E-4
*MATERIAL,NAME=STRING
**  Nylon type 6 general purpose.
*DENSITY
1.07E-4, 
*ELASTIC
1.E6, 
********************
**
** Initial conditions, boundary conditions.
**
********************
*INITIAL CONDITIONS,TYPE=STRESS
**  Tension in strings is 10 lb.
STRINGS,1273.
SPHERE,22.5,22.5
*INITIAL CONDITIONS,TYPE=VELOCITY
**  Initial velocity of ball is 22 fps at a 15 degree angle of attack.
**SPHERE, 1, 68.3
**SPHERE, 2, 0.
**SPHERE, 3,-255.
10000, 1, -68.3
10000, 3,255.
*INITIAL CONDITIONS,TYPE=FLUID PRESSURE
30001,6.
**
*RESTART,WRITE,NUM=1
************************************************************
**
*SURFACE,TYPE=ELEMENT,NAME=SPHERE
SPHERE,SPOS
*SURFACE,TYPE=NODE,NAME=STRINGS
 NCONTACT,
*SURFACE INTERACTION,NAME=SPH_STR
*FRICTION
 0.1,
*RIGID BODY,ELSET=RACKET,REF NODE=10000
**
*STEP
*DYNAMIC,EXPLICIT
,15.E-3
********************
**
** Contact definitions.
**
********************
*CONTACT
*CONTACT INCLUSIONS
 SPHERE,STRINGS
*CONTACT PROPERTY ASSIGNMENT
 SPHERE,STRINGS,SPH_STR
********************
**
** Output requests.
**
********************
*FILE OUTPUT, NUM=4, TIMEMARKS=YES
*NODE FILE,NSET=CAVREF
CVOL,PCAV,
*NODE FILE,NSET=REFNODE
 U,
*OUTPUT,FIELD, NUMBER INTERVAL=20, VARIABLE=PRESELECT
*ELSET, ELSET=QA_TEST
RACKET, STRINGS, SPHERE
*OUTPUT,HISTORY, TIME INTERVAL=75.E-6, VARIABLE=PRESELECT
*INTEGRATED OUTPUT, ELSET=QA_TEST
COORDCOM, MASS
*NSET, NSET=QA_TEST
30016,30021,30022,30023,30028, 10000,
*OUTPUT, FIELD, NUMBER INTERVAL=3
*NODE OUTPUT, NSET=QA_TEST
U,
*END STEP