How can I define individual degree of freedom constraints for 3D structural static PDE problems?

5 vues (au cours des 30 derniers jours)
In the below code, face 1 of the 3D block is fixed. I would like to change the structuralBC such that the constraint on face 1 only prevents motion in the X and Z directions. I can't come across any information on using constraints outside of the standard fixed, free, roller, etc. Is there a way to constrain individual degrees of freedom or is there a better way to approach this problem?
% Define material properties
E = 10.3e6; %Elastic Modulus - psi
nu = 0.33; %Poissons ratio - unitless
% Define loads
f1 = [0;-15;40]; %force in lbf
%Define geometry for use
body = 'Block.stl';
%Create static structureal model, import geomtery
smodel = createpde('structural','static-solid');
importGeometry(smodel,body);
%Plot with face labels
pdegplot(smodel,'FaceLabels','on','FaceAlpha',.5);
%Define properties and boundary conditions
structuralProperties(smodel,'YoungsModulus',E,'PoissonsRatio',nu);
structuralBC(smodel,'Face',1,'Constraint','fixed');
structuralBoundaryLoad(smodel,'Face',3, 'SurfaceTraction',f1);
%Mesh model
msh = generateMesh(smodel,'Hmax',5,'Hmin',1);
%Solve
R = solve(smodel);
%Plot results
pdeplot3D(smodel,'ColorMapData',R.VonMisesStress,'Deformation',R.Displacement);
figure()
pdeplot3D(smodel,'ColorMapData',R.Displacement.uz,'Deformation',R.Displacement);
figure()
pdeplot3D(smodel,'ColorMapData',R.Displacement.uy,'Deformation',R.Displacement);

Réponse acceptée

Ravi Kumar
Ravi Kumar le 11 Fév 2022
You can constrain X and Z displacement individually setting them to zero in one command like this:
structuralBC(smodel,'Face',1,'XDisplacement',0,'ZDisplacement',0,)
Regards,
Ravi
  1 commentaire
xl z
xl z le 25 Nov 2022
How to constrain the deformation of the cylindrical surface of a cylinder to not exceed a cylindrical surface? This case is similar to the restriction of the hole on the piston. The piston can only not deform radially to the outside of the hole. The "roller" constraint in the PDE toolbox seems to fix the radial displacement of the piston surface to 0. However, the actual situation is that the piston is able to deform toward the center.

Connectez-vous pour commenter.

Plus de réponses (0)

Produits


Version

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by