online advertising
Processing math: 100%

Saturday, January 23, 2016

Differential Geometry and Its Application - Exercise 3.2.18

Problem:


Solution:

This is a very long question, and we will tackle this part-by-part. This is for part (a)

Let's start with the formula:

xu=(β(u)+vδ(u))
xv=(δ(u))
xuu=(β(u)+vδ(u))
xuv=(δ(u))
xvv=(0)

The observation is that n=Uxvv=U0=0. So we have got the first equality.

Now we need to compute xu×xv

xu×xv=(β(u)+vδ(u))×δ(u)=β(u)×δ(u)+vδ(u)×δ(u)

So that explain the denominator, as per the hint. Finally we compute m

Uxuu=(β(u)×δ(u)+vδ(u)×δ(u))δ(u)=β(u)×δ(u)δ(u). So we finally also explain the numerator!

Notice the numerator is not exactly the same form I had, but it is the same because it is a circular shift of the scalar triple product.

Now we moved on to part (b), we can parametrize (x,y,xy)=(u,0,0)+v(0,1,u) to make it a ruled surface.

The code for computing the Gaussian curvature is as follow:

syms u
syms v
beta = [u; 0; 0];
delta = [0; 1; u];

beta_u = diff(beta, u);
delta_u = diff(delta, u);

n = -(beta_u.' * cross(delta, delta_u))^2
D = cross(beta_u, delta) + v * cross(delta_u, delta);
d = (D.' * D)^2;

K = simplify(n/d)

We get the answer as 1(u2+v2+1)2.

For part (c) and (d), the Gaussian curvatures for both cone and cylinder are 0, this because either β=0 or δ=0.

For part (e), the helicoid has the parametrization as (vcosu,vsinu,bu)=(0,0,u)+v(cosu,sinu,0), so we use essentially the same code above except

syms b;
beta = [0;0;b*u];
delta= [cos(u);sin(u);0];

So we get the answer as b2(b2+v2)2

For part (f), we will use the ruling patch we found in Exercise 2.1.22, so we simply put in yet another β and δ into the program and get the answer:

a2b2c2(a2b2v2+a2c2(vcosusinu)2+b2c2(vsinu+cosu)2)2

For part (g), in some sense, we have already done with it, for the saddle z=xy is a hyperbolic paraboloid.

For a more general hyperbolic paraboloid, we consider

zc=y2b2x2a2

Now let u=yb+xa and v=ybxa

(a2(uv),b2(u+v),cuv)

Now we obtain the ruled patch (a2u,b2u,0)+v(a2,b2,cu)

Because we wanted the expression to show values in terms of x and y this time, so we modified the program a bit as follow:

% defining the ruled patch
syms a
syms b
syms c
beta =  [ 0.5 * a * u; 0.5 * b * u; 0  ];
delta = [-0.5 * a    ; 0.5 * b    ; u/c];

% back substitute the x, y values

syms x
syms y
simplify(subs(subs(K, u, y/b + x/a), v, y/b - x/a))

So the answer 4a6b6c2(a4b4c2+4a4y2+4b4x2)2

Phew, finally!

No comments:

Post a Comment