*HEADING SUBMERGED CYLINDER UNDEX (KWON & FOX, 1993) ** ** ** ********************** ** *** DEFINE PARTS *** ** ********************** ** *PART, NAME=WATER *INCLUDE, INPUT = submerged_cyl_water.inp *END PART ** *PART, NAME=CYLINDER *INCLUDE, INPUT = submerged_cyl_cylinder.inp *END PART ** *PART, NAME=PULSE *NODE, NSET=PT-MASS 1, 0., 0., 0. *ELEMENT, TYPE=MASS, ELSET=PT-MASS 1, 1 *MASS, ELSET=PT-MASS 1.0E-7, *END PART ** ** ** ****************************** ** *** DEFINE ASSEMBLY DATA *** ** ****************************** ** *ASSEMBLY, NAME=UNDEX ** *INSTANCE, NAME=WATER-1, PART=WATER 0., 0., 0. 0., 0., 0., 1., 0., 0., 90. *END INSTANCE ** *INSTANCE, NAME=CYLINDER-1, PART=CYLINDER 0., 0., -0.5335 0., 0., -0.5335, 0., 0., 0.4665, 4.5 *END INSTANCE ** *INSTANCE, NAME=SHOCK-PULSE, PART=PULSE 0., 0., 0. *END INSTANCE ** ** ** **************************** ** ** SURFACE DEFINITIIONS ** ** ** (ASSEMBLY LEVEL) ** ** **************************** ** *SURFACE, TYPE=ELEMENT, NAME=CYL-FSI, NO THICK CYLINDER-1.ALL-ELEM, SPOS ** *SURFACE, TYPE=ELEMENT, NAME=RAD-MID WATER-1._RAD-MID_S3, S3 WATER-1._RAD-MID_S1, S1 WATER-1._RAD-MID_S2, S2 WATER-1._RAD-MID_S4, S4 *SURFACE, TYPE=ELEMENT, NAME=RAD-ENDL WATER-1._RAD-ENDL_S3, S3 WATER-1._RAD-ENDL_S4, S4 WATER-1._RAD-ENDL_S2, S2 WATER-1._RAD-ENDL_S1, S1 *SURFACE, TYPE=ELEMENT, NAME=RAD-ENDR WATER-1._RAD-ENDR_S3, S3 WATER-1._RAD-ENDR_S4, S4 WATER-1._RAD-ENDR_S2, S2 WATER-1._RAD-ENDR_S1, S1 *SURFACE, TYPE=ELEMENT, NAME=WATER-FSI WATER-1._WATER-FSI_S3, S3 WATER-1._WATER-FSI_S4, S4 WATER-1._WATER-FSI_S2, S2 WATER-1._WATER-FSI_S1, S1 ** ** ** *************************************** ** ** ASSEMBLY LEVEL DEFINITION OF ** ** ** FLUID-STRUCTURE COUPLING VIA *TIE ** ** *************************************** ** *TIE, NAME=FSI, ADJUST=NO, NO ROTATION CYL-FSI, WATER-FSI ** *END ASSEMBLY ** ** ** ********************* ** *** MODEL DATA *** ** ********************* ** *MATERIAL, NAME=ALUMINIUM *DENSITY 2784.5, *ELASTIC 7.56E+10, 0.33 *PLASTIC 3E+08,0. 3.3E+08,1. *DAMPING, ALPHA=750.0 ** ** WATER (C = 1463 METERS/SECOND) ** *MATERIAL, NAME=WATER *ACOUSTIC MEDIUM 2.1404E+09, *DENSITY 1000.0, ** ** ** ********************************** ** *** INCIDENT WAVE DEFINITION *** ** ********************************** ** *INCIDENT WAVE PROPERTY,NAME=SHOCK,TYPE=SPHERE 0.1525, 0., 0., 0., 0., 0. 7.7725, 0., 0. *INCIDENT WAVE FLUID PROPERTY 2.1404e+09, 1000.0 ** ** ** ************************************ ** ** RESTRAIN POINT MASS 1 & 2 DOF ** ** ************************************ ** *BOUNDARY, TYPE=VELOCITY UNDEX.SHOCK-PULSE.PT-MASS, 1, 2, 0.0 ** ** ** ********************************** ** ** EXPLICIT DYNAMIC STEP DATA ** ** ********************************** ** *NSET, NSET=QA_TEST UNDEX.CYLINDER-1.ALL-NODE ********* *STEP, NAME=UNDEX-SHOCK-WAVE INCIDENT WAVE, CYLINDER ALPHA DAMPING = 750 *DYNAMIC, EXPLICIT , 0.008 ** *BULK VISCOSITY 0.06, 1.2 ** *RESTART,WRITE,NUMBERINTERVAL=1,TIME MARKS=NO ** *INCLUDE, INPUT = submerged_cyl_pulse.inp ** ** ** *************************************** ** ** POINT MASS LOADING TO REPRODUCE ** ** ** THE SHOCK PULSE HISTORY CURVE ** ** *************************************** ** *CLOAD, AMPLITUDE=PULSE UNDEX.SHOCK-PULSE.PT-MASS, 3, 1.0e-7 ** ** ** ************************************* ** ** INCIDENT WAVE LOADING AT ** ** ** THE FLUID-STRUCTURE INTERFACE ** ** ************************************* ** *INCIDENT WAVE, PROPERTY=SHOCK, PRESSURE AMPLITUDE=PULSE UNDEX.WATER-FSI, 1. UNDEX.CYL-FSI, 1. ** ** ** ** **************************** ** ** SURFACE IMPEDANCES ** ** ** (RADIATING BOUNDARIES) ** ** **************************** ** *SIMPEDANCE, NON=CIRC UNDEX.RAD-MID, 0.915 *SIMPEDANCE, NON=SPH UNDEX.RAD-ENDL,0.915 UNDEX.RAD-ENDR,0.915 ** ** ** ************************* ** *** OUTPUT REQUESTS *** ** ************************* ** *OUTPUT, FIELD, NUMBERINTERVAL=40 *ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.ALL-ELEM S, LE, PEEQ *NODE OUTPUT, NSET=UNDEX.CYLINDER-1.ALL-NODE U, *NODE OUTPUT, NSET=UNDEX.WATER-1.ALL-NODE POR, *************** *OUTPUT, FIELD, NUM=1 *NODE OUTPUT, NSET=QA_TEST U, ***************** *OUTPUT, HISTORY, TIMEINTERVAL=1.0E-05 *NODE OUTPUT, NSET=UNDEX.SHOCK-PULSE.PT-MASS A3, *NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-R-CTR U,V *NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-L-CTR U,V *NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-B1 U,V *NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-B2 U,V *NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-B3 U,V *NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-B4 U,V *ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.A1 9, LE11, LE22, PEEQ *ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.A2 9, LE11, LE22, PEEQ *ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.B1 9, LE11, LE22, PEEQ *ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.B2 9, LE11, LE22, PEEQ *ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.B3 9, LE11, LE22, PEEQ *ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.C1 9, LE11, LE22, PEEQ *ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.C2 9, LE11, LE22, PEEQ ** *ENERGY OUTPUT,ELSET=UNDEX.CYLINDER-1.ALL-ELEM ALLAE, ALLIE, ALLKE, ALLPD, ALLVD, ALLSE ** *END STEP