← Back to team overview

yade-dev team mailing list archive

[Branch ~yade-pkg/yade/git-trunk] Rev 3903: Add watercolumn example for SPH.

 

------------------------------------------------------------
revno: 3903
committer: Anton Gladky <gladky.anton@xxxxxxxxx>
timestamp: Wed 2014-04-09 16:03:16 +0200
message:
  Add watercolumn example for SPH.
added:
  examples/sph/box.geo
  examples/sph/box.mesh
  examples/sph/cpt/
  examples/sph/cpt/a
  examples/sph/watercolumn.py


--
lp:yade
https://code.launchpad.net/~yade-pkg/yade/git-trunk

Your team Yade developers is subscribed to branch lp:yade.
To unsubscribe from this branch go to https://code.launchpad.net/~yade-pkg/yade/git-trunk/+edit-subscription
=== added file 'examples/sph/box.geo'
--- examples/sph/box.geo	1970-01-01 00:00:00 +0000
+++ examples/sph/box.geo	2014-04-09 14:03:16 +0000
@@ -0,0 +1,27 @@
+acc = 300.0;
+hW = 2000.0;
+thickW = 50.0;
+
+Point(1) = {0,      300,      0, acc};
+Point(2) = {0,      -hW,      0, acc};
+Point(3) = {1000.0, -hW,      0, acc};
+Point(4) = {2000.0, -hW,      0, acc};
+Point(5) = {2300,   -hW+300,  0, acc};
+Point(6) = {2300,   -hW,      0, acc};
+Point(7) = {3700,   -hW,      0, acc};
+Point(8) = {4000,   -hW+300,  0, acc};
+Point(9) = {4000,   300,      0, acc};
+Line(1) = {1, 2};
+Line(2) = {2, 3};
+Line(3) = {3, 4};
+Line(4) = {4, 5};
+Line(5) = {5, 6};
+Line(6) = {6, 7};
+Line(7) = {7, 8};
+Line(8) = {8, 9};
+Line(9) = {9, 1};
+Line Loop(10) = {1, 2, 3, 4, 5, 6, 7, 8, 9};
+Plane Surface(11) = {10};
+Extrude {0, 0, thickW} {
+  Surface{11};
+}

