PDE Matlab Dipolo

Contenido
Introducción.......................................................................................................................................2
Ejemplos...........................................................................................................................................3
Ejemplo 1......................................................................................................................................3
Código PDEModel Matlab 1.......................................................................................................6
Ejemplo 2......................................................................................................................................9
Código PDEModel Matlab 1.....................................................................................................12
Introducción
El electromagnetismo puede ser visto como el estudio de las interacciones entre
cargas en reposo y en movimiento. Así supone el análisis, la síntesis,
interpretación física y aplicación de campos eléctricos y magnéticos.
El electromagnetismo es la rama de la física o ingeniería eléctrica que estudia los
fenómenos eléctricos y magnéticos.
La electrostática es la rama de la física que analiza los efectos mutuos que se
producen entre los cuerpos como consecuencia de sus cargas eléctricas, es decir,
el estudio de las cargas eléctricas en equilibrio. La carga eléctrica es la propiedad
de la materia responsable de los fenómenos electrostáticos, cuyos efectos
aparecen en forma de atracciones y repulsiones entre los cuerpos que la poseen.
En el presente documento se exponen dos ejemplos de Matlab aplicados a la
electrostática y sus gráficos en 2D proporcionados por el mismo programa.
Ejemplos
A continuación de presentan los ejemplos de los programas con sus capturas de
pantalla.
Ejemplo 1
Se colocan dos conductores metálicos circulares sobre un papel secante
empapado en salmuera que sirve como un conductor plano y delgado. El modelo
físico para este problema consiste en la ecuación de Laplace
–∇ · (σ∇V) = 0
para el potencial eléctrico V y estas condiciones de contorno:
V = 1 en el conductor circular izquierdo
V = –1 en el conductor circular derecho
la condición de límite natural de Neumann en los límites exteriores
∂V/∂n = 0
La conductividad es σ = 1.
1.- Definición de espacio en Matlab, para creación de esferas y plano

pderect([-1.2 1.2 -0.6 0.6])
pdecirc(-0.6,0,0.3)
pdecirc(0.6,0,0.3)
Se especifica las condiciones de Dirichlet en las esferas con un peso de 1 y sus

respectivas cargas, positiva (izquierda 1) y negativa (derecha -1) y se procede a
realizar el mallado
Se resuelve el problema con la herramienta y se observa el plano con las cargas
respectivas (positiva y negativa en ambos lados)
Se especifica en la herramienta el contorno de las líneas de campo para que

