在我的理学学士论文中,我必须研究流体注入井周围的孔隙压力和应力分布。我试着用Spyder来解决这一部分,因为我在编码方面几乎没有经验,所以它似乎是我最好的界面。在
尽管我研究了最重要的Python和matplotlib文档,但我无法找到解决问题的方法。在
首先,我的代码是:# -*- coding: utf-8 -*-
"""
Created on Wed Oct 01 10:26:29 2014
@author: Alexander
"""
from mpl_toolkits.mplot3d import Axes3D
from matplotlib import cm
from math import *
from scipy.special import *
import matplotlib.pyplot as plt
import numpy as np
q = 6.0/1000
rhof = 1000
lameu = 11.2*10**9
lame = 8.4*10**9
pi
alpha = 0.8
G = 8.4*10**9
k = 4*10**(-14)
eta = 0.001
t = 10*365*24*3600
kappa = k/eta
print "kappa ist:",kappa
c = ((kappa*(lameu-lame)*(lame+2*G))/((alpha**2)*(