PDE toolbox and load on nodes

Question

haoran liu 2019-12-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/495889-pde-toolbox-and-load-on-nodes

评论： xl z 2024-6-12

1.jpg

I use the PDE Toolbox to analysis a discrete geometry, and I want to put the force vector on a point. But I can't define a point through the syntax "loadedVertex". Is there any other option to define the point? thank you!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Ravi Kumar 2019-12-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/495889-pde-toolbox-and-load-on-nodes#answer_405770

The function addVertex was introduced in R2019b. You might be running an older version of the MATLAB to get this error.

Regards,

Ravi

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

haoran liu 2019-12-12

Thankk you! But when I use the R2019b, I got another error. I can get the VertexID, But I can't put the force on the points.

请先登录，再进行评论。

Answer 2

Ravi Kumar 2019-12-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/495889-pde-toolbox-and-load-on-nodes#answer_406133

OK, in that case can you number the number of verteices in the model:

>> gm

Do you see the new VertexID in the plot you make?

Regards,

Ravi

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

xl z 2024-6-12

Addvertex can not add points on cylindrical surface. Is there a solution to this situation? Thank.

请先登录，再进行评论。

Answer 3

haoran liu 2019-12-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/495889-pde-toolbox-and-load-on-nodes#answer_406232

编辑：haoran liu 2019-12-13

Thank you for your answer. Yes, I can see the new VertexID in my model. But I can't solve it with force.

Regards,

haoran

6 个评论
显示 4更早的评论隐藏 4更早的评论

haoran liu 2019-12-13

在 MATLAB Online 中打开

%%
clear
clc
%%
model = createpde('structural','static-solid');
%% import geometry
gm = importGeometry(model,'sta_simplified_400um.stl');
pdegplot(gm, 'VertexLabels', 'on','FaceAlpha', 0.5)
view(30, 30);
title('KGT STA - Faces')
drawnow
%% mesh
hMax = 1.8; % [mm]
hMin = 0.8; % [mm]
tic
mesh = generateMesh(model, 'Hmax', hMax, 'Hmin', hMin);
elapsedTimeMesh = toc;
numOfElems = length(model.Mesh.Elements);
numOfNodes = length(model.Mesh.Nodes);
numOfRacewayNodes = length(findNodes(model.Mesh, 'region', 'Face', 6));
q = meshQuality(model.Mesh);
elemIDs = find(q < 0.3);
pdemesh(model.Mesh,'FaceAlpha',1)
hold on
pdemesh(model.Mesh.Nodes, model.Mesh.Elements(:, elemIDs), 'FaceColor', 'blue', 'EdgeColor', 'blue')
drawnow
%% Material Props
structuralProperties(model, 'Cell', 1, 'YoungsModulus', 200e9, 'PoissonsRatio', 0.3);
%% Boundary Conditions BC
structuralBC(model, 'Face', 15, 'Constraint', 'fixed');
structuralBC(model, 'Face', 2, 'Constraint', 'fixed');
%% load condition
VertexID = addVertex(gm,'Coordinates',[26.253327790164672,9.003630213068949,6.587610172703529]);
%%
pdegplot(gm, 'VertexLabels', 'on','FaceAlpha', 0.5)
view(30, 30);
title('KGT STA - Faces')
drawnow
%%
structuralBoundaryLoad(model,'Vertex',VertexID, 'Force', [10;10;10]);
%% solver
tic
try 
    result = solve(model);
    elapsedTimeSolver = toc;
    disp(['Solver finished in ' num2str(elapsedTimeSolver) ' [s]'])
   catch e
    elapsedTimeSolver = -1;
    disp(['Solver terminated with error "' e.message '", elapsed time: ' num2str(toc) ' [s]'])
   
    return;
end

Thank you! The model is in the attachment.

Ravi Kumar 2019-12-13

Thank. Now I see the issue. You generated mesh before adding vertex. Do all geometry operation before generating mesh. If you move the addVertex to right below importGeometry, you should not get this error.

Regards,

Ravi

haoran liu 2019-12-14

Thank you! It works

Regards,

haoran

请先登录，再进行评论。

PDE toolbox and load on nodes

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

更多回答（2 个）

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

6 个评论
显示 4更早的评论隐藏 4更早的评论

另请参阅

类别

标签

产品

Community Treasure Hunt

PDE toolbox and load on nodes

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

更多回答（2 个）

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

6 个评论 显示 4更早的评论隐藏 4更早的评论

另请参阅

类别

标签

产品

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

6 个评论
显示 4更早的评论隐藏 4更早的评论