diff --git a/ASMu2DNastran.C b/ASMu2DNastran.C index 6f75f00..4a787ba 100644 --- a/ASMu2DNastran.C +++ b/ASMu2DNastran.C @@ -877,6 +877,23 @@ bool ASMu2DNastran::checkPressSet (int iEl, size_t idx, int iSet) const } +bool ASMu2DNastran::getShellNormals (std::vector& normals) const +{ + normals.reserve(nel); + for (size_t iel = 0; iel < nel; iel++) + if (const IntVec& mnpc = MNPC[iel]; mnpc.size() == 3 || mnpc.size() == 4) + { + const size_t n1 = mnpc.size() - 3; + const size_t n3 = mnpc.size() == 4 ? 3 : 1; + Vec3 V3(coord[mnpc[n3]]-coord[mnpc[0]], coord[mnpc[2]]-coord[mnpc[n1]]); + V3.normalize(); + normals.emplace_back(this->getElementCenter(1+iel),V3); + } + + return !normals.empty(); +} + + #ifdef HAS_ANDES extern "C" void wavgmconstreqn_(const int& iel, const int& lDof, const int& nM, const int& nW, const int* indC, diff --git a/ASMu2DNastran.h b/ASMu2DNastran.h index b7ad385..5081e5d 100644 --- a/ASMu2DNastran.h +++ b/ASMu2DNastran.h @@ -66,6 +66,10 @@ class ASMu2DNastran : public ASMu2DLag //! \param[in] iSet Element set index (1-based) for a pressure load bool checkPressSet(int iEl, size_t idx, int iSet) const; + //! \brief Calculates the shell normal vectors and the element centers + //! \param[out] normals Vector of element-center normal-vector pairs + bool getShellNormals(std::vector& normals) const; + //! \brief Returns an additional geometry to visualize (point masses, etc.). virtual ElementBlock* immersedGeometry(char* name) const; //! \brief Returns an additional geometry to visualize (constraints, etc.). diff --git a/CMakeLists.txt b/CMakeLists.txt index 031e147..0abb1b8 100644 --- a/CMakeLists.txt +++ b/CMakeLists.txt @@ -113,6 +113,7 @@ ifem_add_regression_test( TEST_FILES BeamOnly-dyn.reg BeamOnly.reg + DroppedBox.reg Fine-dyn.reg Fine-mlc.reg Fine-partial.reg diff --git a/SIMAndesShell.C b/SIMAndesShell.C index f5c318f..7a0b2fa 100644 --- a/SIMAndesShell.C +++ b/SIMAndesShell.C @@ -564,6 +564,26 @@ bool SIMAndesShell::writeGlvLoc (std::vector& locfiles, } +bool SIMAndesShell::writeGlvNormal (int& geoBlk, int& nBlock) const +{ + VTF* vtf = this->getVTF(); + if (!vtf || myProblem->hasTractionValues()) return true; + + for (const ASMbase* pch : myModel) + if (const ASMu2DNastran* shl = dynamic_cast(pch); shl) + if (std::vector normals; shl->getShellNormals(normals)) + { + if (msgLevel > 1) + IFEM::cout <<"Writing shell normal vectors"<< std::endl; + + // Write shell normals as discrete point vectors to the VTF-file + return vtf->writeVectors(normals,geoBlk,++nBlock,"Normal vectors",1); + } + + return true; // No shell elements +} + + bool SIMAndesShell::writeGlvG (int& nBlock, const char* inpFile, bool doClear) { if (!this->Parent::writeGlvG(nBlock,inpFile,doClear)) diff --git a/SIMAndesShell.h b/SIMAndesShell.h index a14a307..3a0c032 100644 --- a/SIMAndesShell.h +++ b/SIMAndesShell.h @@ -74,6 +74,11 @@ class SIMAndesShell : public SIMElasticity bool writeGlvLoc(std::vector& locfiles, bool nodal, int& nBlock) const; + //! \brief Writes shell normal vectors to the VTF-file. + //! \param geoBlk Running geometry block counter + //! \param nBlock Running result block counter + bool writeGlvNormal(int& geoBlk, int& nBlock) const; + //! \brief Writes current model geometry to the VTF-file. //! \param nBlock Running result block counter //! \param[in] inpFile File name used to construct the VTF-file name from diff --git a/Test/DroppedBox.reg b/Test/DroppedBox.reg new file mode 100644 index 0000000..360e222 --- /dev/null +++ b/Test/DroppedBox.reg @@ -0,0 +1,1022 @@ +DroppedBox.xinp -dynamic -keep-previous-state + +Input file: DroppedBox.xinp +Using the linear dynamics simulation driver. +Parsing input file DroppedBox.xinp +Parsing +Parsing + Parsing + Reading data file Q4box.nas +Parsing Nastran bulk data file succceeded. +Total number of nodes: 8 +Number of shell elements: 6 +Number of constraint elements: 0 +Number of other elements: 0 +Model extension (diameter): 1.73205 + Reading patch 1 + Parsing +Parsing + Parsing + Gravitation vector: 0 0 -9.81 + Parsing + Sea surface function (expression): wavel=4.0; waveh=0.1; wavespeed=-4.0; tstart=1.0; refx=3.0+wavespeed\*(t-tstart); if(above(x+ux,refx),0.5\*waveh\*sin((x+ux-refx)\*2\*PI/wavel),0.0) + Parsing + Pressure code 1 (hydrostatic): z0 = 0.1 rhow = 1000 +Parsing + alpha1 = 1 alpha2 = 0.01 + beta = 0.25 gamma = 0.5 +Parsing + Parsing + Parsing + Parsing + Sea surface \[-3 -0.5 0.1] to \[2 1.5 0.1], gridSize=0.1 +Parsing input file succeeded. +Equation solver: 2 +Number of Gauss points: 2 +Lagrangian basis functions are used +Problem definition: +Formulation: ANDES shell +Newmark predictor/multicorrector: beta = 0.25 gamma = 0.5 + using constant coefficient matrices + using zero acceleration predictor +Mass-proportional damping (alpha1): 1 +Stiffness-proportional damping (alpha2): 0.01 +Resolving Dirichlet boundary conditions + >>> SAM model summary <<< +Number of elements 6 +Number of nodes 8 +Number of dofs 48 +Number of unknowns 48 + step=1 time=0.05 + Displacement L2-norm : 0.00192022 + Max Z-displacement : 0.00470361 + Velocity L2-norm : 0.0768089 + Max Z-velocity : 0.188144 + Acceleration L2-norm : 3.07236 + Max Z-acceleration : 7.52577 + step=2 time=0.1 + Displacement L2-norm : 0.0094829 + Max Z-displacement : 0.0232283 + Velocity L2-norm : 0.225698 + Max Z-velocity : 0.552846 + Acceleration L2-norm : 2.88322 + Max Z-acceleration : 7.06257 + step=3 time=0.15 + Displacement L2-norm : 0.0241873 + Max Z-displacement : 0.0592465 + Velocity L2-norm : 0.362476 + Max Z-velocity : 0.887883 + Acceleration L2-norm : 2.58791 + Max Z-acceleration : 6.33928 + step=4 time=0.2 + Displacement L2-norm : 0.0452791 + Max Z-displacement : 0.110911 + Velocity L2-norm : 0.481198 + Max Z-velocity : 1.17869 + Acceleration L2-norm : 2.16096 + Max Z-acceleration : 5.29349 + step=5 time=0.25 + Displacement L2-norm : 0.0717048 + Max Z-displacement : 0.17564 + Velocity L2-norm : 0.575828 + Max Z-velocity : 1.41049 + Acceleration L2-norm : 1.62422 + Max Z-acceleration : 3.97875 + step=6 time=0.3 + Displacement L2-norm : 0.102139 + Max Z-displacement : 0.250189 + Velocity L2-norm : 0.641548 + Max Z-velocity : 1.57146 + Acceleration L2-norm : 1.00458 + Max Z-acceleration : 2.46097 + step=7 time=0.35 + Displacement L2-norm : 0.135053 + Max Z-displacement : 0.33081 + Velocity L2-norm : 0.674992 + Max Z-velocity : 1.65339 + Acceleration L2-norm : 0.33318 + Max Z-acceleration : 0.816404 + step=8 time=0.4 + Displacement L2-norm : 0.168788 + Max Z-displacement : 0.413444 + Velocity L2-norm : 0.674417 + Max Z-velocity : 1.65198 + Acceleration L2-norm : 0.356163 + Max Z-acceleration : 0.872705 + step=9 time=0.45 + Displacement L2-norm : 0.201643 + Max Z-displacement : 0.493923 + Velocity L2-norm : 0.639796 + Max Z-velocity : 1.56717 + Acceleration L2-norm : 1.02868 + Max Z-acceleration : 2.52004 + step=10 time=0.5 + Displacement L2-norm : 0.231959 + Max Z-displacement : 0.56818 + Velocity L2-norm : 0.572819 + Max Z-velocity : 1.40311 + Acceleration L2-norm : 1.6504 + Max Z-acceleration : 4.04294 + step=11 time=0.55 + Displacement L2-norm : 0.258199 + Max Z-displacement : 0.632457 + Velocity L2-norm : 0.476812 + Max Z-velocity : 1.16795 + Acceleration L2-norm : 2.18986 + Max Z-acceleration : 5.36432 + step=12 time=0.6 + Displacement L2-norm : 0.279034 + Max Z-displacement : 0.683491 + Velocity L2-norm : 0.356574 + Max Z-velocity : 0.873425 + Acceleration L2-norm : 2.61966 + Max Z-acceleration : 6.41713 + step=13 time=0.65 + Displacement L2-norm : 0.293402 + Max Z-displacement : 0.718685 + Velocity L2-norm : 0.218134 + Max Z-velocity : 0.534317 + Acceleration L2-norm : 2.91795 + Max Z-acceleration : 7.14778 + step=14 time=0.7 + Displacement L2-norm : 0.300566 + Max Z-displacement : 0.736234 + Velocity L2-norm : 0.0684491 + Max Z-velocity : 0.167666 + Acceleration L2-norm : 3.06944 + Max Z-acceleration : 7.51886 + step=15 time=0.75 + Displacement L2-norm : 0.300154 + Max Z-displacement : 0.735224 + Velocity L2-norm : 0.0849424 + Max Z-velocity : 0.208065 + Acceleration L2-norm : 3.06622 + Max Z-acceleration : 7.511 + step=16 time=0.8 + Displacement L2-norm : 0.292173 + Max Z-displacement : 0.715674 + Velocity L2-norm : 0.234303 + Max Z-velocity : 0.573924 + Acceleration L2-norm : 2.90822 + Max Z-acceleration : 7.12397 + step=17 time=0.85 + Displacement L2-norm : 0.277013 + Max Z-displacement : 0.678541 + Velocity L2-norm : 0.372087 + Max Z-velocity : 0.911425 + Acceleration L2-norm : 2.60314 + Max Z-acceleration : 6.37668 + step=18 time=0.9 + Displacement L2-norm : 0.255428 + Max Z-displacement : 0.625668 + Velocity L2-norm : 0.491319 + Max Z-velocity : 1.20348 + Acceleration L2-norm : 2.16613 + Max Z-acceleration : 5.30625 + step=19 time=0.95 + Displacement L2-norm : 0.228496 + Max Z-displacement : 0.559699 + Velocity L2-norm : 0.585949 + Max Z-velocity : 1.43528 + Acceleration L2-norm : 1.61905 + Max Z-acceleration : 3.96614 + step=20 time=1 + Displacement L2-norm : 0.197568 + Max Z-displacement : 0.483942 + Velocity L2-norm : 0.651158 + Max Z-velocity : 1.595 + Acceleration L2-norm : 0.989308 + Max Z-acceleration : 2.4236 + step=21 time=1.05 + Displacement L2-norm : 0.164199 + Max Z-displacement : 0.402205 + Velocity L2-norm : 0.683605 + Max Z-velocity : 1.67448 + Acceleration L2-norm : 0.308569 + Max Z-acceleration : 0.756135 + step=22 time=1.1 + Displacement L2-norm : 0.130069 + Max Z-displacement : 0.318604 + Velocity L2-norm : 0.681597 + Max Z-velocity : 1.66956 + Acceleration L2-norm : 0.388891 + Max Z-acceleration : 0.952881 + step=23 time=1.15 + Displacement L2-norm : 0.0969 + Max Z-displacement : 0.237356 + Velocity L2-norm : 0.645177 + Max Z-velocity : 1.58036 + Acceleration L2-norm : 1.06789 + Max Z-acceleration : 2.61608 + step=24 time=1.2 + Displacement L2-norm : 0.0663674 + Max Z-displacement : 0.162566 + Velocity L2-norm : 0.576127 + Max Z-velocity : 1.41122 + Acceleration L2-norm : 1.69412 + Max Z-acceleration : 4.15002 + step=25 time=1.25 + Displacement L2-norm : 0.0400173 + Max Z-displacement : 0.0980222 + Velocity L2-norm : 0.477877 + Max Z-velocity : 1.17055 + Acceleration L2-norm : 2.23588 + Max Z-acceleration : 5.47706 + step=26 time=1.3 + Displacement L2-norm : 0.019187 + Max Z-displacement : 0.0469984 + Velocity L2-norm : 0.355338 + Max Z-velocity : 0.870397 + Acceleration L2-norm : 2.66568 + Max Z-acceleration : 6.52986 + step=27 time=1.35 + Displacement L2-norm : 0.00493712 + Max Z-displacement : 0.0120935 + Velocity L2-norm : 0.214655 + Max Z-velocity : 0.525796 + Acceleration L2-norm : 2.96163 + Max Z-acceleration : 7.25478 + step=28 time=1.4 + Displacement L2-norm : 0.00200176 + Max Z-displacement : 0.00490335 + Velocity L2-norm : 0.0629 + Max Z-velocity : 0.154073 + Acceleration L2-norm : 3.10858 + Max Z-acceleration : 7.6147 + step=29 time=1.45 + Displacement L2-norm : 0.00126712 + Max Z-displacement : 0.00310386 + Velocity L2-norm : 0.0922854 + Max Z-velocity : 0.226052 + Acceleration L2-norm : 3.09884 + Max Z-acceleration : 7.59085 + step=30 time=1.5 + Displacement L2-norm : 0.00711683 + Max Z-displacement : 0.0174327 + Velocity L2-norm : 0.243073 + Max Z-velocity : 0.595404 + Acceleration L2-norm : 2.93265 + Max Z-acceleration : 7.18377 + step=31 time=1.55 + Displacement L2-norm : 0.0227397 + Max Z-displacement : 0.0557008 + Velocity L2-norm : 0.381843 + Max Z-velocity : 0.93532 + Acceleration L2-norm : 2.61814 + Max Z-acceleration : 6.41339 + step=32 time=1.6 + Displacement L2-norm : 0.0447869 + Max Z-displacement : 0.109732 + Velocity L2-norm : 0.500047 + Max X-velocity : 0.00105501 + Max Z-velocity : 1.22591 + Max y-velocity : 0.0021101 + Acceleration L2-norm : 2.1104 + Max Z-acceleration : 5.21083 + step=33 time=1.65 + Displacement L2-norm : 0.072044 + Max X-displacement : 0.000333931 + Max Z-displacement : 0.176557 + Max y-displacement : 0.00017183 + Velocity L2-norm : 0.590248 + Max X-velocity : 0.0112462 + Max Z-velocity : 1.44709 + Max y-velocity : 0.00264632 + Acceleration L2-norm : 1.50656 + Max Z-acceleration : 3.69988 + step=34 time=1.7 + Displacement L2-norm : 0.103003 + Max X-displacement : 0.00139999 + Max Z-displacement : 0.252418 + Max y-displacement : 0.000232558 + Velocity L2-norm : 0.648208 + Max X-velocity : 0.0316127 + Max Z-velocity : 1.58757 + Max y-velocity : 0.000218367 + Acceleration L2-norm : 0.838713 + Max Z-acceleration : 2.02567 + step=35 time=1.75 + Displacement L2-norm : 0.135954 + Max X-displacement : 0.00350361 + Max Z-displacement : 0.33306 + Max y-displacement : 0.000121025 + Velocity L2-norm : 0.670093 + Max X-velocity : 0.0569959 + Max Z-velocity : 1.64258 + Max y-velocity : 0.00424633 + Acceleration L2-norm : 0.204948 + Max Z-acceleration : 0.175226 + step=36 time=1.8 + Displacement L2-norm : 0.16905 + Max X-displacement : 0.00674742 + Max Z-displacement : 0.414133 + Max y-displacement : 0.000202286 + Velocity L2-norm : 0.654242 + Max X-velocity : 0.077579 + Max Z-velocity : 1.60517 + Max y-velocity : 0.00865253 + Acceleration L2-norm : 0.705039 + Max Z-acceleration : 1.73763 + step=37 time=1.85 + Displacement L2-norm : 0.200457 + Max X-displacement : 0.0108441 + Max Z-displacement : 0.491247 + Max y-displacement : 0.000684053 + Velocity L2-norm : 0.602706 + Max X-velocity : 0.0862895 + Max Z-velocity : 1.47935 + Max y-velocity : 0.0106185 + Acceleration L2-norm : 1.37507 + Max Z-acceleration : 3.37461 + step=38 time=1.9 + Displacement L2-norm : 0.228498 + Max X-displacement : 0.0150606 + Max Z-displacement : 0.560112 + Max y-displacement : 0.00119473 + Velocity L2-norm : 0.519608 + Max X-velocity : 0.082369 + Max Z-velocity : 1.27526 + Max y-velocity : 0.00980882 + Acceleration L2-norm : 1.95278 + Max Z-acceleration : 4.80234 + step=39 time=1.95 + Displacement L2-norm : 0.251744 + Max X-displacement : 0.0187757 + Max Z-displacement : 0.617175 + Max y-displacement : 0.00159368 + Velocity L2-norm : 0.410776 + Max X-velocity : 0.0662347 + Max Z-velocity : 1.00725 + Max y-velocity : 0.00614965 + Acceleration L2-norm : 2.40124 + Max Z-acceleration : 5.91846 + step=40 time=2 + Displacement L2-norm : 0.269088 + Max X-displacement : 0.0214187 + Max Z-displacement : 0.659678 + Max y-displacement : 0.00174587 + Velocity L2-norm : 0.283331 + Max X-velocity : 0.0395505 + Max Z-velocity : 0.692921 + Acceleration L2-norm : 2.6985 + Max Z-acceleration : 6.65793 + step=41 time=2.05 + Displacement L2-norm : 0.279798 + Max X-displacement : 0.0224477 + Max Z-displacement : 0.68582 + Max y-displacement : 0.00160332 + Velocity L2-norm : 0.145231 + Max X-velocity : 0.00731053 + Max Z-velocity : 0.358466 + Max y-velocity : 0.00564044 + Acceleration L2-norm : 2.84049 + Max Z-acceleration : 6.94379 + step=42 time=2.1 + Displacement L2-norm : 0.283537 + Max X-displacement : 0.0213128 + Max Z-displacement : 0.694858 + Max y-displacement : 0.00129115 + Velocity L2-norm : 0.0188073 + Max X-velocity : 0.0470715 + Max Z-velocity : 0.0155382 + Max y-velocity : 0.00684744 + Acceleration L2-norm : 2.82534 + Max Z-acceleration : 6.84968 + step=43 time=2.15 + Displacement L2-norm : 0.280369 + Max X-displacement : 0.0176626 + Max Z-displacement : 0.687035 + Max y-displacement : 0.000981383 + Velocity L2-norm : 0.135975 + Max X-velocity : 0.0989377 + Max Z-velocity : 0.321585 + Max y-velocity : 0.00554394 + Acceleration L2-norm : 2.64638 + Max Z-acceleration : 6.41413 + step=44 time=2.2 + Displacement L2-norm : 0.270737 + Max X-displacement : 0.0116344 + Max Z-displacement : 0.663406 + Max y-displacement : 0.00067086 + Velocity L2-norm : 0.259457 + Max X-velocity : 0.14219 + Max Z-velocity : 0.62359 + Max y-velocity : 0.00687776 + Acceleration L2-norm : 2.32455 + Max Z-acceleration : 5.69121 + step=45 time=2.25 + Displacement L2-norm : 0.255414 + Max X-displacement : 0.00380398 + Max Z-displacement : 0.625725 + Max y-displacement : 0.000204177 + Velocity L2-norm : 0.364684 + Max X-velocity : 0.171026 + Max Z-velocity : 0.883672 + Max y-velocity : 0.0117907 + Acceleration L2-norm : 1.90659 + Max Z-acceleration : 4.71261 + step=46 time=2.3 + Displacement L2-norm : 0.235425 + Max X-displacement : 0.00513769 + Max Z-displacement : 0.576899 + Max y-displacement : 0.00049895 + Velocity L2-norm : 0.447428 + Max X-velocity : 0.186628 + Max Z-velocity : 1.0893 + Max y-velocity : 0.0162864 + Acceleration L2-norm : 1.42341 + Max Z-acceleration : 3.51293 + step=47 time=2.35 + Displacement L2-norm : 0.211969 + Max X-displacement : 0.0145437 + Max Z-displacement : 0.519703 + Max y-displacement : 0.00134922 + Velocity L2-norm : 0.505056 + Max X-velocity : 0.189615 + Max Z-velocity : 1.23253 + Max y-velocity : 0.0177268 + Acceleration L2-norm : 0.906089 + Max Z-acceleration : 2.2204 + step=48 time=2.4 + Displacement L2-norm : 0.186343 + Max X-displacement : 0.0237663 + Max Z-displacement : 0.45697 + Max y-displacement : 0.00219128 + Velocity L2-norm : 0.536392 + Max X-velocity : 0.179288 + Max Z-velocity : 1.3105 + Max y-velocity : 0.0159583 + Acceleration L2-norm : 0.40347 + Max Z-acceleration : 0.970034 + step=49 time=2.45 + Displacement L2-norm : 0.159859 + Max X-displacement : 0.0321609 + Max Z-displacement : 0.391786 + Max y-displacement : 0.00287188 + Velocity L2-norm : 0.541786 + Max X-velocity : 0.1565 + Max Z-velocity : 1.32406 + Max y-velocity : 0.0112686 + Acceleration L2-norm : 0.249599 + Max Z-acceleration : 0.360617 + step=50 time=2.5 + Displacement L2-norm : 0.133782 + Max X-displacement : 0.0391614 + Max Z-displacement : 0.327148 + Max y-displacement : 0.00326168 + Velocity L2-norm : 0.522845 + Max X-velocity : 0.12352 + Max Z-velocity : 1.27709 + Max y-velocity : 0.00432646 + Acceleration L2-norm : 0.652774 + Max Z-acceleration : 1.51858 + step=51 time=2.55 + Displacement L2-norm : 0.109271 + Max X-displacement : 0.0443417 + Max Z-displacement : 0.265834 + Max y-displacement : 0.00327057 + Velocity L2-norm : 0.482106 + Max X-velocity : 0.0876639 + Max Z-velocity : 1.17978 + Max y-velocity : 0.00396962 + Acceleration L2-norm : 1.04905 + Max Z-acceleration : 2.5333 + step=52 time=2.6 + Displacement L2-norm : 0.0873394 + Max X-displacement : 0.0474589 + Max Z-displacement : 0.210325 + Max y-displacement : 0.00285545 + Velocity L2-norm : 0.422732 + Max X-velocity : 0.0536283 + Max Z-velocity : 1.0406 + Max y-velocity : 0.0126353 + Acceleration L2-norm : 1.38456 + Max Z-acceleration : 3.38105 + step=53 time=2.65 + Displacement L2-norm : 0.068816 + Max X-displacement : 0.0484719 + Max Z-displacement : 0.162736 + Max y-displacement : 0.002021 + Velocity L2-norm : 0.348241 + Max X-velocity : 0.02027 + Max Z-velocity : 0.862949 + Max y-velocity : 0.0207431 + Acceleration L2-norm : 1.65031 + Max Z-acceleration : 4.04971 + step=54 time=2.7 + Displacement L2-norm : 0.0543017 + Max X-displacement : 0.0475379 + Max Z-displacement : 0.124785 + Max y-displacement : 0.00081478 + Velocity L2-norm : 0.262358 + Max X-velocity : 0.0368887 + Max Z-velocity : 0.655095 + Max y-velocity : 0.0275061 + Acceleration L2-norm : 1.8447 + Max Z-acceleration : 4.53552 + step=55 time=2.75 + Displacement L2-norm : 0.0441228 + Max X-displacement : 0.0456702 + Max Z-displacement : 0.0984504 + Max y-displacement : 0.000681884 + Velocity L2-norm : 0.169092 + Max X-velocity : 0.0650868 + Max Z-velocity : 0.425545 + Max y-velocity : 0.0323311 + Acceleration L2-norm : 1.96879 + Max Z-acceleration : 4.84002 + step=56 time=2.8 + Displacement L2-norm : 0.0383502 + Max X-displacement : 0.0436525 + Max Z-displacement : 0.0849192 + Max y-displacement : 0.00236089 + Velocity L2-norm : 0.0746819 + Max X-velocity : 0.0827809 + Max Z-velocity : 0.18286 + Max y-velocity : 0.0348303 + Acceleration L2-norm : 2.02503 + Max Z-acceleration : 4.96792 + step=57 time=2.85 + Displacement L2-norm : 0.0370074 + Max X-displacement : 0.0411065 + Max Z-displacement : 0.0837008 + Max y-displacement : 0.00410196 + Velocity L2-norm : 0.0476199 + Max X-velocity : 0.0886973 + Max Z-velocity : 0.099294 + Max y-velocity : 0.0348133 + Acceleration L2-norm : 2.0161 + Max Z-acceleration : 4.95161 + step=58 time=2.9 + Displacement L2-norm : 0.0403031 + Max X-displacement : 0.0384974 + Max Z-displacement : 0.0946952 + Max y-displacement : 0.00578137 + Velocity L2-norm : 0.136127 + Max X-velocity : 0.0828411 + Max Z-velocity : 0.340479 + Max y-velocity : 0.0323632 + Acceleration L2-norm : 1.94445 + Max Z-acceleration : 4.79441 + step=59 time=2.95 + Displacement L2-norm : 0.0484494 + Max X-displacement : 0.0362617 + Max Z-displacement : 0.117402 + Max y-displacement : 0.00729164 + Velocity L2-norm : 0.227542 + Max X-velocity : 0.0669983 + Max Z-velocity : 0.567799 + Max y-velocity : 0.0280482 + Acceleration L2-norm : 1.81128 + Max Z-acceleration : 4.47158 + step=60 time=3 + Displacement L2-norm : 0.0613084 + Max X-displacement : 0.0347899 + Max Z-displacement : 0.150932 + Max y-displacement : 0.00852989 + Velocity L2-norm : 0.311797 + Max X-velocity : 0.0414044 + Max Z-velocity : 0.773416 + Max y-velocity : 0.0214828 + Acceleration L2-norm : 1.62354 + Max Z-acceleration : 4.01635 + step=61 time=3.05 + Displacement L2-norm : 0.0783804 + Max X-displacement : 0.0344168 + Max Z-displacement : 0.19402 + Max y-displacement : 0.00938286 + Velocity L2-norm : 0.385338 + Max X-velocity : 0.00763309 + Max Z-velocity : 0.950072 + Max y-velocity : 0.0126361 + Acceleration L2-norm : 1.37848 + Max Z-acceleration : 3.40438 + step=62 time=3.1 + Displacement L2-norm : 0.0989586 + Max X-displacement : 0.035353 + Max Z-displacement : 0.245059 + Max y-displacement : 0.00976351 + Velocity L2-norm : 0.445265 + Max X-velocity : 0.0324427 + Max Z-velocity : 1.09152 + Max y-velocity : 0.00259045 + Acceleration L2-norm : 1.08113 + Max Z-acceleration : 2.65588 + step=63 time=3.15 + Displacement L2-norm : 0.122217 + Max X-displacement : 0.0376506 + Max Z-displacement : 0.302143 + Max y-displacement : 0.00962782 + Velocity L2-norm : 0.488908 + Max X-velocity : 0.0674802 + Max Z-velocity : 1.19986 + Max y-velocity : 0.00801896 + Acceleration L2-norm : 0.738038 + Max Z-acceleration : 1.78205 + step=64 time=3.2 + Displacement L2-norm : 0.147234 + Max X-displacement : 0.0412148 + Max Z-displacement : 0.363087 + Max y-displacement : 0.0089662 + Velocity L2-norm : 0.513868 + Max X-velocity : 0.101555 + Max Z-velocity : 1.26434 + Max y-velocity : 0.0184481 + Acceleration L2-norm : 0.370305 + Max Z-acceleration : 0.797809 + step=65 time=3.25 + Displacement L2-norm : 0.172997 + Max X-displacement : 0.0458101 + Max Z-displacement : 0.425472 + Max y-displacement : 0.00780804 + Velocity L2-norm : 0.518102 + Max X-velocity : 0.128587 + Max Z-velocity : 1.27738 + Max y-velocity : 0.0278808 + Acceleration L2-norm : 0.228165 + Max Z-acceleration : 0.450938 + step=66 time=3.3 + Displacement L2-norm : 0.198418 + Max X-displacement : 0.0510819 + Max Z-displacement : 0.486703 + Max y-displacement : 0.00622428 + Velocity L2-norm : 0.500083 + Max X-velocity : 0.145634 + Max Z-velocity : 1.23521 + Max y-velocity : 0.0354723 + Acceleration L2-norm : 0.609212 + Max Z-acceleration : 1.53973 + step=67 time=3.35 + Displacement L2-norm : 0.222361 + Max X-displacement : 0.0565912 + Max Z-displacement : 0.544081 + Max y-displacement : 0.00432701 + Velocity L2-norm : 0.458994 + Max X-velocity : 0.150629 + Max Z-velocity : 1.13582 + Max y-velocity : 0.0404214 + Acceleration L2-norm : 1.06221 + Max Z-acceleration : 2.6344 + step=68 time=3.4 + Displacement L2-norm : 0.243672 + Max X-displacement : 0.061861 + Max Z-displacement : 0.5949 + Max y-displacement : 0.00226487 + Velocity L2-norm : 0.394956 + Max X-velocity : 0.14265 + Max Z-velocity : 0.979415 + Max y-velocity : 0.0420671 + Acceleration L2-norm : 1.51084 + Max Z-acceleration : 3.69125 + step=69 time=3.45 + Displacement L2-norm : 0.261237 + Max X-displacement : 0.0664278 + Max Z-displacement : 0.636558 + Max y-displacement : 0.000213258 + Velocity L2-norm : 0.309245 + Max X-velocity : 0.122087 + Max Z-velocity : 0.768982 + Max y-velocity : 0.0399998 + Acceleration L2-norm : 1.92418 + Max Z-acceleration : 4.72669 + step=70 time=3.5 + Displacement L2-norm : 0.274039 + Max X-displacement : 0.0715341 + Max Z-displacement : 0.668338 + Max y-displacement : 0.00164229 + Velocity L2-norm : 0.204466 + Max X-velocity : 0.0906703 + Max Z-velocity : 0.510719 + Max y-velocity : 0.0341761 + Acceleration L2-norm : 2.27335 + Max Z-acceleration : 5.6044 + step=71 time=3.55 + Displacement L2-norm : 0.281221 + Max X-displacement : 0.0750845 + Max Z-displacement : 0.686462 + Max y-displacement : 0.00312179 + Velocity L2-norm : 0.0848009 + Max X-velocity : 0.0513455 + Max Z-velocity : 0.214245 + Max y-velocity : 0.0250039 + Acceleration L2-norm : 2.53035 + Max Z-acceleration : 6.25517 + step=72 time=3.6 + Displacement L2-norm : 0.282149 + Max X-displacement : 0.0764371 + Max Z-displacement : 0.689183 + Max y-displacement : 0.00418338 + Velocity L2-norm : 0.0474828 + Max X-velocity : 0.0147048 + Max Z-velocity : 0.122901 + Max y-velocity : 0.0174612 + Acceleration L2-norm : 2.68383 + Max Z-acceleration : 6.53278 + step=73 time=3.65 + Displacement L2-norm : 0.276468 + Max X-displacement : 0.0751914 + Max Z-displacement : 0.675751 + Max y-displacement : 0.00497552 + Velocity L2-norm : 0.180991 + Max X-velocity : 0.0668113 + Max Z-velocity : 0.446064 + Max y-velocity : 0.0142263 + Acceleration L2-norm : 2.69238 + Max Z-acceleration : 6.52384 + step=74 time=3.7 + Displacement L2-norm : 0.264166 + Max X-displacement : 0.0712946 + Max Z-displacement : 0.646202 + Max y-displacement : 0.00560483 + Velocity L2-norm : 0.31169 + Max X-velocity : 0.114232 + Max Z-velocity : 0.76105 + Max y-velocity : 0.0109465 + Acceleration L2-norm : 2.54078 + Max Z-acceleration : 6.2074 + step=75 time=3.75 + Displacement L2-norm : 0.245616 + Max X-displacement : 0.0652554 + Max Z-displacement : 0.601332 + Max y-displacement : 0.00595048 + Velocity L2-norm : 0.430882 + Max X-velocity : 0.141158 + Max Z-velocity : 1.04758 + Max y-velocity : 0.00288034 + Acceleration L2-norm : 2.23748 + Max Z-acceleration : 5.57734 + step=76 time=3.8 + Displacement L2-norm : 0.221563 + Max X-displacement : 0.0578848 + Max Z-displacement : 0.542793 + Max y-displacement : 0.00583025 + Velocity L2-norm : 0.531541 + Max X-velocity : 0.156544 + Max Z-velocity : 1.29683 + Max y-velocity : 0.00768885 + Acceleration L2-norm : 1.80508 + Max Z-acceleration : 4.50804 + step=77 time=3.85 + Displacement L2-norm : 0.19309 + Max X-displacement : 0.049921 + Max Z-displacement : 0.473183 + Max y-displacement : 0.0052532 + Velocity L2-norm : 0.607483 + Max X-velocity : 0.162008 + Max Z-velocity : 1.48758 + Max y-velocity : 0.0153933 + Acceleration L2-norm : 1.25224 + Max Z-acceleration : 3.12253 + step=78 time=3.9 + Displacement L2-norm : 0.161573 + Max X-displacement : 0.0420117 + Max Z-displacement : 0.39592 + Max y-displacement : 0.00438453 + Velocity L2-norm : 0.653207 + Max X-velocity : 0.154361 + Max Z-velocity : 1.60297 + Max y-velocity : 0.0193537 + Acceleration L2-norm : 0.614588 + Max Z-acceleration : 1.49377 + step=79 time=3.95 + Displacement L2-norm : 0.128632 + Max X-displacement : 0.0347947 + Max Z-displacement : 0.31502 + Max y-displacement : 0.00341422 + Velocity L2-norm : 0.664735 + Max X-velocity : 0.134321 + Max Z-velocity : 1.63304 + Max y-velocity : 0.0194589 + Acceleration L2-norm : 0.23323 + Max Z-acceleration : 0.290877 + step=80 time=4 + Displacement L2-norm : 0.0960444 + Max X-displacement : 0.0288393 + Max Z-displacement : 0.234879 + Max y-displacement : 0.00252356 + Velocity L2-norm : 0.639948 + Max X-velocity : 0.103896 + Max Z-velocity : 1.57257 + Max y-velocity : 0.0161675 + Acceleration L2-norm : 0.891096 + Max Z-acceleration : 2.12845 + step=81 time=4.05 + Displacement L2-norm : 0.0656708 + Max X-displacement : 0.0245887 + Max Z-displacement : 0.160025 + Max y-displacement : 0.00186052 + Velocity L2-norm : 0.578814 + Max X-velocity : 0.0661253 + Max Z-velocity : 1.42161 + Max y-velocity : 0.0103541 + Acceleration L2-norm : 1.59917 + Max Z-acceleration : 3.91076 + step=82 time=4.1 + Displacement L2-norm : 0.0394109 + Max X-displacement : 0.0223179 + Max Z-displacement : 0.0948428 + Max y-displacement : 0.0015242 + Velocity L2-norm : 0.483508 + Max X-velocity : 0.0247079 + Max Z-velocity : 1.18566 + Max y-velocity : 0.00309838 + Acceleration L2-norm : 2.2484 + Max Z-acceleration : 5.5276 + step=83 time=4.15 + Displacement L2-norm : 0.0193672 + Max X-displacement : 0.0216225 + Max Z-displacement : 0.0432033 + Max y-displacement : 0.00134456 + Velocity L2-norm : 0.358396 + Max X-velocity : 0.00311108 + Max Z-velocity : 0.879919 + Max y-velocity : 0.00408878 + Acceleration L2-norm : 2.77791 + Max Z-acceleration : 6.8965 + step=84 time=4.2 + Displacement L2-norm : 0.00913705 + Max X-displacement : 0.0214065 + Max Z-displacement : 0.00819711 + Max y-displacement : 0.00088099 + Velocity L2-norm : 0.209599 + Max X-velocity : 0.00892342 + Max Z-velocity : 0.520328 + Max y-velocity : 0.0144551 + Acceleration L2-norm : 3.18277 + Max Z-acceleration : 7.902 + step=85 time=4.25 + Displacement L2-norm : 0.00916473 + Max X-displacement : 0.0210853 + Max Z-displacement : 0.00784976 + Velocity L2-norm : 0.0463128 + Max X-velocity : 0.0134514 + Max Z-velocity : 0.121542 + Max y-velocity : 0.023907 + Acceleration L2-norm : 3.40939 + Max Z-acceleration : 8.42816 + step=86 time=4.3 + Displacement L2-norm : 0.0087265 + Max X-displacement : 0.0218295 + Max Z-displacement : 0.00351976 + Max y-displacement : 0.00142212 + Velocity L2-norm : 0.127299 + Max X-velocity : 0.0163163 + Max Z-velocity : 0.324586 + Max y-velocity : 0.0298471 + Acceleration L2-norm : 3.43811 + Max Z-acceleration : 8.46128 + step=87 time=4.35 + Displacement L2-norm : 0.0128312 + Max X-displacement : 0.0226678 + Max Z-displacement : 0.0244103 + Max y-displacement : 0.00296077 + Velocity L2-norm : 0.294318 + Max X-velocity : 0.0172158 + Max Z-velocity : 0.73573 + Max y-velocity : 0.031701 + Acceleration L2-norm : 3.26095 + Max Z-acceleration : 7.99091 + step=88 time=4.4 + Displacement L2-norm : 0.0293005 + Max X-displacement : 0.0235001 + Max Z-displacement : 0.0705725 + Max y-displacement : 0.00448433 + Velocity L2-norm : 0.447735 + Max X-velocity : 0.0160768 + Max Z-velocity : 1.11076 + Max y-velocity : 0.0292419 + Acceleration L2-norm : 2.88367 + Max Z-acceleration : 7.10927 + step=89 time=4.45 + Displacement L2-norm : 0.0543065 + Max X-displacement : 0.0242291 + Max Z-displacement : 0.133999 + Max y-displacement : 0.00577953 + Velocity L2-norm : 0.577792 + Max X-velocity : 0.0130822 + Max Z-velocity : 1.42631 + Max y-velocity : 0.0225662 + Acceleration L2-norm : 2.32577 + Max Z-acceleration : 5.78046 + step=90 time=4.5 + Displacement L2-norm : 0.0853813 + Max X-displacement : 0.0246609 + Max Z-displacement : 0.211257 + Max y-displacement : 0.00673266 + Velocity L2-norm : 0.676192 + Max X-velocity : 0.0113703 + Max Z-velocity : 1.66399 + Max y-velocity : 0.015559 + Acceleration L2-norm : 1.61794 + Max Z-acceleration : 4.00782 + step=91 time=4.55 + Displacement L2-norm : 0.120536 + Max X-displacement : 0.0243044 + Max Z-displacement : 0.298144 + Max y-displacement : 0.00750368 + Velocity L2-norm : 0.736539 + Max X-velocity : 0.0337326 + Max Z-velocity : 1.8115 + Max y-velocity : 0.0152824 + Acceleration L2-norm : 0.844113 + Max Z-acceleration : 1.99898 + step=92 time=4.6 + Displacement L2-norm : 0.157691 + Max X-displacement : 0.0225797 + Max Z-displacement : 0.389873 + Max y-displacement : 0.00837104 + Velocity L2-norm : 0.754882 + Max X-velocity : 0.0699516 + Max Z-velocity : 1.85765 + Max y-velocity : 0.0194127 + Acceleration L2-norm : 0.295055 + Max Z-acceleration : 0.22357 + step=93 time=4.65 + Displacement L2-norm : 0.194682 + Max X-displacement : 0.0192512 + Max Z-displacement : 0.48121 + Max y-displacement : 0.00939975 + Velocity L2-norm : 0.729795 + Max X-velocity : 0.104332 + Max Z-velocity : 1.79585 + Max y-velocity : 0.0217357 + Acceleration L2-norm : 0.98625 + Max Z-acceleration : 2.34203 + step=94 time=4.7 + Displacement L2-norm : 0.229359 + Max X-displacement : 0.0143959 + Max Z-displacement : 0.566826 + Max y-displacement : 0.0104719 + Velocity L2-norm : 0.662602 + Max X-velocity : 0.132763 + Max Z-velocity : 1.6288 + Max y-velocity : 0.0211501 + Acceleration L2-norm : 1.78423 + Max Z-acceleration : 4.36327 + step=95 time=4.75 + Displacement L2-norm : 0.259703 + Max X-displacement : 0.00823132 + Max Z-displacement : 0.641691 + Max y-displacement : 0.0114188 + Velocity L2-norm : 0.55731 + Max X-velocity : 0.151692 + Max Z-velocity : 1.3658 + Max y-velocity : 0.0167251 + Acceleration L2-norm : 2.49143 + Max Z-acceleration : 6.1571 + step=96 time=4.8 + Displacement L2-norm : 0.283956 + Max X-displacement : 0.010945 + Max Z-displacement : 0.701391 + Max y-displacement : 0.0120344 + Velocity L2-norm : 0.420477 + Max X-velocity : 0.158631 + Max Z-velocity : 1.02217 + Max y-velocity : 0.0078988 + Acceleration L2-norm : 3.05489 + Max Z-acceleration : 7.5884 + step=97 time=4.85 + Displacement L2-norm : 0.30072 + Max X-displacement : 0.018821 + Max Z-displacement : 0.742358 + Max y-displacement : 0.0120133 + Velocity L2-norm : 0.26188 + Max X-velocity : 0.16515 + Max Z-velocity : 0.625277 + Max y-velocity : 0.00874041 + Acceleration L2-norm : 3.44326 + Max Z-acceleration : 8.63759 + step=98 time=4.9 + Displacement L2-norm : 0.30905 + Max X-displacement : 0.0262575 + Max Z-displacement : 0.762142 + Max y-displacement : 0.011065 + Velocity L2-norm : 0.101142 + Max X-velocity : 0.170243 + Max Z-velocity : 0.203996 + Max y-velocity : 0.029192 + Acceleration L2-norm : 3.61801 + Max Z-acceleration : 9.03182 + step=99 time=4.95 + Displacement L2-norm : 0.308506 + Max X-displacement : 0.0324474 + Max Z-displacement : 0.759689 + Max y-displacement : 0.00920831 + Velocity L2-norm : 0.117184 + Max X-velocity : 0.151616 + Max Z-velocity : 0.272904 + Max y-velocity : 0.045077 + Acceleration L2-norm : 3.58565 + Max Z-acceleration : 8.87702 + step=100 time=5 + Displacement L2-norm : 0.299171 + Max X-displacement : 0.0365336 + Max Z-displacement : 0.735378 + Max y-displacement : 0.00672104 + Velocity L2-norm : 0.277956 + Max X-velocity : 0.111319 + Max Z-velocity : 0.699546 + Max y-velocity : 0.0544148 + Acceleration L2-norm : 3.34887 + Max Z-acceleration : 8.18899 + Time integration completed. diff --git a/Test/DroppedBox.xinp b/Test/DroppedBox.xinp new file mode 100644 index 0000000..9f27519 --- /dev/null +++ b/Test/DroppedBox.xinp @@ -0,0 +1,35 @@ + + + + + Q4box.nas + + + + + + wavel=4.0; waveh=0.1; wavespeed=-4.0; tstart=1.0; + refx=3.0+wavespeed*(t-tstart); + if(above(x+ux,refx),0.5*waveh*sin((x+ux-refx)*2*PI/wavel),0.0) + + + + + + 2 + + + + 0 + + 0.05 + + + + + + + + + + diff --git a/Test/Q4box.nas b/Test/Q4box.nas new file mode 100644 index 0000000..a6d3e61 --- /dev/null +++ b/Test/Q4box.nas @@ -0,0 +1,18 @@ +BEGIN BULK +GRID,1 +GRID,2,,1 +GRID,3,,1,1 +GRID,4,,0,1 +GRID,5,,0,0,1 +GRID,6,,1,0,1 +GRID,7,,1,1,1 +GRID,8,,0,1,1 +CQUAD4,1,1,1,4,3,2 +CQUAD4,2,1,1,2,6,5 +CQUAD4,3,1,2,3,7,6 +CQUAD4,4,1,3,4,8,7 +CQUAD4,5,1,4,1,5,8 +CQUAD4,6,1,5,6,7,8 +PSHELL,1,1,0.01 +MAT1,1,21E10,,0.3,7800 +ENDDATA diff --git a/Test/Robaat-irr.reg b/Test/Robaat-irr.reg new file mode 100644 index 0000000..850d083 --- /dev/null +++ b/Test/Robaat-irr.reg @@ -0,0 +1,499 @@ +Robaat-irr.xinp -dynamic -keep-previous-state -time 1.5 + +Input file: Robaat-irr.xinp +Simulation stop time: 1.5 +Using the linear dynamics simulation driver. +Parsing input file Robaat-irr.xinp +Parsing +Parsing + Parsing + Reading data file Small_Boat.dat + Nastran bulk data starting at line 71 +Parsing Nastran bulk data file succceeded. +Total number of nodes: 518 +Number of shell elements: 504 +Number of constraint elements: 0 +Number of other elements: 3 +Model extension (diameter): 3.5 + Reading patch 1 + Parsing + Topology sets: rib (1,1,5D) + wet (1,2,5D) + Parsing + Parsing +Parsing + Parsing + Gravitation vector: 0 0 -9.81 + Parsing + Sea surface function (spectrum): file = run18-irr482.asc, wave angle = 210, ramp up time = 1 + Parsing + Pressure code 1 (hydrostatic): z0 = 0.1 rhow = 1000 +Parsing + alpha1 = 1 alpha2 = 0.03 + beta = 0.25 gamma = 0.5 +Parsing + Parsing + Parsing + Parsing + Sea surface \[-10 -2.5 0.1] to \[10 2.5 0.1], gridSize=0.5 +Parsing input file succeeded. +Equation solver: 2 +Number of Gauss points: 2 +Lagrangian basis functions are used +Problem definition: +Formulation: ANDES shell +Newmark predictor/multicorrector: beta = 0.25 gamma = 0.5 + using constant coefficient matrices + using zero acceleration predictor +Mass-proportional damping (alpha1): 1 +Stiffness-proportional damping (alpha2): 0.03 +Renumbered 518 nodes. +Resolving Dirichlet boundary conditions + >>> SAM model summary <<< +Number of elements 507 +Number of nodes 518 +Number of dofs 3108 +Detected 492 holes in the element numbers, reordering... +Number of unknowns 3108 +Total element set size: 358 + step=131 time=1.31 + Displacement L2-norm : 0.0449315 + Max X-displacement : 0.00197266 + Max Y-displacement : 0.00691991 + Max Z-displacement : 0.112859 + Max x-displacement : 0.00162378 + Max y-displacement : 0.0019875 + Max z-displacement : 0.00308857 + Velocity L2-norm : 0.106088 + Max X-velocity : 0.0744149 + Max Y-velocity : 0.0270222 + Max Z-velocity : 0.341576 + Max x-velocity : 0.0317302 + Max y-velocity : 0.190009 + Max z-velocity : 0.0101407 + Acceleration L2-norm : 1.13387 + Max X-acceleration : 0.262826 + Max Y-acceleration : 0.122492 + Max Z-acceleration : 3.63258 + Max x-acceleration : 1.04239 + Max y-acceleration : 1.20498 + Max z-acceleration : 0.299877 + step=132 time=1.32 + Displacement L2-norm : 0.0448381 + Max X-displacement : 0.00238757 + Max Y-displacement : 0.00705271 + Max Z-displacement : 0.115851 + Max x-displacement : 0.00129941 + Max y-displacement : 0.00385544 + Max z-displacement : 0.0031867 + Velocity L2-norm : 0.104107 + Max X-velocity : 0.0727882 + Max Y-velocity : 0.0281914 + Max Z-velocity : 0.357206 + Max x-velocity : 0.0339245 + Max y-velocity : 0.183684 + Max z-velocity : 0.0100856 + Acceleration L2-norm : 1.14934 + Max X-acceleration : 0.319886 + Max Y-acceleration : 0.116115 + Max Z-acceleration : 3.88541 + Max x-acceleration : 1.05076 + Max y-acceleration : 1.29707 + Max z-acceleration : 0.310106 + step=133 time=1.33 + Displacement L2-norm : 0.0446592 + Max X-displacement : 0.00278222 + Max Y-displacement : 0.00718745 + Max Z-displacement : 0.118454 + Max x-displacement : 0.000953882 + Max y-displacement : 0.00565305 + Max z-displacement : 0.00328405 + Velocity L2-norm : 0.102607 + Max X-velocity : 0.0706349 + Max Y-velocity : 0.0292389 + Max Z-velocity : 0.370449 + Max x-velocity : 0.0358708 + Max y-velocity : 0.17597 + Max z-velocity : 0.00998249 + Acceleration L2-norm : 1.17954 + Max X-acceleration : 0.355175 + Max Y-acceleration : 0.118692 + Max Z-acceleration : 4.11051 + Max x-acceleration : 0.98385 + Max y-acceleration : 1.47299 + Max z-acceleration : 0.302268 + step=134 time=1.34 + Displacement L2-norm : 0.0443927 + Max X-displacement : 0.00315362 + Max Y-displacement : 0.007324 + Max Z-displacement : 0.120648 + Max x-displacement : 0.000589793 + Max y-displacement : 0.00736686 + Max z-displacement : 0.00338024 + Velocity L2-norm : 0.101703 + Max X-velocity : 0.0679686 + Max Y-velocity : 0.0301512 + Max Z-velocity : 0.381279 + Max x-velocity : 0.0375346 + Max y-velocity : 0.166948 + Max z-velocity : 0.00983673 + Acceleration L2-norm : 1.19281 + Max X-acceleration : 0.40705 + Max Y-acceleration : 0.119365 + Max Z-acceleration : 4.32028 + Max x-acceleration : 0.985346 + Max y-acceleration : 1.55374 + Max z-acceleration : 0.314218 + step=135 time=1.35 + Displacement L2-norm : 0.0440367 + Max X-displacement : 0.00355999 + Max Y-displacement : 0.00759416 + Max Z-displacement : 0.122409 + Max x-displacement : 0.000710039 + Max y-displacement : 0.00898401 + Max z-displacement : 0.00347492 + Velocity L2-norm : 0.101499 + Max X-velocity : 0.064808 + Max Y-velocity : 0.0309181 + Max Z-velocity : 0.390021 + Max x-velocity : 0.0388847 + Max y-velocity : 0.156659 + Max z-velocity : 0.00965048 + Acceleration L2-norm : 1.21895 + Max X-acceleration : 0.43633 + Max Y-acceleration : 0.119011 + Max Z-acceleration : 4.49787 + Max x-acceleration : 0.911791 + Max y-acceleration : 1.71659 + Max z-acceleration : 0.305902 + step=136 time=1.36 + Displacement L2-norm : 0.0435896 + Max X-displacement : 0.00418991 + Max Y-displacement : 0.0079064 + Max Z-displacement : 0.123722 + Max x-displacement : 0.00110104 + Max y-displacement : 0.0104923 + Max z-displacement : 0.00356774 + Velocity L2-norm : 0.10208 + Max X-velocity : 0.061176 + Max Y-velocity : 0.0315296 + Max Z-velocity : 0.396457 + Max x-velocity : 0.0398936 + Max y-velocity : 0.145207 + Max z-velocity : 0.00942779 + Acceleration L2-norm : 1.22754 + Max X-acceleration : 0.48168 + Max Y-acceleration : 0.116964 + Max Z-acceleration : 4.65686 + Max x-acceleration : 0.907981 + Max y-acceleration : 1.78153 + Max z-acceleration : 0.317596 + step=137 time=1.37 + Displacement L2-norm : 0.0430501 + Max X-displacement : 0.00478128 + Max Y-displacement : 0.00822394 + Max Z-displacement : 0.12457 + Max x-displacement : 0.00150028 + Max y-displacement : 0.0118805 + Max z-displacement : 0.00365838 + Velocity L2-norm : 0.103499 + Max X-velocity : 0.0570996 + Max Y-velocity : 0.0319779 + Max Z-velocity : 0.400507 + Max x-velocity : 0.0405393 + Max y-velocity : 0.1327 + Max z-velocity : 0.00917034 + Acceleration L2-norm : 1.24721 + Max X-acceleration : 0.504166 + Max Y-acceleration : 0.11421 + Max Z-acceleration : 4.78087 + Max x-acceleration : 0.829769 + Max y-acceleration : 1.92786 + Max z-acceleration : 0.308872 + step=138 time=1.38 + Displacement L2-norm : 0.0424172 + Max X-displacement : 0.00532983 + Max Y-displacement : 0.0085451 + Max Z-displacement : 0.12494 + Max x-displacement : 0.00190407 + Max y-displacement : 0.0131382 + Max z-displacement : 0.00374655 + Velocity L2-norm : 0.105773 + Max X-velocity : 0.0526098 + Max Y-velocity : 0.0322552 + Max Z-velocity : 0.402211 + Max x-velocity : 0.040805 + Max y-velocity : 0.119171 + Max z-velocity : 0.0088801 + Acceleration L2-norm : 1.249 + Max X-acceleration : 0.542169 + Max Y-acceleration : 0.110642 + Max Z-acceleration : 4.88353 + Max x-acceleration : 0.822935 + Max y-acceleration : 1.97358 + Max z-acceleration : 0.320537 + step=139 time=1.39 + Displacement L2-norm : 0.0416904 + Max X-displacement : 0.00583159 + Max Y-displacement : 0.00886815 + Max Z-displacement : 0.124823 + Max x-displacement : 0.00230856 + Max y-displacement : 0.0142558 + Max z-displacement : 0.00383193 + Velocity L2-norm : 0.108882 + Max X-velocity : 0.0477424 + Max Y-velocity : 0.0323547 + Max Z-velocity : 0.40163 + Max x-velocity : 0.0407563 + Max y-velocity : 0.104775 + Max z-velocity : 0.00855726 + Acceleration L2-norm : 1.25999 + Max X-acceleration : 0.556939 + Max Y-acceleration : 0.11659 + Max Z-acceleration : 4.94804 + Max x-acceleration : 0.829079 + Max y-acceleration : 2.09966 + Max z-acceleration : 0.311662 + step=140 time=1.4 + Displacement L2-norm : 0.0408695 + Max X-displacement : 0.00628299 + Max Y-displacement : 0.00919129 + Max Z-displacement : 0.124212 + Max x-displacement : 0.00270983 + Max y-displacement : 0.015225 + Max z-displacement : 0.00391425 + Velocity L2-norm : 0.112772 + Max X-velocity : 0.0425379 + Max Y-velocity : 0.0322735 + Max Z-velocity : 0.398847 + Max x-velocity : 0.0402662 + Max y-velocity : 0.0895762 + Max z-velocity : 0.00822547 + Acceleration L2-norm : 1.25287 + Max X-acceleration : 0.586495 + Max Y-acceleration : 0.127998 + Max Z-acceleration : 4.98849 + Max x-acceleration : 0.843649 + Max y-acceleration : 2.12198 + Max z-acceleration : 0.322911 + step=141 time=1.41 + Displacement L2-norm : 0.0399548 + Max X-displacement : 0.00668087 + Max Y-displacement : 0.00951272 + Max Z-displacement : 0.123102 + Max x-displacement : 0.0031039 + Max y-displacement : 0.0160387 + Max z-displacement : 0.00399323 + Velocity L2-norm : 0.117355 + Max X-velocity : 0.037037 + Max Y-velocity : 0.0320122 + Max Z-velocity : 0.393942 + Max x-velocity : 0.0394245 + Max y-velocity : 0.0737364 + Max z-velocity : 0.00789351 + Acceleration L2-norm : 1.25376 + Max X-acceleration : 0.593194 + Max Y-acceleration : 0.137813 + Max Z-acceleration : 4.9911 + Max x-acceleration : 0.899013 + Max y-acceleration : 2.22596 + Max z-acceleration : 0.31355 + step=142 time=1.42 + Displacement L2-norm : 0.0389473 + Max X-displacement : 0.00702246 + Max Y-displacement : 0.00983065 + Max Z-displacement : 0.121495 + Max x-displacement : 0.00348682 + Max y-displacement : 0.0166907 + Max z-displacement : 0.00406863 + Velocity L2-norm : 0.122524 + Max X-velocity : 0.0312818 + Max Y-velocity : 0.0315728 + Max Z-velocity : 0.387 + Max x-velocity : 0.0381438 + Max y-velocity : 0.0573203 + Max z-velocity : 0.00753922 + Acceleration L2-norm : 1.23711 + Max X-acceleration : 0.614841 + Max Y-acceleration : 0.148112 + Max Z-acceleration : 4.97026 + Max x-acceleration : 0.912626 + Max y-acceleration : 2.22517 + Max z-acceleration : 0.324495 + step=143 time=1.43 + Displacement L2-norm : 0.0378483 + Max X-displacement : 0.00730545 + Max Y-displacement : 0.0101433 + Max Z-displacement : 0.119391 + Max x-displacement : 0.00385475 + Max y-displacement : 0.0171762 + Max z-displacement : 0.00414024 + Velocity L2-norm : 0.12816 + Max X-velocity : 0.0253156 + Max Y-velocity : 0.0309589 + Max Z-velocity : 0.378963 + Max x-velocity : 0.036529 + Max y-velocity : 0.0404913 + Max z-velocity : 0.00716191 + Acceleration L2-norm : 1.22738 + Max X-acceleration : 0.613653 + Max Y-acceleration : 0.156555 + Max Z-acceleration : 4.91145 + Max x-acceleration : 0.965506 + Max y-acceleration : 2.30666 + Max z-acceleration : 0.314844 + step=144 time=1.44 + Displacement L2-norm : 0.0366598 + Max X-displacement : 0.00752794 + Max Y-displacement : 0.010449 + Max Z-displacement : 0.116796 + Max x-displacement : 0.00420398 + Max y-displacement : 0.0174913 + Max z-displacement : 0.00420784 + Velocity L2-norm : 0.134138 + Max X-velocity : 0.0191834 + Max Y-velocity : 0.0301781 + Max Z-velocity : 0.36983 + Max x-velocity : 0.0345016 + Max y-velocity : 0.0233195 + Max z-velocity : 0.00676572 + Acceleration L2-norm : 1.2009 + Max X-acceleration : 0.627278 + Max Y-acceleration : 0.165011 + Max Z-acceleration : 4.82954 + Max x-acceleration : 0.976246 + Max y-acceleration : 2.28179 + Max z-acceleration : 0.325019 + step=145 time=1.45 + Displacement L2-norm : 0.0353845 + Max X-displacement : 0.00768851 + Max Y-displacement : 0.0107461 + Max Z-displacement : 0.11372 + Max x-displacement : 0.00453097 + Max y-displacement : 0.0176334 + Max z-displacement : 0.00427128 + Velocity L2-norm : 0.140337 + Max X-velocity : 0.0145159 + Max Y-velocity : 0.0292407 + Max Z-velocity : 0.363778 + Max x-velocity : 0.0321768 + Max y-velocity : 0.00597141 + Max z-velocity : 0.00635341 + Acceleration L2-norm : 1.18051 + Max X-acceleration : 0.626286 + Max Y-acceleration : 0.173195 + Max Z-acceleration : 4.7105 + Max x-acceleration : 1.02467 + Max y-acceleration : 2.34014 + Max z-acceleration : 0.314829 + step=146 time=1.46 + Displacement L2-norm : 0.0340254 + Max X-displacement : 0.00778619 + Max Y-displacement : 0.0110331 + Max Z-displacement : 0.110173 + Max x-displacement : 0.00483245 + Max y-displacement : 0.0176013 + Max z-displacement : 0.00433043 + Velocity L2-norm : 0.146639 + Max X-velocity : 0.0169148 + Max Y-velocity : 0.0281583 + Max Z-velocity : 0.381209 + Max x-velocity : 0.0294889 + Max y-velocity : 0.0119412 + Max z-velocity : 0.00593913 + Acceleration L2-norm : 1.1445 + Max X-acceleration : 0.639066 + Max Y-acceleration : 0.180822 + Max Z-acceleration : 4.56946 + Max x-acceleration : 1.03013 + Max y-acceleration : 2.29052 + Max z-acceleration : 0.324275 + step=147 time=1.47 + Displacement L2-norm : 0.0325866 + Max X-displacement : 0.00782045 + Max Y-displacement : 0.0113086 + Max Z-displacement : 0.10617 + Max x-displacement : 0.00510542 + Max y-displacement : 0.0173946 + Max z-displacement : 0.00438519 + Velocity L2-norm : 0.152934 + Max X-velocity : 0.0211789 + Max Y-velocity : 0.0269445 + Max Z-velocity : 0.422707 + Max x-velocity : 0.026562 + Max y-velocity : 0.0293859 + Max z-velocity : 0.00555738 + Acceleration L2-norm : 1.11396 + Max X-acceleration : 0.631807 + Max Y-acceleration : 0.18646 + Max Z-acceleration : 4.39296 + Max x-acceleration : 1.0725 + Max y-acceleration : 2.32565 + Max z-acceleration : 0.313611 + step=148 time=1.48 + Displacement L2-norm : 0.0310726 + Max X-displacement : 0.00779127 + Max Y-displacement : 0.0115714 + Max Z-displacement : 0.101728 + Max x-displacement : 0.00534723 + Max y-displacement : 0.0170145 + Max z-displacement : 0.00443549 + Velocity L2-norm : 0.159119 + Max X-velocity : 0.0272314 + Max Y-velocity : 0.0256128 + Max Z-velocity : 0.46566 + Max x-velocity : 0.0233394 + Max y-velocity : 0.0466379 + Max z-velocity : 0.00516851 + Acceleration L2-norm : 1.06965 + Max X-acceleration : 0.635373 + Max Y-acceleration : 0.191528 + Max Z-acceleration : 4.19752 + Max x-acceleration : 1.07107 + Max y-acceleration : 2.25196 + Max z-acceleration : 0.32254 + step=149 time=1.49 + Displacement L2-norm : 0.0294885 + Max X-displacement : 0.00769904 + Max Y-displacement : 0.0118203 + Max Z-displacement : 0.0968672 + Max x-displacement : 0.00555557 + Max y-displacement : 0.0164631 + Max z-displacement : 0.00448127 + Velocity L2-norm : 0.165102 + Max X-velocity : 0.0331574 + Max Y-velocity : 0.0241807 + Max Z-velocity : 0.506501 + Max x-velocity : 0.0199514 + Max y-velocity : 0.0636527 + Max z-velocity : 0.00477786 + Acceleration L2-norm : 1.03076 + Max X-acceleration : 0.619074 + Max Y-acceleration : 0.19469 + Max Z-acceleration : 3.97064 + Max x-acceleration : 1.10542 + Max y-acceleration : 2.26493 + Max z-acceleration : 0.311082 + step=150 time=1.5 + Displacement L2-norm : 0.0278403 + Max X-displacement : 0.00754463 + Max Y-displacement : 0.0120546 + Max Z-displacement : 0.0916098 + Max x-displacement : 0.00572857 + Max y-displacement : 0.0157435 + Max z-displacement : 0.00452251 + Velocity L2-norm : 0.170797 + Max X-velocity : 0.0389186 + Max Y-velocity : 0.0226689 + Max Z-velocity : 0.544994 + Max x-velocity : 0.0163443 + Max y-velocity : 0.0802556 + Max z-velocity : 0.00438615 + Acceleration L2-norm : 0.980338 + Max X-acceleration : 0.614104 + Max Y-acceleration : 0.197027 + Max Z-acceleration : 3.72807 + Max x-acceleration : 1.09555 + Max y-acceleration : 2.16959 + Max z-acceleration : 0.319503 + Time integration completed. diff --git a/Test/Robaat-irr.xinp b/Test/Robaat-irr.xinp new file mode 100644 index 0000000..7e97e33 --- /dev/null +++ b/Test/Robaat-irr.xinp @@ -0,0 +1,40 @@ + + + + + + Small_Boat.dat + + + 418:449 794:825 + + + 450:517 533:611 737:793 856:945 + + + + + + + run18-irr482.asc + + + + + 2 + + + + 0 + + 0.01 + + + + + + + + + + diff --git a/Test/Robaat.reg b/Test/Robaat.reg new file mode 100644 index 0000000..1047fa8 --- /dev/null +++ b/Test/Robaat.reg @@ -0,0 +1,498 @@ +Robaat.xinp -dynamic -keep-previous-state -time 1.5 + +Input file: Robaat.xinp +Using the linear dynamics simulation driver. +Parsing input file Robaat.xinp +Parsing +Parsing + Parsing + Reading data file Small_Boat.dat + Nastran bulk data starting at line 71 +Parsing Nastran bulk data file succceeded. +Total number of nodes: 518 +Number of shell elements: 504 +Number of constraint elements: 0 +Number of other elements: 3 +Model extension (diameter): 3.5 + Reading patch 1 + Parsing + Topology sets: rib (1,1,5D) + wet (1,2,5D) + Parsing + Parsing +Parsing + Parsing + Gravitation vector: 0 0 -9.81 + Parsing + Sea surface function (expression): g=9.81; Tp=5.0; Hs=1.5; tstart=1.0; omega=2.0\*PI/Tp; wavespeed=g/omega; wavel=wavespeed/omega; refx=10.0-wavespeed\*(t-tstart); if(above(x+ux,refx),0.5\*Hs\*sin((x+ux-refx)/wavel),0.0) + Parsing + Pressure code 1 (hydrostatic): z0 = 0.1 rhow = 1000 +Parsing + alpha1 = 1 alpha2 = 0.03 + beta = 0.25 gamma = 0.5 +Parsing + Parsing + Parsing + Parsing + Sea surface \[-10 -2.5 0.1] to \[10 2.5 0.1], gridSize=0.5 +Parsing input file succeeded. +Equation solver: 2 +Number of Gauss points: 2 +Lagrangian basis functions are used +Problem definition: +Formulation: ANDES shell +Newmark predictor/multicorrector: beta = 0.25 gamma = 0.5 + using constant coefficient matrices + using zero acceleration predictor +Mass-proportional damping (alpha1): 1 +Stiffness-proportional damping (alpha2): 0.03 +Renumbered 518 nodes. +Resolving Dirichlet boundary conditions + >>> SAM model summary <<< +Number of elements 507 +Number of nodes 518 +Number of dofs 3108 +Detected 492 holes in the element numbers, reordering... +Number of unknowns 3108 +Total element set size: 358 + step=131 time=1.31 + Displacement L2-norm : 0.0295227 + Max X-displacement : 0.00429948 + Max Y-displacement : 0.00894124 + Max Z-displacement : 0.0723902 + Max x-displacement : 0.000971693 + Max y-displacement : 0.00033999 + Max z-displacement : 0.00365259 + Velocity L2-norm : 0.145108 + Max X-velocity : 0.107808 + Max Y-velocity : 0.018985 + Max Z-velocity : 0.534831 + Max x-velocity : 0.00591296 + Max y-velocity : 0.249759 + Max z-velocity : 0.00676321 + Acceleration L2-norm : 0.555647 + Max X-acceleration : 0.133943 + Max Y-acceleration : 0.0805338 + Max Z-acceleration : 1.83633 + Max x-acceleration : 0.86293 + Max y-acceleration : 0.972945 + Max z-acceleration : 0.312862 + step=132 time=1.32 + Displacement L2-norm : 0.0290918 + Max X-displacement : 0.00436702 + Max Y-displacement : 0.00910385 + Max Z-displacement : 0.0753266 + Max x-displacement : 0.000910449 + Max y-displacement : 0.00281297 + Max z-displacement : 0.00371759 + Velocity L2-norm : 0.144571 + Max X-velocity : 0.106078 + Max Y-velocity : 0.0187566 + Max Z-velocity : 0.539923 + Max x-velocity : 0.00640996 + Max y-velocity : 0.245654 + Max z-velocity : 0.00659086 + Acceleration L2-norm : 0.552264 + Max X-acceleration : 0.211929 + Max Y-acceleration : 0.0811573 + Max Z-acceleration : 2.0391 + Max x-acceleration : 0.887454 + Max y-acceleration : 1.08756 + Max z-acceleration : 0.323647 + step=133 time=1.33 + Displacement L2-norm : 0.0286777 + Max X-displacement : 0.00494231 + Max Y-displacement : 0.00926173 + Max Z-displacement : 0.0783461 + Max x-displacement : 0.000844195 + Max y-displacement : 0.00523654 + Max z-displacement : 0.00378047 + Velocity L2-norm : 0.143264 + Max X-velocity : 0.103655 + Max Y-velocity : 0.0184516 + Max Z-velocity : 0.541325 + Max x-velocity : 0.00691366 + Max y-velocity : 0.239952 + Max z-velocity : 0.00638696 + Acceleration L2-norm : 0.577416 + Max X-acceleration : 0.272669 + Max Y-acceleration : 0.0819649 + Max Z-acceleration : 2.22385 + Max x-acceleration : 0.844962 + Max y-acceleration : 1.29023 + Max z-acceleration : 0.315538 + step=134 time=1.34 + Displacement L2-norm : 0.0282824 + Max X-displacement : 0.00550242 + Max Y-displacement : 0.00941432 + Max Z-displacement : 0.0811435 + Max x-displacement : 0.000773392 + Max y-displacement : 0.00759473 + Max z-displacement : 0.00384093 + Velocity L2-norm : 0.141204 + Max X-velocity : 0.100558 + Max Y-velocity : 0.0180705 + Max Z-velocity : 0.539064 + Max x-velocity : 0.00731765 + Max y-velocity : 0.232646 + Max z-velocity : 0.00614813 + Acceleration L2-norm : 0.60072 + Max X-acceleration : 0.346811 + Max Y-acceleration : 0.0820953 + Max Z-acceleration : 2.40105 + Max x-acceleration : 0.868125 + Max y-acceleration : 1.39488 + Max z-acceleration : 0.325928 + step=135 time=1.35 + Displacement L2-norm : 0.0279065 + Max X-displacement : 0.00633953 + Max Y-displacement : 0.00956107 + Max Z-displacement : 0.0837013 + Max x-displacement : 0.00069856 + Max y-displacement : 0.00987212 + Max z-displacement : 0.0038987 + Velocity L2-norm : 0.138421 + Max X-velocity : 0.0968119 + Max Y-velocity : 0.0176149 + Max Z-velocity : 0.533201 + Max x-velocity : 0.00771714 + Max y-velocity : 0.223848 + Max z-velocity : 0.00587465 + Acceleration L2-norm : 0.645068 + Max X-acceleration : 0.402386 + Max Y-acceleration : 0.0824857 + Max Z-acceleration : 2.55641 + Max x-acceleration : 0.825178 + Max y-acceleration : 1.58673 + Max z-acceleration : 0.317732 + step=136 time=1.36 + Displacement L2-norm : 0.0275495 + Max X-displacement : 0.00728582 + Max Y-displacement : 0.00970149 + Max Z-displacement : 0.0860037 + Max x-displacement : 0.00062025 + Max y-displacement : 0.0120539 + Max z-displacement : 0.00395354 + Velocity L2-norm : 0.134947 + Max X-velocity : 0.0924459 + Max Y-velocity : 0.0170869 + Max Z-velocity : 0.523803 + Max x-velocity : 0.00801029 + Max y-velocity : 0.213569 + Max z-velocity : 0.00557782 + Acceleration L2-norm : 0.685412 + Max X-acceleration : 0.470811 + Max Y-acceleration : 0.0821758 + Max Z-acceleration : 2.70446 + Max x-acceleration : 0.847349 + Max y-acceleration : 1.67814 + Max z-acceleration : 0.327668 + step=137 time=1.37 + Displacement L2-norm : 0.0272096 + Max X-displacement : 0.00818551 + Max Y-displacement : 0.00983512 + Max Z-displacement : 0.0880362 + Max x-displacement : 0.000539054 + Max y-displacement : 0.0141257 + Max z-displacement : 0.0040052 + Velocity L2-norm : 0.130822 + Max X-velocity : 0.0874911 + Max Y-velocity : 0.0164886 + Max Z-velocity : 0.510946 + Max x-velocity : 0.00829088 + Max y-velocity : 0.201923 + Max z-velocity : 0.00525529 + Acceleration L2-norm : 0.739214 + Max X-acceleration : 0.52015 + Max Y-acceleration : 0.0821434 + Max Z-acceleration : 2.83016 + Max x-acceleration : 0.804004 + Max y-acceleration : 1.85804 + Max z-acceleration : 0.319303 + step=138 time=1.38 + Displacement L2-norm : 0.0268836 + Max X-displacement : 0.00903288 + Max Y-displacement : 0.00996939 + Max Z-displacement : 0.0897859 + Max x-displacement : 0.0004556 + Max y-displacement : 0.0160743 + Max z-displacement : 0.00405347 + Velocity L2-norm : 0.126102 + Max X-velocity : 0.0819839 + Max Y-velocity : 0.015823 + Max Z-velocity : 0.494742 + Max x-velocity : 0.00845831 + Max y-velocity : 0.188939 + Max z-velocity : 0.0049109 + Acceleration L2-norm : 0.786841 + Max X-acceleration : 0.581302 + Max Y-acceleration : 0.0813849 + Max Z-acceleration : 2.94568 + Max x-acceleration : 0.825159 + Max y-acceleration : 1.93271 + Max z-acceleration : 0.328666 + step=139 time=1.39 + Displacement L2-norm : 0.0265675 + Max X-displacement : 0.00982263 + Max Y-displacement : 0.0101236 + Max Z-displacement : 0.0912417 + Max x-displacement : 0.000441197 + Max y-displacement : 0.0178869 + Max z-displacement : 0.00409814 + Velocity L2-norm : 0.120857 + Max X-velocity : 0.0759652 + Max Y-velocity : 0.0150938 + Max Z-velocity : 0.475338 + Max x-velocity : 0.00860722 + Max y-velocity : 0.174762 + Max z-velocity : 0.00454075 + Acceleration L2-norm : 0.840844 + Max X-acceleration : 0.622436 + Max Y-acceleration : 0.0828756 + Max Z-acceleration : 3.03513 + Max x-acceleration : 0.793536 + Max y-acceleration : 2.09484 + Max z-acceleration : 0.320125 + step=140 time=1.4 + Displacement L2-norm : 0.0262564 + Max X-displacement : 0.0105498 + Max Y-displacement : 0.0102706 + Max Z-displacement : 0.0923943 + Max x-displacement : 0.000518314 + Max y-displacement : 0.0195519 + Max z-displacement : 0.00413904 + Velocity L2-norm : 0.11517 + Max X-velocity : 0.0694798 + Max Y-velocity : 0.0143048 + Max Z-velocity : 0.452904 + Max x-velocity : 0.00864093 + Max y-velocity : 0.159445 + Max z-velocity : 0.00415139 + Acceleration L2-norm : 0.887603 + Max X-acceleration : 0.674635 + Max Y-acceleration : 0.0824005 + Max Z-acceleration : 3.11293 + Max x-acceleration : 0.802513 + Max y-acceleration : 2.14778 + Max z-acceleration : 0.328889 + step=141 time=1.41 + Displacement L2-norm : 0.0259448 + Max X-displacement : 0.0112101 + Max Y-displacement : 0.0104094 + Max Z-displacement : 0.0932363 + Max x-displacement : 0.000595587 + Max y-displacement : 0.0210588 + Max z-displacement : 0.00417601 + Velocity L2-norm : 0.109143 + Max X-velocity : 0.0625754 + Max Y-velocity : 0.0134593 + Max Z-velocity : 0.427638 + Max x-velocity : 0.00865445 + Max y-velocity : 0.143151 + Max z-velocity : 0.00374425 + Acceleration L2-norm : 0.93569 + Max X-acceleration : 0.706249 + Max Y-acceleration : 0.0873356 + Max Z-acceleration : 3.16252 + Max x-acceleration : 0.797123 + Max y-acceleration : 2.28827 + Max z-acceleration : 0.320382 + step=142 time=1.42 + Displacement L2-norm : 0.0256266 + Max X-displacement : 0.0117995 + Max Y-displacement : 0.0105395 + Max Z-displacement : 0.0937623 + Max x-displacement : 0.000672398 + Max y-displacement : 0.0223982 + Max z-displacement : 0.00420889 + Velocity L2-norm : 0.102895 + Max X-velocity : 0.0553008 + Max Y-velocity : 0.0125609 + Max Z-velocity : 0.399745 + Max x-velocity : 0.00855361 + Max y-velocity : 0.125949 + Max z-velocity : 0.00332412 + Acceleration L2-norm : 0.976624 + Max X-acceleration : 0.748672 + Max Y-acceleration : 0.0923406 + Max Z-acceleration : 3.20071 + Max x-acceleration : 0.779963 + Max y-acceleration : 2.31687 + Max z-acceleration : 0.328684 + step=143 time=1.43 + Displacement L2-norm : 0.0252958 + Max X-displacement : 0.0123145 + Max Y-displacement : 0.0106604 + Max Z-displacement : 0.0939689 + Max x-displacement : 0.000748137 + Max y-displacement : 0.0235619 + Max z-displacement : 0.00423756 + Velocity L2-norm : 0.0965608 + Max X-velocity : 0.047706 + Max Y-velocity : 0.0116149 + Max Z-velocity : 0.369442 + Max x-velocity : 0.00843232 + Max y-velocity : 0.108005 + Max z-velocity : 0.00288942 + Acceleration L2-norm : 1.01595 + Max X-acceleration : 0.770276 + Max Y-acceleration : 0.0968647 + Max Z-acceleration : 3.21093 + Max x-acceleration : 0.800709 + Max y-acceleration : 2.43473 + Max z-acceleration : 0.319914 + step=144 time=1.44 + Displacement L2-norm : 0.0249462 + Max X-displacement : 0.0127523 + Max Y-displacement : 0.0107716 + Max Z-displacement : 0.0938548 + Max x-displacement : 0.000822194 + Max y-displacement : 0.0245429 + Max z-displacement : 0.0042619 + Velocity L2-norm : 0.0903061 + Max X-velocity : 0.0398428 + Max Y-velocity : 0.0106268 + Max Z-velocity : 0.336958 + Max x-velocity : 0.00819927 + Max y-velocity : 0.0893954 + Max z-velocity : 0.00244534 + Acceleration L2-norm : 1.0482 + Max X-acceleration : 0.802366 + Max Y-acceleration : 0.100755 + Max Z-acceleration : 3.34572 + Max x-acceleration : 0.776433 + Max y-acceleration : 2.43761 + Max z-acceleration : 0.327529 + step=145 time=1.45 + Displacement L2-norm : 0.0245721 + Max X-displacement : 0.0131103 + Max Y-displacement : 0.0108727 + Max Z-displacement : 0.0934203 + Max x-displacement : 0.000893983 + Max y-displacement : 0.0253354 + Max z-displacement : 0.00428183 + Velocity L2-norm : 0.0843211 + Max X-velocity : 0.0317614 + Max Y-velocity : 0.00960174 + Max Z-velocity : 0.302521 + Max x-velocity : 0.00794942 + Max y-velocity : 0.0702875 + Max z-velocity : 0.00199199 + Acceleration L2-norm : 1.07728 + Max X-acceleration : 0.813914 + Max Y-acceleration : 0.104249 + Max Z-acceleration : 3.54169 + Max x-acceleration : 0.803615 + Max y-acceleration : 2.53251 + Max z-acceleration : 0.318631 + step=146 time=1.46 + Displacement L2-norm : 0.0241679 + Max X-displacement : 0.0133867 + Max Y-displacement : 0.0109634 + Max Z-displacement : 0.092668 + Max x-displacement : 0.000962939 + Max y-displacement : 0.0259348 + Max z-displacement : 0.00429728 + Velocity L2-norm : 0.0788311 + Max X-velocity : 0.0235131 + Max Y-velocity : 0.00854504 + Max Z-velocity : 0.266367 + Max x-velocity : 0.00759389 + Max y-velocity : 0.0507562 + Max z-velocity : 0.00153403 + Acceleration L2-norm : 1.09973 + Max X-acceleration : 0.83576 + Max Y-acceleration : 0.10709 + Max Z-acceleration : 3.68899 + Max x-acceleration : 0.778588 + Max y-acceleration : 2.51004 + Max z-acceleration : 0.325707 + step=147 time=1.47 + Displacement L2-norm : 0.0237289 + Max X-displacement : 0.01358 + Max Y-displacement : 0.0110435 + Max Z-displacement : 0.091602 + Max x-displacement : 0.00102853 + Max y-displacement : 0.0263378 + Max z-displacement : 0.00430821 + Velocity L2-norm : 0.0740908 + Max X-velocity : 0.0151485 + Max Y-velocity : 0.00746326 + Max Z-velocity : 0.228735 + Max x-velocity : 0.00722645 + Max y-velocity : 0.0309687 + Max z-velocity : 0.0010715 + Acceleration L2-norm : 1.11772 + Max X-acceleration : 0.837156 + Max Y-acceleration : 0.109267 + Max Z-acceleration : 3.83759 + Max x-acceleration : 0.8055 + Max y-acceleration : 2.58193 + Max z-acceleration : 0.316515 + step=148 time=1.48 + Displacement L2-norm : 0.0232509 + Max X-displacement : 0.0136893 + Max Y-displacement : 0.0111126 + Max Z-displacement : 0.0902284 + Max x-displacement : 0.00109023 + Max y-displacement : 0.0265422 + Max z-displacement : 0.00431458 + Velocity L2-norm : 0.070372 + Max X-velocity : 0.00671743 + Max Y-velocity : 0.00636273 + Max Z-velocity : 0.190307 + Max x-velocity : 0.00675977 + Max y-velocity : 0.0109937 + Max z-velocity : 0.000609366 + Acceleration L2-norm : 1.13029 + Max X-acceleration : 0.849056 + Max Y-acceleration : 0.110838 + Max Z-acceleration : 3.93837 + Max x-acceleration : 0.779748 + Max y-acceleration : 2.53548 + Max z-acceleration : 0.322883 + step=149 time=1.49 + Displacement L2-norm : 0.0227307 + Max X-displacement : 0.0137142 + Max Y-displacement : 0.0111707 + Max Z-displacement : 0.0885547 + Max x-displacement : 0.00114757 + Max y-displacement : 0.0265469 + Max z-displacement : 0.00431641 + Velocity L2-norm : 0.0679332 + Max X-velocity : 0.00492448 + Max Y-velocity : 0.0052509 + Max Z-velocity : 0.182158 + Max x-velocity : 0.00628839 + Max y-velocity : 0.00953113 + Max z-velocity : 0.000151628 + Acceleration L2-norm : 1.13707 + Max X-acceleration : 0.840655 + Max Y-acceleration : 0.111527 + Max Z-acceleration : 4.03907 + Max x-acceleration : 0.806222 + Max y-acceleration : 2.58506 + Max z-acceleration : 0.313296 + step=150 time=1.5 + Displacement L2-norm : 0.0221661 + Max X-displacement : 0.0136549 + Max Y-displacement : 0.0112176 + Max Z-displacement : 0.0865903 + Max x-displacement : 0.00120125 + Max y-displacement : 0.0263521 + Max z-displacement : 0.00431374 + Velocity L2-norm : 0.0669709 + Max X-velocity : 0.0108183 + Max Y-velocity : 0.00505786 + Max Z-velocity : 0.210726 + Max x-velocity : 0.00572492 + Max y-velocity : 0.0294269 + Max z-velocity : 0.00050112 + Acceleration L2-norm : 1.13871 + Max X-acceleration : 0.842544 + Max Y-acceleration : 0.111515 + Max Z-acceleration : 4.08939 + Max x-acceleration : 0.779578 + Max y-acceleration : 2.51378 + Max z-acceleration : 0.318757 + Time integration completed. diff --git a/Test/Robaat.xinp b/Test/Robaat.xinp new file mode 100644 index 0000000..dbb0579 --- /dev/null +++ b/Test/Robaat.xinp @@ -0,0 +1,45 @@ + + + + + + Small_Boat.dat + + + 418:449 794:825 + + + 450:517 533:611 737:793 856:945 + + + + + + + + g=9.81; Tp=5.0; Hs=1.5; tstart=1.0; + omega=2.0*PI/Tp; wavespeed=g/omega; wavel=wavespeed/omega; + refx=10.0-wavespeed*(t-tstart); + if(above(x+ux,refx),0.5*Hs*sin((x+ux-refx)/wavel),0.0) + + + + + + 2 + + + + 0 + + 0.01 + + + + + + + + + + diff --git a/Test/Small_Boat.dat b/Test/Small_Boat.dat new file mode 100644 index 0000000..976ca7c --- /dev/null +++ b/Test/Small_Boat.dat @@ -0,0 +1,1668 @@ +$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ +$* +$* Simcenter v2412.0.0.1700 Translator +$* for Simcenter Nastran version 2412.0 +$* +$* FEM FILE: D:\Projects\dTHOR\SMALL_BOAT\Small_Boat_fem1.fem +$* SIM FILE: D:\Projects\dTHOR\SMALL_BOAT\Small_Boat_sim1.sim +$* ANALYSIS TYPE: Structural +$* SOLUTION NAME: Solution 1 +$* SOLUTION TYPE: SOL 103 Real Eigenvalues +$* +$* SOLVER INPUT FILE: Small_Boat.dat +$* CREATION DATE: 16-Oct-2025 +$* CREATION TIME: 15:58:51 +$* HOSTNAME: SINTEFPC10813 +$* NASTRAN LICENSE: Desktop Bundle +$* +$* UNITS: Meter (newton) +$* ... LENGTH : meter +$* ... TIME : sec +$* ... MASS : kilogram (kg) +$* ... TEMPERATURE : deg Celsius +$* ... FORCE : newton (N) +$* ... THERMAL ENERGY : N-m (joule) +$* +$* IMPORTANT NOTE: +$* This banner was generated by Simcenter and altering this +$* information may compromise the pre and post processing of results +$* +$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ +$* +$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ +$* +$* FILE MANAGEMENT +$* +$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ +$* +$* +$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ +$* +$* EXECUTIVE CONTROL +$* +$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ +$* +ID,NASTRAN,Small_Boat +SOL 103 +CEND +$* +$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ +$* +$* CASE CONTROL +$* +$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ +$* +ECHO = NONE +METHOD = 100 +OUTPUT +DISPLACEMENT(PLOT,REAL) = ALL +SPCFORCES(PLOT,REAL) = ALL +STRESS(PLOT,REAL,VONMISES,CENTER) = ALL +$* Step: Subcase - Normal Modes 1 +SUBCASE 1 + LABEL = Subcase - Normal Modes 1 +$* +$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ +$* +$* BULK DATA +$* +$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ +$* +BEGIN BULK +$* +$* SOLUTION CARDS +$* +$* Modeling Object: Real Eigenvalue - Lanczos1 +EIGRL 100 10 0 7 MASS +$* +$* PARAM CARDS +$* +PARAM K6ROT100.0000 +PARAM OIBULK YES +PARAM OMACHPR YES +PARAM POST -2 +PARAM POSTEXT YES +PARAM UNITSYS N-M +$* +$* GRID CARDS +$* +GRID* 168 0-9.100000000E-01-3.206917725E-01+ +* 4.1868663788E-02 0 +GRID* 172 0-2.647000000E+00-3.929727173E-01+ +* 7.0732658386E-02 0 +GRID* 173 0-1.778500000E+00-4.426176758E-01+ +* 2.4343343735E-02 0 +GRID* 333 0-9.100000000E-01-2.198462090E-01+ +* -3.802750811E-02 0 +GRID* 334 0-9.100000000E-01-1.118771179E-01+ +* -1.080980271E-01 0 +GRID* 342 0-1.778500000E+00-3.330545977E-01+ +* -2.621007141E-02 0 +GRID* 343 0-1.778500000E+00-2.225023031E-01+ +* -7.456509009E-02 0 +GRID* 344 0-1.778500000E+00-1.114081224E-01+ +* -1.216640968E-01 0 +GRID* 351 0-2.647000000E+00-2.948150081E-01+ +* 3.5194804897E-02 0 +GRID* 352 0-2.646999272E+00-1.966304790E-01+ +* -2.690303284E-04 0 +GRID* 353 0-2.647000000E+00-9.834014214E-02+ +* -3.543810082E-02 0 +GRID* 366 0-9.100000000E-010.0000000000E+00+ +* -1.717741394E-01 0 +GRID* 367 0-1.778500000E+000.0000000000E+00+ +* -1.680161133E-01 0 +GRID* 368 0-2.647000000E+000.0000000000E+00+ +* -7.046782684E-02 0 +GRID* 404 0-1.185833214E-01-7.049280204E-02+ +* 3.7061827461E-01 0 +GRID* 405 0-2.240167438E-01-1.353582149E-01+ +* 3.6400349536E-01 0 +GRID* 406 0-3.326482553E-01-1.946931883E-01+ +* 3.5718808886E-01 0 +GRID* 407 0-4.440454704E-01-2.486379000E-01+ +* 3.5019916899E-01 0 +GRID* 408 0-5.578134057E-01-2.973641901E-01+ +* 3.4306149237E-01 0 +GRID* 409 0-6.735961803E-01-3.410709891E-01+ +* 3.3579743442E-01 0 +GRID* 410 0-7.910813111E-01-3.799534616E-01+ +* 3.2842655892E-01 0 +GRID* 412 0-1.778500000E+00-4.929146423E-01+ +* 1.5298835278E-01 0 +GRID* 413 0-1.526714830E+00-5.255740319E-01+ +* 2.9085500486E-01 0 +GRID* 414 0-1.404336362E+00-5.116544273E-01+ +* 2.9640030513E-01 0 +GRID* 415 0-1.282460457E+00-4.938579535E-01+ +* 3.0191671374E-01 0 +GRID* 416 0-1.161219705E+00-4.721450967E-01+ +* 3.0741188733E-01 0 +GRID* 417 0-2.396736508E+00-5.234849023E-01+ +* 2.7722630262E-01 0 +GRID* 418 0-2.275330089E+00-5.356596362E-01+ +* 2.7759690672E-01 0 +GRID* 419 0-2.153562944E+00-5.434516780E-01+ +* 2.7797059296E-01 0 +GRID* 420 0-2.031600585E+00-5.470956326E-01+ +* 2.7834121936E-01 0 +GRID* 421 0-2.647000000E+00-4.385916380E-01+ +* 1.7490880670E-01 0 +GRID* 422 0-2.894620026E+00-4.222064904E-01+ +* 2.8555937770E-01 0 +GRID* 423 0-3.008594558E+00-3.826284583E-01+ +* 2.8851520089E-01 0 +GRID* 424 0-3.119326223E+00-3.347468106E-01+ +* 2.9138693011E-01 0 +GRID* 425 0-3.224988908E+00-2.765749711E-01+ +* 2.9412718537E-01 0 +GRID* 426 0-3.321790552E+00-2.047365308E-01+ +* 2.9663766094E-01 0 +GRID* 427 0-3.400494834E+00-1.137850231E-01+ +* 2.9867878250E-01 0 +GRID* 428 0-9.100000000E-01-4.141855774E-01+ +* 3.2096572876E-01 0 +GRID* 429 0-1.681310654E-020.0000000000E+00+ +* 3.7700320435E-01 0 +GRID* 430 0-1.778500000E+00-5.432116089E-01+ +* 2.8163336182E-01 0 +GRID* 431 0-1.649458130E+00-5.357761230E-01+ +* 2.8529586792E-01 0 +GRID* 432 0-1.040760986E+00-4.464409790E-01+ +* 3.1286700439E-01 0 +GRID* 433 0-2.517585205E+00-5.066571960E-01+ +* 2.7685693359E-01 0 +GRID* 434 0-1.909583984E+00-5.467756958E-01+ +* 2.7871499634E-01 0 +GRID* 435 0-3.435873779E+000.0000000000E+00+ +* 2.9959631348E-01 0 +GRID* 436 0-2.647000000E+00-4.840902100E-01+ +* 2.7913757324E-01 0 +GRID* 437 0-2.778475098E+00-4.548870850E-01+ +* 2.8272592163E-01 0 +GRID* 438 0-1.135342055E-01-7.913442109E-02+ +* 5.2354313497E-01 0 +GRID* 439 0-2.171671113E-01-1.513433925E-01+ +* 5.1023978375E-01 0 +GRID* 440 0-3.252461638E-01-2.168794963E-01+ +* 4.9776184846E-01 0 +GRID* 441 0-4.370366342E-01-2.760367748E-01+ +* 4.8610718779E-01 0 +GRID* 442 0-5.519008605E-01-3.291560720E-01+ +* 4.7526135452E-01 0 +GRID* 443 0-6.693028033E-01-3.765831564E-01+ +* 4.6520484952E-01 0 +GRID* 444 0-7.887895611E-01-4.186636718E-01+ +* 4.5591500263E-01 0 +GRID* 445 0-1.593644414E-020.0000000000E+00+ +* 4.5733231600E-01 0 +GRID* 446 0-1.031615149E+00-4.876963341E-01+ +* 4.3962198373E-01 0 +GRID* 447 0-1.154368817E+00-5.152263807E-01+ +* 4.3257838712E-01 0 +GRID* 448 0-1.278034823E+00-5.385251187E-01+ +* 4.2622428162E-01 0 +GRID* 449 0-1.402411589E+00-5.578718107E-01+ +* 4.2053732858E-01 0 +GRID* 450 0-1.527391767E+00-5.730936770E-01+ +* 4.1557220683E-01 0 +GRID* 451 0-1.652788554E+00-5.846575785E-01+ +* 4.1127069955E-01 0 +GRID* 452 0-1.902988657E+00-5.963690481E-01+ +* 4.0487235730E-01 0 +GRID* 453 0-2.027553574E+00-5.967746829E-01+ +* 4.0278376588E-01 0 +GRID* 454 0-2.152084314E+00-5.934131328E-01+ +* 4.0149876586E-01 0 +GRID* 455 0-2.276456397E+00-5.861917195E-01+ +* 4.0107715959E-01 0 +GRID* 456 0-2.400527925E+00-5.749404671E-01+ +* 4.0157052508E-01 0 +GRID* 457 0-2.524117231E+00-5.593198132E-01+ +* 4.0299355136E-01 0 +GRID* 458 0-2.763084046E+00-5.146299764E-01+ +* 4.0838953309E-01 0 +GRID* 459 0-2.877902886E+00-4.849940914E-01+ +* 4.1223319935E-01 0 +GRID* 460 0-2.990896447E+00-4.491258498E-01+ +* 4.1684189938E-01 0 +GRID* 461 0-3.101282474E+00-4.059983776E-01+ +* 4.2216162484E-01 0 +GRID* 462 0-3.207680556E+00-3.539044085E-01+ +* 4.2809990678E-01 0 +GRID* 463 0-3.307647284E+00-2.904511598E-01+ +* 4.3447065102E-01 0 +GRID* 464 0-3.396275782E+00-2.120879524E-01+ +* 4.4085490154E-01 0 +GRID* 465 0-3.463058944E+00-1.147168171E-01+ +* 4.4623718573E-01 0 +GRID* 466 0-9.100000000E-01-4.556706238E-01+ +* 4.4737646484E-01 0 +GRID* 467 0-1.520623016E-020.0000000000E+00+ +* 5.3766290283E-01 0 +GRID* 468 0-1.778500000E+00-5.924151001E-01+ +* 4.0768661499E-01 0 +GRID* 469 0-3.489314453E+000.0000000000E+00+ +* 4.4854382324E-01 0 +GRID* 470 0-2.647000000E+00-5.389627075E-01+ +* 4.0532687378E-01 0 +GRID* 471 0-2.524417443E+00-1.028339991E-01+ +* -5.516510583E-02 0 +GRID* 472 0-2.400888330E+00-1.063895791E-01+ +* -7.204230505E-02 0 +GRID* 473 0-2.276719834E+00-1.090348247E-01+ +* -8.624674694E-02 0 +GRID* 474 0-2.152175975E+00-1.107644732E-01+ +* -9.801935758E-02 0 +GRID* 475 0-2.027598970E+00-1.116269524E-01+ +* -1.077213544E-01 0 +GRID* 476 0-1.903023512E+00-1.118005647E-01+ +* -1.155492880E-01 0 +GRID* 477 0-1.902934171E+00-2.233195210E-01+ +* -6.960605853E-02 0 +GRID* 478 0-1.902993215E+00-3.348821824E-01+ +* -2.267622020E-02 0 +GRID* 479 0-2.027606445E+00-3.347509958E-01+ +* -1.753144694E-02 0 +GRID* 480 0-2.152303436E+00-3.321628455E-01+ +* -1.077794870E-02 0 +GRID* 481 0-2.276930593E+00-3.269352398E-01+ +* -2.291192729E-03 0 +GRID* 482 0-2.401028627E+00-3.189325666E-01+ +* 8.0632240250E-03 0 +GRID* 483 0-2.524414411E+00-3.080914123E-01+ +* 2.0394715338E-02 0 +GRID* 484 0-2.524382648E+00-2.052516920E-01+ +* -1.757177242E-02 0 +GRID* 485 0-2.276797441E+00-2.176714310E-01+ +* -4.450384965E-02 0 +GRID* 486 0-2.400876262E+00-2.123578564E-01+ +* -3.223722288E-02 0 +GRID* 487 0-2.152211384E+00-2.211794887E-01+ +* -5.468839157E-02 0 +GRID* 488 0-2.027626937E+00-2.229898940E-01+ +* -6.300486483E-02 0 +GRID* 489 0-2.524319164E+00-4.123608983E-01+ +* 5.9139231055E-02 0 +GRID* 490 0-2.400871464E+00-4.272832028E-01+ +* 4.9195251562E-02 0 +GRID* 491 0-2.276869223E+00-4.379978434E-01+ +* 4.0842828813E-02 0 +GRID* 492 0-2.152479898E+00-4.447322751E-01+ +* 3.4274199809E-02 0 +GRID* 493 0-2.027861545E+00-4.475955256E-01+ +* 2.9440442802E-02 0 +GRID* 494 0-1.903159015E+00-4.468752457E-01+ +* 2.6221842229E-02 0 +GRID* 495 0-2.524051387E+000.0000000000E+00+ +* -9.277310677E-02 0 +GRID* 496 0-2.400548634E+000.0000000000E+00+ +* -1.117887530E-01 0 +GRID* 497 0-2.276630521E+000.0000000000E+00+ +* -1.278751542E-01 0 +GRID* 498 0-2.152383815E+000.0000000000E+00+ +* -1.411878052E-01 0 +GRID* 499 0-2.027903529E+000.0000000000E+00+ +* -1.521079164E-01 0 +GRID* 500 0-1.903262046E+000.0000000000E+00+ +* -1.610095274E-01 0 +GRID* 501 0-2.276140823E+00-4.867368554E-01+ +* 1.5880273659E-01 0 +GRID* 502 0-2.522145353E+00-4.593478463E-01+ +* 1.6782742372E-01 0 +GRID* 503 0-2.399090497E+00-4.752590411E-01+ +* 1.6275797168E-01 0 +GRID* 504 0-2.029357317E+00-4.973066341E-01+ +* 1.5368466685E-01 0 +GRID* 505 0-1.905022263E+00-4.968792478E-01+ +* 1.5258574545E-01 0 +GRID* 506 0-2.152933881E+00-4.940132331E-01+ +* 1.5578476223E-01 0 +GRID* 507 0-1.039570849E+00-2.349519047E-01+ +* -4.478524337E-02 0 +GRID* 508 0-1.181581295E+00-2.476882117E-01+ +* -5.060432255E-02 0 +GRID* 509 0-1.298921486E+00-3.012975594E-01+ +* -2.868047084E-02 0 +GRID* 510 0-1.411241031E+00-3.114905408E-01+ +* -2.998372319E-02 0 +GRID* 511 0-1.531516820E+00-3.203045732E-01+ +* -3.023129993E-02 0 +GRID* 512 0-1.654555379E+00-3.275232569E-01+ +* -2.900574375E-02 0 +GRID* 513 0-1.655129495E+00-2.193869160E-01+ +* -7.854565160E-02 0 +GRID* 514 0-1.654882082E+00-1.102215397E-01+ +* -1.263492937E-01 0 +GRID* 515 0-1.531346070E+00-1.089143118E-01+ +* -1.294525926E-01 0 +GRID* 516 0-1.409548928E+00-1.077187091E-01+ +* -1.307169746E-01 0 +GRID* 517 0-1.289656382E+00-1.094864024E-01+ +* -1.286961895E-01 0 +GRID* 518 0-1.166025896E+00-1.190557034E-01+ +* -1.202977911E-01 0 +GRID* 519 0-1.037237710E+00-1.174528873E-01+ +* -1.145400399E-01 0 +GRID* 520 0-1.413605261E+00-2.115454234E-01+ +* -8.119502640E-02 0 +GRID* 521 0-1.532852982E+00-2.154518212E-01+ +* -8.088509457E-02 0 +GRID* 522 0-1.302597189E+00-2.093684946E-01+ +* -7.882911811E-02 0 +GRID* 523 0-1.653242005E+00-4.348821947E-01+ +* 2.3759446460E-02 0 +GRID* 524 0-1.528217504E+00-4.240082589E-01+ +* 2.4413012730E-02 0 +GRID* 525 0-1.403539310E+00-4.097889810E-01+ +* 2.6089347365E-02 0 +GRID* 526 0-1.279299453E+00-3.922471011E-01+ +* 2.8667150861E-02 0 +GRID* 527 0-1.155539522E+00-3.717526900E-01+ +* 3.2320637177E-02 0 +GRID* 528 0-1.032425434E+00-3.477928415E-01+ +* 3.6636159714E-02 0 +GRID* 529 0-1.654494831E+000.0000000000E+00+ +* -1.733722870E-01 0 +GRID* 530 0-1.530434555E+000.0000000000E+00+ +* -1.772491157E-01 0 +GRID* 531 0-1.406337917E+000.0000000000E+00+ +* -1.796996932E-01 0 +GRID* 532 0-1.282220356E+000.0000000000E+00+ +* -1.805668573E-01 0 +GRID* 533 0-1.158102842E+000.0000000000E+00+ +* -1.796853377E-01 0 +GRID* 534 0-1.034015167E+000.0000000000E+00+ +* -1.768588189E-01 0 +GRID* 541 0-4.442239331E-01-1.014285299E-01+ +* -3.390114485E-02 0 +GRID* 542 0-6.757558708E-01-1.061310193E-01+ +* -8.311656569E-02 0 +GRID* 543 0-5.581627983E-01-1.111734527E-01+ +* -5.612178050E-02 0 +GRID* 544 0-7.922778474E-01-1.078989451E-01+ +* -9.827799138E-02 0 +GRID* 545 0-7.954188651E-01-2.065420423E-01+ +* -2.841590444E-02 0 +GRID* 546 0-3.375190142E-01-8.834396004E-02+ +* -8.071639532E-03 0 +GRID* 547 0-6.864128597E-01-1.966486484E-01+ +* -1.245924740E-02 0 +GRID* 548 0-7.951097861E-01-2.919181492E-01+ +* 4.7224853858E-02 0 +GRID* 549 0-6.810435770E-01-2.601742375E-01+ +* 5.3348770423E-02 0 +GRID* 550 0-5.679272551E-01-2.252761375E-01+ +* 5.9910429926E-02 0 +GRID* 551 0-4.558846427E-01-1.871643181E-01+ +* 6.6988721094E-02 0 +GRID* 552 0-3.450631385E-01-1.457065356E-01+ +* 7.4466839161E-02 0 +GRID* 553 0-2.356266164E-01-1.007882787E-01+ +* 8.2339259167E-02 0 +GRID* 554 0-1.277629710E-01-5.227001531E-02+ +* 9.0541820939E-02 0 +GRID* 555 0-7.892795381E-010.0000000000E+00+ +* -1.642844208E-01 0 +GRID* 556 0-6.687915351E-010.0000000000E+00+ +* -1.537081868E-01 0 +GRID* 557 0-5.487038340E-010.0000000000E+00+ +* -1.393007533E-01 0 +GRID* 558 0-4.293386585E-010.0000000000E+00+ +* -1.198224742E-01 0 +GRID* 559 0-3.113797759E-010.0000000000E+00+ +* -9.317236624E-02 0 +GRID* 560 0-1.965971477E-010.0000000000E+00+ +* -5.526003311E-02 0 +GRID* 561 0-9.082152418E-020.0000000000E+00+ +* 2.6979527866E-03 0 +GRID* 562 0-2.169024849E-020.0000000000E+00+ +* 9.9032272339E-02 0 +GRID* 563 0-4.498738278E-01-2.182306399E-01+ +* 2.0797007624E-01 0 +GRID* 564 0-6.772150142E-01-3.008499072E-01+ +* 1.9389405186E-01 0 +GRID* 565 0-7.930104826E-01-3.361525358E-01+ +* 1.8731410007E-01 0 +GRID* 566 0-5.627692732E-01-2.615993389E-01+ +* 2.0080314537E-01 0 +GRID* 567 0-2.297465083E-01-1.184399113E-01+ +* 2.2282681086E-01 0 +GRID* 568 0-1.230420127E-01-6.167263902E-02+ +* 2.3041152280E-01 0 +GRID* 569 0-3.387784423E-01-1.705637539E-01+ +* 2.1532156003E-01 0 +GRID* 570 0-1.403992200E+00-4.606257513E-01+ +* 1.6079912017E-01 0 +GRID* 571 0-1.652257009E+00-4.854452468E-01+ +* 1.5449674359E-01 0 +GRID* 572 0-1.527737268E+00-4.747588628E-01+ +* 1.5728222642E-01 0 +GRID* 573 0-1.280744343E+00-4.430118215E-01+ +* 1.6484182680E-01 0 +GRID* 574 0-1.157822424E+00-4.218735481E-01+ +* 1.6944465677E-01 0 +GRID* 575 0-1.034749003E+00-3.969001023E-01+ +* 1.7473622960E-01 0 +GRID* 576 0-9.100000000E-01-3.677395814E-01+ +* 1.8131639677E-01 0 +GRID* 577 0-1.879513998E-020.0000000000E+00+ +* 2.3800971759E-01 0 +GRID* 578 0-2.855718929E+00-1.108973621E-01+ +* 1.3644174676E-02 0 +GRID* 579 0-2.851622412E+00-2.258945824E-01+ +* 5.1467175688E-02 0 +GRID* 580 0-2.749955650E+00-1.110839026E-01+ +* -1.032564281E-02 0 +GRID* 581 0-3.171136798E+00-1.149028271E-01+ +* 1.0590571899E-01 0 +GRID* 582 0-2.737376801E+00-2.397402063E-01+ +* 3.2132335531E-02 0 +GRID* 583 0-3.055259060E+00-1.409368497E-01+ +* 7.6674544675E-02 0 +GRID* 584 0-2.954223359E+00-1.064534802E-01+ +* 3.7449651755E-02 0 +GRID* 585 0-2.951687772E+00-2.086595169E-01+ +* 6.9634355750E-02 0 +GRID* 586 0-3.267659560E+00-9.437601659E-02+ +* 1.3548982042E-01 0 +GRID* 587 0-3.351693755E+00-1.015148380E-01+ +* 1.6950783496E-01 0 +GRID* 588 0-3.270138065E+00-1.757932859E-01+ +* 1.5483448285E-01 0 +GRID* 589 0-3.175207471E+00-2.322127932E-01+ +* 1.3907494464E-01 0 +GRID* 590 0-3.074324272E+00-2.773719277E-01+ +* 1.2360740687E-01 0 +GRID* 591 0-2.970117898E+00-3.145726637E-01+ +* 1.0887650296E-01 0 +GRID* 592 0-2.863767442E+00-3.455797730E-01+ +* 9.5101384493E-02 0 +GRID* 593 0-2.755931435E+00-3.714863831E-01+ +* 8.2346449495E-02 0 +GRID* 594 0-2.757045092E+000.0000000000E+00+ +* -4.704905400E-02 0 +GRID* 595 0-2.866313396E+000.0000000000E+00+ +* -2.024652313E-02 0 +GRID* 596 0-2.974597700E+000.0000000000E+00+ +* 1.0287551517E-02 0 +GRID* 597 0-3.081639535E+000.0000000000E+00+ +* 4.4921349029E-02 0 +GRID* 598 0-3.187131043E+000.0000000000E+00+ +* 8.4020682565E-02 0 +GRID* 599 0-3.290704953E+00-2.642402345E-14+ +* 1.2794297721E-01 0 +GRID* 600 0-3.391911865E+000.0000000000E+00+ +* 1.7706784058E-01 0 +GRID* 601 0-2.765128149E+00-4.134522101E-01+ +* 1.8189006293E-01 0 +GRID* 602 0-2.877887015E+00-3.838552277E-01+ +* 1.8941477519E-01 0 +GRID* 603 0-2.986651861E+00-3.483294130E-01+ +* 1.9601553396E-01 0 +GRID* 604 0-3.086945736E+00-3.056842468E-01+ +* 1.9695889959E-01 0 +GRID* 608 0-1.778500000E+00-9.980258989E-02+ +* 1.7037169519E-01 0 +GRID* 609 0-1.778500000E+00-2.030772059E-01+ +* 1.6633120794E-01 0 +GRID* 610 0-1.778500000E+00-3.075126497E-01+ +* 1.7125866004E-01 0 +GRID* 611 0-1.778500000E+00-4.026828441E-01+ +* 1.7022350194E-01 0 +GRID* 612 0-1.778500000E+00-3.844844716E-01+ +* 8.3294985777E-02 0 +GRID* 613 0-1.778500000E+00-3.039327920E-01+ +* 6.8264603738E-02 0 +GRID* 614 0-1.778500000E+00-1.901729020E-01+ +* 4.1036474243E-02 0 +GRID* 615 0-1.778500000E+00-9.211497558E-02+ +* -2.956452257E-02 0 +GRID* 616 0-1.778500000E+00-8.930707420E-02+ +* 6.6972431782E-02 0 +GRID* 628 0-1.778500000E+00-4.345692871E-01+ +* 2.8163336182E-01 0 +GRID* 629 0-1.778500000E+00-3.259269653E-01+ +* 2.8163336182E-01 0 +GRID* 630 0-1.778500000E+00-2.172846436E-01+ +* 2.8163336182E-01 0 +GRID* 631 0-1.778500000E+00-1.086423218E-01+ +* 2.8163336182E-01 0 +GRID* 632 0-1.778500000E+000.0000000000E+00+ +* -5.560374451E-02 0 +GRID* 633 0-1.778500000E+000.0000000000E+00+ +* 5.6808624268E-02 0 +GRID* 634 0-1.778500000E+000.0000000000E+00+ +* 1.6922099304E-01 0 +GRID* 642 0-1.778500000E+000.0000000000E+00+ +* 2.8163336182E-01 0 +GRID* 643 0-9.100000000E-01-1.199017526E-01+ +* 2.0202970891E-01 0 +GRID* 644 0-9.100000000E-01-2.413386574E-01+ +* 1.9724298495E-01 0 +GRID* 645 0-9.100000000E-01-1.735977062E-01+ +* 1.6257755685E-02 0 +GRID* 646 0-9.100000000E-01-2.076447803E-01+ +* 8.9052868080E-02 0 +GRID* 647 0-9.100000000E-01-1.046998209E-01+ +* 9.0726615423E-02 0 +GRID* 648 0-9.100000000E-01-9.726152528E-02+ +* -7.983462865E-03 0 +GRID* 649 0-9.100000000E-01-1.380618591E-01+ +* 3.2105260213E-01 0 +GRID* 650 0-9.100000000E-01-2.761237183E-01+ +* 3.2100916545E-01 0 +GRID* 651 0-9.100000000E-010.0000000000E+00+ +* 1.9787849426E-01 0 +GRID* 652 0-9.100000000E-010.0000000000E+00+ +* 7.4660949707E-02 0 +GRID* 653 0-9.100000000E-010.0000000000E+00+ +* -4.855659485E-02 0 +GRID* 654 0-9.100000000E-010.0000000000E+00+ +* 3.2109603882E-01 0 +GRID* 655 0-2.647000000E+00-8.899781882E-02+ +* 1.9257813172E-01 0 +GRID* 656 0-2.647000000E+00-1.813406005E-01+ +* 1.8814579327E-01 0 +GRID* 657 0-2.647000000E+00-2.750687921E-01+ +* 1.9181897042E-01 0 +GRID* 658 0-2.647000000E+00-3.590384023E-01+ +* 1.9022780570E-01 0 +GRID* 659 0-2.647000000E+00-3.427752245E-01+ +* 1.1988070379E-01 0 +GRID* 660 0-2.647000000E+00-2.721700702E-01+ +* 1.1005838971E-01 0 +GRID* 661 0-2.647000000E+00-1.698863967E-01+ +* 9.0490623981E-02 0 +GRID* 662 0-2.647000000E+00-8.254464837E-02+ +* 3.6434906548E-02 0 +GRID* 663 0-2.647000000E+00-8.017132580E-02+ +* 1.1201270128E-01 0 +GRID* 664 0-2.647000000E+00-9.681804199E-02+ +* 2.7921596680E-01 0 +GRID* 665 0-2.647000000E+00-1.936360840E-01+ +* 2.7919636841E-01 0 +GRID* 666 0-2.647000000E+00-2.904541260E-01+ +* 2.7917677002E-01 0 +GRID* 667 0-2.647000000E+00-3.872721680E-01+ +* 2.7915717163E-01 0 +GRID* 668 0-2.647000000E+000.0000000000E+00+ +* 1.6958021164E-02 0 +GRID* 669 0-2.647000000E+000.0000000000E+00+ +* 1.0438386917E-01 0 +GRID* 670 0-2.647000000E+000.0000000000E+00+ +* 1.9180971718E-01 0 +GRID* 671 0-2.647000000E+000.0000000000E+00+ +* 2.7923556519E-01 0 +GRID* 672 0-9.720010698E-01-2.713022689E-01+ +* 3.1707759975E-01 0 +GRID* 673 0-9.666131302E-01-3.533481345E-01+ +* 3.1741875787E-01 0 +GRID* 674 0-9.666067112E-01-1.909428286E-01+ +* 3.1735980928E-01 0 +GRID* 675 0-2.778475098E+00-9.097741699E-02+ +* 2.8272592163E-01 0 +GRID* 676 0-2.778475098E+00-1.819548340E-01+ +* 2.8272592163E-01 0 +GRID* 677 0-2.778475098E+00-2.729322510E-01+ +* 2.8272592163E-01 0 +GRID* 678 0-2.778475098E+00-3.639096680E-01+ +* 2.8272592163E-01 0 +GRID* 679 0-2.517585205E+00-4.053257568E-01+ +* 2.7685693359E-01 0 +GRID* 680 0-2.517585205E+00-3.039943176E-01+ +* 2.7685693359E-01 0 +GRID* 681 0-2.517585205E+00-2.026628784E-01+ +* 2.7685693359E-01 0 +GRID* 682 0-2.517585205E+00-1.013314392E-01+ +* 2.7685693359E-01 0 +GRID* 683 0-1.909583984E+00-1.093551392E-01+ +* 2.7871499634E-01 0 +GRID* 684 0-1.909583984E+00-2.187102783E-01+ +* 2.7871499634E-01 0 +GRID* 685 0-1.909583984E+00-3.280654175E-01+ +* 2.7871499634E-01 0 +GRID* 686 0-1.909583984E+00-4.374205566E-01+ +* 2.7871499634E-01 0 +GRID* 687 0-1.649458130E+00-4.286208984E-01+ +* 2.8529586792E-01 0 +GRID* 688 0-1.649458130E+00-3.214656738E-01+ +* 2.8529586792E-01 0 +GRID* 689 0-1.649458130E+00-2.143104492E-01+ +* 2.8529586792E-01 0 +GRID* 690 0-1.649458130E+00-1.071552246E-01+ +* 2.8529586792E-01 0 +GRID* 691 0-1.040760986E+00-8.928819580E-02+ +* 3.1286700439E-01 0 +GRID* 692 0-1.040760986E+00-1.785763916E-01+ +* 3.1286700439E-01 0 +GRID* 693 0-1.040760986E+00-2.678645874E-01+ +* 3.1286700439E-01 0 +GRID* 694 0-1.040760986E+00-3.571527832E-01+ +* 3.1286700439E-01 0 +GRID* 695 0-7.812894287E-01-2.514958700E-01+ +* 3.3011343384E-01 0 +GRID* 696 0-7.812894287E-01-1.257479350E-01+ +* 3.3011343384E-01 0 +GRID* 697 0-2.778475098E+002.1815882434E-14+ +* 2.8272592163E-01 0 +GRID* 698 0-2.517585205E+000.0000000000E+00+ +* 2.7685693359E-01 0 +GRID* 699 0-1.909583984E+000.0000000000E+00+ +* 2.7871499634E-01 0 +GRID* 700 0-1.649458130E+000.0000000000E+00+ +* 2.8529586792E-01 0 +GRID* 701 0-1.040760986E+000.0000000000E+00+ +* 3.1286700439E-01 0 +GRID* 702 0-7.812894287E-010.0000000000E+00+ +* 3.3011343384E-01 0 +GRID* 703 0-7.812894287E-01-3.772438049E-01+ +* 3.3011343384E-01 0 +GRID* 704 0-9.100000000E-013.2069177246E-01+ +* 4.1868663788E-02 0 +GRID* 705 0-2.647000000E+003.9297271729E-01+ +* 7.0732658386E-02 0 +GRID* 706 0-1.778500000E+004.4261767578E-01+ +* 2.4343343735E-02 0 +GRID* 707 0-9.100000000E-012.1984620895E-01+ +* -3.802750811E-02 0 +GRID* 708 0-9.100000000E-011.1187711792E-01+ +* -1.080980271E-01 0 +GRID* 709 0-1.778500000E+003.3305459772E-01+ +* -2.621007141E-02 0 +GRID* 710 0-1.778500000E+002.2250230306E-01+ +* -7.456509009E-02 0 +GRID* 711 0-1.778500000E+001.1140812238E-01+ +* -1.216640968E-01 0 +GRID* 712 0-2.647000000E+002.9481500809E-01+ +* 3.5194804897E-02 0 +GRID* 713 0-2.646999272E+001.9663047903E-01+ +* -2.690303284E-04 0 +GRID* 714 0-2.647000000E+009.8340142141E-02+ +* -3.543810082E-02 0 +GRID* 715 0-1.185833214E-017.0492802043E-02+ +* 3.7061827461E-01 0 +GRID* 716 0-2.240167438E-011.3535821493E-01+ +* 3.6400349536E-01 0 +GRID* 717 0-3.326482553E-011.9469318832E-01+ +* 3.5718808886E-01 0 +GRID* 718 0-4.440454704E-012.4863790000E-01+ +* 3.5019916899E-01 0 +GRID* 719 0-5.578134057E-012.9736419014E-01+ +* 3.4306149237E-01 0 +GRID* 720 0-6.735961803E-013.4107098911E-01+ +* 3.3579743442E-01 0 +GRID* 721 0-7.910813111E-013.7995346163E-01+ +* 3.2842655892E-01 0 +GRID* 722 0-1.778500000E+004.9291464233E-01+ +* 1.5298835278E-01 0 +GRID* 723 0-1.526714830E+005.2557403191E-01+ +* 2.9085500486E-01 0 +GRID* 724 0-1.404336362E+005.1165442734E-01+ +* 2.9640030513E-01 0 +GRID* 725 0-1.282460457E+004.9385795353E-01+ +* 3.0191671374E-01 0 +GRID* 726 0-1.161219705E+004.7214509670E-01+ +* 3.0741188733E-01 0 +GRID* 727 0-2.396736508E+005.2348490227E-01+ +* 2.7722630262E-01 0 +GRID* 728 0-2.275330089E+005.3565963622E-01+ +* 2.7759690672E-01 0 +GRID* 729 0-2.153562944E+005.4345167801E-01+ +* 2.7797059296E-01 0 +GRID* 730 0-2.031600585E+005.4709563259E-01+ +* 2.7834121936E-01 0 +GRID* 731 0-2.647000000E+004.3859163799E-01+ +* 1.7490880670E-01 0 +GRID* 732 0-2.894620026E+004.2220649036E-01+ +* 2.8555937770E-01 0 +GRID* 733 0-3.008594558E+003.8262845827E-01+ +* 2.8851520089E-01 0 +GRID* 734 0-3.119326223E+003.3474681055E-01+ +* 2.9138693011E-01 0 +GRID* 735 0-3.224988908E+002.7657497107E-01+ +* 2.9412718537E-01 0 +GRID* 736 0-3.321790552E+002.0473653083E-01+ +* 2.9663766094E-01 0 +GRID* 737 0-3.400494834E+001.1378502309E-01+ +* 2.9867878250E-01 0 +GRID* 738 0-9.100000000E-014.1418557739E-01+ +* 3.2096572876E-01 0 +GRID* 739 0-1.778500000E+005.4321160889E-01+ +* 2.8163336182E-01 0 +GRID* 740 0-1.649458130E+005.3577612305E-01+ +* 2.8529586792E-01 0 +GRID* 741 0-1.040760986E+004.4644097900E-01+ +* 3.1286700439E-01 0 +GRID* 742 0-2.517585205E+005.0665719604E-01+ +* 2.7685693359E-01 0 +GRID* 743 0-1.909583984E+005.4677569580E-01+ +* 2.7871499634E-01 0 +GRID* 744 0-2.647000000E+004.8409020996E-01+ +* 2.7913757324E-01 0 +GRID* 745 0-2.778475098E+004.5488708496E-01+ +* 2.8272592163E-01 0 +GRID* 746 0-1.135342055E-017.9134421086E-02+ +* 5.2354313497E-01 0 +GRID* 747 0-2.171671113E-011.5134339249E-01+ +* 5.1023978375E-01 0 +GRID* 748 0-3.252461638E-012.1687949628E-01+ +* 4.9776184846E-01 0 +GRID* 749 0-4.370366342E-012.7603677483E-01+ +* 4.8610718779E-01 0 +GRID* 750 0-5.519008605E-013.2915607200E-01+ +* 4.7526135452E-01 0 +GRID* 751 0-6.693028033E-013.7658315641E-01+ +* 4.6520484952E-01 0 +GRID* 752 0-7.887895611E-014.1866367176E-01+ +* 4.5591500263E-01 0 +GRID* 753 0-1.031615149E+004.8769633406E-01+ +* 4.3962198373E-01 0 +GRID* 754 0-1.154368817E+005.1522638072E-01+ +* 4.3257838712E-01 0 +GRID* 755 0-1.278034823E+005.3852511865E-01+ +* 4.2622428162E-01 0 +GRID* 756 0-1.402411589E+005.5787181066E-01+ +* 4.2053732858E-01 0 +GRID* 757 0-1.527391767E+005.7309367697E-01+ +* 4.1557220683E-01 0 +GRID* 758 0-1.652788554E+005.8465757846E-01+ +* 4.1127069955E-01 0 +GRID* 759 0-1.902988657E+005.9636904809E-01+ +* 4.0487235730E-01 0 +GRID* 760 0-2.027553574E+005.9677468291E-01+ +* 4.0278376588E-01 0 +GRID* 761 0-2.152084314E+005.9341313278E-01+ +* 4.0149876586E-01 0 +GRID* 762 0-2.276456397E+005.8619171946E-01+ +* 4.0107715959E-01 0 +GRID* 763 0-2.400527925E+005.7494046711E-01+ +* 4.0157052508E-01 0 +GRID* 764 0-2.524117231E+005.5931981316E-01+ +* 4.0299355136E-01 0 +GRID* 765 0-2.763084046E+005.1462997645E-01+ +* 4.0838953309E-01 0 +GRID* 766 0-2.877902886E+004.8499409136E-01+ +* 4.1223319935E-01 0 +GRID* 767 0-2.990896447E+004.4912584977E-01+ +* 4.1684189938E-01 0 +GRID* 768 0-3.101282474E+004.0599837758E-01+ +* 4.2216162484E-01 0 +GRID* 769 0-3.207680556E+003.5390440852E-01+ +* 4.2809990678E-01 0 +GRID* 770 0-3.307647284E+002.9045115977E-01+ +* 4.3447065102E-01 0 +GRID* 771 0-3.396275782E+002.1208795245E-01+ +* 4.4085490154E-01 0 +GRID* 772 0-3.463058944E+001.1471681712E-01+ +* 4.4623718573E-01 0 +GRID* 773 0-9.100000000E-014.5567062378E-01+ +* 4.4737646484E-01 0 +GRID* 774 0-1.778500000E+005.9241510010E-01+ +* 4.0768661499E-01 0 +GRID* 775 0-2.647000000E+005.3896270752E-01+ +* 4.0532687378E-01 0 +GRID* 776 0-2.524417443E+001.0283399910E-01+ +* -5.516510583E-02 0 +GRID* 777 0-2.400888330E+001.0638957912E-01+ +* -7.204230505E-02 0 +GRID* 778 0-2.276719834E+001.0903482465E-01+ +* -8.624674694E-02 0 +GRID* 779 0-2.152175975E+001.1076447324E-01+ +* -9.801935758E-02 0 +GRID* 780 0-2.027598970E+001.1162695239E-01+ +* -1.077213544E-01 0 +GRID* 781 0-1.903023512E+001.1180056470E-01+ +* -1.155492880E-01 0 +GRID* 782 0-1.902934171E+002.2331952099E-01+ +* -6.960605853E-02 0 +GRID* 783 0-1.902993215E+003.3488218236E-01+ +* -2.267622020E-02 0 +GRID* 784 0-2.027606445E+003.3475099581E-01+ +* -1.753144694E-02 0 +GRID* 785 0-2.152303436E+003.3216284552E-01+ +* -1.077794870E-02 0 +GRID* 786 0-2.276930593E+003.2693523980E-01+ +* -2.291192729E-03 0 +GRID* 787 0-2.401028627E+003.1893256661E-01+ +* 8.0632240250E-03 0 +GRID* 788 0-2.524414411E+003.0809141229E-01+ +* 2.0394715338E-02 0 +GRID* 789 0-2.524382648E+002.0525169197E-01+ +* -1.757177242E-02 0 +GRID* 790 0-2.276797441E+002.1767143098E-01+ +* -4.450384965E-02 0 +GRID* 791 0-2.400876262E+002.1235785644E-01+ +* -3.223722288E-02 0 +GRID* 792 0-2.152211384E+002.2117948871E-01+ +* -5.468839157E-02 0 +GRID* 793 0-2.027626937E+002.2298989395E-01+ +* -6.300486483E-02 0 +GRID* 794 0-2.524319164E+004.1236089827E-01+ +* 5.9139231055E-02 0 +GRID* 795 0-2.400871464E+004.2728320283E-01+ +* 4.9195251562E-02 0 +GRID* 796 0-2.276869223E+004.3799784345E-01+ +* 4.0842828813E-02 0 +GRID* 797 0-2.152479898E+004.4473227513E-01+ +* 3.4274199809E-02 0 +GRID* 798 0-2.027861545E+004.4759552561E-01+ +* 2.9440442802E-02 0 +GRID* 799 0-1.903159015E+004.4687524571E-01+ +* 2.6221842229E-02 0 +GRID* 800 0-2.276140823E+004.8673685540E-01+ +* 1.5880273659E-01 0 +GRID* 801 0-2.522145353E+004.5934784632E-01+ +* 1.6782742372E-01 0 +GRID* 802 0-2.399090497E+004.7525904110E-01+ +* 1.6275797168E-01 0 +GRID* 803 0-2.029357317E+004.9730663405E-01+ +* 1.5368466685E-01 0 +GRID* 804 0-1.905022263E+004.9687924778E-01+ +* 1.5258574545E-01 0 +GRID* 805 0-2.152933881E+004.9401323314E-01+ +* 1.5578476223E-01 0 +GRID* 806 0-1.039570849E+002.3495190472E-01+ +* -4.478524337E-02 0 +GRID* 807 0-1.181581295E+002.4768821171E-01+ +* -5.060432255E-02 0 +GRID* 808 0-1.298921486E+003.0129755936E-01+ +* -2.868047084E-02 0 +GRID* 809 0-1.411241031E+003.1149054084E-01+ +* -2.998372319E-02 0 +GRID* 810 0-1.531516820E+003.2030457323E-01+ +* -3.023129993E-02 0 +GRID* 811 0-1.654555379E+003.2752325692E-01+ +* -2.900574375E-02 0 +GRID* 812 0-1.655129495E+002.1938691596E-01+ +* -7.854565160E-02 0 +GRID* 813 0-1.654882082E+001.1022153969E-01+ +* -1.263492937E-01 0 +GRID* 814 0-1.531346070E+001.0891431176E-01+ +* -1.294525926E-01 0 +GRID* 815 0-1.409548928E+001.0771870909E-01+ +* -1.307169746E-01 0 +GRID* 816 0-1.289656382E+001.0948640237E-01+ +* -1.286961895E-01 0 +GRID* 817 0-1.166025896E+001.1905570341E-01+ +* -1.202977911E-01 0 +GRID* 818 0-1.037237710E+001.1745288734E-01+ +* -1.145400399E-01 0 +GRID* 819 0-1.413605261E+002.1154542339E-01+ +* -8.119502640E-02 0 +GRID* 820 0-1.532852982E+002.1545182116E-01+ +* -8.088509457E-02 0 +GRID* 821 0-1.302597189E+002.0936849461E-01+ +* -7.882911811E-02 0 +GRID* 822 0-1.653242005E+004.3488219465E-01+ +* 2.3759446460E-02 0 +GRID* 823 0-1.528217504E+004.2400825887E-01+ +* 2.4413012730E-02 0 +GRID* 824 0-1.403539310E+004.0978898103E-01+ +* 2.6089347365E-02 0 +GRID* 825 0-1.279299453E+003.9224710113E-01+ +* 2.8667150861E-02 0 +GRID* 826 0-1.155539522E+003.7175269005E-01+ +* 3.2320637177E-02 0 +GRID* 827 0-1.032425434E+003.4779284146E-01+ +* 3.6636159714E-02 0 +GRID* 828 0-4.442239331E-011.0142852989E-01+ +* -3.390114485E-02 0 +GRID* 829 0-6.757558708E-011.0613101934E-01+ +* -8.311656569E-02 0 +GRID* 830 0-5.581627983E-011.1117345271E-01+ +* -5.612178050E-02 0 +GRID* 831 0-7.922778474E-011.0789894507E-01+ +* -9.827799138E-02 0 +GRID* 832 0-7.954188651E-012.0654204225E-01+ +* -2.841590444E-02 0 +GRID* 833 0-3.375190142E-018.8343960041E-02+ +* -8.071639532E-03 0 +GRID* 834 0-6.864128597E-011.9664864837E-01+ +* -1.245924740E-02 0 +GRID* 835 0-7.951097861E-012.9191814919E-01+ +* 4.7224853858E-02 0 +GRID* 836 0-6.810435770E-012.6017423745E-01+ +* 5.3348770423E-02 0 +GRID* 837 0-5.679272551E-012.2527613752E-01+ +* 5.9910429926E-02 0 +GRID* 838 0-4.558846427E-011.8716431814E-01+ +* 6.6988721094E-02 0 +GRID* 839 0-3.450631385E-011.4570653563E-01+ +* 7.4466839161E-02 0 +GRID* 840 0-2.356266164E-011.0078827866E-01+ +* 8.2339259167E-02 0 +GRID* 841 0-1.277629710E-015.2270015305E-02+ +* 9.0541820939E-02 0 +GRID* 842 0-4.498738278E-012.1823063986E-01+ +* 2.0797007624E-01 0 +GRID* 843 0-6.772150142E-013.0084990723E-01+ +* 1.9389405186E-01 0 +GRID* 844 0-7.930104826E-013.3615253582E-01+ +* 1.8731410007E-01 0 +GRID* 845 0-5.627692732E-012.6159933889E-01+ +* 2.0080314537E-01 0 +GRID* 846 0-2.297465083E-011.1843991128E-01+ +* 2.2282681086E-01 0 +GRID* 847 0-1.230420127E-016.1672639018E-02+ +* 2.3041152280E-01 0 +GRID* 848 0-3.387784423E-011.7056375389E-01+ +* 2.1532156003E-01 0 +GRID* 849 0-1.403992200E+004.6062575131E-01+ +* 1.6079912017E-01 0 +GRID* 850 0-1.652257009E+004.8544524679E-01+ +* 1.5449674359E-01 0 +GRID* 851 0-1.527737268E+004.7475886284E-01+ +* 1.5728222642E-01 0 +GRID* 852 0-1.280744343E+004.4301182150E-01+ +* 1.6484182680E-01 0 +GRID* 853 0-1.157822424E+004.2187354810E-01+ +* 1.6944465677E-01 0 +GRID* 854 0-1.034749003E+003.9690010229E-01+ +* 1.7473622960E-01 0 +GRID* 855 0-9.100000000E-013.6773958143E-01+ +* 1.8131639677E-01 0 +GRID* 856 0-2.855718929E+001.1089736205E-01+ +* 1.3644174676E-02 0 +GRID* 857 0-2.851622412E+002.2589458236E-01+ +* 5.1467175688E-02 0 +GRID* 858 0-2.749955650E+001.1108390260E-01+ +* -1.032564281E-02 0 +GRID* 859 0-3.171136798E+001.1490282712E-01+ +* 1.0590571899E-01 0 +GRID* 860 0-2.737376801E+002.3974020630E-01+ +* 3.2132335531E-02 0 +GRID* 861 0-3.055259060E+001.4093684971E-01+ +* 7.6674544675E-02 0 +GRID* 862 0-2.954223359E+001.0645348016E-01+ +* 3.7449651755E-02 0 +GRID* 863 0-2.951687772E+002.0865951686E-01+ +* 6.9634355750E-02 0 +GRID* 864 0-3.267659560E+009.4376016594E-02+ +* 1.3548982042E-01 0 +GRID* 865 0-3.351693755E+001.0151483798E-01+ +* 1.6950783496E-01 0 +GRID* 866 0-3.270138065E+001.7579328592E-01+ +* 1.5483448285E-01 0 +GRID* 867 0-3.175207471E+002.3221279325E-01+ +* 1.3907494464E-01 0 +GRID* 868 0-3.074324272E+002.7737192769E-01+ +* 1.2360740687E-01 0 +GRID* 869 0-2.970117898E+003.1457266370E-01+ +* 1.0887650296E-01 0 +GRID* 870 0-2.863767442E+003.4557977297E-01+ +* 9.5101384493E-02 0 +GRID* 871 0-2.755931435E+003.7148638311E-01+ +* 8.2346449495E-02 0 +GRID* 872 0-2.765128149E+004.1345221010E-01+ +* 1.8189006293E-01 0 +GRID* 873 0-2.877887015E+003.8385522767E-01+ +* 1.8941477519E-01 0 +GRID* 874 0-2.986651861E+003.4832941299E-01+ +* 1.9601553396E-01 0 +GRID* 875 0-3.086945736E+003.0568424680E-01+ +* 1.9695889959E-01 0 +GRID* 876 0-1.778500000E+009.9802589893E-02+ +* 1.7037169519E-01 0 +GRID* 877 0-1.778500000E+002.0307720591E-01+ +* 1.6633120794E-01 0 +GRID* 878 0-1.778500000E+003.0751264972E-01+ +* 1.7125866004E-01 0 +GRID* 879 0-1.778500000E+004.0268284405E-01+ +* 1.7022350194E-01 0 +GRID* 880 0-1.778500000E+003.8448447161E-01+ +* 8.3294985777E-02 0 +GRID* 881 0-1.778500000E+003.0393279201E-01+ +* 6.8264603738E-02 0 +GRID* 882 0-1.778500000E+001.9017290200E-01+ +* 4.1036474243E-02 0 +GRID* 883 0-1.778500000E+009.2114975578E-02+ +* -2.956452257E-02 0 +GRID* 884 0-1.778500000E+008.9307074201E-02+ +* 6.6972431782E-02 0 +GRID* 885 0-1.778500000E+004.3456928711E-01+ +* 2.8163336182E-01 0 +GRID* 886 0-1.778500000E+003.2592696533E-01+ +* 2.8163336182E-01 0 +GRID* 887 0-1.778500000E+002.1728464355E-01+ +* 2.8163336182E-01 0 +GRID* 888 0-1.778500000E+001.0864232178E-01+ +* 2.8163336182E-01 0 +GRID* 889 0-9.100000000E-011.1990175263E-01+ +* 2.0202970891E-01 0 +GRID* 890 0-9.100000000E-012.4133865743E-01+ +* 1.9724298495E-01 0 +GRID* 891 0-9.100000000E-011.7359770623E-01+ +* 1.6257755685E-02 0 +GRID* 892 0-9.100000000E-012.0764478033E-01+ +* 8.9052868080E-02 0 +GRID* 893 0-9.100000000E-011.0469982088E-01+ +* 9.0726615423E-02 0 +GRID* 894 0-9.100000000E-019.7261525284E-02+ +* -7.983462865E-03 0 +GRID* 895 0-9.100000000E-011.3806185913E-01+ +* 3.2105260213E-01 0 +GRID* 896 0-9.100000000E-012.7612371826E-01+ +* 3.2100916545E-01 0 +GRID* 897 0-2.647000000E+008.8997818822E-02+ +* 1.9257813172E-01 0 +GRID* 898 0-2.647000000E+001.8134060050E-01+ +* 1.8814579327E-01 0 +GRID* 899 0-2.647000000E+002.7506879210E-01+ +* 1.9181897042E-01 0 +GRID* 900 0-2.647000000E+003.5903840228E-01+ +* 1.9022780570E-01 0 +GRID* 901 0-2.647000000E+003.4277522445E-01+ +* 1.1988070379E-01 0 +GRID* 902 0-2.647000000E+002.7217007022E-01+ +* 1.1005838971E-01 0 +GRID* 903 0-2.647000000E+001.6988639666E-01+ +* 9.0490623981E-02 0 +GRID* 904 0-2.647000000E+008.2544648374E-02+ +* 3.6434906548E-02 0 +GRID* 905 0-2.647000000E+008.0171325803E-02+ +* 1.1201270128E-01 0 +GRID* 906 0-2.647000000E+009.6818041992E-02+ +* 2.7921596680E-01 0 +GRID* 907 0-2.647000000E+001.9363608398E-01+ +* 2.7919636841E-01 0 +GRID* 908 0-2.647000000E+002.9045412598E-01+ +* 2.7917677002E-01 0 +GRID* 909 0-2.647000000E+003.8727216797E-01+ +* 2.7915717163E-01 0 +GRID* 910 0-9.720010698E-012.7130226890E-01+ +* 3.1707759975E-01 0 +GRID* 911 0-9.666131302E-013.5334813449E-01+ +* 3.1741875787E-01 0 +GRID* 912 0-9.666067112E-011.9094282864E-01+ +* 3.1735980928E-01 0 +GRID* 913 0-2.778475098E+009.0977416992E-02+ +* 2.8272592163E-01 0 +GRID* 914 0-2.778475098E+001.8195483398E-01+ +* 2.8272592163E-01 0 +GRID* 915 0-2.778475098E+002.7293225098E-01+ +* 2.8272592163E-01 0 +GRID* 916 0-2.778475098E+003.6390966797E-01+ +* 2.8272592163E-01 0 +GRID* 917 0-2.517585205E+004.0532575684E-01+ +* 2.7685693359E-01 0 +GRID* 918 0-2.517585205E+003.0399431763E-01+ +* 2.7685693359E-01 0 +GRID* 919 0-2.517585205E+002.0266287842E-01+ +* 2.7685693359E-01 0 +GRID* 920 0-2.517585205E+001.0133143921E-01+ +* 2.7685693359E-01 0 +GRID* 921 0-1.909583984E+001.0935513916E-01+ +* 2.7871499634E-01 0 +GRID* 922 0-1.909583984E+002.1871027832E-01+ +* 2.7871499634E-01 0 +GRID* 923 0-1.909583984E+003.2806541748E-01+ +* 2.7871499634E-01 0 +GRID* 924 0-1.909583984E+004.3742055664E-01+ +* 2.7871499634E-01 0 +GRID* 925 0-1.649458130E+004.2862089844E-01+ +* 2.8529586792E-01 0 +GRID* 926 0-1.649458130E+003.2146567383E-01+ +* 2.8529586792E-01 0 +GRID* 927 0-1.649458130E+002.1431044922E-01+ +* 2.8529586792E-01 0 +GRID* 928 0-1.649458130E+001.0715522461E-01+ +* 2.8529586792E-01 0 +GRID* 929 0-1.040760986E+008.9288195801E-02+ +* 3.1286700439E-01 0 +GRID* 930 0-1.040760986E+001.7857639160E-01+ +* 3.1286700439E-01 0 +GRID* 931 0-1.040760986E+002.6786458740E-01+ +* 3.1286700439E-01 0 +GRID* 932 0-1.040760986E+003.5715278320E-01+ +* 3.1286700439E-01 0 +GRID* 933 0-7.812894287E-012.5149586995E-01+ +* 3.3011343384E-01 0 +GRID* 934 0-7.812894287E-011.2574793498E-01+ +* 3.3011343384E-01 0 +GRID* 935 0-7.812894287E-013.7724380493E-01+ +* 3.3011343384E-01 0 +$* +$* ELEMENT CARDS +$* +$* Mesh Collector: Thin Shell(Alu8mm) +$* Mesh: 2d_mesh(1) +CQUAD4 450 1 471 353 368 495 +CQUAD4 451 1 471 495 496 472 +CQUAD4 452 1 472 496 497 473 +CQUAD4 453 1 473 497 498 474 +CQUAD4 454 1 474 498 499 475 +CQUAD4 455 1 475 499 500 476 +CQUAD4 456 1 476 500 367 344 +CQUAD4 457 1 476 344 343 477 +CQUAD4 458 1 477 343 342 478 +CQUAD4 459 1 478 342 173 494 +CQUAD4 460 1 478 494 493 479 +CQUAD4 461 1 479 493 492 480 +CQUAD4 462 1 480 492 491 481 +CQUAD4 463 1 481 491 490 482 +CQUAD4 464 1 482 490 489 483 +CQUAD4 465 1 483 489 172 351 +CQUAD4 466 1 483 351 352 484 +CQUAD4 467 1 484 352 353 471 +CQUAD4 468 1 486 482 483 484 +CQUAD4 469 1 486 484 471 472 +CQUAD4 470 1 485 481 482 486 +CQUAD4 471 1 485 486 472 473 +CQUAD4 472 1 487 480 481 485 +CQUAD4 473 1 487 485 473 474 +CQUAD4 474 1 488 475 476 477 +CQUAD4 475 1 488 477 478 479 +CQUAD4 476 1 487 474 475 488 +CQUAD4 477 1 487 488 479 480 +CQUAD4 492 1 507 333 168 528 +CQUAD4 493 1 507 528 527 508 +CQUAD4 494 1 508 527 526 509 +CQUAD4 495 1 509 526 525 510 +CQUAD4 496 1 510 525 524 511 +CQUAD4 497 1 511 524 523 512 +CQUAD4 498 1 512 523 173 342 +CQUAD4 499 1 512 342 343 513 +CQUAD4 500 1 513 343 344 514 +CQUAD4 501 1 514 344 367 529 +CQUAD4 502 1 514 529 530 515 +CQUAD4 503 1 515 530 531 516 +CQUAD4 504 1 516 531 532 517 +CQUAD4 505 1 517 532 533 518 +CQUAD4 506 1 518 533 534 519 +CQUAD4 507 1 519 534 366 334 +CQUAD4 508 1 519 334 333 507 +CQUAD4 509 1 518 519 507 508 +CQUAD4 510 1 508 522 517 518 +CQUAD4 511 1 520 516 517 522 +CQUAD4 512 1 520 522 509 510 +CQUAD4 513 1 521 511 512 513 +CQUAD4 514 1 521 513 514 515 +CQUAD4 515 1 520 510 511 521 +CQUAD4 516 1 520 521 515 516 +CTRIA3 517 1 508 509 522 +CQUAD4 533 1 561 554 553 560 +CQUAD4 534 1 553 546 559 560 +CQUAD4 535 1 541 558 559 546 +CQUAD4 536 1 541 546 552 551 +CQUAD4 537 1 543 542 556 557 +CQUAD4 538 1 543 557 558 541 +CQUAD4 539 1 543 541 551 550 +CQUAD4 540 1 550 547 542 543 +CQUAD4 541 1 544 334 366 555 +CQUAD4 542 1 544 555 556 542 +CQUAD4 543 1 545 547 549 548 +CQUAD4 544 1 545 548 168 333 +CQUAD4 545 1 544 542 547 545 +CQUAD4 546 1 544 545 333 334 +CTRIA3 547 1 561 562 554 +CTRIA3 548 1 553 552 546 +CTRIA3 549 1 550 549 547 +CQUAD4 580 1 580 594 368 353 +CQUAD4 581 1 580 578 595 594 +CQUAD4 582 1 582 579 578 580 +CQUAD4 583 1 582 593 592 579 +CQUAD4 584 1 582 351 172 593 +CQUAD4 585 1 582 580 353 352 +CQUAD4 586 1 578 584 596 595 +CQUAD4 587 1 585 579 592 591 +CQUAD4 588 1 585 584 578 579 +CQUAD4 589 1 584 583 597 596 +CQUAD4 590 1 585 591 590 583 +CQUAD4 591 1 581 583 590 589 +CQUAD4 592 1 581 598 597 583 +CQUAD4 593 1 587 600 599 586 +CQUAD4 594 1 581 586 599 598 +CQUAD4 595 1 581 589 588 586 +CTRIA3 596 1 582 352 351 +CTRIA3 597 1 585 583 584 +CTRIA3 598 1 587 586 588 +$* Mesh Collector: Thin Shell(Alu8mm) +$* Mesh: 2d_mesh(1)_reflected(1) +CQUAD4 856 1 792 779 778 790 90.2262 +CQUAD4 857 1 793 782 781 780 88.9994 +CQUAD4 858 1 793 784 783 782 90.9937 +CQUAD4 859 1 819 809 808 821 89.2459 +CQUAD4 860 1 820 812 811 810 88.5269 +CTRIA3 861 1 807 821 808 47.3587 +CQUAD4 862 1 792 785 784 793 90.6769 +CQUAD4 863 1 817 807 806 818 95.4892 +CQUAD4 864 1 807 817 816 821 74.8168 +CQUAD4 865 1 783 784 798 799 89.0359 +CQUAD4 866 1 784 785 797 798 89.7888 +CQUAD4 867 1 785 786 796 797 90.6974 +CQUAD4 868 1 786 787 795 796 91.7344 +CQUAD4 869 1 787 788 794 795 92.7296 +CQUAD4 870 1 788 712 705 794 93.4147 +CQUAD4 871 1 788 789 713 712 86.6337 +CQUAD4 872 1 789 776 714 713 88.9990 +CQUAD4 873 1 831 708 707 832 89.7237 +CTRIA3 874 1 561 841 562 62.7115 +CTRIA3 875 1 840 833 839 45.8130 +CTRIA3 876 1 837 834 836 39.9805 +CQUAD4 877 1 780 781 500 499 88.8921 +CQUAD4 878 1 783 799 706 709 91.5048 +CQUAD4 879 1 791 789 788 787 90.5961 +CQUAD4 880 1 792 790 786 785 89.7347 +CQUAD4 881 1 792 793 780 779 89.3630 +CQUAD4 882 1 807 808 825 826 74.1584 +CQUAD4 883 1 811 812 710 709 92.5883 +CQUAD4 884 1 818 806 707 708 91.6746 +CQUAD4 885 1 819 821 816 815 87.4982 +CQUAD4 886 1 820 814 813 812 92.8721 +CQUAD4 887 1 561 560 840 841 95.5334 +CQUAD4 888 1 840 560 559 833 66.2273 +CQUAD4 889 1 828 833 559 558 87.4069 +CQUAD4 890 1 828 838 839 833 89.1647 +CQUAD4 891 1 830 557 556 829 84.2589 +CQUAD4 892 1 830 828 558 557 88.9660 +CQUAD4 893 1 812 813 711 710 92.1865 +CQUAD4 894 1 813 529 367 711 91.5433 +CQUAD4 895 1 813 814 530 529 88.6942 +CQUAD4 896 1 814 815 531 530 88.8092 +CQUAD4 897 1 815 816 532 531 89.6178 +CQUAD4 898 1 816 817 533 532 92.1492 +CQUAD4 899 1 817 818 534 533 87.1083 +CQUAD4 900 1 818 708 366 534 87.7693 +CQUAD4 901 1 819 820 810 809 89.4758 +CQUAD4 902 1 819 815 814 820 93.7815 +CQUAD4 903 1 837 830 829 834 63.0482 +CQUAD4 904 1 831 555 366 708 90.5308 +CQUAD4 905 1 808 809 824 825 96.3150 +CQUAD4 906 1 809 810 823 824 90.3229 +CQUAD4 907 1 810 811 822 823 88.3863 +CQUAD4 908 1 811 709 706 822 87.7745 +CQUAD4 909 1 858 714 368 594 82.6171 +CQUAD4 910 1 859 864 866 867 94.7378 +CQUAD4 911 1 806 827 704 707 90.7572 +CQUAD4 912 1 806 807 826 827 90.1839 +CQUAD4 913 1 832 835 836 834 88.2500 +CQUAD4 914 1 832 707 704 835 88.4240 +CQUAD4 915 1 858 594 595 856 90.4983 +CQUAD4 916 1 860 858 856 857 81.3176 +CQUAD4 917 1 860 857 870 871 85.6235 +CQUAD4 918 1 860 871 705 712 67.6207 +CQUAD4 919 1 860 713 714 858 66.1723 +CQUAD4 920 1 856 595 596 862 86.4655 +CQUAD4 921 1 863 869 870 857 93.5554 +CQUAD4 922 1 863 857 856 862 97.4384 +CQUAD4 923 1 862 596 597 861102.1806 +CQUAD4 924 1 863 861 868 869110.5914 +CQUAD4 925 1 859 867 868 861 83.8044 +CQUAD4 926 1 859 861 597 598106.9777 +CQUAD4 927 1 865 864 599 600107.1337 +CQUAD4 928 1 859 598 599 864 74.5177 +CTRIA3 929 1 860 712 713 59.2676 +CQUAD4 930 1 776 495 368 714 91.3682 +CQUAD4 931 1 776 777 496 495 88.6819 +CQUAD4 932 1 777 778 497 496 88.6877 +CQUAD4 933 1 778 779 498 497 88.7740 +CQUAD4 934 1 779 780 499 498 88.8535 +CQUAD4 935 1 781 711 367 500 88.8883 +CQUAD4 936 1 781 782 710 711 91.1963 +CQUAD4 937 1 782 783 709 710 91.2621 +CQUAD4 938 1 791 777 776 789 89.3069 +CQUAD4 939 1 790 791 787 786 90.1820 +CQUAD4 940 1 790 778 777 791 89.7815 +CQUAD4 941 1 830 837 838 828 89.0260 +CQUAD4 942 1 831 829 556 555 91.8018 +CQUAD4 943 1 831 832 834 829 87.9490 +CTRIA3 944 1 863 862 861 58.4179 +CTRIA3 945 1 865 866 864 47.6120 +$* Mesh Collector: Thin Shell(Alu6mm) +$* Mesh: 2d_mesh(2) +CQUAD4 478 2 502 433 436 421 +CQUAD4 479 2 502 421 172 489 +CQUAD4 480 2 503 501 418 417 +CQUAD4 481 2 503 417 433 502 +CQUAD4 482 2 490 491 501 503 +CQUAD4 483 2 490 503 502 489 +CQUAD4 484 2 505 504 493 494 +CQUAD4 485 2 505 494 173 412 +CQUAD4 486 2 434 420 504 505 +CQUAD4 487 2 434 505 412 430 +CQUAD4 488 2 506 501 491 492 +CQUAD4 489 2 506 492 493 504 +CQUAD4 490 2 419 418 501 506 +CQUAD4 491 2 419 506 504 420 +CQUAD4 550 2 565 564 409 410 +CQUAD4 551 2 565 410 428 576 +CQUAD4 552 2 548 549 564 565 +CQUAD4 553 2 548 565 576 168 +CQUAD4 554 2 566 563 407 408 +CQUAD4 555 2 566 408 409 564 +CQUAD4 556 2 550 551 563 566 +CQUAD4 557 2 550 566 564 549 +CQUAD4 558 2 568 567 553 554 +CQUAD4 559 2 568 554 562 577 +CQUAD4 560 2 404 405 567 568 +CQUAD4 561 2 404 568 577 429 +CQUAD4 562 2 569 563 551 552 +CQUAD4 563 2 569 552 553 567 +CQUAD4 564 2 406 407 563 569 +CQUAD4 565 2 406 569 567 405 +CQUAD4 566 2 571 431 430 412 +CQUAD4 567 2 571 412 173 523 +CQUAD4 568 2 572 570 414 413 +CQUAD4 569 2 572 413 431 571 +CQUAD4 570 2 524 525 570 572 +CQUAD4 571 2 524 572 571 523 +CQUAD4 572 2 573 415 414 570 +CQUAD4 573 2 573 570 525 526 +CQUAD4 574 2 574 416 415 573 +CQUAD4 575 2 574 573 526 527 +CQUAD4 576 2 575 576 428 432 +CQUAD4 577 2 575 432 416 574 +CQUAD4 578 2 528 168 576 575 +CQUAD4 579 2 528 575 574 527 +CQUAD4 599 2 601 593 172 421 +CQUAD4 600 2 601 421 436 437 +CQUAD4 601 2 602 592 593 601 +CQUAD4 602 2 602 601 437 422 +CQUAD4 603 2 603 423 424 604 +CQUAD4 604 2 603 604 590 591 +CQUAD4 605 2 602 422 423 603 +CQUAD4 606 2 602 603 591 592 +CQUAD4 607 2 427 435 600 587 +CQUAD4 608 2 427 587 588 426 +CQUAD4 609 2 425 426 588 589 +CQUAD4 610 2 589 604 424 425 +CTRIA3 611 2 589 590 604 +$* Mesh Collector: Thin Shell(Alu6mm) +$* Mesh: 2d_mesh(3) +CQUAD4 418 2 444 466 428 410 +CQUAD4 419 2 444 410 409 443 +CQUAD4 420 2 408 407 441 442 +CQUAD4 421 2 408 442 443 409 +CQUAD4 422 2 440 441 407 406 +CQUAD4 423 2 440 406 405 439 +CQUAD4 424 2 438 439 405 404 +CQUAD4 425 2 445 438 404 429 +CTRIA3 426 2 445 467 438 +CQUAD4 427 2 413 414 449 450 +CQUAD4 428 2 451 468 430 431 +CQUAD4 429 2 451 431 413 450 +CQUAD4 430 2 448 449 414 415 +CQUAD4 431 2 447 448 415 416 +CQUAD4 432 2 432 428 466 446 +CQUAD4 433 2 432 446 447 416 +CQUAD4 434 2 417 418 455 456 +CQUAD4 435 2 457 470 436 433 +CQUAD4 436 2 457 433 417 456 +CQUAD4 437 2 454 455 418 419 +CQUAD4 438 2 453 454 419 420 +CQUAD4 439 2 434 430 468 452 +CQUAD4 440 2 434 452 453 420 +CQUAD4 441 2 465 469 435 427 +CQUAD4 442 2 465 427 426 464 +CQUAD4 443 2 425 424 461 462 +CQUAD4 444 2 463 426 425 462 +CQUAD4 445 2 437 436 470 458 +CQUAD4 446 2 422 437 458 459 +CQUAD4 447 2 460 461 424 423 +CQUAD4 448 2 460 423 422 459 +CTRIA3 449 2 463 464 426 +$* Mesh Collector: Thin Shell(Alu6mm) +$* Mesh: 2d_mesh(2)_reflected(1) +CQUAD4 737 2 845 843 720 719 90.4821 +CQUAD4 738 2 837 845 842 838 89.3831 +CQUAD4 739 2 847 841 840 846 86.5685 +CQUAD4 740 2 847 577 562 841 93.3566 +CQUAD4 741 2 845 719 718 842 89.9260 +CQUAD4 742 2 837 836 843 845 90.8090 +CQUAD4 743 2 715 847 846 716 87.3287 +CQUAD4 744 2 850 822 706 722 90.1815 +CQUAD4 745 2 851 723 724 849 90.3926 +CQUAD4 746 2 851 850 740 723 89.8604 +CQUAD4 747 2 823 851 849 824 91.4097 +CQUAD4 748 2 823 822 850 851 88.7003 +CQUAD4 749 2 852 849 724 725 88.2608 +CQUAD4 750 2 852 825 824 849 91.6518 +CQUAD4 751 2 853 852 725 726 87.3584 +CQUAD4 752 2 853 826 825 852 92.2303 +CQUAD4 753 2 854 741 738 855 93.9159 +CQUAD4 754 2 854 853 726 741 86.1304 +CQUAD4 755 2 827 854 855 704 92.7841 +CQUAD4 756 2 827 826 853 854 87.3009 +CQUAD4 757 2 850 722 739 740 90.5649 +CQUAD4 758 2 801 731 744 742 92.9635 +CQUAD4 759 2 801 794 705 731 88.1800 +CQUAD4 760 2 802 727 728 800 88.5439 +CQUAD4 761 2 802 801 742 727 91.7716 +CQUAD4 762 2 795 802 800 796 90.7966 +CQUAD4 763 2 795 794 801 802 89.2799 +CQUAD4 764 2 804 799 798 803 91.3291 +CQUAD4 765 2 804 722 706 799 88.7242 +CQUAD4 766 2 743 804 803 730 91.8184 +CQUAD4 767 2 743 739 722 804 88.6877 +CQUAD4 768 2 805 797 796 800 90.2273 +CQUAD4 769 2 805 803 798 797 89.4738 +CQUAD4 770 2 729 805 800 728 88.7382 +CQUAD4 771 2 729 730 803 805 90.5301 +CQUAD4 772 2 844 721 720 843 91.2687 +CQUAD4 773 2 844 855 738 721 89.0087 +CQUAD4 774 2 835 844 843 836 91.0163 +CQUAD4 775 2 835 704 855 844 89.0789 +CQUAD4 776 2 737 865 600 435 91.5863 +CQUAD4 777 2 737 736 866 865 79.7243 +CQUAD4 778 2 735 867 866 736 95.7265 +CQUAD4 779 2 867 735 734 875 68.8287 +CQUAD4 780 2 872 731 705 871 86.8258 +CQUAD4 781 2 872 745 744 731 95.3508 +CQUAD4 782 2 873 872 871 870 84.4119 +CQUAD4 783 2 873 732 745 872 96.9169 +CQUAD4 784 2 874 875 734 733 85.7052 +CQUAD4 785 2 874 869 868 875 91.6676 +CQUAD4 786 2 873 874 733 732 84.7861 +CTRIA3 787 2 867 875 868 37.8933 +CQUAD4 788 2 715 429 577 847 92.4567 +CQUAD4 789 2 848 839 838 842 88.2102 +CQUAD4 790 2 848 846 840 839 91.5133 +CQUAD4 791 2 717 848 842 718 88.8266 +CQUAD4 792 2 717 716 846 848 90.6785 +CQUAD4 793 2 873 870 869 874 93.8878 +$* Mesh Collector: Thin Shell(Alu6mm) +$* Mesh: 2d_mesh(3)_reflected(1) +CQUAD4 794 2 752 721 738 773 90.2281 +CQUAD4 795 2 752 751 720 721 89.4283 +CQUAD4 796 2 719 750 749 718 89.7871 +CQUAD4 797 2 719 720 751 750 90.8773 +CQUAD4 798 2 748 717 718 749 86.3666 +CQUAD4 799 2 748 747 716 717 93.5182 +CQUAD4 800 2 746 715 716 747 84.3370 +CQUAD4 801 2 445 429 715 746117.3685 +CTRIA3 802 2 445 746 467 62.6179 +CQUAD4 803 2 723 757 756 724 90.1839 +CQUAD4 804 2 758 740 739 774 91.1655 +CQUAD4 805 2 758 757 723 740 88.5165 +CQUAD4 806 2 755 725 724 756 88.6480 +CQUAD4 807 2 754 726 725 755 87.8094 +CQUAD4 808 2 741 753 773 738 87.1676 +CQUAD4 809 2 741 726 754 753 92.5519 +CQUAD4 810 2 727 763 762 728 89.2602 +CQUAD4 811 2 764 742 744 775 90.0845 +CQUAD4 812 2 764 763 727 742 89.4488 +CQUAD4 813 2 761 729 728 762 87.9403 +CQUAD4 814 2 760 730 729 761 87.1486 +CQUAD4 815 2 743 759 774 739 86.5995 +CQUAD4 816 2 743 730 760 759 92.8943 +CQUAD4 817 2 772 737 435 469 95.7771 +CQUAD4 818 2 772 771 736 737 75.1857 +CQUAD4 819 2 735 769 768 734 71.5032 +CQUAD4 820 2 770 769 735 736108.0548 +CQUAD4 821 2 745 765 775 744 79.8711 +CQUAD4 822 2 732 766 765 745 77.7900 +CQUAD4 823 2 767 733 734 768 76.1424 +CQUAD4 824 2 767 766 732 733102.6901 +CTRIA3 825 2 770 736 771 68.4493 +$* Mesh Collector: Thin Shell(Alu10mm) +$* Mesh: 2d_mesh(5) +CQUAD4 664 2 664 671 697 675 +CQUAD4 665 2 665 664 675 676 +CQUAD4 666 2 666 665 676 677 +CQUAD4 667 2 678 437 436 667 +CQUAD4 668 2 678 667 666 677 +CQUAD4 669 2 682 698 671 664 +CQUAD4 670 2 681 682 664 665 +CQUAD4 671 2 680 681 665 666 +CQUAD4 672 2 667 436 433 679 +CQUAD4 673 2 667 679 680 666 +CQUAD4 674 2 631 642 699 683 +CQUAD4 675 2 631 683 684 630 +CQUAD4 676 2 629 630 684 685 +CQUAD4 677 2 686 434 430 628 +CQUAD4 678 2 686 628 629 685 +CQUAD4 679 2 690 700 642 631 +CQUAD4 680 2 690 631 630 689 +CQUAD4 681 2 688 689 630 629 +CQUAD4 682 2 628 430 431 687 +CQUAD4 683 2 628 687 688 629 +CQUAD4 684 2 649 654 701 691 +CQUAD4 685 2 694 432 428 673 +CQUAD4 686 2 428 650 672 673 +CQUAD4 687 2 694 673 672 693 +CQUAD4 688 2 649 691 692 674 +CQUAD4 689 2 692 693 672 674 +CQUAD4 690 2 649 674 672 650 +CQUAD4 691 2 696 702 654 649 +CQUAD4 692 2 650 428 703 695 +CQUAD4 693 2 650 695 696 649 +$* Mesh Collector: Thin Shell(Alu10mm) +$* Mesh: 2d_mesh(5)_reflected(1) +CQUAD4 826 2 906 913 697 671 87.4570 +CQUAD4 827 2 907 914 913 906 84.9242 +CQUAD4 828 2 908 915 914 907 82.4113 +CQUAD4 829 2 916 909 744 745 79.9276 +CQUAD4 830 2 916 915 908 909100.0724 +CQUAD4 831 2 920 906 671 698 88.0029 +CQUAD4 832 2 919 907 906 920 86.0107 +CQUAD4 833 2 918 908 907 919 84.0280 +CQUAD4 834 2 909 917 742 744 82.0594 +CQUAD4 835 2 909 908 918 917 97.9406 +CQUAD4 836 2 888 921 699 642 90.3115 +CQUAD4 837 2 888 887 922 921 89.6885 +CQUAD4 838 2 886 923 922 887 90.9344 +CQUAD4 839 2 924 885 739 743 91.2458 +CQUAD4 840 2 924 923 886 885 88.7542 +CQUAD4 841 2 928 888 642 700 90.6600 +CQUAD4 842 2 928 927 887 888 89.3400 +CQUAD4 843 2 926 886 887 927 91.9793 +CQUAD4 844 2 885 925 740 739 92.6382 +CQUAD4 845 2 885 886 926 925 87.3618 +CQUAD4 846 2 895 929 701 654 69.5823 +CQUAD4 847 2 932 911 738 741 92.9319 +CQUAD4 848 2 738 911 910 896 42.9973 +CQUAD4 849 2 932 931 910 911 87.0681 +CQUAD4 850 2 895 912 930 929 63.4093 +CQUAD4 851 2 930 912 910 931 80.5491 +CQUAD4 852 2 895 896 910 912 47.0085 +CQUAD4 853 2 934 895 654 702 95.4515 +CQUAD4 854 2 896 933 935 738100.8071 +CQUAD4 855 2 896 895 934 933 79.1929 +$* Mesh Collector: Thin Shell(Alu12mm) +$* Mesh: 2d_mesh(4) +CQUAD4 621 2 608 634 642 631 +CQUAD4 622 2 608 631 630 609 +CQUAD4 623 2 609 630 629 610 +CQUAD4 624 2 610 629 628 611 +CQUAD4 625 2 611 628 430 412 +CQUAD4 626 2 611 412 173 612 +CQUAD4 627 2 612 173 342 613 +CQUAD4 628 2 613 342 343 614 +CQUAD4 629 2 614 343 344 615 +CQUAD4 630 2 615 344 367 632 +CQUAD4 631 2 615 632 633 616 +CQUAD4 632 2 616 633 634 608 +CQUAD4 633 2 614 616 608 609 +CQUAD4 634 2 610 611 612 613 +CQUAD4 635 2 610 613 614 609 +CTRIA3 636 2 614 615 616 +CQUAD4 637 2 649 654 651 643 +CQUAD4 638 2 644 576 428 650 +CQUAD4 639 2 644 650 649 643 +CQUAD4 640 2 168 576 644 646 +CQUAD4 641 2 168 646 645 333 +CQUAD4 642 2 647 643 651 652 +CQUAD4 643 2 645 648 334 333 +CQUAD4 644 2 653 366 334 648 +CQUAD4 645 2 653 648 647 652 +CQUAD4 646 2 647 648 645 646 +CQUAD4 647 2 646 644 643 647 +CQUAD4 648 2 655 670 671 664 +CQUAD4 649 2 655 664 665 656 +CQUAD4 650 2 656 665 666 657 +CQUAD4 651 2 657 666 667 658 +CQUAD4 652 2 658 667 436 421 +CQUAD4 653 2 658 421 172 659 +CQUAD4 654 2 659 172 351 660 +CQUAD4 655 2 660 351 352 661 +CQUAD4 656 2 661 352 353 662 +CQUAD4 657 2 662 353 368 668 +CQUAD4 658 2 662 668 669 663 +CQUAD4 659 2 663 669 670 655 +CQUAD4 660 2 661 663 655 656 +CQUAD4 661 2 657 658 659 660 +CQUAD4 662 2 657 660 661 656 +CTRIA3 663 2 661 662 663 +$* Mesh Collector: Thin Shell(Alu12mm) +$* Mesh: 2d_mesh(4)_reflected(1) +CQUAD4 694 2 895 889 651 654 81.3429 +CQUAD4 695 2 890 896 738 855 81.4832 +CQUAD4 696 2 890 889 895 896103.4410 +CQUAD4 697 2 704 892 890 855 85.9887 +CQUAD4 698 2 704 707 891 892 61.0434 +CQUAD4 699 2 893 652 651 889106.5011 +CQUAD4 700 2 891 707 708 894112.8117 +CQUAD4 701 2 653 894 708 366112.6437 +CQUAD4 702 2 893 892 891 894 93.3779 +CQUAD4 703 2 892 893 889 890106.3668 +CQUAD4 704 2 899 902 901 900 90.9449 +CTRIA3 705 2 882 884 883 50.1739 +CQUAD4 706 2 876 888 642 634 95.2032 +CQUAD4 707 2 876 877 887 888 87.6979 +CQUAD4 708 2 877 878 886 887 80.2742 +CQUAD4 709 2 878 879 885 886 81.1515 +CQUAD4 710 2 879 722 739 885 84.8422 +CQUAD4 711 2 879 880 706 722 91.0102 +CQUAD4 712 2 880 881 709 706124.0301 +CQUAD4 713 2 881 882 710 709 93.6715 +CQUAD4 714 2 882 883 711 710 69.8707 +CQUAD4 715 2 883 632 367 711 86.0467 +CQUAD4 716 2 883 884 633 632104.1186 +CQUAD4 717 2 884 876 634 633102.2887 +CQUAD4 718 2 882 877 876 884 81.4600 +CQUAD4 719 2 878 881 880 879 91.3675 +CQUAD4 720 2 878 877 882 881 85.3080 +CQUAD4 721 2 653 652 893 894 67.3563 +CQUAD4 722 2 899 898 903 902 85.7252 +CTRIA3 723 2 903 905 904 45.2433 +CQUAD4 724 2 897 906 671 670 95.6524 +CQUAD4 725 2 897 898 907 906 87.5903 +CQUAD4 726 2 898 899 908 907 80.0650 +CQUAD4 727 2 899 900 909 908 81.0972 +CQUAD4 728 2 900 731 744 909 83.2857 +CQUAD4 729 2 900 901 705 731 92.1176 +CQUAD4 730 2 901 902 712 705127.6853 +CQUAD4 731 2 902 903 713 712 95.9993 +CQUAD4 732 2 903 904 714 713 74.6653 +CQUAD4 733 2 904 668 368 714 89.1184 +CQUAD4 734 2 904 905 669 668101.4778 +CQUAD4 735 2 905 897 670 669101.6879 +CQUAD4 736 2 903 898 897 905 83.1998 +$* +$* PROPERTY CARDS +$* +$* Property: PSHELL1::Small_Boat_fem1::[1] +PSHELL 1 10.008000 11.000000 10.833333 0.0000 +$* Property: PSHELL2::Small_Boat_fem1::[1] +PSHELL 2 10.010000 1 12.0000 10.833333 0.0000 +$* +$* MATERIAL CARDS +$* +$* Material: Aluminum_A356::Small_Boat_fem1::[1] +MAT1 17.000+10 0.330000 2670.002.1400-5 +$* +CONM2,997,642,,100.0 +CONM2,998,654,,80.0 +CONM2,999,671,,60.0 +ENDDATA diff --git a/Test/run18-irr482.asc b/Test/run18-irr482.asc new file mode 100644 index 0000000..9548c71 --- /dev/null +++ b/Test/run18-irr482.asc @@ -0,0 +1,483 @@ +# i A_i omega_i epsilon_i + 1 1.37461E-04 3.65632E+00 1.54521E+00 + 2 1.42152E-04 3.66792E+00 4.07520E+00 + 3 1.55325E-04 3.69962E+00 5.76430E-01 + 4 1.57619E-04 3.70503E+00 2.99970E+00 + 5 1.66050E-04 3.72464E+00 5.49987E+00 + 6 1.74987E-04 3.74501E+00 1.52616E+00 + 7 1.81704E-04 3.76007E+00 4.98991E+00 + 8 1.95526E-04 3.79048E+00 4.66342E+00 + 9 2.00135E-04 3.80047E+00 5.99118E+00 + 10 2.08346E-04 3.81809E+00 1.78070E+00 + 11 2.18749E-04 3.84014E+00 3.71410E-01 + 12 2.23720E-04 3.85058E+00 3.23852E+00 + 13 2.37507E-04 3.87925E+00 4.28906E+00 + 14 2.43692E-04 3.89199E+00 2.04177E+00 + 15 2.55125E-04 3.91539E+00 3.07633E+00 + 16 2.59595E-04 3.92448E+00 3.16463E+00 + 17 2.65016E-04 3.93549E+00 1.58621E+00 + 18 2.80839E-04 3.96744E+00 1.32207E+00 + 19 2.87923E-04 3.98169E+00 3.18059E+00 + 20 2.94848E-04 3.99560E+00 6.03835E+00 + 21 3.02298E-04 4.01054E+00 5.67005E+00 + 22 3.09363E-04 4.02470E+00 1.06213E+00 + 23 3.26328E-04 4.05871E+00 3.90012E+00 + 24 3.33149E-04 4.07240E+00 2.02779E+00 + 25 3.38049E-04 4.08225E+00 2.99598E+00 + 26 3.46368E-04 4.09900E+00 4.99207E+00 + 27 3.59100E-04 4.12473E+00 2.55239E+00 + 28 3.63414E-04 4.13348E+00 2.34858E+00 + 29 3.79503E-04 4.16627E+00 1.03210E+00 + 30 3.86437E-04 4.18050E+00 1.59647E+00 + 31 3.97739E-04 4.20382E+00 6.23912E-01 + 32 3.97942E-04 4.20424E+00 1.87259E+00 + 33 4.14637E-04 4.23903E+00 4.80595E+00 + 34 4.16342E-04 4.24261E+00 6.09807E+00 + 35 4.29658E-04 4.27071E+00 1.00531E+00 + 36 4.37679E-04 4.28779E+00 4.98638E+00 + 37 4.41787E-04 4.29658E+00 3.40877E+00 + 38 4.55390E-04 4.32590E+00 3.64132E+00 + 39 4.64717E-04 4.34619E+00 5.83264E+00 + 40 4.68329E-04 4.35409E+00 1.72709E+00 + 41 4.76198E-04 4.37137E+00 2.49289E+00 + 42 4.84223E-04 4.38910E+00 1.13665E+00 + 43 4.96600E-04 4.41661E+00 1.78049E+00 + 44 5.03711E-04 4.43249E+00 1.83825E+00 + 45 5.11712E-04 4.45042E+00 2.85382E+00 + 46 5.15820E-04 4.45964E+00 1.69322E+00 + 47 5.29047E-04 4.48937E+00 1.20835E+00 + 48 5.36671E-04 4.50649E+00 3.23773E-01 + 49 5.43953E-04 4.52281E+00 2.35339E+00 + 50 5.49443E-04 4.53509E+00 5.12833E+00 + 51 5.54423E-04 4.54618E+00 5.12701E+00 + 52 5.65917E-04 4.57163E+00 5.48799E+00 + 53 5.74632E-04 4.59072E+00 1.76010E+00 + 54 5.84446E-04 4.61195E+00 4.49547E-01 + 55 5.86646E-04 4.61667E+00 2.62663E+00 + 56 6.01367E-04 4.64776E+00 5.16518E+00 + 57 6.08710E-04 4.66293E+00 5.47161E+00 + 58 6.17931E-04 4.68162E+00 1.27573E+00 + 59 6.27659E-04 4.70090E+00 5.12506E+00 + 60 6.35294E-04 4.71569E+00 4.27014E+00 + 61 6.47928E-04 4.73949E+00 1.37917E+00 + 62 6.56572E-04 4.75528E+00 2.36731E+00 + 63 6.64348E-04 4.76914E+00 1.30147E+00 + 64 6.73421E-04 4.78492E+00 6.05718E+00 + 65 6.90207E-04 4.81297E+00 2.14538E+00 + 66 6.91218E-04 4.81462E+00 5.42836E+00 + 67 7.13072E-04 4.84895E+00 5.38340E+00 + 68 7.23783E-04 4.86497E+00 1.65867E+00 + 69 7.26526E-04 4.86899E+00 4.79254E-01 + 70 7.46204E-04 4.89693E+00 4.75645E+00 + 71 7.54225E-04 4.90789E+00 6.19616E+00 + 72 7.73002E-04 4.93263E+00 4.39635E+00 + 73 7.90806E-04 4.95502E+00 2.63669E+00 + 74 7.99610E-04 4.96575E+00 6.18202E+00 + 75 8.18402E-04 4.98793E+00 5.70800E+00 + 76 8.24144E-04 4.99454E+00 6.21054E+00 + 77 8.40742E-04 5.01319E+00 2.70761E+00 + 78 8.59562E-04 5.03365E+00 2.92114E+00 + 79 8.77957E-04 5.05301E+00 5.32006E+00 + 80 9.05658E-04 5.08116E+00 5.93883E+00 + 81 9.20598E-04 5.09592E+00 3.12112E+00 + 82 9.42942E-04 5.11753E+00 4.55992E+00 + 83 9.45289E-04 5.11978E+00 6.05193E+00 + 84 9.77019E-04 5.14967E+00 5.44215E+00 + 85 9.96447E-04 5.16768E+00 4.13132E+00 + 86 1.00730E-03 5.17767E+00 3.95243E+00 + 87 1.03883E-03 5.20659E+00 6.21214E+00 + 88 1.04609E-03 5.21323E+00 1.04527E+00 + 89 1.06958E-03 5.23486E+00 5.12648E+00 + 90 1.08864E-03 5.25257E+00 6.26396E+00 + 91 1.11460E-03 5.27719E+00 2.18579E-01 + 92 1.13215E-03 5.29431E+00 4.46693E+00 + 93 1.14369E-03 5.30584E+00 4.85334E+00 + 94 1.16604E-03 5.32903E+00 1.88253E+00 + 95 1.18347E-03 5.34817E+00 5.03066E+00 + 96 1.19860E-03 5.36584E+00 3.20959E+00 + 97 1.21492E-03 5.38645E+00 4.06172E+00 + 98 1.22877E-03 5.40576E+00 2.16819E+00 + 99 1.23151E-03 5.40984E+00 4.81483E+00 + 100 1.24600E-03 5.43360E+00 1.29115E+00 + 101 1.25672E-03 5.45477E+00 2.04414E+00 + 102 1.26567E-03 5.47720E+00 3.65557E-02 + 103 1.26785E-03 5.48393E+00 1.50397E+00 + 104 1.27346E-03 5.50662E+00 4.89712E+00 + 105 1.27472E-03 5.51425E+00 3.95288E+00 + 106 1.27660E-03 5.54554E+00 2.72301E+00 + 107 1.27603E-03 5.55952E+00 4.55974E+00 + 108 1.27474E-03 5.57384E+00 5.29081E-01 + 109 1.26998E-03 5.60288E+00 5.53959E-01 + 110 1.26671E-03 5.61675E+00 5.86364E+00 + 111 1.26544E-03 5.62152E+00 4.92497E+00 + 112 1.25507E-03 5.65347E+00 3.58770E+00 + 113 1.24863E-03 5.66964E+00 4.33306E+00 + 114 1.24398E-03 5.68029E+00 4.51927E+00 + 115 1.23756E-03 5.69388E+00 3.34174E+00 + 116 1.22326E-03 5.72108E+00 1.61438E+00 + 117 1.21715E-03 5.73174E+00 5.79840E+00 + 118 1.20540E-03 5.75107E+00 4.59838E-01 + 119 1.18847E-03 5.77687E+00 5.02394E+00 + 120 1.17653E-03 5.79400E+00 4.87449E+00 + 121 1.17121E-03 5.80142E+00 4.28657E+00 + 122 1.14983E-03 5.83008E+00 6.99815E-01 + 123 1.13687E-03 5.84678E+00 3.67408E+00 + 124 1.12791E-03 5.85810E+00 3.30864E+00 + 125 1.11644E-03 5.87236E+00 4.96324E-01 + 126 1.09665E-03 5.89648E+00 5.40507E+00 + 127 1.08066E-03 5.91566E+00 2.55588E+00 + 128 1.07030E-03 5.92798E+00 4.39877E+00 + 129 1.05630E-03 5.94454E+00 5.65357E-01 + 130 1.03192E-03 5.97330E+00 1.12831E+00 + 131 1.01283E-03 5.99586E+00 1.54011E+00 + 132 9.99881E-04 6.01124E+00 3.99102E+00 + 133 9.91167E-04 6.02165E+00 6.21419E-01 + 134 9.78265E-04 6.03716E+00 5.48243E+00 + 135 9.61765E-04 6.05722E+00 4.57185E+00 + 136 9.42502E-04 6.08105E+00 3.40015E+00 + 137 9.25574E-04 6.10243E+00 2.98812E-01 + 138 9.13610E-04 6.11783E+00 2.11175E+00 + 139 9.09079E-04 6.12373E+00 1.13679E+00 + 140 8.92056E-04 6.14630E+00 1.28636E+00 + 141 8.72203E-04 6.17349E+00 2.90021E+00 + 142 8.66445E-04 6.18158E+00 5.58889E+00 + 143 8.47255E-04 6.20926E+00 2.35970E+00 + 144 8.32869E-04 6.23084E+00 2.32936E+00 + 145 8.28000E-04 6.23832E+00 3.01028E+00 + 146 8.13853E-04 6.26063E+00 4.55668E+00 + 147 7.99884E-04 6.28357E+00 6.14150E+00 + 148 7.91325E-04 6.29813E+00 5.46129E+00 + 149 7.85130E-04 6.30892E+00 2.53614E+00 + 150 7.76806E-04 6.32379E+00 8.17081E-01 + 151 7.65385E-04 6.34493E+00 1.51088E+00 + 152 7.58605E-04 6.35792E+00 4.41286E+00 + 153 7.41877E-04 6.39157E+00 2.38365E+00 + 154 7.37056E-04 6.40173E+00 3.89477E+00 + 155 7.25802E-04 6.42634E+00 2.43421E+00 + 156 7.21828E-04 6.43536E+00 7.50626E-01 + 157 7.13370E-04 6.45516E+00 4.69615E+00 + 158 7.05190E-04 6.47518E+00 5.57271E+00 + 159 6.99637E-04 6.48929E+00 3.60242E+00 + 160 6.93252E-04 6.50610E+00 3.79708E+00 + 161 6.85671E-04 6.52690E+00 1.09453E+00 + 162 6.81105E-04 6.53992E+00 1.24687E+00 + 163 6.71341E-04 6.56910E+00 5.25692E+00 + 164 6.68810E-04 6.57698E+00 6.93876E-01 + 165 6.63313E-04 6.59459E+00 8.09162E-01 + 166 6.59044E-04 6.60876E+00 4.08483E+00 + 167 6.52440E-04 6.63159E+00 1.14414E+00 + 168 6.47830E-04 6.64821E+00 3.13254E+00 + 169 6.41532E-04 6.67193E+00 2.02355E+00 + 170 6.35660E-04 6.69515E+00 4.78725E+00 + 171 6.34517E-04 6.69981E+00 3.73709E+00 + 172 6.29546E-04 6.72056E+00 5.42421E+00 + 173 6.23562E-04 6.74673E+00 5.74294E+00 + 174 6.19781E-04 6.76398E+00 2.51438E+00 + 175 6.16330E-04 6.78022E+00 4.04884E+00 + 176 6.14311E-04 6.78994E+00 4.58756E+00 + 177 6.09859E-04 6.81199E+00 3.73295E+00 + 178 6.06431E-04 6.82956E+00 6.55529E-02 + 179 6.01474E-04 6.85589E+00 5.60808E+00 + 180 5.98768E-04 6.87073E+00 1.32815E+00 + 181 5.96413E-04 6.88392E+00 4.11083E+00 + 182 5.91444E-04 6.91259E+00 1.05230E+00 + 183 5.90295E-04 6.91939E+00 4.61349E+00 + 184 5.85489E-04 6.94847E+00 4.01695E+00 + 185 5.83517E-04 6.96071E+00 1.37657E-01 + 186 5.81272E-04 6.97486E+00 2.71196E+00 + 187 5.78289E-04 6.99401E+00 5.89933E+00 + 188 5.76575E-04 7.00520E+00 5.21978E+00 + 189 5.73350E-04 7.02658E+00 5.69745E+00 + 190 5.71560E-04 7.03864E+00 6.03158E+00 + 191 5.68743E-04 7.05789E+00 5.49465E+00 + 192 5.64132E-04 7.09006E+00 1.72021E+00 + 193 5.61588E-04 7.10815E+00 5.27877E+00 + 194 5.59756E-04 7.12132E+00 2.36098E+00 + 195 5.56811E-04 7.14274E+00 1.88392E+00 + 196 5.55120E-04 7.15516E+00 2.59343E+00 + 197 5.52199E-04 7.17683E+00 6.32432E-01 + 198 5.50608E-04 7.18874E+00 5.24285E+00 + 199 5.47185E-04 7.21461E+00 4.32397E+00 + 200 5.46538E-04 7.21953E+00 3.01972E+00 + 201 5.44392E-04 7.23595E+00 5.79657E+00 + 202 5.41576E-04 7.25766E+00 4.23604E+00 + 203 5.38058E-04 7.28504E+00 2.25479E+00 + 204 5.35635E-04 7.30406E+00 3.89270E+00 + 205 5.34983E-04 7.30920E+00 3.90415E+00 + 206 5.31620E-04 7.33583E+00 2.16076E+00 + 207 5.29704E-04 7.35112E+00 4.84800E+00 + 208 5.28338E-04 7.36205E+00 2.29181E+00 + 209 5.25008E-04 7.38882E+00 1.35769E-01 + 210 5.22110E-04 7.41228E+00 2.91500E+00 + 211 5.20403E-04 7.42617E+00 3.11538E+00 + 212 5.18845E-04 7.43888E+00 3.50989E+00 + 213 5.15574E-04 7.46569E+00 1.93503E+00 + 214 5.13630E-04 7.48170E+00 1.77377E+00 + 215 5.12412E-04 7.49175E+00 1.42715E+00 + 216 5.08751E-04 7.52211E+00 5.44456E+00 + 217 5.07234E-04 7.53474E+00 3.62659E+00 + 218 5.05241E-04 7.55140E+00 3.21075E+00 + 219 5.04291E-04 7.55935E+00 4.22473E+00 + 220 5.01578E-04 7.58214E+00 5.41403E+00 + 221 4.98332E-04 7.60954E+00 4.46105E+00 + 222 4.95914E-04 7.63005E+00 2.81659E-01 + 223 4.94881E-04 7.63883E+00 8.86734E-01 + 224 4.93310E-04 7.65222E+00 4.66888E+00 + 225 4.90286E-04 7.67809E+00 4.22378E+00 + 226 4.88239E-04 7.69567E+00 4.16875E-01 + 227 4.86568E-04 7.71007E+00 4.27207E+00 + 228 4.85257E-04 7.72140E+00 4.48032E+00 + 229 4.82792E-04 7.74276E+00 5.15659E+00 + 230 4.79841E-04 7.76845E+00 6.04974E+00 + 231 4.78499E-04 7.78017E+00 1.07457E+00 + 232 4.76574E-04 7.79705E+00 2.88394E+00 + 233 4.73724E-04 7.82212E+00 3.98392E+00 + 234 4.72591E-04 7.83212E+00 2.92302E+00 + 235 4.70785E-04 7.84811E+00 5.63140E+00 + 236 4.67247E-04 7.87958E+00 1.66606E+00 + 237 4.67012E-04 7.88169E+00 1.48535E+00 + 238 4.63353E-04 7.91447E+00 1.05429E+00 + 239 4.61777E-04 7.92865E+00 6.91657E-01 + 240 4.59935E-04 7.94529E+00 1.82552E+00 + 241 4.58013E-04 7.96272E+00 1.49996E+00 + 242 4.56747E-04 7.97423E+00 5.75781E+00 + 243 4.54700E-04 7.99291E+00 3.17904E+00 + 244 4.53004E-04 8.00844E+00 1.26667E+00 + 245 4.50464E-04 8.03179E+00 2.75104E+00 + 246 4.47849E-04 8.05596E+00 2.18412E+00 + 247 4.46235E-04 8.07095E+00 2.77529E+00 + 248 4.44537E-04 8.08676E+00 1.75382E+00 + 249 4.43182E-04 8.09943E+00 5.40320E+00 + 250 4.41273E-04 8.11734E+00 6.10492E+00 + 251 4.39615E-04 8.13294E+00 3.75326E+00 + 252 4.36368E-04 8.16368E+00 5.89759E+00 + 253 4.35098E-04 8.17575E+00 1.95123E+00 + 254 4.32308E-04 8.20241E+00 5.48466E+00 + 255 4.31727E-04 8.20798E+00 6.26640E+00 + 256 4.28606E-04 8.23805E+00 1.78666E+00 + 257 4.28337E-04 8.24065E+00 2.58798E+00 + 258 4.25431E-04 8.26885E+00 3.28910E-01 + 259 4.23836E-04 8.28442E+00 4.15963E+00 + 260 4.21505E-04 8.30726E+00 5.12112E+00 + 261 4.19485E-04 8.32716E+00 1.37436E+00 + 262 4.18393E-04 8.33796E+00 4.03963E+00 + 263 4.16023E-04 8.36149E+00 4.55622E+00 + 264 4.14288E-04 8.37881E+00 4.33559E+00 + 265 4.12225E-04 8.39950E+00 8.52791E-01 + 266 4.11650E-04 8.40529E+00 6.01462E+00 + 267 4.09133E-04 8.43071E+00 2.10893E+00 + 268 4.08160E-04 8.44058E+00 4.38058E+00 + 269 4.05551E-04 8.46717E+00 2.32127E+00 + 270 4.04163E-04 8.48139E+00 5.60517E+00 + 271 4.02418E-04 8.49935E+00 5.44053E+00 + 272 4.00854E-04 8.51550E+00 1.13805E-01 + 273 3.98516E-04 8.53978E+00 1.99889E+00 + 274 3.96317E-04 8.56276E+00 3.52958E+00 + 275 3.96281E-04 8.56313E+00 4.46275E+00 + 276 3.93652E-04 8.59079E+00 1.68615E+00 + 277 3.92150E-04 8.60668E+00 3.99206E+00 + 278 3.90524E-04 8.62395E+00 2.87261E+00 + 279 3.88829E-04 8.64204E+00 2.10570E+00 + 280 3.86593E-04 8.66603E+00 4.03455E+00 + 281 3.85688E-04 8.67579E+00 5.83360E+00 + 282 3.84081E-04 8.69317E+00 1.63197E+00 + 283 3.81520E-04 8.72105E+00 1.25746E+00 + 284 3.80458E-04 8.73267E+00 2.29975E+00 + 285 3.79377E-04 8.74453E+00 7.36735E-01 + 286 3.77355E-04 8.76681E+00 9.58005E-01 + 287 3.75911E-04 8.78282E+00 6.01273E+00 + 288 3.74100E-04 8.80298E+00 5.02066E+00 + 289 3.71648E-04 8.83046E+00 1.85346E+00 + 290 3.71250E-04 8.83494E+00 5.83648E+00 + 291 3.69882E-04 8.85038E+00 3.00296E-01 + 292 3.68076E-04 8.87086E+00 3.30809E+00 + 293 3.65246E-04 8.90318E+00 3.71022E+00 + 294 3.64572E-04 8.91091E+00 3.82945E+00 + 295 3.63031E-04 8.92868E+00 4.00397E+00 + 296 3.61039E-04 8.95176E+00 2.95279E+00 + 297 3.60015E-04 8.96368E+00 5.47311E+00 + 298 3.57888E-04 8.98856E+00 1.12168E+00 + 299 3.57040E-04 8.99854E+00 4.80219E+00 + 300 3.55281E-04 9.01930E+00 1.55111E+00 + 301 3.53801E-04 9.03687E+00 6.08995E-01 + 302 3.51967E-04 9.05876E+00 2.36233E+00 + 303 3.50384E-04 9.07776E+00 6.15409E+00 + 304 3.48715E-04 9.09790E+00 4.88265E+00 + 305 3.47114E-04 9.11733E+00 2.11896E+00 + 306 3.46426E-04 9.12571E+00 1.04092E+00 + 307 3.44736E-04 9.14637E+00 4.54608E+00 + 308 3.42960E-04 9.16821E+00 1.88818E+00 + 309 3.41576E-04 9.18533E+00 4.82717E+00 + 310 3.39777E-04 9.20771E+00 6.06954E+00 + 311 3.39231E-04 9.21452E+00 5.08807E+00 + 312 3.37077E-04 9.24155E+00 1.45455E+00 + 313 3.35719E-04 9.25870E+00 2.20353E+00 + 314 3.34535E-04 9.27371E+00 1.87286E-01 + 315 3.33830E-04 9.28268E+00 3.42058E+00 + 316 3.32302E-04 9.30219E+00 3.68762E+00 + 317 3.29945E-04 9.33252E+00 1.31897E+00 + 318 3.28589E-04 9.35009E+00 5.34963E+00 + 319 3.28221E-04 9.35487E+00 5.37587E+00 + 320 3.26228E-04 9.38087E+00 3.31494E+00 + 321 3.25661E-04 9.38831E+00 4.29571E+00 + 322 3.23900E-04 9.41149E+00 3.94410E+00 + 323 3.21647E-04 9.44137E+00 1.24787E+00 + 324 3.20794E-04 9.45276E+00 7.49383E-01 + 325 3.19050E-04 9.47613E+00 2.87842E+00 + 326 3.17809E-04 9.49288E+00 1.00801E+00 + 327 3.17340E-04 9.49922E+00 2.39949E+00 + 328 3.15768E-04 9.52057E+00 1.56978E-01 + 329 3.14736E-04 9.53466E+00 4.22818E+00 + 330 3.13042E-04 9.55790E+00 1.11766E+00 + 331 3.11708E-04 9.57632E+00 5.43057E+00 + 332 3.10004E-04 9.59997E+00 9.91319E-01 + 333 3.08801E-04 9.61678E+00 4.27092E+00 + 334 3.07857E-04 9.63002E+00 3.85311E+00 + 335 3.06138E-04 9.65426E+00 2.64513E-01 + 336 3.05352E-04 9.66540E+00 2.74719E-01 + 337 3.04082E-04 9.68347E+00 5.90299E+00 + 338 3.02958E-04 9.69955E+00 6.13863E+00 + 339 3.00943E-04 9.72854E+00 2.37575E+00 + 340 3.00679E-04 9.73236E+00 5.31687E-01 + 341 2.99598E-04 9.74803E+00 2.87529E+00 + 342 2.98223E-04 9.76807E+00 4.10590E+00 + 343 2.96327E-04 9.79590E+00 2.79948E+00 + 344 2.95492E-04 9.80821E+00 5.38882E+00 + 345 2.94055E-04 9.82953E+00 5.59734E+00 + 346 2.93368E-04 9.83975E+00 3.01684E+00 + 347 2.92282E-04 9.85598E+00 5.60600E+00 + 348 2.90108E-04 9.88873E+00 5.94655E+00 + 349 2.89002E-04 9.90548E+00 5.56204E+00 + 350 2.88663E-04 9.91064E+00 4.05350E+00 + 351 2.86914E-04 9.93736E+00 2.20510E+00 + 352 2.85953E-04 9.95214E+00 6.10765E+00 + 353 2.84254E-04 9.97840E+00 3.92342E+00 + 354 2.83655E-04 9.98770E+00 4.72616E+00 + 355 2.82596E-04 1.00042E+01 1.40958E+00 + 356 2.81087E-04 1.00279E+01 3.80983E+00 + 357 2.80310E-04 1.00401E+01 4.63850E-01 + 358 2.78930E-04 1.00620E+01 2.45007E+00 + 359 2.78384E-04 1.00707E+01 1.61936E+00 + 360 2.77164E-04 1.00902E+01 6.10311E+00 + 361 2.75841E-04 1.01114E+01 3.55598E+00 + 362 2.74159E-04 1.01386E+01 4.55032E+00 + 363 2.73916E-04 1.01426E+01 5.18451E+00 + 364 2.72724E-04 1.01620E+01 3.90454E+00 + 365 2.71092E-04 1.01888E+01 3.41364E+00 + 366 2.70117E-04 1.02049E+01 1.27777E+00 + 367 2.68719E-04 1.02281E+01 8.53254E-02 + 368 2.68179E-04 1.02371E+01 8.94078E-01 + 369 2.67204E-04 1.02534E+01 1.89443E+00 + 370 2.66091E-04 1.02722E+01 5.86907E-01 + 371 2.65321E-04 1.02852E+01 2.11165E-01 + 372 2.63548E-04 1.03153E+01 3.73390E+00 + 373 2.62946E-04 1.03256E+01 6.27844E+00 + 374 2.61960E-04 1.03426E+01 6.15931E+00 + 375 2.60273E-04 1.03717E+01 5.42831E+00 + 376 2.59650E-04 1.03825E+01 2.68390E+00 + 377 2.58200E-04 1.04079E+01 5.79947E+00 + 378 2.57577E-04 1.04188E+01 3.45486E+00 + 379 2.56830E-04 1.04320E+01 4.94309E+00 + 380 2.55550E-04 1.04547E+01 5.30896E+00 + 381 2.54322E-04 1.04765E+01 5.12806E+00 + 382 2.53293E-04 1.04950E+01 9.16462E-01 + 383 2.52099E-04 1.05165E+01 3.42563E+00 + 384 2.51098E-04 1.05346E+01 1.33480E+00 + 385 2.50234E-04 1.05504E+01 3.57945E-01 + 386 2.49156E-04 1.05701E+01 2.86239E+00 + 387 2.48550E-04 1.05812E+01 4.10999E-01 + 388 2.48048E-04 1.05905E+01 4.47406E+00 + 389 2.46821E-04 1.06132E+01 6.10749E+00 + 390 2.45335E-04 1.06409E+01 5.63394E+00 + 391 2.44893E-04 1.06492E+01 1.26190E+00 + 392 2.43351E-04 1.06783E+01 3.74044E-01 + 393 2.42604E-04 1.06925E+01 2.61487E-01 + 394 2.41763E-04 1.07085E+01 6.86245E-01 + 395 2.40684E-04 1.07291E+01 1.43539E-01 + 396 2.39861E-04 1.07450E+01 6.15214E+00 + 397 2.38732E-04 1.07668E+01 1.27240E+00 + 398 2.38028E-04 1.07805E+01 6.05108E+00 + 399 2.36768E-04 1.08051E+01 2.45062E+00 + 400 2.36366E-04 1.08130E+01 4.42757E+00 + 401 2.35798E-04 1.08242E+01 4.09007E+00 + 402 2.34463E-04 1.08506E+01 5.09218E-01 + 403 2.33672E-04 1.08663E+01 1.27441E+00 + 404 2.32456E-04 1.08907E+01 4.10183E+00 + 405 2.31469E-04 1.09106E+01 1.34251E+00 + 406 2.31083E-04 1.09184E+01 3.09899E+00 + 407 2.30332E-04 1.09336E+01 8.19130E-01 + 408 2.29506E-04 1.09504E+01 2.31715E+00 + 409 2.28382E-04 1.09735E+01 5.08581E+00 + 410 2.27280E-04 1.09962E+01 2.32448E+00 + 411 2.26602E-04 1.10102E+01 5.05798E+00 + 412 2.25897E-04 1.10249E+01 6.62343E-01 + 413 2.25004E-04 1.10436E+01 3.17963E+00 + 414 2.24334E-04 1.10576E+01 1.01444E-01 + 415 2.22991E-04 1.10860E+01 1.28515E+00 + 416 2.22681E-04 1.10926E+01 6.25880E+00 + 417 2.21642E-04 1.11147E+01 4.22276E-01 + 418 2.20778E-04 1.11333E+01 4.13750E+00 + 419 2.20127E-04 1.11473E+01 9.23098E-02 + 420 2.18684E-04 1.11785E+01 1.19735E+00 + 421 2.18434E-04 1.11840E+01 2.96201E+00 + 422 2.17695E-04 1.12001E+01 3.88642E+00 + 423 2.16526E-04 1.12258E+01 6.43547E-02 + 424 2.16069E-04 1.12359E+01 6.23804E+00 + 425 2.14748E-04 1.12652E+01 4.58868E+00 + 426 2.14198E-04 1.12775E+01 4.11392E+00 + 427 2.13000E-04 1.13044E+01 2.45374E+00 + 428 2.12141E-04 1.13238E+01 6.25987E+00 + 429 2.11980E-04 1.13274E+01 1.57282E+00 + 430 2.10626E-04 1.13583E+01 6.27770E+00 + 431 2.10350E-04 1.13646E+01 5.84165E+00 + 432 2.09131E-04 1.13926E+01 4.56388E+00 + 433 2.08568E-04 1.14056E+01 1.11609E+00 + 434 2.07636E-04 1.14273E+01 1.62757E+00 + 435 2.06959E-04 1.14431E+01 6.07793E+00 + 436 2.06161E-04 1.14618E+01 1.32832E+00 + 437 2.05124E-04 1.14863E+01 2.37320E+00 + 438 2.04663E-04 1.14972E+01 5.55134E+00 + 439 2.04020E-04 1.15125E+01 2.96299E+00 + 440 2.03154E-04 1.15333E+01 3.39101E+00 + 441 2.02601E-04 1.15465E+01 1.91734E+00 + 442 2.01750E-04 1.15671E+01 1.52150E+00 + 443 2.01192E-04 1.15806E+01 7.75597E-01 + 444 2.00077E-04 1.16078E+01 3.53651E+00 + 445 1.99283E-04 1.16273E+01 1.51282E+00 + 446 1.98672E-04 1.16423E+01 5.94301E+00 + 447 1.98344E-04 1.16505E+01 1.74139E+00 + 448 1.97326E-04 1.16758E+01 2.71647E+00 + 449 1.96739E-04 1.16904E+01 2.52776E+00 + 450 1.96255E-04 1.17026E+01 1.28167E+00 + 451 1.95071E-04 1.17324E+01 2.16073E+00 + 452 1.94721E-04 1.17413E+01 1.23441E+00 + 453 1.93953E-04 1.17608E+01 4.24845E+00 + 454 1.93142E-04 1.17816E+01 3.40290E+00 + 455 1.92550E-04 1.17968E+01 2.89214E+00 + 456 1.91370E-04 1.18273E+01 3.44644E+00 + 457 1.91108E-04 1.18341E+01 3.68406E+00 + 458 1.90404E-04 1.18525E+01 3.62257E+00 + 459 1.89403E-04 1.18788E+01 5.05069E+00 + 460 1.88768E-04 1.18955E+01 2.07005E+00 + 461 1.88053E-04 1.19145E+01 3.65525E+00 + 462 1.87575E-04 1.19272E+01 3.25702E+00 + 463 1.86663E-04 1.19516E+01 2.27251E-01 + 464 1.86322E-04 1.19608E+01 3.86857E+00 + 465 1.85647E-04 1.19790E+01 2.59955E-01 + 466 1.84980E-04 1.19971E+01 3.95417E+00 + 467 1.84261E-04 1.20167E+01 5.14513E+00 + 468 1.83738E-04 1.20310E+01 1.13382E+00 + 469 1.82700E-04 1.20596E+01 5.10449E-01 + 470 1.82531E-04 1.20642E+01 5.70460E+00 + 471 1.81466E-04 1.20938E+01 6.01886E+00 + 472 1.80751E-04 1.21138E+01 4.14713E+00 + 473 1.80600E-04 1.21181E+01 3.76787E+00 + 474 1.79811E-04 1.21402E+01 3.25689E+00 + 475 1.79105E-04 1.21602E+01 5.25729E+00 + 476 1.78550E-04 1.21760E+01 3.91694E-01 + 477 1.78100E-04 1.21889E+01 4.50881E+00 + 478 1.77426E-04 1.22082E+01 2.38697E+00 + 479 1.76774E-04 1.22270E+01 2.10155E-01 + 480 1.75725E-04 1.22574E+01 6.59004E-01 + 481 1.75504E-04 1.22638E+01 2.84896E+00 + 482 1.74760E-04 1.22856E+01 4.99661E+00 diff --git a/WaveFunc.C b/WaveFunc.C index 43d07c6..6c9467b 100644 --- a/WaveFunc.C +++ b/WaveFunc.C @@ -34,9 +34,8 @@ WaveSpectrum::WaveSpectrum (const char* file, double angle, double g, double t0) { if (temp[0] == '#') continue; std::stringstream str(temp); - int i; Component wc; - str >> i >> wc.A >> wc.omega >> wc.eps; + str >> wc.id >> wc.Ampl >> wc.omega >> wc.eps; wave.emplace_back(wc); } @@ -82,7 +81,7 @@ double WaveSpectrum::evaluate (const Vec3& X) const double res = 0.0; for (const Component& w : wave) - res += w.A*sin(w.omega*t - w.omega*XoG + w.eps); + res += w.Ampl*sin(w.omega*t - w.omega*XoG + w.eps); if (time > 0.0 && t < time) res *= t/time; @@ -118,6 +117,12 @@ void HydroStaticPressure::setParam (const std::string& name, double value) } +/*! + The pressure is assumed in the direction of the outward-directed + normal vector of the shell surface at the evaluation point. + Therefore, it will always have a non-positive value. +*/ + double HydroStaticPressure::evaluate (const Vec3& X) const { #ifdef USE_OPENMP @@ -126,5 +131,5 @@ double HydroStaticPressure::evaluate (const Vec3& X) const const size_t i = 0; #endif const double depth = X.z + dZ[i] - (w_line + zeta(X)); - return depth < 0.0 ? -rhow*grav*depth : 0.0; + return depth < 0.0 ? rhow*grav*depth : 0.0; } diff --git a/WaveFunc.h b/WaveFunc.h index ca7fe97..e5675f8 100644 --- a/WaveFunc.h +++ b/WaveFunc.h @@ -26,7 +26,8 @@ class WaveSpectrum : public RealFunc //! \brief A struct describing a single wave component of a spectrum. struct Component { - double A = 0.0; //!< Amplitude + int id = 0; //!< Component id + double Ampl = 0.0; //!< Amplitude double omega = 0.0; //!< Angular frequency double eps = 0.0; //!< phase shift }; @@ -54,7 +55,7 @@ class WaveSpectrum : public RealFunc void setParam(const std::string& name, double value) override; protected: - //! \brief Evaluates the function. + //! \brief Evaluates the sea elevation function. double evaluate(const Vec3& X) const override; }; @@ -89,7 +90,7 @@ class HydroStaticPressure : public RealFunc void setParam(const std::string& name, double value) override; protected: - //! \brief Evaluates the function. + //! \brief Evaluates the hydrostatic pressure function. double evaluate(const Vec3& X) const override; }; diff --git a/main.C b/main.C index 662efc3..ad074a2 100644 --- a/main.C +++ b/main.C @@ -493,6 +493,10 @@ int main (int argc, char** argv) if (nodalR != 'm' && !model->writeGlvLoc(locfiles,nodalR,geoBlk)) return terminate(12); + // Write surface normal vectors (but only if no surface pressure output) + if (!model->writeGlvNormal(geoBlk,nBlock)) + return terminate(12); + // Write surface pressures, if any if (!model->writeGlvT(iStep,geoBlk,nBlock)) return terminate(13);