Beispiel #1
0
coefStatic[mask1] = 1.0e+6
mask1 = points[:,1] >= +14.999e+3
coefStatic[mask1] = 1.0e+6
mask1 = points[:,2] <= -14.999e+3
coefStatic[mask1] = 1.0e+6


cs = CSCart()
cs._configure()
cs.initialize()

writer = SimpleIOAscii()
writer.inventory.filename = "empty"
writer._configure()

writer.filename("friction_main.spatialdb")
dataOut = {'points': points,
           'coordsys': cs,
           'data_dim': 2,
           'values': [{'name': 'static-coefficient', 
                       'units': 'None',
                       'data': coefStatic},
                      {'name': 'dynamic-coefficient', 
                       'units': 'None',
                       'data': coefDyn},
                      {'name': 'slip-weakening-parameter', 
                       'units': 'm',
                       'data': d0},
                      {'name': 'cohesion', 
                       'units': 'Pa',
                       'data': cohesion},
Szz = density * gacc * z

Syy = mask1 * (b22 * (Szz + Pf) - 0 * Pf) + ~mask1 * Szz
Sxx = mask1 * (b33 * (Szz + Pf) - 0 * Pf) + ~mask1 * Szz
Sxy = mask1 * (-b23 * (Szz + Pf)) + ~mask1 * 0.0

Syz = 0.0 * z
Sxz = 0.0 * z

writer = SimpleIOAscii()
writer.inventory.filename = "empty"
writer._configure()

# ----------------------------------------------------------------------
# Initial stresses in domain
writer.filename("%s_initial_stress.spatialdb" % sim)
dataOut = {
    'points':
    points,
    'coordsys':
    cs,
    'data_dim':
    1,
    'values': [
        {
            'name': 'stress-xx',
            'units': 'Pa',
            'data': Sxx
        },
        {
            'name': 'stress-yy',
Syy = mask1*(b22*(Szz+Pf)-0*Pf) + ~mask1*Szz
Sxx = mask1*(b33*(Szz+Pf)-0*Pf) + ~mask1*Szz
Sxy = mask1*(-b23*(Szz+Pf)) + ~mask1*0.0

Syz = 0.0*z
Sxz = 0.0*z

writer = SimpleIOAscii()
writer.inventory.filename = "empty"
writer._configure()


# ----------------------------------------------------------------------
# Initial stresses in domain
writer.filename("%s_initial_stress.spatialdb" % sim)
dataOut = {'points': points,
           'coordsys': cs,
           'data_dim': 1,
           'values': [{'name': 'stress-xx', 
                       'units': 'Pa',
                       'data': Sxx},
                      {'name': 'stress-yy', 
                       'units': 'Pa',
                       'data': Syy},
                      {'name': 'stress-zz', 
                       'units': 'Pa',
                       'data': Szz},
                      {'name': 'stress-xy', 
                       'units': 'Pa',
                       'data': Sxy},
coefStatic = data[:, 9]
coefDyn = data[:, 10]
d0 = data[:, 11]
cohesion = data[:, 12]
weakTime = data[:, 13]

cs = CSCart()
cs._configure()
cs.initialize()

writer = SimpleIOAscii()
writer.inventory.filename = "empty"
writer._configure()

writer.filename("%s_tractions.spatialdb" % sim)
dataOut = {
    'points':
    points,
    'coordsys':
    cs,
    'data_dim':
    2,
    'values': [
        {
            'name': 'traction-shear-leftlateral',
            'units': 'Pa',
            'data': tractionShearHoriz
        },
        {
            'name': 'traction-shear-updip',
coefStatic = data[:,9]
coefDyn = data[:,10]
d0 = data[:,11]
cohesion = data[:,12]
weakTime = data[:,13]


cs = CSCart()
cs._configure()
cs.initialize()

writer = SimpleIOAscii()
writer.inventory.filename = "empty"
writer._configure()

writer.filename("%s_tractions.spatialdb" % sim)
dataOut = {'points': points,
           'coordsys': cs,
           'data_dim': 2,
           'values': [{'name': 'traction-shear-leftlateral', 
                       'units': 'Pa',
                       'data': tractionShearHoriz},
                      {'name': 'traction-shear-updip',
                       'units': 'Pa',
                       'data': tractionShearVert},
                      {'name': 'traction-normal',
                       'units': 'Pa',
                       'data': tractionNormal},
                      ],
           }