使用新版opensees建立两跨的多跨简支梁桥时，在建立桥梁整体模型时，因为要研究梁体碰撞，所以对于梁与梁之间的缝隙处应设置碰撞单元，请你给我一个示例建模（两跨梁体均设置4个单元即可，加上缝隙处的碰撞单元以及两端与桥台之间的碰撞单元，碰撞单元选用hertz-damp）

以下是一个使用新版OpenSees建立两跨的多跨简支梁桥，并设置梁与梁之间的碰撞单元的示例建模代码：

# 导入OpenSees库
import openseespy.opensees as ops

# 清空模型
ops.wipe()

# 定义材料参数
E = 30000  # 弹性模量
nu = 0.2   # 泊松比
rho = 2500 # 密度

# 定义截面参数
width = 0.3   # 梁宽度
height = 0.5  # 梁高度

# 定义节点坐标
ops.node(1, 0.0, 0.0)
ops.node(2, 4.0, 0.0)
ops.node(3, 8.0, 0.0)
ops.node(4, 12.0, 0.0)
ops.node(5, 16.0, 0.0)

# 定义材料
ops.nDMaterial('ElasticIsotropic', 1, E, nu, rho)

# 定义截面
ops.section('Fiber', 1)
ops.patch('rect', 1, width/2, height/2, -width/2, -height/2)

# 创建梁单元
ops.beamIntegration('Legendre', 1, 'Lobatto', 1, 3)
ops.element('dispBeamColumn', 1, 1, 2, 1, '-cMass', 1)

# 设置碰撞单元
ops.nDMaterial('Hertz', 2, 1e9, 0.1, 1e6)
ops.element('twoNodeLink', 2, 2, 3, '-mat', 2, '-dir', 1)
ops.element('twoNodeLink', 3, 3, 4, '-mat', 2, '-dir', 1)
ops.element('twoNodeLink', 4, 4, 5, '-mat', 2, '-dir', 1)

# 设置边界条件
ops.fix(1, 1)
ops.fix(5, 1)

# 施加荷载
ops.timeSeries('Linear', 1)
ops.pattern('Plain', 1, 1)
ops.load(2, 0.0, -10.0)

# 定义分析类型
ops.constraints('Transformation')
ops.numberer('RCM')
ops.system('BandGeneral')
ops.test('NormUnbalance', 1e-6, 100)
ops.algorithm('Newton')
ops.integrator('LoadControl', 0.1)
ops.analysis('Static')

# 执行分析
ops.analyze(10)

# 输出结果
ops.printModel()
ops.printNodeResults()

在这个示例中，我们首先定义了材料参数、截面参数和节点坐标。然后创建了梁单元，并使用Hertz-Damp模型设置了梁与梁之间的碰撞单元。接下来设置了边界条件和施加荷载，并定义了分析类型。最后执行分析并输出结果。

请注意，这只是一个示例建模，你可能需要根据具体的要求进行调整和扩展。同时，确保你已经安装了OpenSees并正确配置了Python环境，以便能够运行上述代码。