aparezcan en el plano y se observa el desplazamiento
En este ejemplo se observó un dipolo con cargas equipotenciales en las esferas,
en la placa las condiciones de Neuman se consideraron en 0 por lo que las líneas
de campo son tangenciales a esta misma placa.
Las superficies equipotenciales son aquellas en las que el potencial toma un valor
constante. Por ejemplo, las superficies equipotenciales creadas por cargas
puntuales son esferas concéntricas centradas en la carga, como se deduce de la
definición de potencial (r = cte).
Las líneas de campo eléctrico son, en cada punto, perpendiculares a las
superficies equipotenciales y se dirigen hacia donde el potencial disminuye.
El trabajo para desplazar una carga entre dos puntos de una misma superficie
equipotencial es nulo.
Dos superficies equipotenciales no se pueden cortar.
Código PDEModel Matlab 1
El siguiente Código fue creado a través de la herramienta de PDEModel que
representa todas las figuras y valores asignados mostrados anteriormente, para
abrirlo se abre un nuevo PDEModel y se abre el archivo PDEModel.m (Que es el
nombre que se asignó al proyecto)
% This script is written and read by pdetool and should NOT be edited.
% There are two recommended alternatives:
% 1) Export the required variables from pdetool and create a MATLAB
script
% to perform operations on these.
% 2) Define the problem completely using a MATLAB script. See
% http://www.mathworks.com/help/pde/examples/index.html for examples
% of this approach.
function pdemodel
[pde_fig,ax]=pdeinit;
pdetool('appl_cb',8);
set(ax,'DataAspectRatio',[1.5 1 1]);
set(ax,'PlotBoxAspectRatio',[1 1 1]);
set(ax,'XLim',[-1.5 1.5]);
set(ax,'YLim',[-1 1]);
set(ax,'XTickMode','auto');
set(ax,'YTickMode','auto');
% Geometry description:
pderect([-1.2 1.2 -0.59999999999999998 0.59999999999999998],'R1');
pdecirc(-0.59999999999999998,0,0.29999999999999999,'C1');
pdecirc(0.59999999999999998,0,0.29999999999999999,'C2');
set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String','R1-
(C1+C2)')
% Boundary conditions:
pdetool('changemode',0)
% esfera derecha con carga negativa -1 y peso o valor w de 1 con
condición de Dirichlet
pdesetbd(12,...
'dir',...
1,...
'1',...
'-1')
pdesetbd(11,...
'dir',...
1,...
'1',...
'-1')
pdesetbd(10,...
'dir',...
1,...
'1',...
'-1')
pdesetbd(9,...
'dir',...
1,...
'1',...
'-1')
% esfera izquierda con carga positiva 1 y peso o valor w de 1 con
condición de Dirichlet
pdesetbd(8,...
'dir',...
1,...
'1',...
'1')
pdesetbd(7,...
'dir',...
1,...
'1',...
'1')
pdesetbd(6,...
'dir',...
1,...
'1',...
'1')
pdesetbd(5,...
'dir',...
1,...
'1',...
'1')
% Condición de Neuman en el plano exterior rectangular con valores de 0
pdesetbd(4,...
'neu',...
1,...
'0',...
'0')
pdesetbd(3,...
'neu',...
1,...
'0',...
'0')
pdesetbd(2,...
'neu',...
1,...
'0',...
'0')
pdesetbd(1,...
'neu',...
1,...
'0',...
'0')
% Mesh generation:
setappdata(pde_fig,'Hgrad',1.3);
setappdata(pde_fig,'refinemethod','regular');
setappdata(pde_fig,'jiggle',char('on','mean',''));
setappdata(pde_fig,'MesherVersion','preR2013a');
pdetool('initmesh')
pdetool('refine')
pdetool('jiggle')
% PDE coefficients:
pdeseteq(1,...
'1.0',...
'0.0',...
'0',...
'1.0',...
'0:10',...
'0.0',...
'0.0',...
'[0 100]')
setappdata(pde_fig,'currparam',...
['1.0';...
'0 '])
% Solve parameters:
setappdata(pde_fig,'solveparam',...
char('0','1000','10','pdeadworst',...
'0.5','longest','0','1E-4','','fixed','Inf'))
% Plotflags and user data strings:

setappdata(pde_fig,'plotflags',[1 1 1 1 2 1 1 1 0 0 0 1 1 0 0 1 0 1]);
setappdata(pde_fig,'colstring','');
setappdata(pde_fig,'arrowstring','');
setappdata(pde_fig,'deformstring','');
setappdata(pde_fig,'heightstring','');
% Solve PDE:
pdetool('solve')
Ejemplo 2
Se tienen 2 placas cargadas positiva y negativamente, y entre estas 2 placas se
coloca una esfera descargada para ver su comportamiento al estar entre las líneas
de campo de una placa cargada positivamente y una negativa.
1.- Definición de espacio de placas y esfera en Matlab

clc;
pderect([-1.2 1.2 -1 1])
pderect([-0.8 -0.3 -0.6 0.6])
pderect([0.8 0.3 -0.6 0.6])
pdecirc(0,0,0.2)
2.- Se especifican las condiciones de Dirichlet para las placas y la esfera

descargada, una positiva (4) y la otra negativa (-3)
a) Placa Izquierda (positiva)
b) Plaza Derecha (Negativa)

2.- Se realiza el mallado de para la representación del problema
Se resuelve el modelo por medio del software y se presenta el siguiente gráfico

donde se observa el resultado de las cargas en ambas placas y la esfera
descargada.
Como podemos ver la esfera se comporta como positiva de un lado y negativa de
otro, atrae las líneas de campo del lado izquierdo, donde está posicionada la placa
positiva, y a su vez las conduce hacia el lado derecho, que es donde esta
posicionada la placa de carga negativa, funcionando como un conductor de carga
entre las placas.
En las placas conductoras paralelas, las líneas del campo eléctrico son
perpendiculares a las placas y las líneas equipotenciales son paralelas a las
placas.
Podemos ver sobre y debajo de la esfera que las líneas de campo de una placa a
la otra permanecen perpendiculares a la placa, saliendo en línea recta de la placa
y llegar a la otra, y mientras más se acercan a la esfera, esta las atrae y
distorsiona su trayectoria al entrar a la esfera y salir de la misma.
Código PDEModel Matlab 1
El siguiente Código fue creado a través de la herramienta de PDEModel que
representa todas las figuras y valores asignados mostrados anteriormente, para
abrirlo se abre un nuevo PDEModel y se abre el archivo PlacasCargadas.m (Que
es el nombre que se asignó al proyecto)
% This script is written and read by pdetool and should NOT be edited.
% There are two recommended alternatives:
% 1) Export the required variables from pdetool and create a MATLAB
script
% to perform operations on these.
% 2) Define the problem completely using a MATLAB script. See
% http://www.mathworks.com/help/pde/examples/index.html for examples
% of this approach.
function pdemodel
[pde_fig,ax]=pdeinit;
pdetool('appl_cb',8);
set(ax,'DataAspectRatio',[1 1 1]);
set(ax,'PlotBoxAspectRatio',[2.0150114722054191 1 1]);
set(ax,'XLimMode','auto');
set(ax,'YLim',[-1 1]);
set(ax,'XTickMode','auto');
set(ax,'YTickMode','auto');
% Geometry description:
pderect([-1.2 1.2 -1 1],'R1');
pderect([-0.80000000000000004 -0.29999999999999999 -0.59999999999999998
0.59999999999999998],'R2');
pderect([0.80000000000000004 0.29999999999999999 -0.59999999999999998
0.59999999999999998],'R3');
pdecirc(0,0,0.20000000000000001,'C1');
set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String','R1-
(R2+R3+C1)')
% Boundary conditions:
pdetool('changemode',0)
pdesetbd(16,...
'dir',...
1,...
'1',...
'0')
pdesetbd(15,...
'dir',...
1,...
'1',...
'0')
pdesetbd(14,...
'dir',...
1,...
'1',...
'0')
pdesetbd(13,...
'dir',...
1,...
'1',...
'0')
pdesetbd(12,...
'dir',...
1,...
'1',...
'-3')
pdesetbd(11,...
'dir',...
1,...
'1',...
'4')
pdesetbd(10,...
'dir',...
1,...
'1',...
'-3')
pdesetbd(9,...
'dir',...
1,...
'1',...
'-3')
pdesetbd(8,...
'dir',...
1,...
'1',...
'4')
pdesetbd(7,...
'dir',...
1,...
'1',...
'4')
pdesetbd(6,...
'dir',...
1,...
'1',...
'4')
pdesetbd(5,...
'neu',...
1,...
'2',...
'1')
pdesetbd(4,...
'neu',...
1,...
'2',...
'1')
pdesetbd(3,...
'dir',...
1,...
'1',...
'-3')
pdesetbd(2,...
'neu',...
1,...
'2',...
'1')
pdesetbd(1,...
'neu',...
1,...
'2',...
'1')
% Mesh generation:
setappdata(pde_fig,'Hgrad',1.3);
setappdata(pde_fig,'refinemethod','regular');
setappdata(pde_fig,'jiggle',char('on','mean',''));
setappdata(pde_fig,'MesherVersion','preR2013a');
pdetool('initmesh')
pdetool('refine')
pdetool('jiggle')
% PDE coefficients:
pdeseteq(1,...
'1',...
'0.0',...
'1.0',...
'1.0',...
'0:10',...
'0.0',...
'0.0',...
'[0 100]')
setappdata(pde_fig,'currparam',...
['1 ';...
'1.0'])
% Solve parameters:
setappdata(pde_fig,'solveparam',...
char('0','1392','10','pdeadworst',...
'0.5','longest','0','1E-4','','fixed','Inf'))
% Plotflags and user data strings:

setappdata(pde_fig,'plotflags',[1 1 1 1 2 1 1 1 0 0 0 1 1 0 0 1 0 1]);
setappdata(pde_fig,'colstring','');
setappdata(pde_fig,'arrowstring','');
setappdata(pde_fig,'deformstring','');
setappdata(pde_fig,'heightstring','');
% Solve PDE:
pdetool('solve')

PDE Matlab Dipolo

Cargado por

Copyright:

Formatos disponibles

PDE Matlab Dipolo

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

PDE Matlab Dipolo

Cargado por

Copyright:

Formatos disponibles

Contenido

1.- Definición de espacio en Matlab, para creación de esferas y plano

Se especifica las condiciones de Dirichlet en las esferas con un peso de 1 y sus

Se especifica en la herramienta el contorno de las líneas de campo para que

% Plotflags and user data strings:

1.- Definición de espacio de placas y esfera en Matlab

2.- Se especifican las condiciones de Dirichlet para las placas y la esfera

b) Plaza Derecha (Negativa)

Se resuelve el modelo por medio del software y se presenta el siguiente gráfico

% Plotflags and user data strings:

También podría gustarte