=== added file 'examples/sph/box.mesh'
--- examples/sph/box.mesh	1970-01-01 00:00:00 +0000
+++ examples/sph/box.mesh	2014-04-09 14:03:16 +0000
@@ -0,0 +1,1273 @@
+ MeshVersionFormatted 2
+ Dimension
+ 3
+ Vertices
+ 388
+                    0                       300                         0      1
+                    0                     -2000                         0      2
+                 1000                     -2000                         0      3
+                 2000                     -2000                         0      4
+                 2300                     -1700                         0      5
+                 2300                     -2000                         0      6
+                 3700                     -2000                         0      7
+                 4000                     -1700                         0      8
+                 4000                       300                         0      9
+                    0                       300                        50      10
+                    0                     -2000                        50      11
+                 1000                     -2000                        50      15
+                 2000                     -2000                        50      19
+                 2300                     -1700                        50      23
+                 2300                     -2000                        50      27
+                 3700                     -2000                        50      31
+                 4000                     -1700                        50      35
+                 4000                       300                        50      39
+                    0           12.500000000903                         0      1
+                    0          -274.99999999796                         0      1
+                    0          -562.49999999681                         0      1
+                    0          -849.99999999633                         0      1
+                    0          -1137.4999999966                         0      1
+                    0          -1424.9999999969                         0      1
+                    0          -1712.4999999971                         0      1
+      249.99999999967                     -2000                         0      2
+      499.99999999976                     -2000                         0      2
+      749.99999999991                     -2000                         0      2
+                 1250                     -2000                         0      3
+                 1500                     -2000                         0      3
+                 1750                     -2000                         0      3
+      2150.0000000009          -1849.9999999991                         0      4
+                 2580                     -2000                         0      6
+                 2860                     -2000                         0      6
+                 3140                     -2000                         0      6
+                 3420                     -2000                         0      6
+                 3850                     -1850                         0      7
+                 4000          -1414.2857142857                         0      8
+                 4000          -1128.5714285715                         0      8
+                 4000           -842.8571428572                         0      8
+                 4000           -557.1428571428                         0      8
+                 4000          -271.42857142853                         0      8
+                 4000           14.285714285733                         0      8
+      3714.2857142857                       300                         0      9
+      3428.5714285715                       300                         0      9
+      3142.8571428572                       300                         0      9
+      2857.1428571429                       300                         0      9
+      2571.4285714287                       300                         0      9
+      2285.7142857144                       300                         0      9
+      2000.0000000001                       300                         0      9
+      1714.2857142856                       300                         0      9
+      1428.5714285713                       300                         0      9
+      1142.8571428571                       300                         0      9
+       857.1428571428                       300                         0      9
+      571.42857142853                       300                         0      9
+      285.71428571427                       300                         0      9
+                    0           12.500000000903                        50      13
+                    0          -274.99999999796                        50      13
+                    0          -562.49999999681                        50      13
+                    0          -849.99999999633                        50      13
+                    0          -1137.4999999966                        50      13
+                    0          -1424.9999999969                        50      13
+                    0          -1712.4999999971                        50      13
+      249.99999999967                     -2000                        50      14
+      499.99999999976                     -2000                        50      14
+      749.99999999991                     -2000                        50      14
+                 1250                     -2000                        50      15
+                 1500                     -2000                        50      15
+                 1750                     -2000                        50      15
+      2150.0000000009          -1849.9999999991                        50      16
+                 2580                     -2000                        50      18
+                 2860                     -2000                        50      18
+                 3140                     -2000                        50      18
+                 3420                     -2000                        50      18
+                 3850                     -1850                        50      19
+                 4000          -1414.2857142857                        50      20
+                 4000          -1128.5714285715                        50      20
+                 4000           -842.8571428572                        50      20
+                 4000           -557.1428571428                        50      20
+                 4000          -271.42857142853                        50      20
+                 4000           14.285714285733                        50      20
+      3714.2857142857                       300                        50      21
+      3428.5714285715                       300                        50      21
+      3142.8571428572                       300                        50      21
+      2857.1428571429                       300                        50      21
+      2571.4285714287                       300                        50      21
+      2285.7142857144                       300                        50      21
+      2000.0000000001                       300                        50      21
+      1714.2857142856                       300                        50      21
+      1428.5714285713                       300                        50      21
+      1142.8571428571                       300                        50      21
+       857.1428571428                       300                        50      21
+      571.42857142853                       300                        50      21
+      285.71428571427                       300                        50      21
+      1529.4446589135          -846.19882489058                         0      11
+      2948.8647910372          -797.27139693469                         0      11
+      768.85477021509          -1191.4331778603                         0      11
+      2225.3179510661          -995.96925999185                         0      11
+      882.98506830098          -423.62886416129                         0      11
+      2357.4759469348          -385.72537641386                         0      11
+      3334.1697178436          -1295.4261086936                         0      11
+      1753.5944220141           -1401.251792038                         0      11
+      2785.9051377217          -1352.6096802644                         0      11
+      3340.5444050499          -333.78170779744                         0      11
+      1781.7192577219          -312.00179782218                         0      11
+      1239.0164902535          -1479.6066037279                         0      11
+      3406.8595304718          -826.55892132953                         0      11
+      498.05326042224          -798.11345014194                         0      11
+      2877.3940074457          -158.86269305865                         0      11
+      1312.9094429575          -188.08933892446                         0      11
+      421.26022679498          -1568.3121263798                         0      11
+       479.3644383792           -179.9473566853                         0      11
+       1079.615368893           -851.4864133385                         0      11
+      2094.5747160518          -1373.0348484798                         0      11
+      1914.3929728395          -699.98358083334                         0      11
+      885.37357867263          -1614.8958107276                         0      11
+      2512.5374604002          -723.15889087768                         0      11
+      3033.6727876963          -1618.0590145742                         0      11
+      3535.4401799374          -1637.6777328536                         0      11
+      2120.8436812794          -73.272527652081                         0      11
+      1870.1705689431          -1065.5436320433                         0      11
+      1470.6887651328          -1197.6555917815                         0      11
+      3101.6438544384          -1176.8301328472                         0      11
+      2640.5965260745          -1681.7980989853                         0      11
+      376.50369589741          -1138.3073016349                         0      11
+      934.26001492933          -130.97450694126                         0      11
+      2446.7543465715          -1314.9110827347                         0      11
+      1519.1695021153          -1704.2615789078                         0      11
+        1321.38510647          -490.93211495581                         0      11
+      1856.8766695228          -1713.8800482108                         0      11
+      3629.2731642966          -119.36920619177                         0      11
+       3579.263470053          -1141.5591856709                         0      11
+      2705.0374641419          -1026.2932625922                         0      11
+      3603.2172344318          -594.03118802079                         0      11
+       3110.049443106           -478.9767977651                         0      11
+      2571.4746621187          -46.609234895912                         0      11
+      1135.8000096194          -1171.5142581261                         0      11
+      3187.7055426079           -56.57788934156                         0      11
+      2713.4240465413          -521.10101130376                         0      11
+      1609.7287827223          -27.757030646951                         0      11
+      2218.2808542131          -611.73080288984                         0      11
+      339.77434389417          -501.56273219575                         0      11
+      1628.4095443371           -583.8378231765                         0      11
+      1127.7531780478          -1736.4024554212                         0      11
+      783.45663487986          -911.36971711706                         0      11
+      632.05380854947          -1754.5359385757                         0      11
+      3708.2839439231          -1405.4255247641                         0      11
+      2454.7868603542          -1539.9042746893                         0      11
+      2091.9665388526          -311.17123097013                         0      11
+      2082.1217200227          -1635.0000625566                         0      11
+      3273.4354982145          -1742.8740271813                         0      11
+      3706.0339596851          -848.60652082524                         0      11
+      722.36683469073           42.963518414442                         0      11
+      262.94130347572          -1733.5288300386                         0      11
+       3778.058225824           -1628.058225824                         0      11
+      622.80394092887          -467.13330345548                         0      11
+      1004.3260710173          -1346.6259737243                         0      11
+      3129.5375392764          -925.21562551409                         0      11
+      1920.0539321953           81.716682003892                         0      11
+      264.10502811154           33.186453018978                         0      11
+      3335.9949786522          -564.15567518578                         0      11
+      2971.1394624533          -1381.1740234635                         0      11
+      1351.8207521926          -1026.0372727017                         0      11
+       215.1945587744             -921.44343644                         0      11
+      3777.9039666689          -389.01253121797                         0      11
+      696.60118099605          -1445.3382867397                         0      11
+      728.35605944582          -231.74410256578                         0      11
+      3319.3115215333          -1073.9684694836                         0      11
+      1530.8304268418          -320.52362110518                         0      11
+      1285.7142857142           76.103882369648                         0      11
+      2474.4274814706          -1015.6649119318                         0      11
+      2951.0287882931           71.331937890859                         0      11
+       1520.501381672          -1425.3304293798                         0      11
+       253.6973135114          -1374.8468483467                         0      11
+       1110.952305133          -319.40091038483                         0      11
+      1350.6011578823          -1811.1729631562                         0      11
+      2338.1660764839           65.729833865481                         0      11
+      2858.5107349988          -1760.0160943389                         0      11
+      1627.2275441016          -1039.7316169914                         0      11
+      3665.7543823172          -1786.4216673679                         0      11
+      1923.8707999441          -1275.6821326923                         0      11
+      1085.6123252695          -587.51500034014                         0      11
+      3438.4243515027           68.682495248415                         0      11
+      605.58286398704          -1059.6338500857                         0      11
+       2609.508999673          -311.68836401073                         0      11
+      245.98314002489          -259.25014738251                         0      11
+      2058.4956784593          -885.34925903342                         0      11
+      1764.0325246399          -866.84558719798                         0      11
+       3110.001253574          -276.20095824932                         0      11
+      2265.1599402191          -1488.4548701852                         0      11
+      2468.3827255104           -1797.189068538                         0      11
+      1682.7777746686          -1801.9522820407                         0      11
+      874.99999999995          -1825.0191987092                         0      11
+      2826.2893615091          -1545.2191605369                         0      11
+      2072.1121891632          -1138.5574595677                         0      11
+      798.90642055128           -666.5740256143                         0      11
+      3363.9794589634          -1515.8092628208                         0      11
+      2310.0022045934          -817.77285230863                         0      11
+       956.7723965576          -1038.1425165412                         0      11
+      203.80770423961          -694.60096326135                         0      11
+      1332.7146260126           -731.0012415672                         0      11
+      2736.2963545966          -809.01301810765                         0      11
+       1912.917229461          -1489.1145205076                         0      11
+      2848.8056048724          -1187.6725098473                         0      11
+       1910.311811976          -490.06137012571                         0      11
+      2755.2480385675           99.259138194832                         0      11
+        3137.71838514          -698.10133794125                         0      11
+      3769.8773198013           73.021229696858                         0      11
+      2504.8730852043          -516.86975916373                         0      11
+      1707.6394965756          -1589.2984418475                         0      11
+      1967.7996779093          -1839.7760221533                         0      11
+      2322.1948729117          -158.21196510554                         0      11
+      520.41667523366          -1296.3119318412                         0      11
+      2648.8039626677          -1229.9564880484                         0      11
+       1694.342246968          -1234.1991991544                         0      11
+      2142.5260679531           138.52442710411                         0      11
+      3504.2273541441          -1399.1795629606                         0      11
+       1015.165916226           86.098665248372                         0      11
+      3835.4108901534          -138.50067297114                         0      11
+      476.63184692998           107.35456591656                         0      11
+      3046.5330873028          -1818.6891070206                         0      11
+      3816.8158456071          -1224.5822677841                         0      11
+      2911.6976013938          -987.04932430719                         0      11
+      1288.2093595642          -1287.1762444129                         0      11
+      2633.2953778388          -1455.0724689842                         0      11
+      2909.6711019268          -606.11308583545                         0      11
+      3503.5967245883          -1828.5311691038                         0      11
+        3819.95043295          -657.02581520156                         0      11
+      443.03408776235          -1816.8508568957                         0      11
+      1922.7072374645          -140.38078415628                         0      11
+      2263.0577707571          -1221.0987388152                         0      11
+      165.61538782086          -1558.4628703367                         0      11
+       3820.422655069          -1037.2353091418                         0      11
+      3441.4789343952          -149.50414065605                         0      11
+      1068.5780858665          -1557.7249410525                         0      11
+      3169.8062650757          -1459.8592624457                         0      11
+      124.99999999984          -1856.2499999986                         0      11
+      1144.2671355981           -77.68362293098                         0      11
+      3530.4249014986          -366.64277666758                         0      11
+      2888.3414753778          -392.15715170384                         0      11
+      1352.7439897812          -1626.7393704996                         0      11
+      1529.4446589135          -846.19882489058                        50      58
+      2948.8647910372          -797.27139693469                        50      58
+      768.85477021509          -1191.4331778603                        50      58
+      2225.3179510661          -995.96925999185                        50      58
+      882.98506830098          -423.62886416129                        50      58
+      2357.4759469348          -385.72537641386                        50      58
+      3334.1697178436          -1295.4261086936                        50      58
+      1753.5944220141           -1401.251792038                        50      58
+      2785.9051377217          -1352.6096802644                        50      58
+      3340.5444050499          -333.78170779744                        50      58
+      1781.7192577219          -312.00179782218                        50      58
+      1239.0164902535          -1479.6066037279                        50      58
+      3406.8595304718          -826.55892132953                        50      58
+      498.05326042224          -798.11345014194                        50      58
+      2877.3940074457          -158.86269305865                        50      58
+      1312.9094429575          -188.08933892446                        50      58
+      421.26022679498          -1568.3121263798                        50      58
+       479.3644383792           -179.9473566853                        50      58
+       1079.615368893           -851.4864133385                        50      58
+      2094.5747160518          -1373.0348484798                        50      58
+      1914.3929728395          -699.98358083334                        50      58
+      885.37357867263          -1614.8958107276                        50      58
+      2512.5374604002          -723.15889087768                        50      58
+      3033.6727876963          -1618.0590145742                        50      58
+      3535.4401799374          -1637.6777328536                        50      58
+      2120.8436812794          -73.272527652081                        50      58
+      1870.1705689431          -1065.5436320433                        50      58
+      1470.6887651328          -1197.6555917815                        50      58
+      3101.6438544384          -1176.8301328472                        50      58
+      2640.5965260745          -1681.7980989853                        50      58
+      376.50369589741          -1138.3073016349                        50      58
+      934.26001492933          -130.97450694126                        50      58
+      2446.7543465715          -1314.9110827347                        50      58
+      1519.1695021153          -1704.2615789078                        50      58
+        1321.38510647          -490.93211495581                        50      58
+      1856.8766695228          -1713.8800482108                        50      58
+      3629.2731642966          -119.36920619177                        50      58
+       3579.263470053          -1141.5591856709                        50      58
+      2705.0374641419          -1026.2932625922                        50      58
+      3603.2172344318          -594.03118802079                        50      58
+       3110.049443106           -478.9767977651                        50      58
+      2571.4746621187          -46.609234895912                        50      58
+      1135.8000096194          -1171.5142581261                        50      58
+      3187.7055426079           -56.57788934156                        50      58
+      2713.4240465413          -521.10101130376                        50      58
+      1609.7287827223          -27.757030646951                        50      58
+      2218.2808542131          -611.73080288984                        50      58
+      339.77434389417          -501.56273219575                        50      58
+      1628.4095443371           -583.8378231765                        50      58
+      1127.7531780478          -1736.4024554212                        50      58
+      783.45663487986          -911.36971711706                        50      58
+      632.05380854947          -1754.5359385757                        50      58
+      3708.2839439231          -1405.4255247641                        50      58
+      2454.7868603542          -1539.9042746893                        50      58
+      2091.9665388526          -311.17123097013                        50      58
+      2082.1217200227          -1635.0000625566                        50      58
+      3273.4354982145          -1742.8740271813                        50      58
+      3706.0339596851          -848.60652082524                        50      58
+      722.36683469073           42.963518414442                        50      58
+      262.94130347572          -1733.5288300386                        50      58
+       3778.058225824           -1628.058225824                        50      58
+      622.80394092887          -467.13330345548                        50      58
+      1004.3260710173          -1346.6259737243                        50      58
+      3129.5375392764          -925.21562551409                        50      58
+      1920.0539321953           81.716682003892                        50      58
+      264.10502811154           33.186453018978                        50      58
+      3335.9949786522          -564.15567518578                        50      58
+      2971.1394624533          -1381.1740234635                        50      58
+      1351.8207521926          -1026.0372727017                        50      58
+       215.1945587744             -921.44343644                        50      58
+      3777.9039666689          -389.01253121797                        50      58
+      696.60118099605          -1445.3382867397                        50      58
+      728.35605944582          -231.74410256578                        50      58
+      3319.3115215333          -1073.9684694836                        50      58
+      1530.8304268418          -320.52362110518                        50      58
+      1285.7142857142           76.103882369648                        50      58
+      2474.4274814706          -1015.6649119318                        50      58
+      2951.0287882931           71.331937890859                        50      58
+       1520.501381672          -1425.3304293798                        50      58
+       253.6973135114          -1374.8468483467                        50      58
+       1110.952305133          -319.40091038483                        50      58
+      1350.6011578823          -1811.1729631562                        50      58
+      2338.1660764839           65.729833865481                        50      58
+      2858.5107349988          -1760.0160943389                        50      58
+      1627.2275441016          -1039.7316169914                        50      58
+      3665.7543823172          -1786.4216673679                        50      58
+      1923.8707999441          -1275.6821326923                        50      58
+      1085.6123252695          -587.51500034014                        50      58
+      3438.4243515027           68.682495248415                        50      58
+      605.58286398704          -1059.6338500857                        50      58
+       2609.508999673          -311.68836401073                        50      58
+      245.98314002489          -259.25014738251                        50      58
+      2058.4956784593          -885.34925903342                        50      58
+      1764.0325246399          -866.84558719798                        50      58
+       3110.001253574          -276.20095824932                        50      58
+      2265.1599402191          -1488.4548701852                        50      58
+      2468.3827255104           -1797.189068538                        50      58
+      1682.7777746686          -1801.9522820407                        50      58
+      874.99999999995          -1825.0191987092                        50      58
+      2826.2893615091          -1545.2191605369                        50      58
+      2072.1121891632          -1138.5574595677                        50      58
+      798.90642055128           -666.5740256143                        50      58
+      3363.9794589634          -1515.8092628208                        50      58
+      2310.0022045934          -817.77285230863                        50      58
+       956.7723965576          -1038.1425165412                        50      58
+      203.80770423961          -694.60096326135                        50      58
+      1332.7146260126           -731.0012415672                        50      58
+      2736.2963545966          -809.01301810765                        50      58
+       1912.917229461          -1489.1145205076                        50      58
+      2848.8056048724          -1187.6725098473                        50      58
+       1910.311811976          -490.06137012571                        50      58
+      2755.2480385675           99.259138194832                        50      58
+        3137.71838514          -698.10133794125                        50      58
+      3769.8773198013           73.021229696858                        50      58
+      2504.8730852043          -516.86975916373                        50      58
+      1707.6394965756          -1589.2984418475                        50      58
+      1967.7996779093          -1839.7760221533                        50      58
+      2322.1948729117          -158.21196510554                        50      58
+      520.41667523366          -1296.3119318412                        50      58
+      2648.8039626677          -1229.9564880484                        50      58
+       1694.342246968          -1234.1991991544                        50      58
+      2142.5260679531           138.52442710411                        50      58
+      3504.2273541441          -1399.1795629606                        50      58
+       1015.165916226           86.098665248372                        50      58
+      3835.4108901534          -138.50067297114                        50      58
+      476.63184692998           107.35456591656                        50      58
+      3046.5330873028          -1818.6891070206                        50      58
+      3816.8158456071          -1224.5822677841                        50      58
+      2911.6976013938          -987.04932430719                        50      58
+      1288.2093595642          -1287.1762444129                        50      58
+      2633.2953778388          -1455.0724689842                        50      58
+      2909.6711019268          -606.11308583545                        50      58
+      3503.5967245883          -1828.5311691038                        50      58
+        3819.95043295          -657.02581520156                        50      58
+      443.03408776235          -1816.8508568957                        50      58
+      1922.7072374645          -140.38078415628                        50      58
+      2263.0577707571          -1221.0987388152                        50      58
+      165.61538782086          -1558.4628703367                        50      58
+       3820.422655069          -1037.2353091418                        50      58
+      3441.4789343952          -149.50414065605                        50      58
+      1068.5780858665          -1557.7249410525                        50      58
+      3169.8062650757          -1459.8592624457                        50      58
+      124.99999999984          -1856.2499999986                        50      58
+      1144.2671355981           -77.68362293098                        50      58
+      3530.4249014986          -366.64277666758                        50      58
+      2888.3414753778          -392.15715170384                        50      58
+      1352.7439897812          -1626.7393704996                        50      58
+ Edges
+ 103
+ 1 19 1
+ 19 20 1
+ 20 21 1
+ 21 22 1
+ 22 23 1
+ 23 24 1
+ 24 25 1
+ 25 2 1
+ 2 26 2
+ 26 27 2
+ 27 28 2
+ 28 3 2
+ 3 29 3
+ 29 30 3
+ 30 31 3
+ 31 4 3
+ 4 32 4
+ 32 5 4
+ 5 6 5
+ 6 33 6
+ 33 34 6
+ 34 35 6
+ 35 36 6
+ 36 7 6
+ 7 37 7
+ 37 8 7
+ 8 38 8
+ 38 39 8
+ 39 40 8
+ 40 41 8
+ 41 42 8
+ 42 43 8
+ 43 9 8
+ 9 44 9
+ 44 45 9
+ 45 46 9
+ 46 47 9
+ 47 48 9
+ 48 49 9
+ 49 50 9
+ 50 51 9
+ 51 52 9
+ 52 53 9
+ 53 54 9
+ 54 55 9
+ 55 56 9
+ 56 1 9
+ 10 57 13
+ 57 58 13
+ 58 59 13
+ 59 60 13
+ 60 61 13
+ 61 62 13
+ 62 63 13
+ 63 11 13
+ 11 64 14
+ 64 65 14
+ 65 66 14
+ 66 12 14
+ 12 67 15
+ 67 68 15
+ 68 69 15
+ 69 13 15
+ 13 70 16
+ 70 14 16
+ 14 15 17
+ 15 71 18
+ 71 72 18
+ 72 73 18
+ 73 74 18
+ 74 16 18
+ 16 75 19
+ 75 17 19
+ 17 76 20
+ 76 77 20
+ 77 78 20
+ 78 79 20
+ 79 80 20
+ 80 81 20
+ 81 18 20
+ 18 82 21
+ 82 83 21
+ 83 84 21
+ 84 85 21
+ 85 86 21
+ 86 87 21
+ 87 88 21
+ 88 89 21
+ 89 90 21
+ 90 91 21
+ 91 92 21
+ 92 93 21
+ 93 94 21
+ 94 10 21
+ 1 10 23
+ 2 11 24
+ 3 12 28
+ 4 13 32
+ 5 14 36
+ 6 15 40
+ 7 16 44
+ 8 17 48
+ 9 18 52
+ Triangles
+ 772
+ 162 236 123 11
+ 123 236 101 11
+ 33 178 124 11
+ 99 175 126 11
+ 23 125 164 11
+ 110 129 169 11
+ 44 183 208 11
+ 23 174 125 11
+ 115 143 188 11
+ 1 160 56 11
+ 1 19 160 11
+ 184 108 125 11
+ 127 171 231 11
+ 107 134 161 11
+ 21 142 186 11
+ 112 142 156 11
+ 115 187 141 11
+ 148 127 190 11
+ 125 108 164 11
+ 101 217 132 11
+ 99 126 167 11
+ 141 187 198 11
+ 5 191 148 11
+ 95 143 201 11
+ 51 52 140 11
+ 131 208 183 11
+ 30 192 128 11
+ 171 98 231 11
+ 107 132 152 11
+ 101 132 168 11
+ 30 128 176 11
+ 20 21 186 11
+ 127 231 190 11
+ 107 168 132 11
+ 52 170 140 11
+ 44 45 183 11
+ 129 201 143 11
+ 95 188 143 11
+ 211 130 31 11
+ 101 168 123 11
+ 46 138 183 11
+ 100 149 141 11
+ 33 124 191 11
+ 107 152 134 11
+ 112 156 167 11
+ 117 202 139 11
+ 137 199 157 11
+ 110 175 129 11
+ 33 34 178 11
+ 129 143 169 11
+ 41 165 228 11
+ 130 192 31 11
+ 199 97 157 11
+ 109 206 136 11
+ 116 166 146 11
+ 110 140 170 11
+ 115 141 205 11
+ 8 38 155 11
+ 134 228 165 11
+ 48 177 136 11
+ 115 205 143 11
+ 111 146 166 11
+ 105 230 140 11
+ 139 202 226 11
+ 45 46 183 11
+ 19 186 160 11
+ 107 207 158 11
+ 122 173 215 11
+ 124 148 191 11
+ 117 209 141 11
+ 127 225 214 11
+ 127 148 225 11
+ 19 20 186 11
+ 7 180 227 11
+ 161 207 107 11
+ 5 150 32 11
+ 46 172 138 11
+ 41 42 165 11
+ 134 239 161 11
+ 104 135 161 11
+ 100 141 209 11
+ 122 179 163 11
+ 117 139 209 11
+ 141 149 205 11
+ 3 29 144 11
+ 109 136 185 11
+ 51 140 159 11
+ 108 156 142 11
+ 105 140 169 11
+ 95 163 179 11
+ 113 196 145 11
+ 4 211 31 11
+ 102 215 173 11
+ 104 138 189 11
+ 113 137 163 11
+ 108 142 200 11
+ 100 212 149 11
+ 132 217 147 11
+ 104 189 135 11
+ 110 169 140 11
+ 109 189 138 11
+ 109 172 206 11
+ 109 138 172 11
+ 22 23 164 11
+ 127 214 171 11
+ 5 6 191 11
+ 234 138 104 11
+ 23 24 174 11
+ 117 141 198 11
+ 113 182 196 11
+ 32 150 211 11
+ 140 230 159 11
+ 48 136 206 11
+ 54 153 218 11
+ 126 175 238 11
+ 38 147 155 11
+ 236 151 197 11
+ 98 198 187 11
+ 100 185 212 11
+ 136 212 185 11
+ 54 55 153 11
+ 27 28 146 11
+ 48 49 177 11
+ 97 184 213 11
+ 105 143 205 11
+ 116 157 166 11
+ 113 199 137 11
+ 7 227 36 11
+ 97 166 157 11
+ 118 151 236 11
+ 126 218 153 11
+ 108 145 196 11
+ 112 186 142 11
+ 99 167 156 11
+ 21 200 142 11
+ 133 171 214 11
+ 116 193 144 11
+ 106 224 157 11
+ 104 161 239 11
+ 123 168 158 11
+ 97 213 166 11
+ 137 157 224 11
+ 3 144 193 11
+ 119 155 147 11
+ 120 149 212 11
+ 105 169 143 11
+ 8 155 37 11
+ 108 196 156 11
+ 46 47 172 11
+ 131 234 239 11
+ 108 184 145 11
+ 35 36 151 11
+ 131 239 165 11
+ 50 51 159 11
+ 114 190 231 11
+ 122 163 224 11
+ 116 144 235 11
+ 29 176 144 11
+ 106 173 224 11
+ 108 200 164 11
+ 137 224 163 11
+ 111 229 146 11
+ 130 211 150 11
+ 111 232 154 11
+ 113 201 182 11
+ 106 144 241 11
+ 96 226 202 11
+ 28 193 146 11
+ 103 214 225 11
+ 133 202 171 11
+ 30 31 192 11
+ 106 241 173 11
+ 122 215 179 11
+ 130 203 210 11
+ 116 146 193 11
+ 121 195 187 11
+ 106 235 144 11
+ 29 30 176 11
+ 111 154 229 11
+ 117 171 202 11
+ 130 150 203 11
+ 5 148 190 11
+ 126 153 167 11
+ 97 145 184 11
+ 112 167 153 11
+ 111 174 232 11
+ 103 162 204 11
+ 129 182 201 11
+ 26 229 154 11
+ 97 199 145 11
+ 113 145 199 11
+ 9 208 43 11
+ 9 44 208 11
+ 119 147 217 11
+ 52 53 170 11
+ 37 155 180 11
+ 26 27 229 11
+ 5 190 150 11
+ 107 158 168 11
+ 6 33 191 11
+ 112 160 186 11
+ 123 204 162 11
+ 112 153 220 11
+ 114 203 150 11
+ 109 185 240 11
+ 98 187 195 11
+ 105 149 230 11
+ 96 223 158 11
+ 123 158 223 11
+ 38 222 147 11
+ 120 230 149 11
+ 7 37 180 11
+ 119 197 151 11
+ 113 163 201 11
+ 139 240 185 11
+ 120 216 159 11
+ 2 237 25 11
+ 112 220 160 11
+ 110 238 175 11
+ 103 194 162 11
+ 27 146 229 11
+ 53 218 170 11
+ 105 205 149 11
+ 122 224 173 11
+ 50 159 216 11
+ 102 210 203 11
+ 114 150 190 11
+ 119 151 227 11
+ 121 187 188 11
+ 34 35 221 11
+ 102 173 210 11
+ 119 180 155 11
+ 133 223 202 11
+ 21 22 200 11
+ 118 221 151 11
+ 144 176 241 11
+ 132 147 222 11
+ 128 210 173 11
+ 55 56 220 11
+ 170 218 238 11
+ 109 240 189 11
+ 129 175 182 11
+ 115 188 187 11
+ 121 181 195 11
+ 38 39 222 11
+ 99 156 196 11
+ 125 213 184 11
+ 56 160 220 11
+ 35 151 221 11
+ 3 193 28 11
+ 99 196 182 11
+ 40 228 152 11
+ 139 185 209 11
+ 99 182 175 11
+ 40 152 233 11
+ 26 154 237 11
+ 100 209 185 11
+ 118 162 194 11
+ 136 177 212 11
+ 22 164 200 11
+ 47 48 206 11
+ 132 233 152 11
+ 128 173 241 11
+ 124 225 148 11
+ 124 178 194 11
+ 134 152 228 11
+ 111 213 174 11
+ 40 41 228 11
+ 39 40 233 11
+ 36 227 151 11
+ 135 207 161 11
+ 95 179 188 11
+ 114 181 203 11
+ 135 189 240 11
+ 125 174 213 11
+ 96 202 223 11
+ 138 234 183 11
+ 96 158 207 11
+ 34 221 178 11
+ 98 195 231 11
+ 55 220 153 11
+ 116 235 157 11
+ 24 25 232 11
+ 121 188 179 11
+ 49 50 216 11
+ 42 43 219 11
+ 4 32 211 11
+ 53 54 218 11
+ 98 171 198 11
+ 95 201 163 11
+ 106 157 235 11
+ 120 159 230 11
+ 114 195 181 11
+ 25 154 232 11
+ 117 198 171 11
+ 25 237 154 11
+ 118 194 178 11
+ 124 194 225 11
+ 2 26 237 11
+ 47 206 172 11
+ 130 210 192 11
+ 131 165 219 11
+ 103 225 194 11
+ 111 166 213 11
+ 42 219 165 11
+ 102 203 181 11
+ 104 239 234 11
+ 134 165 239 11
+ 128 192 210 11
+ 118 236 162 11
+ 120 212 177 11
+ 133 204 223 11
+ 103 204 214 11
+ 110 170 238 11
+ 120 177 216 11
+ 121 179 215 11
+ 133 214 204 11
+ 102 181 215 11
+ 118 178 221 11
+ 49 216 177 11
+ 24 232 174 11
+ 126 238 218 11
+ 121 215 181 11
+ 128 241 176 11
+ 96 207 226 11
+ 131 183 234 11
+ 119 217 197 11
+ 101 197 217 11
+ 135 226 207 11
+ 114 231 195 11
+ 101 236 197 11
+ 131 219 208 11
+ 119 227 180 11
+ 43 208 219 11
+ 123 223 204 11
+ 132 222 233 11
+ 139 226 240 11
+ 135 240 226 11
+ 39 233 222 11
+ 19 20 57 25
+ 57 20 58 25
+ 60 23 61 25
+ 60 22 23 25
+ 61 24 62 25
+ 61 23 24 25
+ 10 19 57 25
+ 10 1 19 25
+ 20 59 58 25
+ 20 21 59 25
+ 59 22 60 25
+ 59 21 22 25
+ 25 62 24 25
+ 25 63 62 25
+ 11 25 2 25
+ 11 63 25 25
+ 11 2 64 29
+ 2 26 64 29
+ 26 65 64 29
+ 27 65 26 29
+ 65 28 66 29
+ 65 27 28 29
+ 66 3 12 29
+ 66 28 3 29
+ 12 3 67 33
+ 3 29 67 33
+ 68 30 69 33
+ 30 31 69 33
+ 69 31 13 33
+ 31 4 13 33
+ 67 30 68 33
+ 67 29 30 33
+ 13 32 70 37
+ 13 4 32 37
+ 32 14 70 37
+ 32 5 14 37
+ 14 5 15 41
+ 5 6 15 41
+ 6 71 15 45
+ 6 33 71 45
+ 33 34 71 45
+ 71 34 72 45
+ 34 73 72 45
+ 34 35 73 45
+ 36 73 35 45
+ 74 73 36 45
+ 16 36 7 45
+ 16 74 36 45
+ 7 37 16 49
+ 16 37 75 49
+ 75 37 17 49
+ 37 8 17 49
+ 17 38 76 53
+ 17 8 38 53
+ 38 39 76 53
+ 76 39 77 53
+ 81 43 18 53
+ 43 9 18 53
+ 77 40 78 53
+ 77 39 40 53
+ 40 79 78 53
+ 41 79 40 53
+ 79 42 80 53
+ 79 41 42 53
+ 80 43 81 53
+ 42 43 80 53
+ 18 9 82 57
+ 9 44 82 57
+ 44 45 82 57
+ 82 45 83 57
+ 83 45 84 57
+ 45 46 84 57
+ 46 47 84 57
+ 84 47 85 57
+ 49 50 88 57
+ 49 88 87 57
+ 51 90 89 57
+ 51 52 90 57
+ 54 55 92 57
+ 92 55 93 57
+ 93 56 94 57
+ 93 55 56 57
+ 48 85 47 57
+ 86 85 48 57
+ 49 86 48 57
+ 87 86 49 57
+ 89 50 51 57
+ 89 88 50 57
+ 52 91 90 57
+ 52 53 91 57
+ 91 54 92 57
+ 91 53 54 57
+ 94 1 10 57
+ 94 56 1 57
+ 309 383 270 58
+ 270 383 248 58
+ 71 325 271 58
+ 246 322 273 58
+ 61 272 311 58
+ 257 276 316 58
+ 82 330 355 58
+ 61 321 272 58
+ 262 290 335 58
+ 10 307 94 58
+ 10 57 307 58
+ 331 255 272 58
+ 274 318 378 58
+ 254 281 308 58
+ 59 289 333 58
+ 259 289 303 58
+ 262 334 288 58
+ 295 274 337 58
+ 272 255 311 58
+ 248 364 279 58
+ 246 273 314 58
+ 288 334 345 58
+ 14 338 295 58
+ 242 290 348 58
+ 89 90 287 58
+ 278 355 330 58
+ 68 339 275 58
+ 318 245 378 58
+ 254 279 299 58
+ 248 279 315 58
+ 68 275 323 58
+ 58 59 333 58
+ 274 378 337 58
+ 254 315 279 58
+ 90 317 287 58
+ 82 83 330 58
+ 276 348 290 58
+ 242 335 290 58
+ 358 277 69 58
+ 248 315 270 58
+ 84 285 330 58
+ 247 296 288 58
+ 71 271 338 58
+ 254 299 281 58
+ 259 303 314 58
+ 264 349 286 58
+ 284 346 304 58
+ 257 322 276 58
+ 71 72 325 58
+ 276 290 316 58
+ 79 312 375 58
+ 277 339 69 58
+ 346 244 304 58
+ 256 353 283 58
+ 263 313 293 58
+ 257 287 317 58
+ 262 288 352 58
+ 17 76 302 58
+ 281 375 312 58
+ 86 324 283 58
+ 262 352 290 58
+ 258 293 313 58
+ 252 377 287 58
+ 286 349 373 58
+ 83 84 330 58
+ 57 333 307 58
+ 254 354 305 58
+ 269 320 362 58
+ 271 295 338 58
+ 264 356 288 58
+ 274 372 361 58
+ 274 295 372 58
+ 57 58 333 58
+ 16 327 374 58
+ 308 354 254 58
+ 14 297 70 58
+ 84 319 285 58
+ 79 80 312 58
+ 281 386 308 58
+ 251 282 308 58
+ 247 288 356 58
+ 269 326 310 58
+ 264 286 356 58
+ 288 296 352 58
+ 12 67 291 58
+ 256 283 332 58
+ 89 287 306 58
+ 255 303 289 58
+ 252 287 316 58
+ 242 310 326 58
+ 260 343 292 58
+ 13 358 69 58
+ 249 362 320 58
+ 251 285 336 58
+ 260 284 310 58
+ 255 289 347 58
+ 247 359 296 58
+ 279 364 294 58
+ 251 336 282 58
+ 257 316 287 58
+ 256 336 285 58
+ 256 319 353 58
+ 256 285 319 58
+ 60 61 311 58
+ 274 361 318 58
+ 14 15 338 58
+ 381 285 251 58
+ 61 62 321 58
+ 264 288 345 58
+ 260 329 343 58
+ 70 297 358 58
+ 287 377 306 58
+ 86 283 353 58
+ 92 300 365 58
+ 273 322 385 58
+ 76 294 302 58
+ 383 298 344 58
+ 245 345 334 58
+ 247 332 359 58
+ 283 359 332 58
+ 92 93 300 58
+ 65 66 293 58
+ 86 87 324 58
+ 244 331 360 58
+ 252 290 352 58
+ 263 304 313 58
+ 260 346 284 58
+ 16 374 74 58
+ 244 313 304 58
+ 265 298 383 58
+ 273 365 300 58
+ 255 292 343 58
+ 259 333 289 58
+ 246 314 303 58
+ 59 347 289 58
+ 280 318 361 58
+ 263 340 291 58
+ 253 371 304 58
+ 251 308 386 58
+ 270 315 305 58
+ 244 360 313 58
+ 284 304 371 58
+ 12 291 340 58
+ 266 302 294 58
+ 267 296 359 58
+ 252 316 290 58
+ 17 302 75 58
+ 255 343 303 58
+ 84 85 319 58
+ 278 381 386 58
+ 255 331 292 58
+ 73 74 298 58
+ 278 386 312 58
+ 88 89 306 58
+ 261 337 378 58
+ 269 310 371 58
+ 263 291 382 58
+ 67 323 291 58
+ 253 320 371 58
+ 255 347 311 58
+ 284 371 310 58
+ 258 376 293 58
+ 277 358 297 58
+ 258 379 301 58
+ 260 348 329 58
+ 253 291 388 58
+ 243 373 349 58
+ 66 340 293 58
+ 250 361 372 58
+ 280 349 318 58
+ 68 69 339 58
+ 253 388 320 58
+ 269 362 326 58
+ 277 350 357 58
+ 263 293 340 58
+ 268 342 334 58
+ 253 382 291 58
+ 67 68 323 58
+ 258 301 376 58
+ 264 318 349 58
+ 277 297 350 58
+ 14 295 337 58
+ 273 300 314 58
+ 244 292 331 58
+ 259 314 300 58
+ 258 321 379 58
+ 250 309 351 58
+ 276 329 348 58
+ 64 376 301 58
+ 244 346 292 58
+ 260 292 346 58
+ 18 355 81 58
+ 18 82 355 58
+ 266 294 364 58
+ 90 91 317 58
+ 75 302 327 58
+ 64 65 376 58
+ 14 337 297 58
+ 254 305 315 58
+ 15 71 338 58
+ 259 307 333 58
+ 270 351 309 58
+ 259 300 367 58
+ 261 350 297 58
+ 256 332 387 58
+ 245 334 342 58
+ 252 296 377 58
+ 243 370 305 58
+ 270 305 370 58
+ 76 369 294 58
+ 267 377 296 58
+ 16 75 327 58
+ 266 344 298 58
+ 260 310 348 58
+ 286 387 332 58
+ 267 363 306 58
+ 11 384 63 58
+ 259 367 307 58
+ 257 385 322 58
+ 250 341 309 58
+ 65 293 376 58
+ 91 365 317 58
+ 252 352 296 58
+ 269 371 320 58
+ 88 306 363 58
+ 249 357 350 58
+ 261 297 337 58
+ 266 298 374 58
+ 268 334 335 58
+ 72 73 368 58
+ 249 320 357 58
+ 266 327 302 58
+ 280 370 349 58
+ 59 60 347 58
+ 265 368 298 58
+ 291 323 388 58
+ 279 294 369 58
+ 275 357 320 58
+ 93 94 367 58
+ 317 365 385 58
+ 256 387 336 58
+ 276 322 329 58
+ 262 335 334 58
+ 268 328 342 58
+ 76 77 369 58
+ 246 303 343 58
+ 272 360 331 58
+ 94 307 367 58
+ 73 298 368 58
+ 12 340 66 58
+ 246 343 329 58
+ 78 375 299 58
+ 286 332 356 58
+ 246 329 322 58
+ 78 299 380 58
+ 64 301 384 58
+ 247 356 332 58
+ 265 309 341 58
+ 283 324 359 58
+ 60 311 347 58
+ 85 86 353 58
+ 279 380 299 58
+ 275 320 388 58
+ 271 372 295 58
+ 271 325 341 58
+ 281 299 375 58
+ 258 360 321 58
+ 78 79 375 58
+ 77 78 380 58
+ 74 374 298 58
+ 282 354 308 58
+ 242 326 335 58
+ 261 328 350 58
+ 282 336 387 58
+ 272 321 360 58
+ 243 349 370 58
+ 285 381 330 58
+ 243 305 354 58
+ 72 368 325 58
+ 245 342 378 58
+ 93 367 300 58
+ 263 382 304 58
+ 62 63 379 58
+ 268 335 326 58
+ 87 88 363 58
+ 80 81 366 58
+ 13 70 358 58
+ 91 92 365 58
+ 245 318 345 58
+ 242 348 310 58
+ 253 304 382 58
+ 267 306 377 58
+ 261 342 328 58
+ 63 301 379 58
+ 264 345 318 58
+ 63 384 301 58
+ 265 341 325 58
+ 271 341 372 58
+ 11 64 384 58
+ 85 353 319 58
+ 277 357 339 58
+ 278 312 366 58
+ 250 372 341 58
+ 258 313 360 58
+ 80 366 312 58
+ 249 350 328 58
+ 251 386 381 58
+ 281 312 386 58
+ 275 339 357 58
+ 265 383 309 58
+ 267 359 324 58
+ 280 351 370 58
+ 250 351 361 58
+ 257 317 385 58
+ 267 324 363 58
+ 268 326 362 58
+ 280 361 351 58
+ 249 328 362 58
+ 265 325 368 58
+ 87 363 324 58
+ 62 379 321 58
+ 273 385 365 58
+ 268 362 328 58
+ 275 388 323 58
+ 243 354 373 58
+ 278 330 381 58
+ 266 364 344 58
+ 248 344 364 58
+ 282 373 354 58
+ 261 378 342 58
+ 248 383 344 58
+ 278 366 355 58
+ 266 374 327 58
+ 81 355 366 58
+ 270 370 351 58
+ 279 369 380 58
+ 286 373 387 58
+ 282 387 373 58
+ 77 380 369 58
+ End

