The OpenSees Community

Hello everyone, I am using pvd recorder in OpenSeesPy and trying to post-process the results using Paraview. However,When I try to open and apply the pvd file in Paraview, an error message appears:
ERROR: In C:\glr\builds\paraview\paraviewci\build\superbuild\paraview\src\VTK\IO\XML\vtkXMLDataReader.cxx, line 410 vtkXMLUnstructuredGridReader (000001B56822CBC0): Cannot read point data array "Velocity" from PointData in piece 0. The data array in the element may be too short.

Please help me solve this error. Here's my OpenSeesPy code:
from openseespy.opensees import *
import matplotlib.pyplot as plt
import time
import os
wipe()
model('basic', '-ndm', 2, '-ndf' , 2)
E = 15005714.286
nu = 0.3
rho = 2020
nDMaterial('ElasticIsotropic', 2000, E, nu, rho)
nx = 2
ny = 100
e1 = 1
n1 = 1
eleArgs = [1, 'PlaneStrain', 2000]
points = [1, 0.0, 0.0,
2, 0.2, 0.0,
3, 0.2, 10.0,
4, 0.0, 10.0]
block2D(nx, ny, e1, n1,'quad', *eleArgs, *points)
# -------- Soil B.C ---------------
for i in range(ny+1):
equalDOF(3*i+1, 3*i+3, 1,2)
#------------ Load Pattern -------
timeSeries('Linear',705)
pattern('Plain',703, 705)
load(1,0,1)
load(2,0,1)
# ===Create an output directory====
filename = 'soilcolumn'
if not os.path.exists(filename):
os.makedirs(filename)
recorder('PVD', filename, 'vel','eleResponse', 'stresses')
system("UmfPack")
numberer("RCM")
constraints("Transformation")
integrator("Newmark", 0.5, 0.25)
algorithm("Newton")
test('EnergyIncr',1e-8, 200)
analysis("Transient")
analyze(8000,1e-4)
print("finish analyze:0 ~ 0.8s")
[https://ibb.co/nmz5FKV][https://ibb.co/xz99DTK]

Hello everyone, as I know we can use eleLoad -type -beamUniform to apply a distributed load on beam element. If I want to apply a uniformly distributed load on quad element, could I use the equalDOF command to connect an elastic BeamColumn element and a quad element, and use the eleLoad command to apply distributed load to transfer load to quad element ?