#[ 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]
