flac中修正剑桥本构模型源码中文注解

#include "usercam.h"

#include <math.h>

static const unsigned long mShearNow = 0x01; 定义静态全局常量 （0x08 表示退格）

static const unsigned long mShearPast = 0x02;

static const double dC1D3 = 1.0 / 3.0;

static const double dC2D3 = 2.0 / 3.0;

static const double dC4D3 = 4.0 / 3.0;

static UserCamClayModel usercamclaymodel(true);

力 ；基准压力 =1pa

dMV_L(0.0) 在基准压力下，正常固结线的比容 =3.3 2, dMV(0.0) 初始比容 v0,

dEV(0.0) 累计总体积应变, dEV_P(0.0) 累计塑性体积应变, dMP(0.0) 初始压力（当前的平

均有效应力）, dMQ(0.0) 当前的平均偏应力{

}；

const char *UserCamClayModel::Keyword(void) const { return("usercam-clay"); }

const char *UserCamClayModel::Name(void) const { return("User-Cam-Clay"); }

const char **UserCamClayModel::Properties(void) const {

static const char *strKey[] = {与命令流交换时的形式

"bulk", "shear", "bulk_bound", "poisson", "kappa", "lambda",

const char **UserCamClayModel::States(void) const {

static const char *strKey[] = {

"shear-n", "shear-p", 0

};

return(strKey);

}

double UserCamClayModel::GetProperty(unsigned ul) const {

switch (ul) {

case 1: return(dBulk); 弹性体积模量

case 2: return(dShear); 剪切模量 G=250·

case 4: return(dPoisson); 泊松比

case 5: return(dKappa); 弹性膨胀线的斜率

case 6: return(dLambda); 正常固结线的斜率

case 7: return(dMM); 摩擦常数 M

case 14: return(dMP); 初始压力

case 15: return(dMQ); 当前的平均偏应力

}

return(0.0);

}

void UserCamClayModel::SetProperty(unsigned ul,const double &dVal) {

switch (ul) {

case 1: dBulk = dVal; break; 弹性体积模量 K

case 2: dShear = dVal; break; 剪切模量 G=250·

case 3: dBulkB = dVal; break;

case 4: dPoisson = dVal; break; 泊松比

case 5: dKappa = dVal; break; 弹性膨胀线的斜率

case 6: dLambda = dVal; break; 正常固结线的斜率

case 7: dMM = dVal; break; 摩擦常数 M

case 14: dMP = dVal; break; 初始压力

case 15: dMQ = dVal; break; 当前的平均偏应力

}

}

if (str) return(str);

UserCamClayModel *em = (UserCamClayModel *)m;

dBulk = em->dBulk; 弹性体积模量 K ->是分量运算符，为成员选择

dBulkB = em->dBulkB;

dPoisson = em->dPoisson; 泊松比

dKappa = em->dKappa; 弹性膨胀线的斜率

dLambda = em->dLambda; 正常固结线的斜率

dMPC = em->dMPC; 先期固结压力

dEV_P = em->dEV_P; 累计塑性体积应变

return(0);

}

// Initialize mean pressure