Error on new uniaxialmaterial

[aylsworth]

Hi,

Need help on my new uniaxialmaterial. I already compiled and obtained a dll file. However, when I loaded it in Opensees, it says:

WARNING: numeric analysis returns 1 -- Umfpackgenlinsolver::solve
WARNING NewtonRaphson::solveCurrentStep() - the LinearSysofEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at load factor 0.01
Domain::update - domain failed in update
Opensees > analyze failed, returned: -3 error flag -3


The code I used was the Mattest.tcl which is included in the command language manual. I am trying to model a material with degrading bilinear hysteretic rule.

Thanks for your help!

Patrick



# matTest.tcl: SDOF truss to test uniaxial material models
# Units: kip, in
# MHS, Sept 1999
# email: mhscott@ce.berkeley.edu

wipe
model BasicBuilder -ndm 1 -ndf 1
# Define nodes
node 1 0.0
node 2 1.0
# Fix node 1
fix 1 1

# Define uniaxialMaterial
#uniaxialMaterial DegradingBilinear mattag Ky Dy Fy Ku B0 dmaxT dmaxC
uniaxialMaterial DegradingBilinear 1 [expr 3.50/0.002] 0.006 3.50 [expr (5.0-3.5)/(0.01-0.006)] 0.4 0.10 -0.10;

#uniaxialMaterial ElasticPPcpp 1 3000.0 0.03;

# Define truss element with unit area
# tag ndI ndJ A matTag
element truss 3 1 2 1.0 1
set dt 0.010 ;# Increment between data points
set filename pattern1.txt ;# Filename containing data points
set factor 0.006 ;# Factor applied to data values
# Read displacement pattern from file
# Note, any pattern type can be used here: Linear, Path, Sine, etc.
pattern Plain 1 "Series -dt $dt -filePath $filename -factor $factor" {
# Set reference displacement value
# node dof value
sp 2 1 1.0
}

# Impose monotonic displacements
#pattern Plain 2 "Linear -factor $factor" {
# sp 2 1 1.0
#}
# Record nodal displacements (same as strains since truss length is 1.0)
#recorder Node -file truss.out disp -load -node 2 -dof 1
recorder Node -file truss.txt -time -node 2 -dof 1 disp;
# Record truss force (same as stress since truss area is 1.0)
recorder Element -file force.txt -time -ele 3 localForce
system UmfPack
constraints Penalty 1.0e12 1.0e12
# Set increment in load factor used for integration
# Does not have to be the same as dt used to read in displacement pattern
set dl $dt
integrator LoadControl $dl 1 $dl $dl
test NormDispIncr 1.0e-6 10
algorithm Newton
numberer RCM
analysis Static
analyze 10000

[ismailqeshta]

I think the load factor and number of steps need to be corrected. The first step to identify the problem is to make sure that their multiplication is always equal to 1. This can, at least, ensure that you obtain the targeted total applied load.

[aylsworth]

ismailqeshta wrote:
> I think the load factor and number of steps need to be corrected. The first
> step to identify the problem is to make sure that their multiplication is
> always equal to 1. This can, at least, ensure that you obtain the targeted
> total applied load.


Thanks for your reply! Apparently the error was on my code. Thank you!

Patrick