请参考相关opensees建模代码,建立以下opensees模型,圆形桥墩总长37米,上5米为直径1.7m,下32米为直径1.9米,使用钢筋HRB335,材料C30混凝土,保护层厚度10cm,需考虑纤维截面划分、钢筋插入及截面组装, 编写的代码请每行注释得仔细一些,包括每个位置的数字的含义,请注意,这是圆形桥墩,纤维截面划分请使用圆形相关指令
时间: 2024-04-09 10:28:18 浏览: 155
以下是根据您提供的信息编写的Opensees圆形桥墩模型代码。每行都有详细注释,包括每个位置数字的含义。
```tcl
# 清空之前的模型
wipe
# 定义材料属性
uniaxialMaterial Steel 1 300000 200000 0.01 # 钢筋材料,参数分别为:标号、强度、弹性模量、屈服应变
uniaxialMaterial Concrete 2 -30 -0.002 -0.01 0.002 # 混凝土材料,参数分别为:标号、压力强度、压力应变、拉伸强度、拉伸应变、断裂应变
# 定义截面属性
section Aggregator 1 {
patch circle 0.0 0.0 0.85 100 # 上部截面圆形区域定义,参数分别为:圆心x坐标、圆心y坐标、半径、纤维数量
patch circle 0.0 0.0 0.95 300 # 下部截面圆形区域定义,参数同上
}
# 定义钢筋属性
set num_bars 8 # 钢筋数量
set bar_diameter 1.2 # 钢筋直径
set bar_spacing [expr (2 * 3.14159 * 0.85) / $num_bars] # 钢筋间距,根据上部截面计算
# 定义截面组装
section Fiber 2 {
patch circle 0.0 0.0 0.85 100 # 上部截面圆形区域定义,参数同上
layer straight $bar_diameter $num_bars $bar_spacing 0.0 # 上部截面钢筋插入,参数分别为:钢筋直径、钢筋数量、钢筋间距、钢筋起始深度
patch circle 0.0 0.0 0.95 300 # 下部截面圆形区域定义,参数同上
}
# 定义桥墩几何属性
set L_upper 5.0 # 上部长度
set L_lower 32.0 # 下部长度
set D_upper 1.7 # 上部直径
set D_lower 1.9 # 下部直径
set cover 0.1 # 保护层厚度
# 定义桥墩节点和单元
node 1 0.0 0.0 # 桥墩底部节点
node 2 [expr $L_upper] 0.0 # 上部截面节点
node 3 [expr $L_upper + $L_lower] 0.0 # 桥墩顶部节点
element zeroLength SectionAggregator 1 1 2 # 上部截面单元
element zeroLength SectionAggregator 2 2 3 # 下部截面单元
# 施加边界条件
fix 1 1 1 1 # 底部节点固定
# 施加加载
pattern Plain 1 Linear {
load 3 0 -1000000 # 桥墩顶部节点施加向下的垂直力
}
# 定义分析参数
set delta_t 1.0 # 时间步长
set T_total 10.0 # 总分析时间
# 执行动力分析
integrator Newmark $delta_t $delta_t/2.0 # Newmark积分器,参数分别为:时间步长、加速度积分系数
analysis Transient # 动力分析
analyze $T_total # 执行动力分析
```
以上是根据您提供的要求编写的Opensees圆形桥墩模型代码。每行都有详细的注释,包括每个位置数字的含义。请根据您的实际需求进行适当修改和调整。
阅读全文