Unable to resolve the name results.NodalSolution.

5 次查看(过去 30 天)
I am using R2020b. While solving a PDE with myufunction for updating the boundary condition at each iteration, I keep getting the error message 'Unable to resolve the name results.NodalSolution'. Any help is welcome.
  3 个评论
Shreya George
Shreya George 2021-1-26
model=createpde;
gm=importGeometry(model,'tankdimension3d.stl');
pdegplot(gm,'FaceLabels','on','FaceAlpha',0.5)
phi=zeros(25,205);
for i=1:25
applyBoundaryCondition(model,'dirichlet','Face',6,'u',phi(i,:));
applyBoundaryCondition(model,'dirichlet','Face',1,'u',0);
applyBoundaryCondition(model,'neumann','Face',2:5,'g',0.00484);
myufunction = @(location,state)phi(i,:);
applyBoundaryCondition(model,'dirichlet','Face',6,...
'u',myufunction,...
'Vectorized','on');
specifyCoefficients(model,'m',0,'d',0,'c',1,'a',0,'f',0);
generateMesh(model,'Hmax',1);
result= solvepde(model);
u=results.NodalSolution;
k1=0.2*0.01505;
phi(i+1)=phi(i)+k1;
end
pdeplot3D(model,'ColorMapData',u)
Shreya George
Shreya George 2021-1-26
This was the code for the stl file
x=[-5
5
5
-5
-5
5
5
-5];
y=[-2
-2
2
2
-2
-2
2
2];
z=[-5
-5
-5
-5
0
0
0
0];
v = [x(:) y(:) z(:)];
V=v*1000;
faces=zeros(12,3);
faces(1,:)=[1, 2, 3];
faces(2,:)=[1, 3, 4];
faces(3,:)=[1, 2, 6];
faces(4,:)=[1, 6 ,5];
faces(5,:)=[2, 3, 7];
faces(6,:)=[2, 7, 6];
faces(7,:)=[3, 4, 8];
faces(8,:)=[3, 8, 7];
faces(9,:)=[4, 5, 8];
faces(10,:)=[4, 1, 5];
faces(11,:)=[5, 6 7];
faces(12,:)=[5, 7, 8];
T=[1 2 3
1 3 4
1 2 6
1 6 5
2 3 7
2 7 6
3 4 8
3 8 7
4 5 8
4 1 5
5 6 7
5 7 8];
TR=triangulation(T,v);
a = 1;
shp = alphaShape(x,y,z,a)
[faces, v] = boundaryFacets(shp);
fv.faces = faces;
fv.vertices = v;
trisurf(T,x,y,z);
stlwrite(TR,'tankdimension3d.stl')

请先登录,再进行评论。

采纳的回答

Ravi Kumar
Ravi Kumar 2021-1-26
You have a typo. It should be u=result.NodalSolution; not u=results.NodalSolution;
Regards,
Ravi

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Partial Differential Equation Toolbox 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by