Un cable sujeto por dos puntos de una superficie cilíndrica
En la figura, se muestra la forma de un cable suspendido de dos puntos A y B situados en una superficie cilíndrica de radio R.
Un punto P sobre una superficie cilíndrica de radio R tiene coordenadas (φ, z), donde φ es el ángulo que forma con el eje X la línea que une el origen con la proyección de dicho punto sobre el plano XY. z es la altura del punto P sobre el plano XY. Por ejemplo, las coordenadas de A(-π/6,0) y las de B(π/6,1).

Una curva sobre la superficie cilíndrica viene descrita por la función z=z(φ)
La longitud ds de un arco diferencial de curva situada en el punto φ,z
La longitud de un cable suspendido de los puntos A (-φ0, 0) y B (φ0, h) es
La energía potencial de un pequeño arco de longitud ds del cable localizado en el punto (φ,z) es dm·gz. Su masa es dm=λ·ds, λ la densidad lineal en kg/m
La energía potencial del cable es
Tenemos que determinar la ecuación de la curva z=z(φ) que haga que la energía potencial del cable sea mínima, con la condición de que el cable tenga una longitud fija L.
Un problema similar isoperimétrico estudiado en la página titulada 'Cálculo de variaciones'. Definimos una función F dependiente de un parámetro k
Aplicando la ecuación de Euler-Lagrange a la función F que no depende de φ
Esta es la ecuación diferencial de un cable suppendido de dos puntos de una superficie cilíndrica. Elevando al cuadrado despejamos dz/dφ
Integramos
Haciendo el cambio de variable u=z+k/(λg), la integral es inmediata
La ecuación de la catenaria depende de tres parámetros a, b y c que se determinan a partir de las coordenadas de los puntos A (-φ0, 0), B (φ0, h) y de la longitud L del cable
La longitud del cable es
-
Cable suspendido de los puntos A y B
Para obtener la última ecuación hay que expresar el seno hiperbólico como diferencia de dos exponenciales
Habiendo eliminado el parámetro c del sistema de dos ecuaciones, nos quedan dos ecuaciones con dos incógnitas a y b
Dividiendo
Elevando al cuadrado la segunda ecuación y empleando la relación cosh2x-sinh2x=1, obtenemos una ecuación transcendente en a
Calculamos el parámetro a, luego b y por último c, de una de las dos ecuaciones del sistema
Resumiendo, para determinar a, b y c empleamos las relaciones
Ejemplo
Sea un cable de longitud L=3 m, suspendido de los puntos A(-π/6,0) y B(π/6, 1) situados en una superficie cilíndrica de radio R=0.5 m
La mayor dificultad estriba en determinar el valor del parámetro a. Para tener una idea de dónde se encuentra la raíz de la primera ecuación, representamos las funciones
Una raíz es y=0. Para que exista una raíz y≠0, la pendiente de la recta ha de ser mayor que la pendiente de la función f(y) en el origen
2Rφ0 es la longitud del arco de circunferencia entre -φ0 y φ0
L=3; %longitud R=0.5; %radio phi_0=pi/6; %posición h=1; hold on f=@(y) 2*R*sinh(phi_0*y); g=@(y) sqrt(L^2-h^2)*y; fplot(f,[0,10]) fplot(g,[0,10]) raiz=fzero(@(y) f(y)-g(y), [0.1,10]); fprintf('El parámetro a=%1.4f\n',1/raiz) line([raiz,raiz],[0, f(raiz)],'lineStyle','--') hold off ylim([0,40]) grid on xlabel('y') ylabel('f(y), g(y)') title('Raíz de la ecuación transcendente')
Vemos que la intersección se produce en y≈7, (a≈1/7). Calculamos la raíz con la función
El parámetro a=0.1421
Conocido el parámetro a calculamos b y c y representamos la catenaria
L=3; %longitud R=0.5; %radio phi_0=pi/6; %posición de A h=1; %altura de B %parámetros a,b,c f=@(y) 2*R*sinh(phi_0*y)-sqrt(L^2-h^2)*y; a=1/fzero(f, [0.1,10]); b=a*log((L+h)/(L-h))/2; c=-a*R*cosh((-phi_0+b)/a); %cilindro r=R*ones(30,1); phi=linspace(0,2*pi,30); [r,phi]=meshgrid(r,phi); x=r.*cos(phi); y=r.*sin(phi); z=repmat(linspace(-1.1,1.1,30),30,1); hold on surfl(x,y,z); shading interp colormap(gray); %catenaria X=@(t) R*cos(t); Y=@(t) R*sin(t); Z=@(t) a*R*cosh((t+b)/a)+c; fplot3(X,Y,Z, [-phi_0,phi_0],'color','r') plot3(R*cos(-phi_0),R*sin(-phi_0),0, 'co','markersize',3,'markerfacecolor','c') plot3(R*cos(phi_0),R*sin(phi_0), h, 'co','markersize',3,'markerfacecolor','c') %ejes fplot3(@(t) R*cos(t), @(t) R*sin(t), @(t) t-1*t, [0,2*pi], 'color','k') line([0,R+0.25],[0,0],[0,0], 'color','k') line([0,0],[0,R+0.25],[0,0], 'color','k') line([0,0],[0,0],[-1,1], 'color','k') hold off axis equal xlabel('x'); ylabel('y'); zlabel('z') title('Superficie cilíndrica') view(70,40)
Catenaria simétrica
Cuando h=0, los puntos A y B están a la misma altura, la catenaria es simétrica y las fórmulas se reducen
Las ecuaciones paramétricas de la catenaria son
L=3; %longitud R=0.5; %radio phi_0=pi/6; %posición de A %parámetros a,c f=@(y) 2*R*sinh(phi_0*y)-L*y; a=1/fzero(f, [0.1,10]); c=-a*R*cosh(-phi_0/a); %cilindro r=R*ones(30,1); phi=linspace(0,2*pi,30); [r,phi]=meshgrid(r,phi); x=r.*cos(phi); y=r.*sin(phi); z=repmat(linspace(-1.5,0.1,30),30,1); hold on surfl(x,y,z); shading interp colormap(gray); %catenaria X=@(t) R*cos(t); Y=@(t) R*sin(t); Z=@(t) a*R*cosh(t/a)+c; fplot3(X,Y,Z, [-phi_0,phi_0],'color','r') plot3(R*cos(-phi_0),R*sin(-phi_0),0, 'co','markersize',3,'markerfacecolor','c') plot3(R*cos(phi_0),R*sin(phi_0), 0, 'co','markersize',3,'markerfacecolor','c') %ejes fplot3(@(t) R*cos(t), @(t) R*sin(t), @(t) t-1*t, [0,2*pi], 'color','k') line([0,R+0.25],[0,0],[0,0], 'color','k') %eje X line([0,0],[0,R+0.25],[0,0], 'color','k') %eje Y line([0,0],[0,0],[-1,0.5], 'color','k') %eje Z hold off axis equal xlabel('x'); ylabel('y'); zlabel('z') title('Superficie cilíndrica') view(70,40)
Referencias
E. G. M. de Lacerda, H. L. Carrión. A catenary-like cable confined in a circular cylinder. Revista Mexicana de Física, E 21 010211 1–10. January–June 2024