int SteelZ02::setTrialStrain (double strain, double strainRate) { // Reset history variables to last converged state TminStrain = CminStrain; TmaxStrain = CmaxStrain; TloadingState = CloadingState; TloopPathState = CloopPathState; for (int i=0; i<SIZE; i++) { TreverseTopStrain[i] = CreverseTopStrain[i]; TreverseTopStress[i] = CreverseTopStress[i]; TreverseBottomStrain[i] = CreverseBottomStrain[i]; TreverseBottomStress[i] = CreverseBottomStress[i]; } TreverseTopNum = CreverseTopNum; TreverseBottomNum = CreverseBottomNum; // Set trial strain Tstrain = strain; // Determine change in strain from last converged state double dStrain = Tstrain - Cstrain; // Calculate the trial state given the trial strain // if (fabs(dStrain) > DBL_EPSILON) if (fabs(dStrain) > 1e-10) determineTrialState (dStrain); ttStrain = dStrain; return 0; }
int Bond_SP01::setTrial (double strain, double &stress, double &tangent, double strainRate) { // Reset history variables to last converged state TRSlip = CRSlip; // Return slip TRLoad = CRLoad; // Return load TRSlope = CRSlope; // Return slope TmaxHSlip = CmaxHSlip; // Maximum slip in tension TminHSlip = CminHSlip; // Maximum slip in compression Tloading = Cloading; // Loading flag TYieldFlag = CYieldFlag; // Yield flag // Set trial Strain Tslip = strain; // Determine change in Strain from last converged state double dslip; dslip = Tslip - Cslip; // Calculate the trial state given the trial Strain determineTrialState (Tslip, dslip); stress = Tload; tangent = Ttangent; return 0; }
int TendonL01::setTrial (double strain, double &stress, double &tangent, double strainRate) { // Reset history variables to last converged state TminStrain = CminStrain; TmaxStrain = CmaxStrain; TloadingState = CloadingState; TloopPathState = CloopPathState; for (int i=0; i<SIZE1; i++) { TreverseTopStrain[i] = CreverseTopStrain[i]; TreverseTopStress[i] = CreverseTopStress[i]; TreverseBottomStrain[i] = CreverseBottomStrain[i]; TreverseBottomStress[i] = CreverseBottomStress[i]; } TreverseTopNum = CreverseTopNum; TreverseBottomNum = CreverseBottomNum; // Set trial strain Tstrain = strain + epsp; // Determine change in strain from last converged state double dStrain; dStrain = Tstrain - Cstrain; // Calculate the trial state given the trial strain // if (fabs(dStrain) > DBL_EPSILON) determineTrialState (dStrain); stress = Tstress; tangent = Ttangent; return 0; }
int Steel01::setTrialStrain (double strain, double strainRate) { // Reset history variables to last converged state TminStrain = CminStrain; TmaxStrain = CmaxStrain; TshiftP = CshiftP; TshiftN = CshiftN; Tloading = Cloading; Tstrain = Cstrain; Tstress = Cstress; Ttangent = Ctangent; // Determine change in strain from last converged state double dStrain = strain - Cstrain; if (fabs(dStrain) > DBL_EPSILON) { // Set trial strain Tstrain = strain; // Calculate the trial state given the trial strain determineTrialState (dStrain); } return 0; }
int SteelMP::setTrialStrain (double strain, double strainRate) { // Reset history variables to last converged state TminStrain = CminStrain; TmaxStrain = CmaxStrain; Tloading = Cloading; TYieldStrain = CYieldStrain; TYieldStress = CYieldStress; TReverStrain = CReverStrain; TReverStress = CReverStress; TPlasticExcursion = CPlasticExcursion; // Set trial strain Tstrain = strain; // Determine change in strain from last converged state double dStrain = Tstrain - Cstrain; // Calculate the trial state given the trial strain if (fabs(dStrain) > DBL_EPSILON) determineTrialState (dStrain); return 0; }
int DoddRestr::setTrial (double strain, double &stress, double &tangent, double strainRate) { // Reset history variables to last converged state Te_so = Ce_so; Tf_so = Cf_so; Tyield1 = Cyield1; Tregion = Cregion; Tpoint11= Cpoint11; Tpoint12= Cpoint12; Tpoint13= Cpoint13; Tpoint21= Cpoint21; Tpoint22= Cpoint22; Tpoint23= Cpoint23; Tpoint31= Cpoint31; Tpoint32= Cpoint32; Tpoint33= Cpoint33; Tpoint41= Cpoint41; Tpoint42= Cpoint42; Tpoint43= Cpoint43; Tpoint51= Cpoint51; Tpoint52= Cpoint52; Tpoint53= Cpoint53; Tep_o1 = Cep_o1; Tep_o2 = Cep_o2; Tep_M = Cep_M; Tfps_so = Cfps_so; Thist11 = Chist11; Thist12 = Chist12; Tpoint61= Cpoint61; Tpoint62= Cpoint62; Tpoint63= Cpoint63; Tsim1 = Csim1; TDam = CDam; Tstrain = Cstrain; Tstress = Cstress; Ttangent = Ctangent; // Determine change in strain from last converged state double dStrain = strain - Cstrain; // if (fabs(dStrain) > DBL_EPSILON) { // Set trial strain Tstrain = strain; // Calculate the trial state given the trial strain determineTrialState (dStrain); // } stress = Tstress; tangent = Ttangent; return 0; }