I believe that reliability analysis for version 2.5.0 has a bug. I have a code which works with 2.4.6 correctly but with current version of 2.5.0 gives me zero for displacement sensitivity with respect to E.
The code is as follows (simple cantilever beam with elastic section):
#------------------------------------------------------------
#simple model
#Units: N & m
set ee 200000000000.0;
set II 1e-05;
model BasicBuilder -ndm 2 -ndf 3;
reliability;
node 1 0.0 0.0;
node 2 1.0 0.0;
fix 1 1 1 1;
section Elastic 1 $ee 0.012 $II;
geomTransf Linear 1;
#element dispBeamColumnWithSensitivity $eleTag $iNode $jNode $numIntgrPts $secTag $transfTag <integration method>
element dispBeamColumnWithSensitivity 2 1 2 5 1 1;
parameter 1 -element 2 E
timeSeries Linear 1;
pattern Plain 1 1 {
load 2 0.0 1.0 0.0;
}
recorder Node -file enddisp.txt -node 2 -dof 2 disp;
recorder Node -file ddmdisp.txt -node 2 -dof 2 "sensitivity 1";
system BandGeneral;
test NormDispIncr 1.0e-12 50 1;
constraints Transformation;
integrator LoadControl 1;
algorithm Newton;
numberer RCM;
sensitivityIntegrator -static
sensitivityAlgorithm -computeAtEachStep
analysis Static
analyze 1;
wipe;
wipeReliability;
#------------------------------------------------------------
Reliability bug
Moderators: silvia, selimgunay, Moderators
-
- Posts: 16
- Joined: Thu Jun 08, 2006 3:17 pm
- Location: University of California at San Diego
- Contact:
-
- Posts: 16
- Joined: Thu Jun 08, 2006 3:17 pm
- Location: University of California at San Diego
- Contact:
Re: Reliability bug
I am wondering if OpenSees confirms this bug? Or maybe using new version I missed something?