[Noti-MNS] Re: otra consulta

[Mario Storti]

>>>>> On Fri, 2 Dec 2005 10:29:13 -0300, 
>>>>>      gustavo bellino <gushbellino en gmail.com> said:

> Mario estamos teniendo problemas con el ejercicio 3 de la guia dos.
> Por un lado nos tira unos warnings:

> warning: matrix singular to machine precision, rcond = 0
> warning: near line 59, column 3:

>>>> q = A \ b;


> warning: attempting to find minimum norm solution
> warning: near line 59, column 3:

>>>> q = A \ b;

> Y por otro lado la solucion que nos da no es logica.
> Te agradeceriamos si nos podes dar una ayuda o por lo menos fijarte si
> estamos bien encaminados.
> A la viga la dividimos en 4 elementos con 6 nodos.

> --
> Atentamente
> Gustavo E. Bellino

Estan poniendo mal las fijaciones, al poner el elemento diagonal en
1. 

Tambien guarda que los elmentos estaban numerados en sentido horario,
en vez de antihorario (que es lo usual).  Eso no trae error en este
caso, ya que solo hace que todos los volumenes sean negativos y
equivale a multiplicar todo el sistema de ecuaciones por -1. Guarda
que si seria error si algunos los numeraran en sentido horario y otros
antihorario.

Al final estan usando `mesh'. Aqui me parece que no tiene demasiado
sentido, no?

Les mando el parche. 

Mario

diff -c /u/mstorti/TEX/METNUMSIM/g2ej3.m.\~1\~ /u/mstorti/TEX/METNUMSIM/g2ej3.m
*** /u/mstorti/TEX/METNUMSIM/g2ej3.m.~1~	Fri Dec  2 12:09:11 2005
--- /u/mstorti/TEX/METNUMSIM/g2ej3.m	Fri Dec  2 12:16:41 2005
***************
*** 1,6 ****
  
  		x=[0 0; 0 0.1; 0.5 0.1; 0.5 0; 1 0.1; 1 0]; %colocamos las coordenadas de 6 nodos
! 	icone=[1 2 3; 3 4 1; 4 3 5; 4 5 6];	%coordenadas de los elementos
  	fijacion = [1 1 0; 1 2 0; 2 1 0; 2 2 0];	
  	%fijacion se entiende asi: nodo campo valor, ejemplo nodo 1 del campo 2 (que es v, campo 1 es u) y desplazamiento         	%cero.
  	
--- 1,6 ----
  
  		x=[0 0; 0 0.1; 0.5 0.1; 0.5 0; 1 0.1; 1 0]; %colocamos las coordenadas de 6 nodos
! 		icone=[1 3 2; 3 1 4; 4 5 3; 4 6 5];	%coordenadas de los elementos
  	fijacion = [1 1 0; 1 2 0; 2 1 0; 2 2 0];	
  	%fijacion se entiende asi: nodo campo valor, ejemplo nodo 1 del campo 2 (que es v, campo 1 es u) y desplazamiento         	%cero.
  	
***************
*** 50,64 ****
  	%"condicion de contorno"
  	%fijacion = [1 1 0; 1 2 0; 2 1 0; 2 2 0];
  	for i=1:nc
! 		A(fijacion(i,1)*2-2+fijacion(i,2),:)=0;				
! 		A(fijacion(i,1)*2-2+fijacion(i,2),fijacion(i,1))=1;	
! 		b(fijacion(i,1)*2-2+fijacion(i,2))=fijacion(i,3);		
  	end
  	A;
  	b;
  	q=A\b;
  	
! 	mesh(reshape(q(:),3,4));
  	
  			
  
--- 50,65 ----
  	%"condicion de contorno"
  	%fijacion = [1 1 0; 1 2 0; 2 1 0; 2 2 0];
  	for i=1:nc
! 	  eq = fijacion(i,1)*2-2+fijacion(i,2);
! 	  A(eq,:)=0;				
! 	  A(eq,eq)=1;	
! 	  b(eq)=fijacion(i,3);		
  	end
  	A;
  	b;
  	q=A\b;
  	
! #	mesh(reshape(q(:),3,4));

-- 
-------------------------
Mario Alberto Storti     [cel. +54-342-156144983]
CIMEC (INTEC/CONICET-UNL), Guemes 3450 - 3000 Santa Fe, Argentina
Tel: +54-342-4511594 (ext 1015), Tel/Fax: +54-342-4511169
e-mail: mstorti at intec dot unl dot edu dot ar
http://www.cimec.org.ar/mstorti
-------------------------