Warning: Explicit integral could not be found, not taking integral

1 visualización (últimos 30 días)
May u please find, where i am doing wrong? It is not working...
Warning: Explicit integral could not be found.
It is giving me this error message. It will find firstly Y then , with the result of Y it will find X value...This will be a loop. I wrote where the error is happening .U can see on below %%%%%%Should i take the integral numerically ?
clc;clear all
a =10;
b =10;
h =0.1;
E =210*10^9;
v =0.3;
q=100;
Dp = E*h^3/(12*(1-v^2));
%%%%------------
pi=3.1415926;
syms x
X1 = 1-cos(2*pi*x/a);
p=X1
for i=1:2
dXdxx=diff(p,x,2);
dXdxxxx=diff(p,x,4);
b1= simplify(int(p*dXdxxxx,0,a));
b2= 2*(int(p*dXdxx,0,a));
b3=int(p*p,0,a);
bv= int(p*(q/Dp),0,a)
k1=b1;
k2=b2;
k3=b3;
kv=bv;
syms y %d1 d2 d3 d4 h1 h2 h3 hv
denk1='k1*Y+k2*D2Y+k3*D4Y-kv=0'
bs1='Y(0)=0,DY(0)=0,Y(10)=0,Y(10)=0,DY(10)=0','y';
Y1=dsolve(denk1,bs1,'y');
m=Y1
dYdyy = diff(m,y,2);
dYdyyyy = diff(m,y,4);
%%%%error is happening in there. It is not taking integral of function. And also function is including k1,k2,k3... insteading of showing me the integers of k1 k2 k3... Because k1 k2 k3 they are already integers...But i am seeing symbolic .For example dYdyyyy , it is taking diferantial with symbolicly. When i copy and paste the result on a below line, it is giving me integer values of k1 k2 k3....in dYdyyyy(for example)
d1 = int(m*m,0,10);
d2 = 2*(int(m*dYdyy,0,b));
d3 = int(m*dYdyyyy,0,b);
dv = int(m*(q/Dp),0,b);
h1=d1;
h2=d2;
h3=d3;
hv=dv;
A=dsolve('h1*D4X+h2*D2X+h3*D4X-hv=0','X(0)=0','DX(0)=0','X(10)=0','DX(10)=0','x');
p=A
end
Xson=p
Yson=m
function w=fin(x,y)
w=p*m
end

Respuesta aceptada

Walter Roberson
Walter Roberson el 12 de Mzo. de 2014
denk1 = subs(sym('k1*Y+k2*D2Y+k3*D4Y-kv=0'));
  1 comentario
Ilke
Ilke el 13 de Mzo. de 2014
Thanks a lot.I wrote like that. Because i know b1, b2 , b3 , bv values. I put this values into the k1 k2 k3 kv. Now it is working :) denk1 = subs(sym('k1*Y+k2*D2Y+k3*D4Y-kv=0'),{k1,k2,k3,kv},{b1,b2,b3,bv})

Iniciar sesión para comentar.

Más respuestas (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by