Paraview error

[twsleeper]

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]