def main():
truss = Domain(3,3,'test03')
truss.add_mat(Material1D.Material_Elastic,1,2.e5)
truss.add_sec(Section.Elastic_Truss,1,truss.mat[1],100.)
L = 1000.0
mass = [200.,200.,200.]
truss.add_node(1,(0.0,0.,0.),mass)
truss.add_node(2,(L ,0.,0.),mass)
truss.add_refpoint(1,(1.,1.,0.))
truss.add_ele(Element.Truss3D,1,truss.sec[1],(truss.node[1],truss.node[2]),truss.refp[1])
F = -4.e4*np.linspace(0.1,1,10)
truss.add_load(1,truss.node[2],[1],[F])
truss.add_cons(1,truss.node[1],[1,2,3],[0.0,0.0,0.0])
truss.add_cons(2,truss.node[2],[2,3],[0.0,0.0])
truss.build_model()
truss.apply_cons(0)
analyse = Analysis.Analysis_Eigen(truss)
analyse.execute(1)
analyse = Analysis.Analysis_Linear_Static(truss)
analyse.execute(10)
return truss
def main():
truss = Domain(2,2,'test01')
truss.add_mat(Material1D.Material_Elastic,1,2.e5)
truss.add_sec(Section.Elastic_Truss,1,truss.mat[1],100.)
L = 1000.0
mass = [200.,200.,0.]
truss.add_node(1,(0.0,0.,0.),mass)
truss.add_node(2,(L ,0.,0.),mass)
truss.add_ele(Element.Truss2D,1,truss.sec[1],(truss.node[1],truss.node[2]))
F = -4.e4*np.linspace(0.1,1,10)
truss.add_load(1,truss.node[2],[1],[F])
truss.add_cons(1,truss.node[1],[1,2],[0.0,0.0])
truss.add_cons(2,truss.node[2],[2],[0.0])
truss.build_model()
truss.apply_cons(0)
analyse = Analysis.Analysis_Eigen(truss)
analyse.execute(1)
analyse = Analysis.Analysis_Linear_Static(truss)
analyse.execute(10)
print 'T = ',truss.T[0],'s.'
print 'Element Disp = ',truss.ele[1].local_deformation[0],'mm.'
return truss