=== added directory 'examples/sph/cpt'
=== added file 'examples/sph/cpt/a'
=== added file 'examples/sph/watercolumn.py'
--- examples/sph/watercolumn.py	1970-01-01 00:00:00 +0000
+++ examples/sph/watercolumn.py	2014-04-09 14:03:16 +0000
@@ -0,0 +1,63 @@
+#!/usr/bin/env python
+# encoding: utf-8
+
+from yade import utils, plot, qt, ymport
+o = Omega()
+
+# Physical parameters
+fr = 0.0;
+rho=1000.0
+
+k = 5.0
+mu = 0.001
+tc = 0.0001; en = 0.7; et = 0.7;
+vel = 0.05
+Rad = 15.0e-3
+o.dt = 0.0002
+
+
+scaleF =  0.001
+
+# Add material
+mat1 = O.materials.append(ViscElMat(frictionAngle=fr,density=rho, SPHmode=True,mu=mu,tc=tc, en=en, et=et))
+id1 = O.bodies.append(ymport.gmsh("box.mesh", scale=scaleF, material=mat1, color=(1,0,0), mask = 1, wire=True))
+
+d = 15.0*scaleF
+
+print d
+idSpheres = O.bodies.append(
+  pack.regularHexa(
+    pack.inAlignedBox(
+    (0,       -2000.0*scaleF,  0.0),
+    (1000*scaleF, 200*scaleF,  50.0*scaleF)),
+    radius=d,gap=0.01*d, material=mat1, mask=1, color=(0,1,1)))
+
+
+print len(idSpheres)
+
+# Add engines
+o.engines = [
+  ForceResetter(),
+  InsertionSortCollider([Bo1_Sphere_Aabb(label='is2aabb'),Bo1_Facet_Aabb(label='is3aabb')]),
+  InteractionLoop(
+    [Ig2_Sphere_Sphere_ScGeom(label='ss2sc'),Ig2_Facet_Sphere_ScGeom()],
+    [Ip2_ViscElMat_ViscElMat_ViscElPhys()],
+    [Law2_ScGeom_ViscElPhys_Basic()],
+  ),
+  NewtonIntegrator(damping=0.05,gravity=[0,-9.81,0]),
+  SPHEngine(mask=1, k=k, rho0 = rho),
+  VTKRecorder(iterPeriod=100,fileName='./cpt/spheres-', recorders=['spheres','velocity','colors','intr','ids','mask','materialId','stress']),
+  VTKRecorder(iterPeriod=100,fileName='./cpt/facet-',   recorders=['facets'],label='VTK_box2'),
+  PyRunner(command='addPlotData()',iterPeriod=50,dead=False),
+]
+
+def addPlotData(): 
+  plot.addData(t=O.time, Ekin=utils.kineticEnergy())
+
+plot.plots={'t':('Ekin')}; plot.plot()
+
+O.run(1, True)
+
+qt.View()
+#O.run(10000, True)
+#plot.saveGnuplot('sim-data_0.05')