Ondas y campos vectoriales en un arco (Grupo 32)

Trabajo realizado por estudiantes
Título	Ondas y campos vectoriales en un arco. Grupo 32
Asignatura	Teoría de Campos
Curso	2025-26
Autores	Isabella Eugenia Acosta Montoya Macarena Gil Andrés Guillermo Polo Toledo Marta de la Quintana Zubiría Liam O'Hea Kith
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

1 Dibujo del Mallado

Mallado del arco realizado a través de Matlab:

Dibujo del mallado del arco mediante MatLab.

% Paso de muestreo
h = 0.1;

% Valor de u (radios)
u = 1:h:2;

% Calculamos cuántos puntos caben aproximadamente
puntos = round(pi/h) + 1;

% Vector v (ángulos) de 0 a pi
v = linspace(0, pi, puntos);

% Matrices de coordenadas polares
[rho, theta] = meshgrid(u, v);

% Conversión a Cartesianas
xx = rho .* cos(theta);
yy = rho .* sin(theta);


plot(xx, yy, 'r'); 

hold on % Mantenemos el gráfico para añadir las líneas radiales

plot(xx', yy', 'r');

hold off

% Ajustes finales visuales
axis equal
axis([-3, 3, -0.5, 2.5]);
title('Arco II');

4 Dibujo del campo de vectores en los puntos del mallado del sólido

Texto en negrita

Representamos el campo de vectores de deslizamiento en una mitad de anillo en dos dimensiones. Nos dan el campo [math]\ltmi\gtu\lt/mi\gt\ltmo\gt(\lt/mo\gt\ltmi\gtρ\lt/mi\gt\ltmo\gt,\lt/mo\gt\ltmi\gtθ\lt/mi\gt\ltmo\gt)\lt/mo\gt\ltmo\gt=\lt/mo\gt\ltmfrac\gt\ltmn\gt1\lt/mn\gt\ltmn\gt5\lt/mn\gt\lt/mfrac\gt\ltmfenced\gt\ltmrow\gt\ltmi\gtρ\lt/mi\gt\ltmo\gt-\lt/mo\gt\ltmn\gt1\lt/mn\gt\lt/mrow\gt\lt/mfenced\gt\ltmi\gtρ\lt/mi\gt\ltmsub\gt\ltmi\gte\lt/mi\gt\ltmi\gtρ\lt/mi\gt\lt/msub\gt[/math]

6 Divergencia de u

Se calculará la divergencia siguiendo la siguiente fórmula:

[math] \nabla ·\overrightarrow{u}= \frac{1}{\rho}\left( \frac{\partial }{\partial \rho}\left( \rho u_{\rho} \right) + \frac{\partial }{\partial \theta}\left( u_{\theta} \right) \right)=\frac{3\rho-2}{5} [/math]
La divergencia aumenta en función de [math] \rho [/math].

Cálculo y dibujo de la divergencia de [math]\overrightarrow{u}[/math] a través de matlab:

Dibujo de la divergencia en el arco mediante MatLab.

%creamos los intervalos de las variables
w=10; 
p=80;
u=linspace(1,2,w);
v=linspace(0,pi,p); 

%mallado
[U,V] = meshgrid(u,v);
X = U.*cos(V);
Y = U.*sin(V);

%divergencia
DIVu = (3.*U - 2)./5;

%gráfica
surf(X,Y,DIVu)
view(2)
axis equal;
colorbar;
title('Divergencia en el arco');
axis([-3, 3, -0.5, 2.5]);

7 Rotacional de u

Se calculará el rotacional siguiendo la siguiente fórmula:
[math] \nabla \times \overrightarrow{u}=\frac{1}{\rho} \begin{vmatrix} \overrightarrow{e_{\rho}} & \overrightarrow{e_{\theta}} & \overrightarrow{e_{z}}\\ \frac{\partial }{\partial \rho}& \frac{\partial }{\partial \theta}& \frac{\partial }{\partial z} \\ u_{\rho}& \rho u_{\theta} & u_{z} \end{vmatrix} [/math]

9 Tensiones tangenciales respecto al plano ortogonal a [math]\vec e_\rho[/math]

En este apartado se calcularán las tensiones tangenciales, especificamente con [math]\left | \sigma\cdot \vec e_{\rho}-(\vec e_{\rho}\cdot \sigma \cdot\vec e_{\rho}) \vec e_{\rho} \right |[/math].

Haciendo referencia a los datos adquiridos en el apartado anterior, comenzamos:

[math]|σ·\vec e_ρ-(\vec e_ρ·σ·\vec e_ρ)·\vec e_ρ| = \left | \begin{pmatrix} \ \frac{sen(θ)}{2\rho} & \frac{cos(θ)}{2\rho} & 0 \\ \frac{cos(θ)}{2\rho} & \frac{3sen(θ)}{2\rho} & 0 \\ 0 & 0 & \frac{sen(θ)}{2\rho} \end{pmatrix}·\begin{pmatrix} 1\\0\\0 \end{pmatrix} - \frac{sen(θ)}{2\rho}·\begin{pmatrix} 1\\0\\0 \end{pmatrix}\right | = \left | \begin{pmatrix} 0\\ \frac{cos(θ)}{2\rho}\\0 \end{pmatrix}\right| [/math]

10 Tensiones tangenciales respecto al plano ortogonal a [math]\frac{1}{\rho}[/math][math]\vec e_\theta[/math]

En este apartado se calcularán las tensiones tangenciales, especificamente con [math]\left | \sigma\cdot \vec e_{\rho}-(\vec e_{\theta}\cdot \sigma \cdot\vec e_{\theta}) \vec e_{\theta} \right |[/math].

11 Masa aproximando la integral numéricamente

La densidad de la placa viene dada por la función [math] d(\rho,\theta)=1 + e^{\rho^2 \cos\theta} [/math] A partir de este dato se va a calcular la masa aproximando la integral numéricamente. La coordenada [math] \rho [/math] conservará el radio del arco y [math] \theta [/math] es un ángulo que queda comprendido entre 0 y [math] \pi [/math], por tanto: [math] (\rho,\theta)=\left[ 1;2 \right]\times \left[ 0,\pi \right] [/math]

                                                                                                                                          Masa del arco calculada numéricamente con Matlab (Método Trapecio):

{{matlab|codigo= % MATLAB Code para calcular la masa % Densidad: d(rho, theta) = 1 + exp(rho^2 * cos(theta)) % Región: 1 <= rho <= 2, 0 <= theta <= pi % Parámetros y Límites de integración rho_min = 1; rho_max = 2; theta_min = 0; theta_max = pi; % Subintervalos N_rho = 50; N_theta = 100; % Tamaño del paso (h) para cada variable h_rho = (rho_max - rho_min) / N_rho; h_theta = (theta_max - theta_min) / N_theta; % Mallado de la superficie rho_vec = linspace(rho_min, rho_max, N_rho + 1); theta_vec = linspace(theta_min, theta_max, N_theta + 1); [RHO, THETA] = meshgrid(rho_vec, theta_vec); % 4. Funcion densidad F = RHO .* (1 + exp(RHO.^2 .* cos(THETA))); % Cálculo de los Pesos de la Regla del Trapecio % Pesos para la dirección rho w_rho = ones(N_rho + 1, 1); w_rho(1) = 0.5;