function [curvsin_,curvcon_,h,coordcurv_,coord_int]=curv_todo_punto(tipo_interp,matsin,matcon,coord,n_elem,nbarra,n_puntos_barra); puntos_barra=[1;n_puntos_barra]; [fil,col]=size(matsin); coord_int=[];sdf=[];auxiliar=0;coordcurv_=[]; a=0;b=0; if tipo_interp==1 for j=1:1:nbarra a=puntos_barra(j)+a;b=b+puntos_barra(j+1); for i=1:1:col modosin((j-1)*n_elem+j:j*(n_elem+1),i)=interp1(coord(a:b),... matsin(a:b,i),coord(a):(coord(b)-coord(a))/n_elem:coord(b),'cubic')'; modocon((j-1)*n_elem+j:j*(n_elem+1),i)=interp1(coord(a:b),... matcon(a:b,i),coord(a):(coord(b)-coord(a))/n_elem:coord(b),'cubic')'; end coord_inta=[coord(a):(coord(b)-coord(a))/n_elem:coord(b)]; coord_int=[coord_int;coord_inta']; auxiliar=coord_inta(length(coord_inta)); h(j)=coord_inta(2)-coord_inta(1); end for i=1:1:col z=1; for j=1:1:nbarra for k=((j-1)*n_elem+j+1):1:j*n_elem+(j-1) curvsin_(z,i)=(1/(h(j)^2))*(modosin(k-1,i)-2*modosin(k,i)+modosin(k+1,i)); curvcon_(z,i)=(1/(h(j)^2))*(modocon(k-1,i)-2*modocon(k,i)+modocon(k+1,i)); z=z+1; end end end p=1; for j=1:1:nbarra for k=((j-1)*n_elem+j+1):1:j*n_elem+(j-1) coordcurv_(p,1)=coord_int(k); p=p+1; end end end if tipo_interp==2 for j=1:1:nbarra a=puntos_barra(j)+a;b=b+puntos_barra(j+1); for i=1:1:col modosin((j-1)*n_elem+j:j*(n_elem+1),i)=interp1(coord(a:b),... matsin(a:b,i),coord(a):(coord(b)-coord(a))/n_elem:coord(b),'spline')'; modocon((j-1)*n_elem+j:j*(n_elem+1),i)=interp1(coord(a:b),... matcon(a:b,i),coord(a):(coord(b)-coord(a))/n_elem:coord(b),'spline')'; end coord_inta=[coord(a):(coord(b)-coord(a))/n_elem:coord(b)]; coord_int=[coord_int;coord_inta']; auxiliar=coord_inta(length(coord_inta)); h(j)=coord_inta(2)-coord_inta(1); end for i=1:1:col z=1; for j=1:1:nbarra for k=((j-1)*n_elem+j+1):1:j*n_elem+(j-1) curvsin_(z,i)=(1/(h(j)^2))*(modosin(k-1,i)-2*modosin(k,i)+modosin(k+1,i)); curvcon_(z,i)=(1/(h(j)^2))*(modocon(k-1,i)-2*modocon(k,i)+modocon(k+1,i)); z=z+1; end end end p=1; for j=1:1:nbarra for k=((j-1)*n_elem+j+1):1:j*n_elem+(j-1) coordcurv_(p,1)=coord_int(k); p=p+1; end end end