#// imperfection to plate nodes

pi=3.1415926

pImp0=PB/200

p=mdb.models['Model-1'].parts['plate']

for n in range(len(p.nodes)):

x = p.nodes[n].coordinates[0]

y = p.nodes[n].coordinates[1]

offset = p.nodes[n].coordinates[2] + pImp0*cos(pi*x/PB)*cos(pi*y/PL)

p.editNode(nodes=p.nodes[n:n+1],offset3=offset,projectToGeometry=OFF)

#// imperfection to web nodes

pi=3.1415926

pImp0=WH/200

p = mdb.models['Model-1'].parts['FlgStiff']

for n in range(len(p.nodes)):

x = p.nodes[n].coordinates[0]

y = p.nodes[n].coordinates[1]

offset = p.nodes[n].coordinates[2] + pImp0*cos(pi*x/PB)*cos(pi*y/WH)

p.editNode(nodes=p.nodes[n:n+1],offset3=offset,projectToGeometry=OFF)