cv.back = [False, False, False] m.velocity = Link(iob, "out") m.prescribed_velocity = Link(cv, "value_of_constraint") m.location_prescribed_velocity = Link(cv, "location_of_constraint") m.rel_tol = 0.01 m.expansion_coefficient = 0. m.bulk_modulus = 1000. m.shear_modulus = 1. m.plastic_stress = 0. m.friction_parameter = 0. m.dilatancy_parameter = 0. m.shear_length = m.shear_modulus * 0.9 ug = UpdateGeometry(debug=dbg) ug.domain = Link(dom, "domain") ug.displacement = Link(m, "displacement") vis = WriteVTK(debug=dbg) vis.t = Link(sq) vis.data0 = Link(m, "plastic_stress") vis.data1 = Link(m, "velocity") vis.data2 = Link(m, "stress") vis.dt = 0.01 vis.filename = os.path.join(WORKDIR, "temp.vtu") s = Simulation([sq, cv, m, vis], debug=True) s.writeXML() s.run()
m.velocity=Link(iob,"out") m.prescribed_velocity=Link(cv,"value_of_constraint") m.location_prescribed_velocity=Link(cv,"location_of_constraint") m.rel_tol=0.01 m.expansion_coefficient= 0. m.bulk_modulus=1000. m.shear_modulus=1. m.plastic_stress=0. m.friction_parameter=0. m.dilatancy_parameter=0. m.shear_length=m.shear_modulus*0.9 ug=UpdateGeometry(debug=dbg) ug.domain=Link(dom,"domain") ug.displacement=Link(m,"displacement") vis=WriteVTK(debug=dbg) vis.t=Link(sq) vis.data0=Link(m,"plastic_stress") vis.data1=Link(m,"velocity") vis.data2=Link(m,"stress") vis.dt=0.01 vis.filename=os.path.join(WORKDIR,"temp.vtu") s=Simulation([sq,cv,m,vis],debug=True) s.writeXML() s.run()