####################################################
# 
# --- Command log file ---
# PLAXIS 3D 2025.1
# Tutorial - 03
# Loading of a suction pile
# Note: Use the PLAXIS Commands runner to execute this
# 
####################################################
_setproperties "Title" "Loading of a suction pile"
_initializerectangular SoilContour -30 0 30 30
_borehole 0 0
_soillayer 0
_setsoillayerlevel Borehole_1 1 -30
_set Borehole_1.Head 1
_soilmat "Identification" "Clay" "SoilModel" "Mohr-Coulomb" "DrainageType" "Undrained B" "Colour" 15262369 "gammaUnsat" 20 "gammaSat" 20 "ERef" 1000  "nu" 0.35 "EInc" 1000 "sURef" 1 "sUInc" 4 "TensionCutOff" False "K0PrimaryIsK0Secondary" True "K0Determination" "Manual" "K0Primary" 0.5  
_soilmat "Identification" "Interface" "SoilModel" "Mohr-Coulomb" "DrainageType" "Undrained B"  "Colour" 10676870 "gammaUnsat" 20 "gammaSat" 20 "ERef" 1000 "nu" 0.35 "EInc" 1000  "sURef" 1 "sUInc" 4  "K0Determination" "Manual" "TensionCutOff" False "K0PrimaryIsK0Secondary" True "K0Primary" 0.5  
_set Soillayer_1.Soil.Material Clay
_gotostructures
_polycurve 2.5 0 0 1 0 0 0 1 0
_add Polycurve_1
_set Polycurve_1.Segments[0].SegmentType "Arc"
_set Polycurve_1.Segments[0].ArcProperties.RelativeStartAngle1 90
_set Polycurve_1.Segments[0].ArcProperties.Radius 2.5
_set Polycurve_1.Segments[0].ArcProperties.CentralAngle 180
_extrude (Polycurve_1) 0 0 -10
_posinterface Surface_1
_neginterface Surface_1
_close Polycurve_1
_surface Polycurve_1
_neginterface Surface_2
_set (Surface_2.NegativeInterface.MaterialMode Surface_1.NegativeInterface.MaterialMode Surface_1.PositiveInterface.MaterialMode) "Custom"
_set (Surface_2.NegativeInterface.Material Surface_1.NegativeInterface.Material Surface_1.PositiveInterface.Material) Interface
_rigidbody (Surface_2 Surface_1)
_set (Surface_1.RigidBody.Xref Surface_2.RigidBody.Xref) 2.5
_set (Surface_1.RigidBody.Yref Surface_2.RigidBody.Yref) 0
_set (Surface_1.RigidBody.Zref Surface_2.RigidBody.Zref) -7
_set (Surface_1.RigidBody.TranslationConditiony Surface_2.RigidBody.TranslationConditiony) "Displacement"
_set (Surface_1.RigidBody.RotationConditionx Surface_2.RigidBody.RotationConditionx) "Rotation"
_set (Surface_1.RigidBody.RotationConditionz Surface_2.RigidBody.RotationConditionz) "Rotation"
_polycurve 7.5 0 0 1 0 0 0 1 0
_add Polycurve_2
_set Polycurve_2.Segments[0].SegmentType "Arc"
_set Polycurve_2.Segments[0].ArcProperties.RelativeStartAngle1 90
_set Polycurve_2.Segments[0].ArcProperties.Radius 7.5
_set Polycurve_2.Segments[0].ArcProperties.CentralAngle 180
_close Polycurve_2
_extrude (Polycurve_2) 0 0 -15
_delete Geometry (Polycurve_1 Polycurve_2)
_gotomesh
_set (Surface_1_1.CoarsenessFactor Surface_3_2.CoarsenessFactor Surface_3_1.CoarsenessFactor Surface_2_1.CoarsenessFactor) 0.25
_mesh 0.05 256 True 1.2 0.005
_gotostages
_phase InitialPhase
_set Model.CurrentPhase Phase_1
_set Phase_1.Identification "Install pile"
_set Phase_1.Deform.IgnoreUndrainedBehaviour True
_activate RigidBodies Phase_1
_activate Interfaces Phase_1
_phase Phase_1
_set Model.CurrentPhase Phase_2
_set Phase_2.Identification "Load pile 30 degrees"
_set Phase_2.Deform.ResetDisplacementsToZero True
_set Phase_2.Solver "Pardiso (multicore direct)"
_set Phase_2.Deform.UseDefaultIterationParams False
_set Phase_2.Deform.MaxLoadFractionPerStep 0.1
_set (RigidBody_1_1.Fx RigidBody_2_1.Fx) Phase_2 1949
_set (RigidBody_1_1.Fz RigidBody_2_1.Fz) Phase_2 1125
_set Model.CurrentPhase Phase_1
_phase Phase_1
_set Model.CurrentPhase Phase_3
_set Phase_3.Identification "Load pile 40 degrees"
_set Phase_3.Deform.ResetDisplacementsToZero True
_set Phase_3.Solver "Pardiso (multicore direct)"
_set Phase_3.Deform.UseDefaultIterationParams False
_set Phase_3.Deform.MaxLoadFractionPerStep 0.1
_set (RigidBody_1_1.Fx RigidBody_2_1.Fx) Phase_3 1724
_set (RigidBody_1_1.Fz RigidBody_2_1.Fz) Phase_3 1447
_set Model.CurrentPhase Phase_1
_phase Phase_1
_set Model.CurrentPhase Phase_4
_set Phase_4.Identification "Load pile 50 degrees"
_set Phase_4.Deform.ResetDisplacementsToZero True
_set Phase_4.Solver "Pardiso (multicore direct)"
_set Phase_4.Deform.UseDefaultIterationParams False
_set Phase_4.Deform.MaxLoadFractionPerStep 0.1
_set (RigidBody_1_1.Fx RigidBody_2_1.Fx) Phase_4 1447
_set (RigidBody_1_1.Fz RigidBody_2_1.Fz) Phase_4 1724
_set Model.CurrentPhase Phase_1
_phase Phase_1
_set Model.CurrentPhase Phase_5
_set Phase_5.Identification "Load pile 60 degrees"
_set Phase_5.Deform.ResetDisplacementsToZero True
_set Phase_5.Solver "Pardiso (multicore direct)"
_set Phase_5.Deform.UseDefaultIterationParams False
_set Phase_5.Deform.MaxLoadFractionPerStep 0.1
_set (RigidBody_1_1.Fx RigidBody_2_1.Fx) Phase_5 1125
_set (RigidBody_1_1.Fz RigidBody_2_1.Fz) Phase_5 1949