请参考相关opensees建模代码,建立以下opensees模型,圆形桥墩总长37米,上5米为直径1.7m,下32米为直径1.9米,使用钢筋HRB335,材料C30混凝土,保护层厚度10cm,需考虑纤维截面划分及组装, 编写的代码请每行注释得仔细一些,包括每个位置的数字的含义
时间: 2024-04-09 18:28:17 浏览: 193
下面是根据您提供的信息编写的基于OpenSees的圆形桥墩模型代码。请注意,注释中提供了每个位置数字的含义。
```tcl
# 清空之前的模型
wipe
# 定义材料属性
uniaxialMaterial Steel 1 3000 200000 0.01 # 钢筋材料,参数分别为:标号、强度、弹性模量、屈服应变
uniaxialMaterial Concrete 2 -30 -0.004 -0.01 0.002 # 混凝土材料,参数分别为:标号、压力强度、压力应变、拉伸强度、拉伸应变、断裂应变
# 定义截面属性
section Fiber 1 {
patch rect -0.85 -0.85 0.85 0.85 100 100 # 上部截面矩形区域定义,参数分别为:左下角x坐标、左下角y坐标、右上角x坐标、右上角y坐标、横向纤维数量、纵向纤维数量
patch rect -0.95 -0.95 0.95 0.95 100 300 # 下部截面矩形区域定义,参数同上
}
# 定义截面组装
section Aggregator 2 {
fiberSection 1 # 上部截面
fiberSection 1 # 下部截面
}
# 定义桥墩几何属性
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圆形桥墩模型代码,其中包含了每个位置数字的含义的详细注释。请根据您的实际需求进行适当修改和调整。
阅读全文