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