bug

For posts concerning the documentation, errors, ommissions, general comments, etc.

Moderators: silvia, selimgunay, Moderators

Post Reply
civilmanjun
Posts: 32
Joined: Mon Nov 19, 2007 12:47 am
Location: China

bug

Post by civilmanjun » Thu Dec 13, 2007 2:34 am

Q1: I want define axial force of truss as performance function, as follows

performanceFunction 1 "671000+{rec_element_1_localForce_1}"

limit-state function number: 1
child killed: segmentation violation

Q2:I'm in trouble with following problem, please help me, thanks.
I loaded an example3b.tcl from http://www.inrisk.ubc.ca/software1.htm.
but when I run example3b.tcl in opensees version1.7.4 and 1.7.5 ,the following message occur:

WARNING could not create uniaxialMaterial SmoothSteel01
while executing
"uniaxialMaterial SmoothSteel01 1 $sigmaY $E $b $smoothgama $smootheta"
<file "exampke3b.tcl" line 19>


#example3b.tcl
set load 200000.0
set length 2000.0
set E 200000.0
set sigmaY 400.0
set A 400.0
set b 0.05

set smoothgamma 0.9
set smootheta 4.0

##### Run reference analysis and compute DDM gradients #####
model basic -ndm 2 -ndf 2
reliability

node 1 0.0 0.0
node 2 $length 0.0

#uniaxialMaterial Steel01 1 $sigmaY $E $b
uniaxialMaterial SmoothSteel01 1 $sigmaY $E $b $smoothgamma $smootheta

element truss 1 1 2 $A 1

fix 1 1 1
fix 2 0 1

pattern Plain 1 Linear {
load 2 $load 0.0
}

randomVariablePositioner 1 -createRV3 lognormal $sigmaY [expr 0.1*$sigmaY] -element 1 -material sigmaY
randomVariablePositioner 2 -createRV3 lognormal $load [expr 0.1*$load] -loadPattern 1 -loadAtNode 2 1


constraints Plain
numberer RCM
test NormUnbalance 1.0e-5 25 0
integrator LoadControl 0.01 1 0.01 0.01
algorithm Newton
system BandGeneral
analysis Static

recorder Node response.out disp -time -node 2 -dof 1

performanceFunction 1 "4.2 - {u_2_1}"

# CREATE NECESSARY RELIABILITY ANALYSIS TOOLS
probabilityTransformation Nataf -print 0
reliabilityConvergenceCheck Standard -e1 1.0e-3 -e2 1.0e-3 -print 1
gFunEvaluator OpenSees -analyze 100
gradGEvaluator FiniteDifference -pert 1000
searchDirection iHLRF
meritFunctionCheck AdkZhang -multi 2.0 -add 10.0 -factor 0.5
stepSizeRule Armijo -maxNum 10 -base 0.5 -initial 1.0 2 -print 0
startPoint Mean
findDesignPoint StepSearch -maxNumIter 100 -printAllPointsY allPointsY.out
rootFinding Secant -maxIter 50 -tol 1.0e-5 -maxStepLength 1.0

runFORMAnalysis Example3b.out






#auxiliary file
#aux_3.tcl
integrator LoadControl 0.05 3 0.05 0.05
analyze 16

integrator LoadControl 0.02 3 0.02 0.02
analyze 10
learn more

silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia » Thu Dec 13, 2007 6:21 am

please contact the authors of that example.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104

civilmanjun
Posts: 32
Joined: Mon Nov 19, 2007 12:47 am
Location: China

Post by civilmanjun » Thu Dec 13, 2007 4:29 pm

ok, thanks, silvia
I'll contact author of that example.

but for Q1, I think it's a bug.

Q1: I want define axial force of truss as performance function, as follows

performanceFunction 1 "671000+{rec_element_1_localForce_1}"

limit-state function number: 1
child killed: segmentation violation
learn more

rjaeger
Posts: 102
Joined: Thu Aug 31, 2006 9:57 pm
Location: UC Davis

Post by rjaeger » Thu Dec 13, 2007 7:53 pm

If I remember correctly, Frank said he wanted everyone to email him with segmentation faults (from OpenSees days 2007).

civilmanjun
Posts: 32
Joined: Mon Nov 19, 2007 12:47 am
Location: China

Post by civilmanjun » Thu Dec 13, 2007 10:20 pm

rjaeger, thanks for your advice, but I don't know how to contact with Frank, Do you have his/her email? thank you very much.
learn more

silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia » Fri Dec 14, 2007 6:19 am

opensees AT berkeley DOT edu
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104

Post Reply