yade-dev team mailing list archive
-
yade-dev team
-
Mailing list archive
-
Message #10703
[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')