PP constraint

[nbastias]

Dear STKO team,

I am performing a dynamic analysis of a slope considering pore pressure generation through SSPquadUP elements. Given the geometry and subsequent meshing, it is not possible to define the PP constraint under the surface. Attached is a figure of the geometry and the approximate location where the PP constraint should be (I also attach .txt with the coordinates of the constraint). Is there any way to address this issue without having to define subsections in the geometry?, this because triangular elements would be generated. I also attach the .scd, thank you in advance.

Best regards

[STKO Team]

If you don't want to cut the geometry you cannot select that sub-surface to assign the condition.
In this case you can insert a customCommand in which you can write a Tcl script to loop over all nodes and based on their spatial location you can decide whether to assign the condition or not.

[nbastias]

I understand, can you please help me a little with the general structure of this customCommand?

[STKO Team]

A little hint:

Code: Select all

foreach ntag [getNodeTags] {
   set npos [nodeCoord $ntag]
   set nx [expr double([lindex $npos 0])]
   set ny [expr double([lindex $npos 1])]
   set nz [expr double([lindex $npos 2])]
   puts "Note $ntag @ ($nx $ny $nz)"
   # here you can put an if-condition checking the coordinates to assign an sp 
}