Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from fenics import *
- from mshr import *
- from math import pi, sin, cos, sqrt
- # Parameters
- R = 1.1
- r = 0.4
- t = 10
- x = R*cos(float(t) / 180 * pi)
- y = 0
- z = R*sin(t)
- # Create geometry
- s1 = Sphere(Point(0, 0, 0), 1)
- s2 = Sphere(Point(x, y, z), r)
- b1 = Box(Point(-2, -2, -0.03), Point(2, 2, 0.03))
- geometry = s1 - s2 - b1
- # Create mesh
- mesh = generate_mesh(geometry, 32)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement