谁能解释下cam clay模型中 dBa dBb dBoa dAlam 的含义？

zhaoyang107

谁能解释下cam clay模型中 dBa dBb dBoa   dAlam 的含义？

见以下代码


ps->iworking[Pind] = 1;
      ps->mState |= mShearNow;
      double dSa    = 6.0 * dShear * dQVal;
      double dSc    = dM2 * (2.0 * dPVal - dMPC);
      double dSb    = dBulk * dSc;
      double dBa    = dSa*dSa + dM2 * dSb*dSb;
      double dBb    = dSa * dQVal + 0.5 * dSb * dSc;
      double dBc    = dFS;
      double dAlam  = 0.0;
      double dAlam1 = 0.0;
      if (dBa != 0.0) {
        double dBoa   = dBb / dBa;
        double dVal   = dBoa*dBoa - dBc / dBa;
        if (dVal < 0.0) return("Cam-clay:  Yield envelope cannot be reached");
        double dVal1   = sqrt(dVal);
        dAlam  = dBoa + dVal1;
        dAlam1 = dBoa - dVal1;
      } else {
        if (dBb != 0.0) {
           dAlam  = 0.5 * dBc / dBb;
           dAlam1 = dAlam;

zhaoyang107

明白了，这个要求根dBa dBb dBoa 是求根公式的系数