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}, ], }