← Back to branch summary

~esys-p-dev/esys-particle/esys-foam-dem

« back to all changes in this revision

Viewing changes to trunk/example/params.py

  • Committer: Feng Chen
  • Date: 2012-06-16 03:15:16 UTC
  • Revision ID: fchen3@gmail.com-20120616031516-0m2hw2zyhdkhu1uv
add 2d sheetpile problem files to repository

Show diffs side-by-side

added added

removed removed

Lines of Context:
trunk/example/params.py
 
1
"""
 
2
"""
 
3
import sys
 
4
import math
 
5
import fileinput 
 
6
import time
 
7
# from esys.lsm import *
 
8
# from esys.lsm.util import *
 
9
# from esys.lsm.geometry import *
 
10
from LsmPy import *
 
11
from FoundationPy import Vec3, BoundingBox
 
12
from GeometryPy import *
 
13
 
 
14
N_total = 5000 #int(sys.argv[1]) #500000
 
15
N_int = 10 #int(sys.argv[2]) #1000
 
16
ndf = 10 #int(sys.argv[3])
 
17
 
 
18
npx=1 #int(sys.argv[4])
 
19
npy=1 #int(sys.argv[5])
 
20
npz=1 #int(sys.argv[6])
 
21
 
 
22
foamCaseName="2dshpl" #str(sys.argv[7])
 
23
 
 
24
dem_fvm_combined_incr = 1
 
25
 
 
26
timeStepChkIncrDEM= ndf*dem_fvm_combined_incr
 
27
 
 
28
minPoint=Vec3(-200,-200,-200)
 
29
maxPoint=Vec3(200,200,200)
 
30
 
 
31
dt = 5.0E-6
 
32
 
 
33
rad=0.5
 
34
density=1.0
 
35
vol=4.0/3.0*math.pi*rad*rad*rad
 
36
mass_p=vol*density
 
37
 
 
38
kn = (3.77E+07) #*2
 
39
 
 
40
gravityConst=Vec3(0.0,-9806.65,0.0)
 
41
 
 
42
rhos = density
 
43
rhof = 1000.0/2650.0*density
 
44
grav = -gravityConst[1]
 
45
fbg = (rhos - rhof)*grav*vol
 
46
delta_fbg = fbg/(2.0*kn)
 
47
 
 
48
geoFileName="shpl100.geo" #"geo_2dshpl.geo" #"shpl100.geo"
 
49
mshName="2dshplWall"
 
50
mshFilename="test_mesh2d.msh"
 
51
fileNamePrefix = "bingle_data"
 
52
 
 
53
foamLibPathName="/home/fchen/ESyS-Particle-2.0/icoFoamESyS/Make/linuxGccDPOpt/"
 
54
foamLibFileName="libIcoFoam.so"
 
55
particleLevel=1
 
56
fixParticle=0
 
57
 
 
58
dampingViscosity=2.0
 
59
 
 
60
sx = 0.0
 
61
sy = 0.0
 
62
wx = 20.0
 
63
wy = 10.0
 
64
zt = 1.0
 
65
zb = 0.0
 
66
swx = sx + wx
 
67
swy = sy + wy
 
68
s2wx = sx + 2.0*swx
 
69
s2wy = sy + 2.0*swy