Categoría discusión:TC22/23 - Historial de revisiones

Leyrecg en 18:55 9 dic 2022

2022-12-09T18:55:18Z

Leyrecg en 18:50 9 dic 2022

2022-12-09T18:50:36Z

Leyrecg en 18:49 9 dic 2022

2022-12-09T18:49:50Z

Leyrecg en 18:14 9 dic 2022

2022-12-09T18:14:56Z

Leyrecg: Sección nueva: /* CAMPOS ESCALARES Y VECTORIALES EN ELASTICIDAD, EN UN CUARTO DE CORONA CIRCULAR - GRUPO 3C */

2022-12-09T18:10:43Z

Sección nueva: ‎CAMPOS ESCALARES Y VECTORIALES EN ELASTICIDAD, EN UN CUARTO DE CORONA CIRCULAR - GRUPO 3C

Mostrar los cambios

Matias.robon: Página blanqueada

2022-12-08T17:26:41Z

Página blanqueada

Matias.robon: Sección nueva: /* Modelización del comportamiento de una placa bidimensional */

2022-12-08T17:26:01Z

Sección nueva: ‎Modelización del comportamiento de una placa bidimensional

Matias.robon: Página blanqueada

2022-12-08T17:20:50Z

Página blanqueada

Matias.robon: Sección nueva: /* Modelización del comportamiento de una placa bidimensional */

2022-12-08T17:19:44Z

Sección nueva: ‎Modelización del comportamiento de una placa bidimensional

Adriansanchez: Página blanqueada

2022-12-05T19:54:02Z

Página blanqueada

← Revisión anterior		Revisión del 18:55 9 dic 2022
Línea 1:		Línea 1:
		+	<inputbox>
		+	type=create
		+	width=100
		+	break=no
		+	buttonlabel=Create new article
		+	default=(== CAMPOS ESCALARES Y VECTORIALES EN ELASTICIDAD, EN UN CUARTO DE CORONA CIRCULAR - GRUPO 3C ==)
		+	</inputbox>
	== CAMPOS ESCALARES Y VECTORIALES EN ELASTICIDAD, EN UN CUARTO DE CORONA CIRCULAR - GRUPO 3C ==		== CAMPOS ESCALARES Y VECTORIALES EN ELASTICIDAD, EN UN CUARTO DE CORONA CIRCULAR - GRUPO 3C ==

← Revisión anterior		Revisión del 18:50 9 dic 2022
Línea 144:		Línea 144:
	[[Archivo:Apartado11Mallado.png\|500px\|left\|Código Matlab]]		[[Archivo:Apartado11Mallado.png\|500px\|left\|Código Matlab]]
	[[Archivo:FotoMallado11.png\|700px\|Mallado resultante]]		[[Archivo:FotoMallado11.png\|700px\|Mallado resultante]]
−
−
−
−	~~[[Categoría:Teoría de Campos]]~~
−	~~[[Categoría:TC22/23]]~~

← Revisión anterior		Revisión del 18:49 9 dic 2022
Línea 149:		Línea 149:
	[[Categoría:Teoría de Campos]]		[[Categoría:Teoría de Campos]]
	[[Categoría:TC22/23]]		[[Categoría:TC22/23]]
−
−	~~== CAMPOS ESCALARES Y VECTORIALES EN ELASTICIDAD, EN UN CUARTO DE CORONA CIRCULAR - GRUPO 3C ==~~
−
−	{{ TrabajoED \| Deformaciones de una placa plana. Grupo 3-C \| [[:Categoría: Teoría de Campos\|Teoría de Campos]]\|[[:Categoría:TC22/23\|2022-23]] \| Leyre Castañeda Gato, Natalia Cano Martín, Miguel Ángel De Gregorio Ávila, Juan Marquez Alba, Jaime San Vicente Lara}}
−	~~'''CAMPOS ESCALARES Y VECTORIALES EN ELASTICIDAD, EN UN CUARTO DE CORONA CIRCULAR - GRUPO 3C'''~~
−
−	Para este artículo, nos disponemos a analizar y a visualizar, los distintos efectos que fijan los campos escalares y vectoriales sobre un sólido. Para ello, consideraremos una placa plana que ocupa un cuarto de anillo circular, centrado en el origen, comprendido entre los radios 1 y 2, estando en el plano y≥x.
−	~~Supondremos que tenemos dos cantidades físicas.~~
−	Una de ellas, la temperatura <math> T(x,y) = log ((x-2)^2 + y^2 + 1) </math>, y la otra serán los desplazamientos, definidos por: <math> \vec u(ρ,θ) = (log(ρ)/2) · sen(2θ- \Pi/2) </math>, siendo (x,y) la posición de cada punto, y <math> \vec r_0(ρ,θ) </math>, el vector de posición previo a la deformación.
−
−
−	~~== '''''MALLADO DE LOS PUNTOS INTERIORES DEL SÓLIDO''''' ==~~
−
−	Para poder representar de forma correcta el mallado de la placa plana, vamos a proceder a parametrizar la superficie usando coordenadas cilíndricas. Nos apoyaremos en Matlab para este, y posteriores visualizaciones de gráficas.
−	~~Elegimos un paso de muestreo de 2/10 o 0.2 para x e y, y representaremos la placa en un rectángulo de [-3,3]x[-1,3]. Insertamos también una imagen del mallado que resulta.~~
−
−	~~[[Archivo:Apartado1Mallado1.png\|500px\|left\|Código Matlab]]~~
−	~~[[Archivo:FotoMallado1.png\|500px\|Mallado resultante]]~~
−
−	~~== '''''CURVAS DE NIVEL DE LA TEMPERATURA''''' ==~~
−
−	A continuación, veremos las curvas de nivel del campo escalar de temperaturas, dado por la función T definida al inicio, que serán representadas, usando el comando "contour", y veremos sus máximos. Como vemos, usaremos una gama cromática de colores que nos indican las zonas cálidas (amarillos) y las zonas más frías (azules). Como podemos ver en la imagen, disponemos de dos gráficas. La situada a la izquierda representa el campo escalar, mientras que la situada a la derecha, representa las curvas del nivel del campo.
−
−	~~[[Archivo:Apartado2Mallado.png\|600px\|left\|Código Matlab]]~~
−	~~[[Archivo:FotoMallado2.png\|600px\|Mallado resultante]]~~
−
−
−	~~== '''''GRADIENTE DE LA TEMPERATURA''''' ==~~
−
−	~~En este apartado, calcularemos el gradiente que representa la variación de temperatura. Derivaremos la temperatura T(x,y) respecto de cada variable~~
−	~~Hemos de apreciar que el gradiente de T, ∇(T), es ortogonal a las curvas de nivel.~~
−
−	~~[[Archivo:Apartado3Mallado.png\|600px\|left\|Código Matlab]]~~
−	~~[[Archivo:FotoMallado3.png\|600px\|Mallado resultante]]~~
−
−
−	~~== '''''CAMPO DE VECTORES''''' ==~~
−
−	~~Mostraremos, en este apartado, el campo de vectores en los puntos del mallado de la placa.~~
−	~~Aunque nos piden los resultados en t=0, nuestra función no depende de t. Mostraremos igualmente el campo de desplazamientos.~~
−
−	~~[[Archivo:Apartado4Mallado.png\|600px\|left\|Código Matlab]]~~
−	~~[[Archivo:FotoMallado4.png\|600px\|Mallado resultante]]~~
−
−
−	~~== '''''SOLIDO ANTES Y DESPUÉS DEL DESPLAZAMIENTO''''' ==~~
−
−	El desplazamiento del sólido dado por el campo de vectores <math> \vec u </math> . Lo representamos antes y después de que se produzca, así como una comparación entre ambos usando el comando "subplot".
−	Comparando las gráficas, podemos comprobar que la zona central superior de la placa es la que más se desplaza y por el contrario, en la parte inferior por los laterales de la placa el desplazamiento es mínimo.
−
−
−	~~[[Archivo:Apartado5desplazamiento.png\|600px\|left]]~~
−	~~[[Archivo:foto5desplazamientos.jpeg\|600px]]~~
−
−
−	~~== '''''DIVERGENCIA DEL CAMPO''''' ==~~
−
−	~~Primero calculamos analíticamente la divergencia del campo ū (en t = 0), cuando es máxima, mínima y nula.~~
−	~~La divergencia es una medida del cambio de volumen local debido al desplazamiento y lo podemos observar en esta gráfica:~~
−
−	~~[[Archivo:Apartado6Mallado.png\|600px\|left]]~~
−	~~[[Archivo:FotoMallado6.png\|600px]]~~
−
−	~~Además, aportamos los datos que nos arroja Matlab en cuanto al máximo y mínimo de la divergencia, siendo el máximo = 0.3444 y el mínimo = -0.2984.~~
−
−	~~== '''''ROTACIONAL DEL CAMPO''''' ==~~
−	~~Calculo del módulo del rotacional del campo ū (en t = 0).~~
−	~~Aplicamos la fórmula y tenemos que el rotacional es~~
−	~~<math> \mathbb{\nabla \times \vec u} = {1 \over ρ} \;~~
−	~~\begin{vmatrix}~~
−	~~{\vec e_ρ} & {ρ· \vec e_θ} & {\vec e_z}\\~~
−	~~{\partial \over \partial ρ} & {\partial \over \partial θ} & {\partial \over \partial z}\\~~
−	~~{-log(ρ)*cos(2θ) \over 2} & {0} & {0}\\~~
−	~~\end{vmatrix} </math> = <math> -{1 \over ρ}·{\partial \over \partial θ}[{-log(ρ)*cos(2θ) \over 2}]·{\vec e_z} </math>~~
−
−	~~Siendo el resultado final = <math> {-log(ρ)sin(2θ) \over ρ}·{\vec e_z} </math>, y si hacemos el módulo, nos queda, <math> \|{\nabla \times \vec u}\| = {log(ρ)sin(2θ) \over ρ} </math>~~
−
−	~~== '''''ECUACIÓN DEL TENSOR DE TENSIONES σ CON LOS COEFICIENTES DE LAMÉ''''' ==~~
−	~~Tras haber realizado los cálculos analíticos, y suponiendo que los coeficientes de Lamé, "λ" y "μ", son iguales a 1, hemos llegado a la matriz resultante:~~
−
−	~~<math> \mathbb{σ} = \;~~
−	~~\begin{pmatrix}~~
−	~~{-cos(2θ)(ln(ρ)+3) \over 2ρ} & {ln(ρ)sen(2θ) \over ρ} & {0}\\~~
−	~~{ln(ρ)sen(2θ) \over ρ} & {-cos(2θ)(3*ln(ρ)+2) \over 2ρ} & {0}\\~~
−	~~{0} & {0} & {-cos(2θ)*(ln(ρ)+1) \over 2ρ} \\~~
−	~~\end{pmatrix} </math>~~
−
−	Una vez obtenida dicha matriz σ, hemos de calcular las tensiones normales en la dirección <math> \vec i, \vec j, \vec k, </math>, coincidiendo los resultados de elementos de la diagonal principal con los calculados.
−	Es decir, ''<math> \vec i · σ · \vec i </math>'' coincide con el elemento σ(1,1) de la matriz; así como ''<math> \vec j · σ · \vec j </math>'' coincide con el elemento σ(2,2) de la matriz y ''<math> \vec k · σ · \vec k </math>'' coincide con el elemento σ(3,3) de la matriz.
−	Mostramos a continuación el código de Matlab, así como las 3 gráficas, habiendo usado el comando subplot. Hemos añadido, como apartado extra, la dirección de las tensiones según los ejes, pudiendo apreciarse en la 3 imagen de este apartado.
−
−	~~[[Archivo:Apartado8Mallado_.png\|300px\|left]]~~
−	~~[[Archivo:FotoMallado8.png\|900px]]~~
−	~~[[Archivo:FotoMallado8_.png\|900px]]~~
−
−	~~== '''''TENSIONES TANGENCIALES RESPECTO AL PLANO ORTOGONAL''''' ==~~
−	Esta vez, calcularemos las tensiones tangenciales respecto al plano ortogonal <math> \vec i </math>, siendo la fórmula a aplicar: <math> \ \left\| \sigma ·\vec i -(\vec i·\sigma·\vec i)\vec i \right\| </math>.
−	~~Aunque nos piden mostrar el valor de t=0, como ya hemos comentado anteriormente, nuestra función no depende de t.~~
−	~~Vemos primero las fórmulas, y después, la codificicación en Matlab, seguido de las gráficas.~~
−
−	~~<math> \ \left\| \sigma ·\vec i -(\vec i·\sigma·\vec i)\vec i \right\| </math> =~~
−	~~<math> \mathbb\; \begin{pmatrix}~~
−	~~{-cos(2θ)(ln(ρ)+3) \over 2ρ} & {ln(ρ)sen(2θ) \over ρ} & {0}\\~~
−	~~{ln(ρ)sen(2θ) \over ρ} & {-cos(2θ)(3*ln(ρ)+2) \over 2ρ} & {0}\\~~
−	~~{0} & {0} & {-cos(2θ)*(ln(ρ)+1) \over 2ρ} \\~~
−	~~\end{pmatrix}~~
−	~~· \begin{pmatrix}~~
−	~~1 \\~~
−	~~0 \\~~
−	~~0 \\~~
−	~~\end{pmatrix} - (\frac{-cos(2θ)*(ln(ρ)+3)}{2ρ}) · \begin{pmatrix}~~
−	~~1 \\~~
−	~~0 \\~~
−	~~0 \\~~
−	~~\end{pmatrix} = \begin{pmatrix}~~
−	~~0 \\~~
−	~~ln(ρ)*sin(2θ) \over ρ \\~~
−	~~0 \\~~
−	~~\end{pmatrix} </math>~~
−
−	~~[[Archivo:Apartado9Mallado_.png\|500px\|left]]~~
−	~~[[Archivo:FotoMallado9_.png\|500px]]~~
−
−	~~== '''''TENSIÓN DE VON MISES''''' ==~~
−
−	~~La tensión de Von Mises se define por la siguiente formula:~~
−	~~[[Archivo:Tension_VonMises.PNG \|400px\|miniatura\|centro\|Fórmula de Von Mises]]~~
−
−	donde \(\ σ_{1}\), \(\ σ_{2}\), \(\ σ_{3}\) son los autovalores de σ, también conocidos como tensores principales. La tensión de Von Mises es una magnitud escalar que es un indicador para saber cuando un material inicia un comportamiento plástico y no elástico puro.
−
−	~~Vamos a proceder a representar dicha tensión con la ayuda del software MatLab.~~
−	~~Podremos observar gracias a la imagen que se adjunta a continuación que en el punto (0,2) se alcanza el mayor valor de la tensión de Von Mises.~~
−
−	~~[[Archivo:Apartado10Mallado10.png\|500px\|left]]~~
−	~~[[Archivo:FotoMallado10.png\|500px]]~~
−
−	~~== '''''CAMPO DE FUERZAS QUE ACTÚA SOBRE LA PLACA''''' ==~~
−
−	~~Vamos a proceder a calcular <math> \vec F </math> usando la ecuación de Lamé equivalente:~~
−	~~[[Archivo:11.21.PNG \|400px\|miniaturadeimagen\|centro\|Ecuación de Lamé]]~~
−
−	Donde <math> λ </math> y <math> μ </math> son conocidos como los coeficientes de Lamé que dependen de las propiedades del material. En este caso vamos a considerar como 1 ambas constantes. En la grafica siguiente podemos apreciar como actúa la fuerza sobre el solido.
−
−
−	~~[[Archivo:Apartado11Mallado.png\|500px\|left\|Código Matlab]]~~
−	~~[[Archivo:FotoMallado11.png\|700px\|Mallado resultante]]~~

← Revisión anterior		Revisión del 18:14 9 dic 2022
Línea 295:		Línea 295:
	[[Archivo:Apartado11Mallado.png\|500px\|left\|Código Matlab]]		[[Archivo:Apartado11Mallado.png\|500px\|left\|Código Matlab]]
	[[Archivo:FotoMallado11.png\|700px\|Mallado resultante]]		[[Archivo:FotoMallado11.png\|700px\|Mallado resultante]]
−
−
−
−	~~[[Categoría:Teoría de Campos]]~~
−	~~[[Categoría:TC22/23]]~~

← Revisión anterior		Revisión del 17:26 8 dic 2022
Línea 1:		Línea 1:
−	~~== Modelización del comportamiento de una placa bidimensional ==~~

−	{{ TrabajoED \| Modelización del comportamiento de una placa bidimensional. Grupo 1-C \| [[:Categoría:Teoría de Campos\|Teoría de Campos]]\|[[:Categoría:TC22/23\|2022-23]] \| Mario Andres Silva Peña Loza, Matias Rodriguez Obon, Javier Nogales Sanchez, Sofia Benito Gomez, Miguel Angel Abad Robles }}
−
−	~~[[Categoría:Teoría de Campos]]~~
−	~~[[Categoría:TC22/23]]~~

← Revisión anterior		Revisión del 17:20 8 dic 2022
Línea 1:		Línea 1:
−	~~== Modelización del comportamiento de una placa bidimensional ==~~

−	{{ TrabajoED \| Modelización del comportamiento de una placa bidimensional. Grupo 1-C \| [[:Categoría:Teoría de Campos\|Teoría de Campos]]\|[[:Categoría:TC22/23\|2022-23]] \| Mario Andres Silva Peña Loza, Matias Rodriguez Obon, Javier Nogales Sanchez, Sofia Benito Gomez, Miguel Angel Abad Robles }}
−
−	~~[[Categoría:Teoría de Campos]]~~
−	~~[[Categoría:TC22/23]]~~

← Revisión anterior	Revisión del 17:19 8 dic 2022
Línea 1:	Línea 1:
	+	== Modelización del comportamiento de una placa bidimensional ==

	+	{{ TrabajoED \| Modelización del comportamiento de una placa bidimensional. Grupo 1-C \| [[:Categoría:Teoría de Campos\|Teoría de Campos]]\|[[:Categoría:TC22/23\|2022-23]] \| Mario Andres Silva Peña Loza, Matias Rodriguez Obon, Javier Nogales Sanchez, Sofia Benito Gomez, Miguel Angel Abad Robles }}
	+
	+	[[Categoría:Teoría de Campos]]
	+	[[Categoría:TC22/23]]

