hey everybody,
here i have modelled a cantilever beam using 3d elements as per opensees example 8 but still it says it is failing kindly check my code and let me know my mistake please
# ----------------------------
# Start of model generation
# ----------------------------
# Create ModelBuilder with 3 dimensions and 6 DOF/node
model basic -ndm 3 -ndf 3
# create the material
nDMaterial ElasticIsotropic 1 100 0.25 1.27
# Define geometry
# ---------------
# define some parameters
set eleArgs "1"
set element stdBrick
#set element BbarBrick
set nx 6
set ny 2
set nz 2
set nn [expr ($nz+1)*($nx+1)*($ny+1) ]
# mesh generation
block3D $nx $ny $nz 1 1 $element $eleArgs {
1 0 0 0
2 10 0 0
3 10 0 2
4 0 0 2
5 0 2 0
6 10 2 0
7 10 2 2
8 0 2 2
}
recorder Node -file Node.out -time -node 7 -dof 2 disp
recorder plot Node.out CenterNodeDisp 625 10 625 450 -columns 1 2
recorder display ShakingBeam 0 0 300 300 -wipe
prp -100 100 120.5
vup 0 1 0
display 1 0 1
# --------------------------
set load -0.1
# Constant point load
pattern Plain 1 Linear {
load $nn 0.0 $load 0.0
}
# boundary conditions
fixX 0.0 1 1 1
# --------------------------------------------------------------------
# Start of static analysis (creation of the analysis & analysis itself)
# --------------------------------------------------------------------
# Load control with variable load steps
# init Jd min max
integrator LoadControl 1.0 1 1.0 10.0
# Convergence test
# tolerance maxIter displayCode
test NormUnbalance 1.0e-10 20 1
# Solution algorithm
algorithm Newton
# DOF numberer
numberer RCM
# Cosntraint handler
constraints Plain
# System of equations solver
system ProfileSPD
# Analysis for gravity load
analysis Static
# Perform the analysis
analyze 5
# --------------------------
# End of static analysis
# --------------------------
# ----------------------------
cantilever modelled with 3d elements
Moderators: silvia, selimgunay, Moderators
-
- Posts: 10
- Joined: Wed Feb 03, 2016 5:34 am
- Location: CEPT UNIVERSITY
Re: cantilever modelled with 3d elements
the block command nodes are labelled incorrectly and as a result generate a negative Jacobian
try:
block3D $nx $ny $nz 1 1 $element $eleArgs {
1 0 0 0
2 10 0 0
3 10 2 0
4 0 2 0
5 0 0 2
6 10 0 2
7 10 2 2
8 0 2 2
}
instead
try:
block3D $nx $ny $nz 1 1 $element $eleArgs {
1 0 0 0
2 10 0 0
3 10 2 0
4 0 2 0
5 0 0 2
6 10 0 2
7 10 2 2
8 0 2 2
}
instead
-
- Posts: 10
- Joined: Wed Feb 03, 2016 5:34 am
- Location: CEPT UNIVERSITY
Re: cantilever modelled with 3d elements
perfect thankyou very much