hi, everyone, I assembly BeamColumnJoint Element to the 3-story frame.
But in this case you see no errors but I get a period above 6 seconds. The period of central Line model is just 0.5 seconds.
# 1601 1701 1801 1901 2001
# -->04 ||========== ==||============||============||============||
# || || || || || 3600
# || 1101 || 1201 || 1301 || 1401 || 1501
# -->03 ||========== ==||============||============||============||
# || || || || ||
# 602 || || || ||
# 603 601 || 701 || 801 || 901 || 1001
# -->02 ||============ ||============||============||============||
# 604 || || || ||
# || || || || || 3600 # 节点区编号参见6节点，逆时针编号
# || 101 || 201 || 301 || 401 || 501
# ==== ==== ==== ==== ====
# |<--- 5000-- ->| <--- 5000 --->|<--- 5000-- ->|<---5000 --->|
####################################程序开始#################################################
model BasicBuilder -ndm 2 -ndf 3
#=====================定义几何参数========================
set numbay 4; # 跨数
set numstory 3; # 层数
set baywidth 5000.0; # 跨度
set storyheight 3600.0; # 层高
#框架柱截面数据
set Bc 500; #Bc -- width of section 柱平面外宽度
set Hc 500; #Hc -- height of section 柱平面内宽度
set cov 30 ; #covc -- thickness of cover concreat
# 框架梁截面数据
set Bb 250
set Hb 500
# 节点几何数据 JointB
# -----------------------
# | |
# | |
# | - | JointH Hb=梁高
# | |
# | Hc=柱宽 |
# -----------------------
set JointB [expr $Hc]
set JointH [expr $Hb]
set JointD [expr $Bc]
set JointV [expr $JointB*$JointH*$JointD]
#======================创建节点===========================
# 定义几何中心点坐标
set X1 0
set X2 [expr $baywidth]
set X3 [expr $baywidth*2]
set X4 [expr $baywidth*3]
set X5 [expr $baywidth*4]
set Y1 0
set Y2 [expr $storyheight]
set Y3 [expr $storyheight*2]
set Y4 [expr $storyheight*3]
# 定义各层节点坐标，从第一层至第三层每个几何中线点周围分布四个节点，编号顺序为逆时针，为 X*100+1 至 X*100+4
# 标高+-0.000
node 101 $X1 $Y1
node 201 $X2 $Y1
node 301 $X3 $Y1
node 401 $X4 $Y1
node 501 $X5 $Y1
# 一层
node 601 [expr $X1+$JointB/2.0] [expr $Y2]
node 602 [expr $X1] [expr $Y2+$JointH/2.0]
node 603 [expr $X1-$JointB/2.0] [expr $Y2]
node 604 [expr $X1] [expr $Y2-$JointH/2.0]
node 701 [expr $X2+$JointB/2.0] [expr $Y2]
node 702 [expr $X2] [expr $Y2+$JointH/2.0]
node 703 [expr $X2-$JointB/2.0] [expr $Y2]
node 704 [expr $X2] [expr $Y2-$JointH/2.0]
node 801 [expr $X3+$JointB/2.0] [expr $Y2]
node 802 [expr $X3] [expr $Y2+$JointH/2.0]
node 803 [expr $X3-$JointB/2.0] [expr $Y2]
node 804 [expr $X3] [expr $Y2-$JointH/2.0]
node 901 [expr $X4+$JointB/2.0] [expr $Y2]
node 902 [expr $X4] [expr $Y2+$JointH/2.0]
node 903 [expr $X4-$JointB/2.0] [expr $Y2]
node 904 [expr $X4] [expr $Y2-$JointH/2.0]
node 1001 [expr $X5+$JointB/2.0] [expr $Y2]
node 1002 [expr $X5] [expr $Y2+$JointH/2.0]
node 1003 [expr $X5-$JointB/2.0] [expr $Y2]
node 1004 [expr $X5] [expr $Y2-$JointH/2.0]
# 二层
node 1101 [expr $X1+$JointB/2.0] [expr $Y3]
node 1102 [expr $X1] [expr $Y3+$JointH/2.0]
node 1103 [expr $X1-$JointB/2.0] [expr $Y3]
node 1104 [expr $X1] [expr $Y3-$JointH/2.0]
node 1201 [expr $X2+$JointB/2.0] [expr $Y3]
node 1202 [expr $X2] [expr $Y3+$JointH/2.0]
node 1203 [expr $X2-$JointB/2.0] [expr $Y3]
node 1204 [expr $X2] [expr $Y3-$JointH/2.0]
node 1301 [expr $X3+$JointB/2.0] [expr $Y3]
node 1302 [expr $X3] [expr $Y3+$JointH/2.0]
node 1303 [expr $X3-$JointB/2.0] [expr $Y3]
node 1304 [expr $X3] [expr $Y3-$JointH/2.0]
node 1401 [expr $X4+$JointB/2.0] [expr $Y3]
node 1402 [expr $X4] [expr $Y3+$JointH/2.0]
node 1403 [expr $X4-$JointB/2.0] [expr $Y3]
node 1404 [expr $X4] [expr $Y3-$JointH/2.0]
node 1501 [expr $X5+$JointB/2.0] [expr $Y3]
node 1502 [expr $X5] [expr $Y3+$JointH/2.0]
node 1503 [expr $X5-$JointB/2.0] [expr $Y3]
node 1504 [expr $X5] [expr $Y3-$JointH/2.0]
# 三层
node 1601 [expr $X1+$JointB/2.0] [expr $Y4]
node 1602 [expr $X1] [expr $Y4+$JointH/2.0]
node 1603 [expr $X1-$JointB/2.0] [expr $Y4]
node 1604 [expr $X1] [expr $Y4-$JointH/2.0]
node 1701 [expr $X2+$JointB/2.0] [expr $Y4]
node 1702 [expr $X2] [expr $Y4+$JointH/2.0]
node 1703 [expr $X2-$JointB/2.0] [expr $Y4]
node 1704 [expr $X2] [expr $Y4-$JointH/2.0]
node 1801 [expr $X3+$JointB/2.0] [expr $Y4]
node 1802 [expr $X3] [expr $Y4+$JointH/2.0]
node 1803 [expr $X3-$JointB/2.0] [expr $Y4]
node 1804 [expr $X3] [expr $Y4-$JointH/2.0]
node 1901 [expr $X4+$JointB/2.0] [expr $Y4]
node 1902 [expr $X4] [expr $Y4+$JointH/2.0]
node 1903 [expr $X4-$JointB/2.0] [expr $Y4]
node 1904 [expr $X4] [expr $Y4-$JointH/2.0]
node 2001 [expr $X5+$JointB/2.0] [expr $Y4]
node 2002 [expr $X5] [expr $Y4+$JointH/2.0]
node 2003 [expr $X5-$JointB/2.0] [expr $Y4]
node 2004 [expr $X5] [expr $Y4-$JointH/2.0]
#========================施加约束===========================
fixY 0 1 1 1 # 对平面内所有Y=0的节点施加约束
#========================定义材料===========================
# 柱-保护层非核心混凝土材料细节
set CUnconfFc -26.1
set CUnconfEc -0.002
set CUnconfFcu -11.6
set CUnconfEcu -0.004
# 柱-核心混凝土材料细节
set CCconfFc -28.5
set CConfEc -0.0022
set CConfFcu -20.3
set CConfEcu -0.0112
#==============================柱混凝土=======================
# CONCRETE manual tag f'c ec0 f'cu ecu
uniaxialMaterial Concrete01 10 $CUnconfFc $CUnconfEc $CUnconfFcu $CUnconfEcu; # 柱保护层混凝土
uniaxialMaterial Concrete01 20 $CCconfFc $CConfEc $CConfFcu $CConfEcu; # 柱核心混凝土
# 梁-保护层非核心混凝土细节
set BUnconfFc -26.1
set BUnconfEc -0.002
set BUnconfFcu -11.6
set BUnconfEcu -0.004
# 梁-核心混凝土细节
set BCconfFc -28.5
set BConfEc -0.0022
set BConfFcu -20.3
set BConfEcu -0.0112
#================================梁混凝土====================
# CONCRETE manual tag f'c ec0 f'cu ecu
uniaxialMaterial Concrete01 40 $BUnconfFc $BUnconfEc $BUnconfFcu $BUnconfEcu; # 梁保护层混凝土
uniaxialMaterial Concrete01 50 $BCconfFc $BConfEc $BConfFcu $BConfEcu; # 梁核心混凝土
# 钢筋
set Fy 388
set Es 2.0e5
set Hratio 0.001
# tag fy E0 b
# uniaxialMaterial Steel01 3 $fy $E $b; # harden modle
uniaxialMaterial Steel02 30 $Fy $Es $Hratio 18.5 0.925 0.15; # harden modle
#========================定义非线性柱单元截面===========================
#定义梁柱配筋直径
set Ac16 201; #Ac16- area of individual steal bar with the diameter of 16mm
set Ac18 254;
#
# |----------- H -----------|
# j----------------------------i — —
# | o o | ┆ -- cover
# | | ┆
# | o y o | ┆
# | ←――― + | B
# | o ┆ o | ┆
# | ↓ Z | ┆
# | o o | ┆ -- cover
# k ----------------------------l — —
#划分纤维柱网格定位坐标数据 some viriables from parameters
set c_yPOUT [expr $Hc/2] ; #P POSITIVE 正值
set c_zPOUT [expr $Bc/2]
set c_yNOUT [expr -$Hc/2] ; #N NEGITIVE 负值
set c_zNOUT [expr -$Bc/2]
set c_yPIN [expr $Hc/2-$cov]
set c_zPIN [expr $Bc/2-$cov]
set c_yNIN [expr $cov-$Hc/2]
set c_zNIN [expr $cov-$Bc/2]
#定义框架柱纤维截面
section Fiber 1 {
# creat the concreat cover fibers
# nfIJ nfJK yI zI yJ zJ yK zK yL zL
patch quad 10 10 1 $c_yNIN $c_zNOUT $c_yPIN $c_zNOUT $c_yPIN $c_zNIN $c_yNIN $c_zNIN
patch quad 10 10 1 $c_yNIN $c_zPIN $c_yPIN $c_zPIN $c_yPIN $c_zPOUT $c_yNIN $c_zPOUT
patch quad 10 2 2 $c_yPIN $c_zNOUT $c_yPOUT $c_zNOUT $c_yPOUT $c_zPOUT $c_yPIN $c_zPOUT
patch quad 10 2 2 $c_yNOUT $c_zNOUT $c_yNIN $c_zNOUT $c_yNIN $c_zPOUT $c_yNOUT $c_zPOUT
# creat the concreat core fibers
patch quad 20 10 2 $c_yNIN $c_zNIN $c_yPIN $c_zNIN $c_yPIN $c_zPIN $c_yNIN $c_zPIN
# creat the reinforced fibers
# areaBar yStat zStart yEnd zEnd
#平面内钢筋
layer straight 30 4 $Ac16 $c_yNIN $c_zNIN $c_yNIN $c_zPIN
layer straight 30 4 $Ac16 $c_yPIN $c_zNIN $c_yPIN $c_zPIN
#平面外钢筋
layer straight 30 2 $Ac16 73.33 $c_zNIN 73.33 $c_zPIN
layer straight 30 2 $Ac16 -73.33 $c_zNIN -73.33 $c_zPIN
} #####
#========================定义非线性梁单元截面===========================
# j ——————————i ———
# | o o o o | ┆—— cover
# | ∧y | ┆
# | ┆ | ┆
# | ┆ | ┆
# | Z ┆ | ┆
# | <------+ | H
# | | ┆
# | | ┆
# | | ┆
# | | ┆
# | o o o o | ┆—— cover
# k ——————————l
# |------- B -------|
# 划分纤维梁网格定位坐标数据 some viriables from parameters
set b_yPOUT [expr $Hb/2]; #P POSITIVE
set b_zPOUT [expr $Bb/2]
set b_yNOUT [expr -$Hb/2]; #N NEGITIVE
set b_zNOUT [expr -$Bb/2]
set b_yPIN [expr $Hb/2-$cov]
set b_zPIN [expr $Bb/2-$cov]
set b_yNIN [expr $cov-$Hb/2]
set b_zNIN [expr $cov-$Bb/2]
#####
section Fiber 2 {
# creat the concrete fibers
# nfIJ nfJK yI zI yJ zJ yK zK yL zL
patch quad 40 10 1 $b_yNIN $b_zNOUT $b_yPIN $b_zNOUT $b_yPIN $b_zNIN $b_yNIN $b_zNIN
patch quad 40 10 1 $b_yNIN $b_zPIN $b_yPIN $b_zPIN $b_yPIN $b_zPOUT $b_yNIN $b_zPOUT
patch quad 40 2 2 $b_yPIN $b_zNOUT $b_yPOUT $b_zNOUT $b_yPOUT $b_zPOUT $b_yPIN $b_zPOUT
patch quad 40 2 2 $b_yNOUT $b_zNOUT $b_yNIN $b_zNOUT $b_yNIN $b_zPOUT $b_yNOUT $b_zPOUT
# creat the concreat core fibers
patch quad 50 10 2 $b_yNIN $b_zNIN $b_yPIN $b_zNIN $b_yPIN $b_zPIN $b_yNIN $b_zPIN
# creat the reinforced fibers
#numBar areaBar yStat zStart yEnd zEnd
layer straight 30 2 $Ac16 $b_yPIN $b_zNIN $b_yPIN $b_zPIN
layer straight 30 1 $Ac18 $b_yPIN 0.0 $b_yPIN 0.0
layer straight 30 3 $Ac16 $b_yNIN $b_zNIN $b_yNIN $b_zPIN
} #####
#========================定义梁柱单元===========================
# 框架柱单元
#element beamWithHinges $eleTag $iNode $jNode $secTagI $Lpi $secTagJ $Lpj $E $A $Iz $transfTag <-mass $massDens> <-iter $maxIters $tol>
set np 5
set colLp [expr $Hc*1.0]
set colE 3.356e4
set colA [expr $Hc*$Bc]
set colIz [expr $Bc*$Hc*$Hc*$Hc/12.0]
geomTransf PDelta 1
# 定义首层柱单元# 柱编号从1开始#############
element beamWithHinges 1 101 604 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 2 201 704 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 3 301 804 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 4 401 904 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 5 501 1004 1 $colLp 1 $colLp $colE $colA $colIz 1
# 定义二层柱单元
element beamWithHinges 6 602 1104 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 7 702 1204 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 8 802 1304 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 9 902 1404 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 10 1002 1504 1 $colLp 1 $colLp $colE $colA $colIz 1
# 定义三层柱单元
element beamWithHinges 11 1102 1604 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 12 1202 1704 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 13 1302 1804 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 14 1402 1904 1 $colLp 1 $colLp $colE $colA $colIz 1
element beamWithHinges 15 1502 2004 1 $colLp 1 $colLp $colE $colA $colIz 1
# 框架梁单元
set beamLp [expr $Hb*1.0]
set beamE 36744.3798
set beamA [expr $Hb*$Bb]
set beamIz [expr $Bb*$Hb*$Hb*$Hb/12.0]
geomTransf Linear 2
# 定义首层梁单元#梁编号从101开始#############
element beamWithHinges 101 601 703 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
element beamWithHinges 102 701 803 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
element beamWithHinges 103 801 903 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
element beamWithHinges 104 901 1003 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
# 定义首层梁单元
element beamWithHinges 105 1101 1203 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
element beamWithHinges 106 1201 1303 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
element beamWithHinges 107 1301 1403 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
element beamWithHinges 108 1401 1503 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
# 定义首层梁单元
element beamWithHinges 109 1601 1703 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
element beamWithHinges 110 1701 1803 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
element beamWithHinges 111 1801 1903 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
element beamWithHinges 112 1901 2003 2 $beamLp 2 $beamLp $beamE $beamA $beamIz 2
# 节点单元 #节点单元编号从1001开始
# 载入PR1.tcl文件
source PR1.tcl
#========================施加节点质量=======================
set p6 186662.5; # point gravity load in node 6 10 11 15 # 单位：N 牛顿
set p7 278050; # point gravity load in node 7 8 9 12 13 14
set p16 151487.5; # point gravity load in node 16 20
set p17 223750; # point gravity load in node 17 18 19
set g 9800
mass 601 [expr ($p6)/$g/3] [expr ($p6)/$g/3] 0.0
mass 602 [expr ($p6)/$g/3] [expr ($p6)/$g/3] 0.0
mass 603 0.0 0.0 0.0
mass 604 [expr ($p6)/$g/3] [expr ($p6)/$g/3] 0.0
mass 701 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 702 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 703 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 704 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 801 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 802 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 803 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 804 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 901 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 902 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 903 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 904 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1001 0.0 0.0 0.0
mass 1002 [expr ($p6)/$g/3] [expr ($p6)/$g/3] 0.0
mass 1003 [expr ($p6)/$g/3] [expr ($p6)/$g/3] 0.0
mass 1004 [expr ($p6)/$g/3] [expr ($p6)/$g/3] 0.0
mass 1101 [expr ($p6)/$g/3] [expr ($p6)/$g/3] 0.0
mass 1102 [expr ($p6)/$g/3] [expr ($p6)/$g/3] 0.0
mass 1103 0.0 0.0 0.0
mass 1104 [expr ($p6)/$g/4] [expr ($p6)/$g/4] 0.0
mass 1201 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1202 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1203 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1204 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1301 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1302 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1303 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1304 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1401 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1402 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1403 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1404 [expr ($p7)/$g/4] [expr ($p7)/$g/4] 0.0
mass 1501 0.0 0.0 0.0
mass 1502 [expr ($p6)/$g/3] [expr ($p6)/$g/3] 0.0
mass 1503 [expr ($p6)/$g/3] [expr ($p6)/$g/3] 0.0
mass 1504 [expr ($p6)/$g/3] [expr ($p6)/$g/3] 0.0
mass 1601 [expr ($p16)/$g/2] [expr ($p16)/$g/2] 0.0
mass 1602 0.0 0.0 0.0
mass 1603 0.0 0.0 0.0
mass 1604 [expr ($p16)/$g/2] [expr ($p16)/$g/2] 0.0
mass 1701 [expr ($p17)/$g/3] [expr ($p17)/$g/3] 0.0
mass 1702 0.0 0.0 0.0
mass 1703 [expr ($p17)/$g/3] [expr ($p17)/$g/3] 0.0
mass 1704 [expr ($p17)/$g/3] [expr ($p17)/$g/3] 0.0
mass 1801 [expr ($p17)/$g/3] [expr ($p17)/$g/3] 0.0
mass 1802 0.0 0.0 0.0
mass 1803 [expr ($p17)/$g/3] [expr ($p17)/$g/3] 0.0
mass 1804 [expr ($p17)/$g/3] [expr ($p17)/$g/3] 0.0
mass 1901 [expr ($p17)/$g/3] [expr ($p17)/$g/3] 0.0
mass 1902 0.0 0.0 0.0
mass 1903 [expr ($p17)/$g/3] [expr ($p17)/$g/3] 0.0
mass 1904 [expr ($p17)/$g/3] [expr ($p17)/$g/3] 0.0
mass 2001 0.0 0.0 0.0
mass 2002 0.0 0.0 0.0
mass 2003 [expr ($p16)/$g/2] [expr ($p16)/$g/2] 0.0
mass 2004 [expr ($p16)/$g/2] [expr ($p16)/$g/2] 0.0
#########输出周期##################
set pi 3.1415926
set lambda [eigen -generalized -fullGenLapack 3];
set lambda1 [lindex $lambda 0];
set lambda2 [lindex $lambda 1];
set lambda3 [lindex $lambda 2];
set omega1 [expr pow($lambda1,0.5)]
set omega2 [expr pow($lambda2,0.5)]
set omega3 [expr pow($lambda3,0.5)]
set T1 [expr 2*$pi/$omega1]
set T2 [expr 2*$pi/$omega2]
set T3 [expr 2*$pi/$omega3]
set fileT [open T.txt w]
puts -nonewline $fileT "$T1 $T2 $T3"
close $fileT
puts "omega1=$omega1"
puts "omega3=$omega3"
puts " Period output "
puts " T1=$T1 "
puts " T2=$T2 "
puts " T3=$T3 "
puts "==========================="
RC frame including BeamColumnJoint Element
Moderators: silvia, selimgunay, Moderators
-
civilape_anhui
- Posts: 9
- Joined: Tue Mar 29, 2011 3:34 am
- Location: China
-
civilape_anhui
- Posts: 9
- Joined: Tue Mar 29, 2011 3:34 am
- Location: China
Re: RC frame including BeamColumnJoint Element
## unit name ---- PR1
set fName "PR1";
source procUniaxialPinching.tcl
#======================================梁钢筋滑移材料的细节======================================
#=========================details for the material models of bar slip of the beam===============
# set bs_fc [expr -$BUnconfFc]; set bs_fs $Fy; set bs_es $Es; set bs_fsu 543; set bs_dbar 16; set bs_esh [expr $Hratio*$Es];
set bs_fc [expr -$BUnconfFc]; set bs_fs $Fy; set bs_es $Es; set bs_fsu 543; set bs_dbar 16; set bs_esh [expr $Hratio*$Es];
set bs_wid $Bc; set bs_dep $Hb;
set bsT_nbars 3; set bsB_nbars 3;
set bs_ljoint $Hc;
#======================================柱钢筋滑移材料的细节=========================================
#========================details for the material models of bar slip of the column=================
# set cs_fc [expr -$CUnconfFc]; set cs_fs $Fy; set cs_es $Es; set cs_fsu 543; set cs_dbar 16; set cs_esh [expr $Hratio*$Es];
set cs_fc [expr -$CUnconfFc]; set cs_fs $Fy; set cs_es $Es; set cs_fsu 543; set cs_dbar 16; set cs_esh [expr $Hratio*$Es];
set cs_wid $Bc; set cs_dep $Hc;
set cs_nbars 2;
set cs_ljoint $Hb;
#==========================================交界面剪切材料=============================================
#command: uniaxialMaterial Elastic tag? E?
uniaxialMaterial Elastic 1 10000000000.0
#=======================================定义钢筋滑移BarSlip材料=======================================
#==============for beam bottom
set matID1 21
set matID2 22
uniaxialMaterial BarSlip $matID1 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsB_nbars $bs_wid $bs_dep strong beamBot
uniaxialMaterial BarSlip $matID2 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsB_nbars $bs_wid $bs_dep strong beamBot
## %%%%%%%%%%%%%%% equivalent statement can be made in other way
#uniaxialMaterial BarSlip $matID1 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsB_nbars $bs_wid $bs_dep 1.0 strong beamBot damage
#uniaxialMaterial BarSlip $matID2 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsB_nbars $bs_wid $bs_dep 1.0 strong beamBot damage
#===============for beam top
set matID3 31
set matID4 32
uniaxialMaterial BarSlip $matID3 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsT_nbars $bs_wid $bs_dep strong beamTop
uniaxialMaterial BarSlip $matID4 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsT_nbars $bs_wid $bs_dep strong beamTop
## %%%%%%%%%%%%%%% equivalent statement can be made in other way
#uniaxialMaterial BarSlip $matID3 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsT_nbars $bs_wid $bs_dep 1.0 strong beamTop damage
#uniaxialMaterial BarSlip $matID4 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsT_nbars $bs_wid $bs_dep 1.0 strong beamTop damage
#===============for columns
set matID5 41
set matID6 42
set matID7 43
set matID8 44
uniaxialMaterial BarSlip $matID5 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep strong column
uniaxialMaterial BarSlip $matID6 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep strong column
uniaxialMaterial BarSlip $matID7 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep strong column
uniaxialMaterial BarSlip $matID8 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep strong column
## %%%%%%%%%%%%%%% equivalent statement can be made in other way
#uniaxialMaterial BarSlip $matID5 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep 1.0 strong column damage
#uniaxialMaterial BarSlip $matID6 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep 1.0 strong column damage
#uniaxialMaterial BarSlip $matID7 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep 1.0 strong column damage
#uniaxialMaterial BarSlip $matID8 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep 1.0 strong column damage
#==============================================定义剪切块材料============================================================
#=========================================material for shear panel======================================================
## Positive/Negative envelope Stress
set p1 2.1932; set p2 4.0872; set p3 4.4862; set p4 [expr $p3*1e-3];
#set p1 4.1932; set p2 6.0872; set p3 6.4862; set p4 [expr $p3*1e-3];
## stress1 stress2 stress3 stress4
set pEnvStrsp [list [expr $p1*$JointV] [expr $p2*$JointV] [expr $p3*$JointV] [expr $p4*$JointV]]
set nEnvStrsp [list [expr -$p1*$JointV] [expr -$p2*$JointV] [expr -$p3*$JointV] [expr -$p4*$JointV]]
## Positive/Negative envelope Strain
## strain1 strain2 strain3 strain4
set pEnvStnsp [list 0.0002 0.004465 0.0131 0.0269]
set nEnvStnsp [list -0.0002 -0.004465 -0.0131 -0.0269]
## Ratio of maximum deformation at which reloading begins
## Pos_env. Neg_env.
set rDispsp [list 0.25 0.25]
## Ratio of envelope force (corresponding to maximum deformation) at which reloading begins
### Pos_env. Neg_env.
set rForcesp [list 0.15 0.15]
## Ratio of monotonic strength developed upon unloading
### Pos_env. Neg_env.
set uForcesp [list 0.0 0.0]
## Coefficients for Unloading Stiffness degradation
## gammaK1 gammaK2 gammaK3 gammaK4 gammaKLimit
set gammaKsp [list 1.13364492409642 0.0 0.10111033064469 0.0 0.91652498468618]
#set gammaKsp [list 0.0 0.0 0.0 0.0 0.0]
#### Coefficients for Reloading Stiffness degradation
### gammaD1 gammaD2 gammaD3 gammaD4 gammaDLimit
set gammaDsp [list 0.12 0.0 0.23 0.0 0.95]
#set gammaDsp [list 0.0 0.0 0.0 0.0 0.0]
#### Coefficients for Strength degradation
### gammaF1 gammaF2 gammaF3 gammaF4 gammaFLimit
set gammaFsp [list 1.11 0.0 0.319 0.0 0.125]
#set gammaFsp [list 0.0 0.0 0.0 0.0 0.0]
set gammaEsp 10.0
uniaxialMaterial Pinching4 5 [lindex $pEnvStrsp 0] [lindex $pEnvStnsp 0] \
[lindex $pEnvStrsp 1] [lindex $pEnvStnsp 1] [lindex $pEnvStrsp 2] \
[lindex $pEnvStnsp 2] [lindex $pEnvStrsp 3] [lindex $pEnvStnsp 3] \
[lindex $nEnvStrsp 0] [lindex $nEnvStnsp 0] \
[lindex $nEnvStrsp 1] [lindex $nEnvStnsp 1] [lindex $nEnvStrsp 2] \
[lindex $nEnvStnsp 2] [lindex $nEnvStrsp 3] [lindex $nEnvStnsp 3] \
[lindex $rDispsp 0] [lindex $rForcesp 0] [lindex $uForcesp 0] \
[lindex $rDispsp 1] [lindex $rForcesp 1] [lindex $uForcesp 1] \
[lindex $gammaKsp 0] [lindex $gammaKsp 1] [lindex $gammaKsp 2] [lindex $gammaKsp 3] [lindex $gammaKsp 4] \
[lindex $gammaDsp 0] [lindex $gammaDsp 1] [lindex $gammaDsp 2] [lindex $gammaDsp 3] [lindex $gammaDsp 4] \
[lindex $gammaFsp 0] [lindex $gammaFsp 1] [lindex $gammaFsp 2] [lindex $gammaFsp 3] [lindex $gammaFsp 4] \
$gammaEsp energy
puts "uniaxialMaterial Pinching4 5 [lindex $pEnvStrsp 0] [lindex $pEnvStnsp 0] \
[lindex $pEnvStrsp 1] [lindex $pEnvStnsp 1] [lindex $pEnvStrsp 2] \
[lindex $pEnvStnsp 2] [lindex $pEnvStrsp 3] [lindex $pEnvStnsp 3] \
[lindex $nEnvStrsp 0] [lindex $nEnvStnsp 0] \
[lindex $nEnvStrsp 1] [lindex $nEnvStnsp 1] [lindex $nEnvStrsp 2] \
[lindex $nEnvStnsp 2] [lindex $nEnvStrsp 3] [lindex $nEnvStnsp 3] \
[lindex $rDispsp 0] [lindex $rForcesp 0] [lindex $uForcesp 0] \
[lindex $rDispsp 1] [lindex $rForcesp 1] [lindex $uForcesp 1] \
[lindex $gammaKsp 0] [lindex $gammaKsp 1] [lindex $gammaKsp 2] [lindex $gammaKsp 3] [lindex $gammaKsp 4] \
[lindex $gammaDsp 0] [lindex $gammaDsp 1] [lindex $gammaDsp 2] [lindex $gammaDsp 3] [lindex $gammaDsp 4] \
[lindex $gammaFsp 0] [lindex $gammaFsp 1] [lindex $gammaFsp 2] [lindex $gammaFsp 3] [lindex $gammaFsp 4] \
$gammaEsp energy"
#return
#====================================定义节点单元=================================================
##element BeamColumnJoint tag? iNode? jNode? kNode? lNode? matTag1? matTag2? matTag3? matTag4?
## matTag5? matTag6? matTag7? matTag8? matTag9? matTag10? matTag11? matTag12? matTag13?
## <element Height factor?> <element Width factor?>
## please note: the four nodes are in anticlockwise direction around the element
## requires material tags for all 13 different components within the element.
## the first 12 being that of spring and the last of the shear panel
# 首层梁柱节点单元#单元编号从1001开始
element beamColumnJoint 1001 604 601 602 603 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1002 704 701 702 703 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1003 804 801 802 803 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1004 904 901 902 903 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1005 1004 1001 1002 1003 41 42 1 21 31 1 43 44 1 22 32 1 5
# 二层梁柱节点单元
element beamColumnJoint 1006 1104 1101 1102 1103 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1007 1204 1201 1202 1203 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1008 1304 1301 1302 1303 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1009 1404 1401 1402 1403 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1010 1504 1501 1502 1503 41 42 1 21 31 1 43 44 1 22 32 1 5
# 三层梁柱节点单元
element beamColumnJoint 1011 1604 1601 1602 1603 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1012 1704 1701 1702 1703 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1013 1804 1801 1802 1803 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1014 1904 1901 1902 1903 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1015 2004 2001 2002 2003 41 42 1 21 31 1 43 44 1 22 32 1 5
set fName "PR1";
source procUniaxialPinching.tcl
#======================================梁钢筋滑移材料的细节======================================
#=========================details for the material models of bar slip of the beam===============
# set bs_fc [expr -$BUnconfFc]; set bs_fs $Fy; set bs_es $Es; set bs_fsu 543; set bs_dbar 16; set bs_esh [expr $Hratio*$Es];
set bs_fc [expr -$BUnconfFc]; set bs_fs $Fy; set bs_es $Es; set bs_fsu 543; set bs_dbar 16; set bs_esh [expr $Hratio*$Es];
set bs_wid $Bc; set bs_dep $Hb;
set bsT_nbars 3; set bsB_nbars 3;
set bs_ljoint $Hc;
#======================================柱钢筋滑移材料的细节=========================================
#========================details for the material models of bar slip of the column=================
# set cs_fc [expr -$CUnconfFc]; set cs_fs $Fy; set cs_es $Es; set cs_fsu 543; set cs_dbar 16; set cs_esh [expr $Hratio*$Es];
set cs_fc [expr -$CUnconfFc]; set cs_fs $Fy; set cs_es $Es; set cs_fsu 543; set cs_dbar 16; set cs_esh [expr $Hratio*$Es];
set cs_wid $Bc; set cs_dep $Hc;
set cs_nbars 2;
set cs_ljoint $Hb;
#==========================================交界面剪切材料=============================================
#command: uniaxialMaterial Elastic tag? E?
uniaxialMaterial Elastic 1 10000000000.0
#=======================================定义钢筋滑移BarSlip材料=======================================
#==============for beam bottom
set matID1 21
set matID2 22
uniaxialMaterial BarSlip $matID1 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsB_nbars $bs_wid $bs_dep strong beamBot
uniaxialMaterial BarSlip $matID2 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsB_nbars $bs_wid $bs_dep strong beamBot
## %%%%%%%%%%%%%%% equivalent statement can be made in other way
#uniaxialMaterial BarSlip $matID1 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsB_nbars $bs_wid $bs_dep 1.0 strong beamBot damage
#uniaxialMaterial BarSlip $matID2 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsB_nbars $bs_wid $bs_dep 1.0 strong beamBot damage
#===============for beam top
set matID3 31
set matID4 32
uniaxialMaterial BarSlip $matID3 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsT_nbars $bs_wid $bs_dep strong beamTop
uniaxialMaterial BarSlip $matID4 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsT_nbars $bs_wid $bs_dep strong beamTop
## %%%%%%%%%%%%%%% equivalent statement can be made in other way
#uniaxialMaterial BarSlip $matID3 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsT_nbars $bs_wid $bs_dep 1.0 strong beamTop damage
#uniaxialMaterial BarSlip $matID4 $bs_fc $bs_fs $bs_es $bs_fsu $bs_esh $bs_dbar $bs_ljoint $bsT_nbars $bs_wid $bs_dep 1.0 strong beamTop damage
#===============for columns
set matID5 41
set matID6 42
set matID7 43
set matID8 44
uniaxialMaterial BarSlip $matID5 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep strong column
uniaxialMaterial BarSlip $matID6 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep strong column
uniaxialMaterial BarSlip $matID7 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep strong column
uniaxialMaterial BarSlip $matID8 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep strong column
## %%%%%%%%%%%%%%% equivalent statement can be made in other way
#uniaxialMaterial BarSlip $matID5 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep 1.0 strong column damage
#uniaxialMaterial BarSlip $matID6 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep 1.0 strong column damage
#uniaxialMaterial BarSlip $matID7 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep 1.0 strong column damage
#uniaxialMaterial BarSlip $matID8 $cs_fc $cs_fs $cs_es $cs_fsu $cs_esh $cs_dbar $cs_ljoint $cs_nbars $cs_wid $cs_dep 1.0 strong column damage
#==============================================定义剪切块材料============================================================
#=========================================material for shear panel======================================================
## Positive/Negative envelope Stress
set p1 2.1932; set p2 4.0872; set p3 4.4862; set p4 [expr $p3*1e-3];
#set p1 4.1932; set p2 6.0872; set p3 6.4862; set p4 [expr $p3*1e-3];
## stress1 stress2 stress3 stress4
set pEnvStrsp [list [expr $p1*$JointV] [expr $p2*$JointV] [expr $p3*$JointV] [expr $p4*$JointV]]
set nEnvStrsp [list [expr -$p1*$JointV] [expr -$p2*$JointV] [expr -$p3*$JointV] [expr -$p4*$JointV]]
## Positive/Negative envelope Strain
## strain1 strain2 strain3 strain4
set pEnvStnsp [list 0.0002 0.004465 0.0131 0.0269]
set nEnvStnsp [list -0.0002 -0.004465 -0.0131 -0.0269]
## Ratio of maximum deformation at which reloading begins
## Pos_env. Neg_env.
set rDispsp [list 0.25 0.25]
## Ratio of envelope force (corresponding to maximum deformation) at which reloading begins
### Pos_env. Neg_env.
set rForcesp [list 0.15 0.15]
## Ratio of monotonic strength developed upon unloading
### Pos_env. Neg_env.
set uForcesp [list 0.0 0.0]
## Coefficients for Unloading Stiffness degradation
## gammaK1 gammaK2 gammaK3 gammaK4 gammaKLimit
set gammaKsp [list 1.13364492409642 0.0 0.10111033064469 0.0 0.91652498468618]
#set gammaKsp [list 0.0 0.0 0.0 0.0 0.0]
#### Coefficients for Reloading Stiffness degradation
### gammaD1 gammaD2 gammaD3 gammaD4 gammaDLimit
set gammaDsp [list 0.12 0.0 0.23 0.0 0.95]
#set gammaDsp [list 0.0 0.0 0.0 0.0 0.0]
#### Coefficients for Strength degradation
### gammaF1 gammaF2 gammaF3 gammaF4 gammaFLimit
set gammaFsp [list 1.11 0.0 0.319 0.0 0.125]
#set gammaFsp [list 0.0 0.0 0.0 0.0 0.0]
set gammaEsp 10.0
uniaxialMaterial Pinching4 5 [lindex $pEnvStrsp 0] [lindex $pEnvStnsp 0] \
[lindex $pEnvStrsp 1] [lindex $pEnvStnsp 1] [lindex $pEnvStrsp 2] \
[lindex $pEnvStnsp 2] [lindex $pEnvStrsp 3] [lindex $pEnvStnsp 3] \
[lindex $nEnvStrsp 0] [lindex $nEnvStnsp 0] \
[lindex $nEnvStrsp 1] [lindex $nEnvStnsp 1] [lindex $nEnvStrsp 2] \
[lindex $nEnvStnsp 2] [lindex $nEnvStrsp 3] [lindex $nEnvStnsp 3] \
[lindex $rDispsp 0] [lindex $rForcesp 0] [lindex $uForcesp 0] \
[lindex $rDispsp 1] [lindex $rForcesp 1] [lindex $uForcesp 1] \
[lindex $gammaKsp 0] [lindex $gammaKsp 1] [lindex $gammaKsp 2] [lindex $gammaKsp 3] [lindex $gammaKsp 4] \
[lindex $gammaDsp 0] [lindex $gammaDsp 1] [lindex $gammaDsp 2] [lindex $gammaDsp 3] [lindex $gammaDsp 4] \
[lindex $gammaFsp 0] [lindex $gammaFsp 1] [lindex $gammaFsp 2] [lindex $gammaFsp 3] [lindex $gammaFsp 4] \
$gammaEsp energy
puts "uniaxialMaterial Pinching4 5 [lindex $pEnvStrsp 0] [lindex $pEnvStnsp 0] \
[lindex $pEnvStrsp 1] [lindex $pEnvStnsp 1] [lindex $pEnvStrsp 2] \
[lindex $pEnvStnsp 2] [lindex $pEnvStrsp 3] [lindex $pEnvStnsp 3] \
[lindex $nEnvStrsp 0] [lindex $nEnvStnsp 0] \
[lindex $nEnvStrsp 1] [lindex $nEnvStnsp 1] [lindex $nEnvStrsp 2] \
[lindex $nEnvStnsp 2] [lindex $nEnvStrsp 3] [lindex $nEnvStnsp 3] \
[lindex $rDispsp 0] [lindex $rForcesp 0] [lindex $uForcesp 0] \
[lindex $rDispsp 1] [lindex $rForcesp 1] [lindex $uForcesp 1] \
[lindex $gammaKsp 0] [lindex $gammaKsp 1] [lindex $gammaKsp 2] [lindex $gammaKsp 3] [lindex $gammaKsp 4] \
[lindex $gammaDsp 0] [lindex $gammaDsp 1] [lindex $gammaDsp 2] [lindex $gammaDsp 3] [lindex $gammaDsp 4] \
[lindex $gammaFsp 0] [lindex $gammaFsp 1] [lindex $gammaFsp 2] [lindex $gammaFsp 3] [lindex $gammaFsp 4] \
$gammaEsp energy"
#return
#====================================定义节点单元=================================================
##element BeamColumnJoint tag? iNode? jNode? kNode? lNode? matTag1? matTag2? matTag3? matTag4?
## matTag5? matTag6? matTag7? matTag8? matTag9? matTag10? matTag11? matTag12? matTag13?
## <element Height factor?> <element Width factor?>
## please note: the four nodes are in anticlockwise direction around the element
## requires material tags for all 13 different components within the element.
## the first 12 being that of spring and the last of the shear panel
# 首层梁柱节点单元#单元编号从1001开始
element beamColumnJoint 1001 604 601 602 603 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1002 704 701 702 703 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1003 804 801 802 803 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1004 904 901 902 903 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1005 1004 1001 1002 1003 41 42 1 21 31 1 43 44 1 22 32 1 5
# 二层梁柱节点单元
element beamColumnJoint 1006 1104 1101 1102 1103 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1007 1204 1201 1202 1203 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1008 1304 1301 1302 1303 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1009 1404 1401 1402 1403 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1010 1504 1501 1502 1503 41 42 1 21 31 1 43 44 1 22 32 1 5
# 三层梁柱节点单元
element beamColumnJoint 1011 1604 1601 1602 1603 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1012 1704 1701 1702 1703 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1013 1804 1801 1802 1803 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1014 1904 1901 1902 1903 41 42 1 21 31 1 43 44 1 22 32 1 5
element beamColumnJoint 1015 2004 2001 2002 2003 41 42 1 21 31 1 43 44 1 22 32 1 5
-
civilape_anhui
- Posts: 9
- Joined: Tue Mar 29, 2011 3:34 am
- Location: China
Re: RC frame including BeamColumnJoint Element
#####################################################################################################
# #
# procUniaxialPinching.tcl #
# procedure for activating the pinching material given its parameters in the form of list #
# created NM (nmitra@u.washington.edu) dated : Feb 2002 #
#####################################################################################################
proc procUniaxialPinching { materialTag pEnvelopeStress nEnvelopeStress pEnvelopeStrain nEnvelopeStrain rDisp rForce uForce gammaK gammaD gammaF gammaE damage} {
# add material - command: uniaxialMaterial ...... paramaters as shown
#uniaxialMaterial Pinching4 tag
#### stress1P strain1P stress2P strain2P stress3P strain3P stress4P strain4P
#### stress1N strain1N stress2N strain2N stress3N strain3N stress4N strain4N
#### rDispP rForceP uForceP rDispN rForceN uForceN
#### gammaK1 gammaK2 gammaK3 gammaK4 gammaKLimit
#### gammaD1 gammaD2 gammaD3 gammaD4 gammaDLimit
#### gammaF1 gammaF2 gammaF3 gammaF4 gammaFLimit gammaE $damage
uniaxialMaterial Pinching4 $materialTag [lindex $pEnvelopeStress 0] [lindex $pEnvelopeStrain 0] [lindex $pEnvelopeStress 1] [lindex $pEnvelopeStrain 1] [lindex $pEnvelopeStress 2] [lindex $pEnvelopeStrain 2] [lindex $pEnvelopeStress 3] [lindex $pEnvelopeStrain 3] [lindex $nEnvelopeStress 0] [lindex $nEnvelopeStrain 0] [lindex $nEnvelopeStress 1] [lindex $nEnvelopeStrain 1] [lindex $nEnvelopeStress 2] [lindex $nEnvelopeStrain 2] [lindex $nEnvelopeStress 3] [lindex $nEnvelopeStrain 3] [lindex $rDisp 0] [lindex $rForce 0] [lindex $uForce 0] [lindex $rDisp 1] [lindex $rForce 1] [lindex $uForce 1] [lindex $gammaK 0] [lindex $gammaK 1] [lindex $gammaK 2] [lindex $gammaK 3] [lindex $gammaK 4] [lindex $gammaD 0] [lindex $gammaD 1] [lindex $gammaD 2] [lindex $gammaD 3] [lindex $gammaD 4] [lindex $gammaF 0] [lindex $gammaF 1] [lindex $gammaF 2] [lindex $gammaF 3] [lindex $gammaF 4] $gammaE $damage
}
# #
# procUniaxialPinching.tcl #
# procedure for activating the pinching material given its parameters in the form of list #
# created NM (nmitra@u.washington.edu) dated : Feb 2002 #
#####################################################################################################
proc procUniaxialPinching { materialTag pEnvelopeStress nEnvelopeStress pEnvelopeStrain nEnvelopeStrain rDisp rForce uForce gammaK gammaD gammaF gammaE damage} {
# add material - command: uniaxialMaterial ...... paramaters as shown
#uniaxialMaterial Pinching4 tag
#### stress1P strain1P stress2P strain2P stress3P strain3P stress4P strain4P
#### stress1N strain1N stress2N strain2N stress3N strain3N stress4N strain4N
#### rDispP rForceP uForceP rDispN rForceN uForceN
#### gammaK1 gammaK2 gammaK3 gammaK4 gammaKLimit
#### gammaD1 gammaD2 gammaD3 gammaD4 gammaDLimit
#### gammaF1 gammaF2 gammaF3 gammaF4 gammaFLimit gammaE $damage
uniaxialMaterial Pinching4 $materialTag [lindex $pEnvelopeStress 0] [lindex $pEnvelopeStrain 0] [lindex $pEnvelopeStress 1] [lindex $pEnvelopeStrain 1] [lindex $pEnvelopeStress 2] [lindex $pEnvelopeStrain 2] [lindex $pEnvelopeStress 3] [lindex $pEnvelopeStrain 3] [lindex $nEnvelopeStress 0] [lindex $nEnvelopeStrain 0] [lindex $nEnvelopeStress 1] [lindex $nEnvelopeStrain 1] [lindex $nEnvelopeStress 2] [lindex $nEnvelopeStrain 2] [lindex $nEnvelopeStress 3] [lindex $nEnvelopeStrain 3] [lindex $rDisp 0] [lindex $rForce 0] [lindex $uForce 0] [lindex $rDisp 1] [lindex $rForce 1] [lindex $uForce 1] [lindex $gammaK 0] [lindex $gammaK 1] [lindex $gammaK 2] [lindex $gammaK 3] [lindex $gammaK 4] [lindex $gammaD 0] [lindex $gammaD 1] [lindex $gammaD 2] [lindex $gammaD 3] [lindex $gammaD 4] [lindex $gammaF 0] [lindex $gammaF 1] [lindex $gammaF 2] [lindex $gammaF 3] [lindex $gammaF 4] $gammaE $damage
}
-
civilape_anhui
- Posts: 9
- Joined: Tue Mar 29, 2011 3:34 am
- Location: China
Re: RC frame including BeamColumnJoint Element
who can tell me what's the problem??
the eigenvalue analysis must have some problems
the eigenvalue analysis must have some problems
Re: RC frame including BeamColumnJoint Element
I can not help you with this big model. If you make small portal frame I can take a look and try to help you.
-
civilape_anhui
- Posts: 9
- Joined: Tue Mar 29, 2011 3:34 am
- Location: China
Re: RC frame including BeamColumnJoint Element
vesna wrote:
> I can not help you with this big model. If you make small portal frame I
> can take a look and try to help you.
hi vesna
Could you give me your email ？
> I can not help you with this big model. If you make small portal frame I
> can take a look and try to help you.
hi vesna
Could you give me your email ？
Re: RC frame including BeamColumnJoint Element
Hi, I have also had issues while using this element. Were you able to solve this? Please offer some tips if you did.