Main Content

evaluateVonMisesStress

Evaluate von Mises stress for dynamic structural analysis problem

Description

example

vmStress = evaluateVonMisesStress(structuralresults) evaluates von Mises stress at nodal locations for all time- or frequency-steps.

Examples

collapse all

Evaluate the von Mises stress in a beam under a harmonic excitation.

Create a transient dynamic model for a 3-D problem.

structuralmodel = createpde("structural","transient-solid");

Create the geometry and include it in the model. Plot the geometry.

gm = multicuboid(0.06,0.005,0.01);
structuralmodel.Geometry = gm;
pdegplot(structuralmodel,"FaceLabels","on","FaceAlpha",0.5)
view(50,20)

Figure contains an axes object. The axes object contains 3 objects of type quiver, patch, line.

Specify Young's modulus, Poisson's ratio, and the mass density of the material.

structuralProperties(structuralmodel,"YoungsModulus",210E9, ...
                                     "PoissonsRatio",0.3, ...
                                     "MassDensity",7800);

Fix one end of the beam.

structuralBC(structuralmodel,"Face",5,"Constraint","fixed");

Apply a sinusoidal displacement along the y-direction on the end opposite the fixed end of the beam.

structuralBC(structuralmodel,"Face",3, ...
                             "YDisplacement",1E-4, ...
                             "Frequency",50);

Generate a mesh.

generateMesh(structuralmodel,"Hmax",0.01);

Specify the zero initial displacement and velocity.

structuralIC(structuralmodel,"Displacement",[0;0;0],"Velocity",[0;0;0]);

Solve the model.

tlist = 0:0.002:0.2;
structuralresults = solve(structuralmodel,tlist);

Evaluate the von Mises stress in the beam.

vmStress = evaluateVonMisesStress(structuralresults);

Plot the von Mises stress for the last time-step.

figure
pdeplot3D(structuralmodel,"ColorMapData",vmStress(:,end))
title("von Mises Stress in the Beam for the Last Time-Step")

Input Arguments

collapse all

Solution of a dynamic structural analysis problem, specified as a TransientStructuralResults or FrequencyStructuralResults object. Create structuralresults by using the solve function.

Example: structuralresults = solve(structuralmodel,tlist)

Output Arguments

collapse all

Von Mises Stress at the nodes, returned as a matrix. The rows of the matrix contain the values of von Mises stress at nodal locations, while the columns correspond to the time or frequency steps.

Version History

Introduced in R2018a