The OpenSees Community

Hello
I am a new user of Openseepy. I try to create 10*20m a soil model (10m in vertical direction and 20 m in the horizontal direction). I try to use 9_4_QuadUP with an element size of 0.5m by using For loop in Python. Nevertheless, it is too complicated for me. There is an example called "freeFieldEffective.py" that works fine if there is one element in the horizontal direction. Is there any reference or code that can help me create my model?

Thank you in advance

Thank you Prof. Scott.

Hello

I am using Openseespy and I have been trying to visualize my model but it seems the elements that I have been created will not appear in the outcome of my model. Then, I decided to create a much simpler model to see the result, it did not work in simpler model as well. As it is evident in the code below in the outcome window nodes are shown but not the element. I would really appreciate if someone helped me.

import openseespymac.opensees as op
import openseespy.postprocessing.Get_Rendering as oplt


op.wipe()

nNodeY = 3
# create ModelBuilder (with two-dimensions and 3 DOF/node)
op.model("BasicBuilder", "-ndm",2, "-ndf",3)

# create nodes & add to Domain - command: node nodeId xCrd yCrd
op.node(1, 0.0, 0.0)
op.node(2, 4.0, 0.0)
op.node(3, 4.0, 4.0)
op.node(4, 0.0, 4.0)

op.fix(1, 0, 1, 0)
op.fix(3, 0, 1, 0)

for i in range(1, (nNodeY - 2), 6):
op.equalDOF(i, i + 2, 1, 2)

op.model('basic', '-ndm', 2, '-ndf', 2)
op.node(5, 2.0, 0.0)
op.node(6, 4.0, 2.0)
op.node(7, 2.0, 4.0)
op.node(8, 0.0, 2.0)
op.node(9, 2.0, 2.0)

op.nDMaterial('PressureDependMultiYield', 1, 2, 1.4, 7.5e4, 2.0e5, 33, 0.1, 80.0, 0.5, 27, 0.07, 0.4, 2, 10,
0.01, 1, 20, 0.7, 0.9, 0.02, 0.7, 101)

op.element('9_4_QuadUP', 1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 1, 5.24E6, 1.0, 1.0, 1.0)

oplt.plot_model("Nodes", "elements")