Opensees 2.4.3 + shellMITC4: bug?

[alz]

I have a problem with a very simple model: one shell element. Gravity Analysis only.
I do the analysis and it works, but when I type "exit" to save the recorder and quit OpenSees shows this error message:

"OpenSees(669,0x7fff7bcb5310) malloc: *** error for object 0x7feb4a616930: pointer being freed was not allocated
*** set a breakpoint in malloc_error_break to debug
Abort trap: 6"

Here all the screen shot. Same problem both with Mac Os and Windows 8 (64 bit version). The problem is only for the multi-layer shell element.
What could I do?

__________________________________________________________________________________

OpenSees -- Open System For Earthquake Engineering Simulation
Pacific Earthquake Engineering Research Center -- 2.4.3 (rev 5665)

(c) Copyright 1999-2013 The Regents of the University of California
All Rights Reserved
(Copyright and Disclaimer @ http://www.berkeley.edu/OpenSees/copyright.html)


OpenSees > cd /Users/Alz/Downloads
OpenSees > source modnew.tcl
Modello OK
Gravity Analysis ok...
OpenSees > exit
OpenSees(780,0x7fff7bcb5310) malloc: *** error for object 0x7ff280f05480: pointer being freed was not allocated
*** set a breakpoint in malloc_error_break to debug
Abort trap: 6
logout

[mhscott]

Would you mind e-mailing me your script, or posting it here?

[alz]

You can use the example in this page ( http://www.luxinzheng.net/download/Open ... enSEES.htm ) with OpenSEES 2.4.3

Otherwise this model is simpler.
Thank you!!

____________

# units: N, mm, sec
wipe;
logFile log.txt
file mkdir Muro;
model basic -ndm 3 -ndf 6


nDMaterial PlaneStressUserMaterial 10 40 7 30 1 -5 -0.0021 -0.006 0.00004 0.1
nDMaterial PlateFromPlaneStress 1 10 11700

uniaxialMaterial Steel01 20 630 200000 0.001

nDMaterial PlateRebar 9 20 90
nDMaterial PlateRebar 11 20 0


section LayeredShell 1 6 1 17.5 1 17.5 9 0.0755190541728316 11 0.0755190541728316 1 17.5 1 17.5



# Modello
# Nodi
node 1 0 0 0
node 2 200 0 0
node 17 0 0 200
node 18 200 0 200

#Element
element shellMITC4 1 1 2 18 17 1


# Incastri
fix 1 1 1 1 1 1 1
fix 2 1 1 1 1 1 1



recorder Node -file Muro/Reac.out -time -node 1 2 -dof 1 2 3 reaction
recorder Node -file Muro/Disp.out -time -node 18 -dof 1 2 3 disp


pattern Plain 1 Linear {
load 17 0 0 -11250 0 0 0
load 18 0 0 -11250 0 0 0


}

puts "Modello OK"




constraints Plain
numberer RCM
system BandGeneral
test NormDispIncr 1.0e-8 200 ;
#algorithm BFGS -count 100
algorithm Newton;
integrator LoadControl 0.1;
analysis Static
analyze 10;

puts "Gravity Analysis ok..."
loadConst -time 0.0;

[xiell]

I think it is caused by stress concentration, you can try to subdivide this element, and then apply the vertical load on several top nodes in the central zone.