← Revisión anterior		Revisión del 19:54 5 dic 2022
Línea 1:		Línea 1:
−	{{ TrabajoED \| Flujo de Couette. Grupo 7C \| [[:Categoría:Teoría de Campos\|Teoría de Campos]]\|[[:Categoría:TC22/23\|2022-23]] \| Adrián Sánchez Gregorio, Marta Isabel Trigueros Díaz de Villafranca, Guillermo Navarro Moya, Lucía Hernando Molinero}}</nowiki></pre>

−
−	~~==INTRODUCCIÓN==~~
−
−	Este artículo está realizado por alumnos de la asignatura Teoría de Campos de Ingeniería Civil y Territorial con el objetivo de estudiar distintos campos escalares y vectoriales, así como su visualización en Matlab.
−
−	~~Consideraremos un fluido incompresible a través de dos planos paralelos y horizontales, en lo que podríamos asemejar a un canal, de tal forma que el plan superior se mueva con velocidad constante.~~
−
−	~~Definimos fluido incompresible como aquel que se opone a ser comprimido y cuya densidad es constante a lo largo del tiempo.~~
−
−
−	~~==MALLADO==~~
−
−	~~Nuestra superficie de trabajo está contenida en el rectángulo <math>(y, z) \in[0,8] \times[0,1]</math>~~
−
−
−	~~{{matlab\|codigo=~~
−	~~y=0:0.1:8;~~
−	~~z=0:0.1:1; %intervalos definidos en enunciado~~
−	~~[ey,ez]=meshgrid(y,z);~~
−	~~figure(1)~~
−	~~[ey,ez]=meshgrid([0:0.1:8],[0:0.1:1]);~~
−	~~mesh(ey,ez,0*ey);~~
−	~~grid on;~~
−	~~axis([0,8,-1,2]);~~
−	~~view(2)}}~~
−
−	~~[[Archivo:Mallado 7 2223.png\|miniaturadeimagen\|centro]]~~