2D cylinder -Force control but not linear with time

[badinhle]

Dear all,

I am a newbie using OpenSees. My target is a static 2D model for a concrete cylinder 100x200 mm. It is tested under force control, but not linear with time (as below).
I try many times using " timeSeries Path " to access Time and Force files. But the analysis can not run, or run but wrong with weird results.
- Some Errors under my efforts:
Error Recorder Node/Envolope #TimeSeries must equal # dof - Ignoring TimeSeries option
Warning PathTimeSeries :: PathTimeSeries <> - file containing data points for path and time do not contain some number of points
My questions are:
1) Is it possible for OpenSees to build my model with nonlinear Force-time relationship?
2) Due to my incomplete understanding about how OpenSees works, please tell me what I did wrong with my code? How to improve it?

Thank you very much for your reading.
Cheers,



# Times.txt (s):
0
1
2
# Forces.txt (factor):
0
0.4
1

#
# CYLINDER - 2D WITH FIBER SECTION
# Unit: kN, mm

# Remove existing model
wipe

# ------------------------------
# Definition of model
# ------------------------------
model basic -ndm 2 -ndf 3
file mkdir DinhData;

# ------------------------------
# Definition of node
# ------------------------------
node 1 0 0
node 2 0 200

# ------------------------------
# Definition of boundary
# ------------------------------
fix 1 1 1 1
fix 2 1 0 1

# ------------------------------
# Define materials
# ------------------------------
set IDconc 1; # material ID tag -- confined core concrete
uniaxialMaterial Concrete02 $IDconc -0.082 -0.0027 -0.081 -0.0029 0.5 0.006 30

# ------------------------------
# Definition of elements
# ------------------------------

set SecTag 2; # set tag for symmetric section

section fiberSec $SecTag {
#patch circ $matTag $numSubdivCirc $numSubdivRad $yCenter $zCenter $intRad $extRad $startAng $endAng
patch circ 1 16 16 0 0 0 50 0 360
}

set TransfTag 1; # associate a tag to column transformation
set TransfType Linear; # options, Linear PDelta Corotational
geomTransf $TransfType $TransfTag;

set numIntgrPts 5;
element nonlinearBeamColumn 1 1 2 $numIntgrPts $SecTag $TransfTag;

# ------------------------------
# Define loads
# ------------------------------

set tsTag 1
#timeSeries Path $tsTag -dt 1 -filePath Forces.txt
timeSeries Path $tsTag -fileTime Times.txt -filePath Forces.txt


# Create a Plain load pattern with a linear TimeSeries: command pattern Plain $tag $timeSeriesTag { $loads }
pattern Plain 1 $tsTag {
# Create the nodal load - command: load nodeID xForce yForce
load 2 0 -648 0
}


# ------------------------------
# Start of analysis generation
# ------------------------------
# Create the system of equation, a SPD using a band storage scheme
system BandSPD

# Create the DOF numberer, the reverse Cuthill-McKee algorithm
numberer RCM

# Create the constraint handler, a Plain handler is used as homo constraints
constraints Plain

# Create the integration scheme, the LoadControl scheme using steps of 1.0
integrator LoadControl 1

# Create the solution algorithm, a Linear algorithm is created
algorithm Linear

# create the analysis object
analysis Static



# ------------------------------
# Start of recorder generation
# ------------------------------
# create a Recorder object for the nodal displacements at node 2
recorder Node -file DinhData/FreeNode.out -timeSeries $tsTag -node 2 -dof 2 disp

# ------------------------------
# Finally perform the analysis
# ------------------------------
# Perform the analysis
analyze 3


# Print the current state at node 4 and at all elements
puts "node 2 displacement: [nodeDisp 2]"
print node 2
print ele

#
puts "Analysis is done"
wipe analysis

[badinhle]

Good morning all,

May I ask your help, please?