Wrench.slv

from extras.mesh import trigen

# domain
domain.planeStress(1.)
domain.gravity(0,-1,0)
db.SQLite("bench")

# materials
material.elastic(1,200e3,0.33,0.)

# mesh
trigen.file='wrench'
trigen.info=1
trigen.iDxf()
trigen.mesh(2,20.)

# loadcase
lc.define(1)
trigen.lineLoad(1,2,2,-2.,-2.)

# analysis
analysis.static()
reorder.rCM()
soe.band()
analysis.run(1,1)

#results
print "uy = ",node.data(1)['disp'][1]
db.store("wrench")
db.exportToVtk("wrench")