#[ 0,0 ], #[ 1,0 ], #[ 0,1 ], #[ 1,1 ], #[ 0,2 ], #[ 1,2 ] #]) #elements = np.array([ #[ node a, node b, node c, property set number] #[ 0, 1, 3, 0 ], #[ 0, 3, 2, 1 ], #[ 2, 3, 5, 0 ], #[ 2, 5, 4, 0 ] #]) nodes, elements, groups = readmsh('paracircle.msh') # properties properties = np.array([ #[ k ] [ 2e1 ], [ 4e1 ] ]) loads = 10*np.ones((len(nodes),2)) for i in np.arange(len(nodes)): loads[i,0]=i #if nodes[i,0]**2+nodes[i,1]**2>=1: #loads[i,1]=0 # contour conditions #loads = np.array([ ##[ node number, fz]
#[ 0,0 ], #[ 1,0 ], #[ 0,1 ], #[ 1,1 ], #[ 0,2 ], #[ 1,2 ] #]) #elements = np.array([ #[ node a, node b, node c, property set number] #[ 0, 1, 3, 0 ], #[ 0, 3, 2, 1 ], #[ 2, 3, 5, 0 ], #[ 2, 5, 4, 0 ] #]) nodes, elements, groups = readmsh('rectangle.msh') # material and section properties properties = np.array([ #[ k ] [ 2e1 ], [ 4e1 ] ]) loads = 10*np.ones((len(nodes),2)) for i in np.arange(len(loads)): loads[i,0]=i if nodes[i,0]**2+nodes[i,1]**2>=1: loads[i,1]=0 # contour conditions #loads = np.array([ ##[ node number, fz]
#[ 0,0 ], #[ 1,0 ], #[ 0,1 ], #[ 1,1 ], #[ 0,2 ], #[ 1,2 ] #]) #elements = np.array([ #[ node a, node b, node c, property set number] #[ 0, 1, 3, 0 ], #[ 0, 3, 2, 1 ], #[ 2, 3, 5, 0 ], #[ 2, 5, 4, 0 ] #]) nodes, elements, groups = readmsh('heart.msh') # material and section properties properties = np.array([ #[ k ] [ 2e1 ], [ 4e1 ] ]) loads = -10*np.ones((len(nodes),2)) for i in np.arange(len(nodes)): loads[i,0]=i #if nodes[i,0]**2+nodes[i,1]**2>=1: #loads[i,1]=0 # contour conditions #loads = np.array([ ##[ node number, fz]