$ vmchk1 DIAG 8,15 SOL 103 Include 'SSSALTERDIR:vmchka.alt' $ checking alter CEND TITLE = VIRTUAL MASS OF A BOX, 2 elements per edge vmchk1 subtitle = alter for vm wetted element checking inserted, label = variation with unit cube, internal fluid, echo = unsort $ METHO = 400 $ reduce cost of run Subcase 1 Super=1 mfluid = 300 Subcase 11 Super=0 Mfluid = 400 BEGIN BULK $ ADD A SPHERICAL COORDINATE SYSTEM WHOSE ORIGIN IS AT A POINT WHERE $ ALL WETTED SURFACES CAN BE SEEN. IF THERE IS ALREADY ONE AVAILABLE, $ USE IT. PARAM VMSPHER 1 CORD2S 1 1. +C1 +C1 1. DAREA 100 10101 3 1.0 FREQ 3 0.0001 .043 .086 .159 .2079 100.0 RLOAD2 2 100 200 TABDMP1 4 , 1.0 0.1 2.0 0.1 ENDT TABLED1 200 , 10. 1.0 100. 1.0 ENDT $EGRID 1 -.5 -.5 -.5 $EGRID 2 +.5 -.5 -.5 $EGRID 3 +.5 +.5 -.5 $EGRID 4 -.5 +.5 -.5 $EGRID 5 -.5 -.5 +.5 $EGRID 6 +.5 -.5 +.5 $EGRID 7 +.5 +.5 +.5 $EGRID 8 -.5 +.5 +.5 $GRIDG 1 2 -1 -4 -8 +G1 $+G1 2 -5 GRID 10000 -.5 -.5 -.5 GRID 10001 -.5 0.0 -.5 GRID 10002 -.5 .5 -.5 GRID 10100 -.5 -.5 0.0 GRID 10101 -.5 0.0 0.0 GRID 10102 -.5 .5 0.0 GRID 10200 -.5 -.5 .5 GRID 10201 -.5 0.0 .5 GRID 10202 -.5 .5 .5 $GRIDG 2 2 -2 -6 -7 +G2 $+G2 2 -3 GRID 20000 .5 -.5 -.5 GRID 20001 .5 -.5 0.0 GRID 20002 .5 -.5 .5 GRID 20100 .5 0.0 -.5 GRID 20101 .5 0.0 0.0 GRID 20102 .5 0.0 .5 GRID 20200 .5 .5 -.5 GRID 20201 .5 .5 0.0 GRID 20202 .5 .5 .5 $GRIDG 3 2 -1 -5 -6 +G3 $+G3 2 -2 GRID 30100 0.0 -.5 -.5 GRID 30101 0.0 -.5 0.0 GRID 30102 0.0 -.5 .5 $GRIDG 4 2 -3 -7 -8 +G4 $+G4 2 -4 GRID 40100 0.0 .5 -.5 GRID 40101 0.0 .5 0.0 GRID 40102 0.0 .5 .5 $GRIDG 5 2 -1 -2 -3 +G5 $+G5 2 -4 GRID 50101 0.0 0.0 -.5 $EQUIV AUTO $CGEN QUAD4 1001 100 1 CQUAD4 1001 100 10000 10001 10101 10100 CQUAD4 1002 100 10001 10002 10102 10101 CQUAD4 1003 100 10100 10101 10201 10200 CQUAD4 1004 100 10101 10102 10202 10201 $CGEN QUAD4 2001 100 2 CQUAD4 2001 100 20000 20001 20101 20100 CQUAD4 2002 100 20001 20002 20102 20101 CQUAD4 2003 100 20100 20101 20201 20200 CQUAD4 2004 100 20101 20102 20202 20201 $CGEN QUAD4 3001 100 3 CQUAD4 3001 100 10000 10100 30101 30100 CQUAD4 3002 100 10100 10200 30102 30101 CQUAD4 3003 100 30100 30101 20001 20000 CQUAD4 3004 100 30101 30102 20002 20001 $CGEN QUAD4 4001 100 4 CQUAD4 4001 100 20200 20201 40101 40100 CQUAD4 4002 100 20201 20202 40102 40101 CQUAD4 4003 100 40100 40101 10102 10002 CQUAD4 4004 100 40101 40102 10202 10102 $CGEN QUAD4 5001 100 5 CQUAD4 5001 100 10000 30100 50101 10001 $CQUAD4 5002 100 30100 20000 20100 50101 $ reverse the order of connected nodes to get error on wetted side $ element 5002 is now in error for VM CQUAD4 5002 100 30100 50101 20100 20000 CQUAD4 5003 100 10001 50101 40100 10002 CQUAD4 5004 100 50101 20100 20200 40100 PSHELL 100 100 .001 100 100 MAT1 100 1.+6 .3 10. $SPCG 200 5 123456 A C MFLUID 400 0 +.5 1. 350 N N ELIST 350 1001 THRU 1004 $ simulate some sign errors $ELIST 350 2001 THRU 2004 $ elements 2001 through 2004 are now in error for vm ELIST, 350 -2001 THRU -2004 ELIST 350 3001 THRU 3004 ELIST 350 4001 THRU 4004 $ELIST 350 5001 THRU 5004 $ make a hole on bottom surface, 5003 ELIST, 350 5001 THRU 5002, 5004 EIGRL 400 20 PARAM PRGPST NO PARAM GRDPNT 0 $ boundary points for se GRID, 100000, 2, -.5 -.5 -.5 GRID, 100010, 2, -.5 .5 -.5 GRID, 200000, 2, .5 -.5 -.5 GRID, 201000, 2, .5 .5 -.5 Cord2r, 2, 0, 10. 0. 0. 10. 0. 1. , 11. 0. 0. $ displace it laterally from other model begin super=1 $ fluid in part superelement $ the VM details of this model are correct. PARAM, VMSPHER 2 $ put in center of wetted elements CORD2S 2 0 10. 0. 0. 10. 0. 1. , 11. 0. 0. CQUAD4 1001 100 10000 10001 10101 10100 CQUAD4 2001 100 20000 20001 20101 20100 CQUAD4 3001 100 10000 10100 20001 20000 CQUAD4 4001 100 20100 20101 10101 10001 CQUAD4 5001 100 10000 20000 20100 10001 EIGRL 400 10 ELIST 350 1001 2001 3001 4001 5001 Cord2r, 1, 0, 10. 0. 0. 10. 0. 1. , 11. 0. 0. $ displace it laterally from other model GRID, 10000, 1, -.5 -.5 -.5 GRID, 10001, 1, -.5 .5 -.5 GRID, 10100, 1, -.5 -.5 .5 GRID, 10101, 1, -.5 .5 .5 GRID, 20000, 1, .5 -.5 -.5 GRID, 20001, 1, .5 -.5 .5 GRID, 20100, 1, .5 .5 -.5 GRID, 20101, 1, .5 .5 .5 MAT1 100 1.+6 .3 10. MFLUID 300 0 +.5 1. 350 N N PSHELL 100 100 .001 100 100 ENDDATA