# encoding: utf-8
from scipy.optimize import fsolve
from numpy import *
beita = 1
gama = 1
#以下为Mohr-Columb准则计算
# def f(x):
# OH = (0.33/(1-0.33)+beita)*(14.12-0.73*4.68)+0.73*4.68
# Oh = (0.33/(1-0.33)+gama)*(14.12-0.73*4.68)+0.73*4.68
# Or = 0.946*x+0.25272
# Osita = 3*OH-Oh-x*0.683552239-1.481
# Oz = 12.639+0.66*(OH-Oh)+0.316447761*x
# x = float(x)
# return Osita-3.4164-(Or-3.4164)*4.4835+10.29
# result1 = fsolve(f, 7)
# print '摩尔-库伦准则计算结果:', result1*10/6/9.8, 'g/cm^3'
# 以下为Druker-Prager准则计算
def g(y):
OH = (0.33/(1-0.33)+beita)*(14.12-0.73*4.68)+0.73*4.68
Oh = (0.33/(1-0.33)+gama)*(14.12-0.73*4.68)+0.73*4.68
Or = 0.946*y+0.25272
Osita = 3*OH-Oh-y*0.683552239-1.481
Oz = 12.639+0.66*(OH-Oh)+0.316447761*y
y = float(y)
return (math.sqrt(1/6*((Or-Osita)**2&
Python解非线性方程
最新推荐文章于 2024-07-29 00:58:39 发布
fsolve函数在寻找非线性方程的根时局限性明显,仅能求解附近根,全局根的搜索仍面临难题,对解方程技术的需求持续存在。
摘要由CSDN通过智能技术生成