solvepde(model)
Здравствуйте!
Помогите пожалуйста решить проблему!
При решении системы дифференциальных уравнений в частных производных
model = createpde(3);
syms u1(x,y) u2(x,y) u3(x,y) ;
pdeeq = [2*(r_radius*r_radius-x*x)*diff(u1,x,y)-x*diff(u1,y)-x*sqrt((r_radius*r_radius-x*x))*diff(u2,x,y)+x*x*diff(u2,y,y)-sqrt((r_radius*r_radius-x*x))*diff(u2,y)+4*(r_radius*r_radius-x*x)*diff(u3,x,x)+2*r_radius*r_radius*(1+nu_coef)*diff(u3,y,y)-2*x*diff(u3,x) + 2*prEh_coef*(1+nu_coef)/(nu_coef*len); sqrt((r_radius*r_radius-x*x))*diff(u1,y,y)+sqrt((r_radius*r_radius-x*x))*diff(u2,x,y)+2*sqrt((r_radius*r_radius-x*x))*diff(u3,x,y); 2*(r_radius*r_radius-x*x)/r_radius*diff(u1,x)-2*x*sqrt((r_radius*r_radius-x*x))/r_radius*diff(u2,x)+r_radius*nu_coef*diff(u3,y) - (R_radius-y)/len*prEh_coef-(3-nu_coef*nu_coef)*prEh_coef];
symCoeffs = pdeCoefficients(pdeeq,[u1 u2 u3],'Symbolic',true);
coeffs = pdeCoefficients(pdeeq,[u1 u2 u3]);
specifyCoefficients(model,'m',coeffs.m,'d',coeffs.d, ...
'c',coeffs.c,'a',coeffs.a,'f',@fcoeffunction);
rmix = [0 0 0];
qmix = [0 0 0; 0 0 0; 0 0 0];
gmix = [0 0 -nu_coef*prEh_coef_2];
applyBoundaryCondition(model,'neumann','Edge',1,'g',[0 prEh_coef_2*(1-nu_coef)/nu_coef 0],'q',[0 0 0; 0 0 0; 0 0 0]);
applyBoundaryCondition(model,'dirichlet','Edge',2,'r',[prEh_coef 0 0],'h',[1 0 0; 0 1 0; 0 0 0]);
applyBoundaryCondition(model,'mixed','Edge',3,'h',@myufun,'r',rmix,'g',gmix,'q',qmix);
applyBoundaryCondition(model,'dirichlet','Edge',4,'r',[-prEh_coef 0 0],'h',[1 0 0; 0 1 0; 0 0 0]);
res = solvepde(model);
возникает следующая ошибка:
Error using formGlobalMatrix
Coefficient evaluation function "makeCoefficient/coefficientFunction" returned a complex matrix but previously returned a real matrix. This function must consistently return the same type of matrix.(Пользуюсь пробной версией R2021b)
Подскажите, пожалуйста, что может быть причиной данной ошибки. Заранее благодарю!
Комментарии
Thank you for sharing this great post, I am very surviv io impressed with your post, the information given is very meticulous and understandable. I will regularly follow your next post.