Is it possible to calculate the reaction forces in a PDEModel if i use enforced displacement as load?

4 次查看（过去 30 天）
Baturay Yalvac
Baturay Yalvac 2018-12-11
评论： Baturay Yalvac 2019-1-13
Hello,
I am analysing cantilever-Beams using pdemodels. (I don't use structuralmodel) As BC i have to DC. First one is for the fixation of the beam and the second one applies enfroced displacement. I would like to calculate the reaction force, which would cause to this displacement. Is it possible to do this with pde Toolbox. Do you have any ideas?
Thank You

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

回答（1 个）

Ravi Kumar
Ravi Kumar 2018-12-11
You can evaluate reaction at the boundary using evaluateReaction function.
  4 个评论
显示 2更早的评论
Baturay Yalvac
Baturay Yalvac 2018-12-12
Hello,
In my project i am loading and fixing the 3D beams not with faces but with solids. Two cubes are in collision with the beam. One Cube-Beam collision gives the faces to fix the geometry and the other Cube-Beam collision gives the faces for enforced displacement or load conditions. I only could do that with using the coefficient like h u r g etc. Thats way i use the general model. I have alread calculated the Stresstensor like this:
....there is bit more code
if ~isempty(p)
[cgradx,cgrady,cgradz] = evaluateCGradient(result,p(:,1),p(:,2),p(:,3),[1 2 3]);
else
[cgradx,cgrady,cgradz] = evaluateCGradient(result);
end
sxx = cgradx(:,1);
sxy = cgradx(:,2);
syy = cgrady(:,2);
syz = cgrady(:,3);
szx = cgradz(:,1);
szz = cgradz(:,3);
stresstensor(:,1) = sxx;
stresstensor(:,2) = syy;
stresstensor(:,3) = szz;
stresstensor(:,4) = sxy;
stresstensor(:,5) = szx;
stresstensor(:,6) = syz;
And then i write a small function to do the nodal sum of the stresses
function [F] = reactionForce_BY(nodeIDs,stresstensor,width,thickness) ;
id_max = size(nodeIDs,2) ;
F = zeros(3,1) ;
Area = width*thickness ;
for ii = 1:id_max
nnumber = nodeIDs(1,ii) ;
F(1,1) = stresstensor(nnumber,1)+F(1,1);
F(2,1) = stresstensor(nnumber,2)+F(2,1);
F(3,1) =stresstensor(nnumber,3)+F(3,1);
end
F = F*(Area);
I want to learn the reaction force, that i need to realize the enforced displacement. So should i take the nodes only in the forcely displaced volume or the nodes in the whole beam? Or something else ?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Structural Mechanics 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by