Hi,kamyarsa wrote: ↑Thu Mar 09, 2023 1:13 pm Hi, I want to use PDMY02 material. I have one problem regarding using it. I tried to follow the instruction to write the code for the material, I just wonder why changing the void ratio(e) does not change any results.
I even model a single 2d element and simulate a triaxial test with this material, still changing void ratio does not change anything!!!
Can someone explain why this happens? Am I doing something wrong?
# nDMaterial('PressureDependMultiYield02', matTag, nd, rho, refShearModul, refBulkModul, frictionAng, peakShearStra, refPress, pressDependCoe, PTAng, contrac[0], contrac[2], dilat[0], dilat[2], noYieldSurf=20.0, *yieldSurf=[], contrac[1]=5.0, dilat[1]=3.0, *liquefac=[1.0,0.0],e=0.6, *params=[0.9, 0.02, 0.7, 101.0], c=0.1)
ops.nDMaterial('PressureDependMultiYield02', matTag, 2 , rho, refShearModul, refBulkModul, frictionAng, peakShearStra, refPress, pressDependCoe, PTAng, 0.067 , 0.23 , 0.06 , 0.27, 20, 5.0, 3.0, *[1.0,0.0], 0.6, *[0.9, 0.02, 0.7, 101.0], 0.1 )
did you use Standard Single Point (SSP) fully coupled (UP) elements? (for ex. SSPquadUP and SSPBrickUP)
As far as I know, these elements also have a void ratio. If you forget to update e for these elements, OpenSees probably overwrite the void ratio because elements are defined after materials.