Skip to content

Issue with IntialStateAnalysis command #153

Description

@bahugunaashish

I was trying to use the _IntialStateAnalysis for my model in Openseespy. To check the applicability I tried to run the example given at https://opensees.berkeley.edu/wiki/index.php/InitialStateAnalysisWrapper

the kernel keeps getting crashed with Windows fatal exception: code 0xc0000374.
following is the conversion from tcl to py
"
import openseespy.opensees as op
op.model ('BasicBuilder', '-ndm', 3, '-ndf', 3)

##create the nodes
op.node (1, 1.0,0.0, 0.0)
op.node (2, 1.0,1.0, 0.0)
op.node (3, 0.0,1.0, 0.0)
op.node (4, 0.0,0.0, 0.0)
op.node (5, 1.0,0.0, 1.0)
op.node (6, 1.0,1.0, 1.0)
op.node (7, 0.0,1.0, 1.0)
op.node (8, 0.0,0.0, 1.0)

##boundary conditions
op.fix (1, 1, 1, 1)
op.fix (2, 1, 1, 1)
op.fix (3, 1, 1, 1)
op.fix (4, 1, 1, 1)
op.fix (5, 1, 1, 0)
op.fix (6, 1, 1, 0)
op.fix (7, 1, 1, 0)
op.fix (8, 1, 1, 0)

##define main material object
op.nDMaterial ('DruckerPrager', 1, 2777.8, 925, 5., 0.398, 0.398, 0, 0, 0, 0, 0, 1, 1.7)
##define material wrapper
op.nDMaterial ('InitialStateAnalysisWrapper', 2, 1, 3)

##create the element (NOTE: the material tag associated with this element is that of the wrapper)
op.element ('SSPbrick', 1, 1, 2, 3, 4, 5, 6, 7, 8, 2, 0.0, 0.0, -17.0)

##create the pre-gravity recorders
step = 0.1

op.recorder ('Node', '-time', '-file', 'Gdisp.out', '-dT', step, '-nodeRange', 5, 8, '-dof', 1, 2, 3, 'disp')
op.recorder ('Element', '-ele', 1, '-time', '-file', 'Gstress.out', '-dT', step, 'stress')
op.recorder ('Element', '-ele', 1, '-time', '-file', 'Gstrain.out', '-dT', step, 'strain')

##create the gravity analysis
op.integrator ('LoadControl', 0.1)
op.numberer ('RCM')
op.system ('SparseGeneral')
op.constraints ('Transformation')
op.test ('NormDispIncr', 1e-5, 40, 1)
op.algorithm ('Newton')
op.analysis ('Static')

##turn on the initial state analysis feature
op.InitialStateAnalysis ('on')

##analyze four steps
op.analyze (4)

##turn off the initial state analysis feature
op.InitialStateAnalysis ('off')

##create post-gravity recorders
op.recorder ('Node', '-time', '-file', 'disp.out', '-dT', step, '-nodeRange', 5, 8, '-dof', 1, 2, 3, 'disp')
op.recorder ('Element', '-ele', 1, '-time', '-file', 'stress.out', '-dT', step, 'stress')
op.recorder ('Element', '-ele', 1, '-time', '-file', 'strain.out', '-dT', step, 'strain')

##analyze for three steps, should have non-zero stress and strain with zero displacement
op.analyze (3)

op.wipe()

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions