Modelling reconstituted clay sample

[TrishReddy]

Hello,
I am new to the Openseespy community I started learning openseespy by modeling the reconstitution of clay. I implemented elasticIsotropic material ('AxiSymmetric'). I choose quad element. But when I am unable to create quad element, I am getting the error as Dead kernel. I am attaching the code below. Kindly, correct my mistakes.

Code: Select all

from openseespy.opensees import *
import openseespy.postprocessing.Get_Rendering as opsplt
import numpy as np
import matplotlib.pyplot as plt
import math
wipe()
nDMaterial('ElasticIsotropic',1,3058,0.453,0)
model('basic','-ndm',2,'-ndf',2)
N = 50
n = 300
L = 25
l = 150
o = 1
for i in range(0,N+1):
    for j in range(0,n+1):
        o = o+1
        node(o,i*L/N,j*l/n,0)
        
 for a in range(2,15052+301,301):
    fix(a,1,1,0)
for b in range(3,303,1):
    fix(b,1,0,0)
for c in range(15053,15352,1):
    fix(c,1,0,0)
    
# element  definition
p =1
for w in range(2,15052+1,1):
    p = p+1
    element('quad',p,w,w+301,w+302,w+1,1,1,1,0,0,0)

I also have another question after I am successful in the model creation. I need to impose a vertical load on the top and I have 8-time steps, in which the load is kept constant during the time step and increased for the next time step. Can I know what the commands I need to use to achieve this are?

Thank you for the support.

[mhscott]

There's a few things going on here. You need to give a type string as input:
https://openseespydoc.readthedocs.io/en ... /quad.html
Run the script from Python command line instead of in notebook and you'll see the error.
Also, you fix three DOFs in your fix commands, but there's only two DOFs/node in your model.

[mhscott]

Start with a simple patch test like this: https://portwooddigital.com/2021/04/04/ ... atch-test/