Contents     Functions         Previous Next     PDF Index



Computation of stresses and energies for given deformations.


Result = fe_stress('Command',MODEL,DEF)
  ...  = fe_stress('Command',node,elt,pl,il, ...)
  ...  = fe_stress( ... ,mode,mdof)


You can display stresses and energies directly using fecom ColorDataEner commands and use fe_stress to analyze results numerically. MODEL can be specified by four input arguments node, elt, pl and il (those used by fe_mk, see also section 7.1 and following), a data structure with fields .Node, .Elt, .pl, .il, or a database wrapper with those fields.

The deformations DEF can be specified using two arguments: mode and associated DOF definition vector mdof or a structure array with fields .def and .DOF.

Ener [m,k]ElementSelection

Element energy computation. For a given shape, the levels of strain and kinetic energy in different elements give an indication of how much influence the modification of the element properties may have on the global system response. This knowledge is a useful analysis tool to determine regions that may need to be updated in a FE model. Accepted command options are

Obsolete options are

The strain and kinetic energies of an element are defined by

φTKelementφ   and     Ekinetice=

For complex frequency responses, one integrates the response over one cycle, which corresponds to summing the energies of the real and imaginary parts and using a factor 1/4 rather than 1/2.


feplot allows the visualization of these energies using a color coding. You should compute energies once, then select how it is displayed. Energy computation clearly require material and element properties to be defined in InitModel.

The earlier high level commands fecom ColorDataK or ColorDataM don't store the result and thus tend to lead to the need to recompute energies multiple times. The preferred strategy is illustrated below.

% Computing, storing and displaying energy data
 demosdt('LoadGartFe'); % load model,def 
 cf=feplot(model,def);cf.sel='eltname quad4';fecom ch7
 % Compute energy and store in Stack
 Ek=fe_stress('ener -MatDes 1 -curve',model,def)
 % Color is energy density by element
 feplot('ColorDataElt  -dens -ColorBarTitle "Ener Dens"',Ek);
 % Color by group of elements
 cf.sel={'eltname quad4', ... % Just the plates
   'ColorDataElt -ColorBarTitle "ener" -bygroup -edgealpha .1', ...
   Ek}; % Data with no need to recompute
 fecom(cf,'ColorScale One Off Tight') % Default color scaling for energies

Accepted ColorDataElt options are

The color animation mode is set to ScaleColorOne.


out=fe_stress('stress CritFcn Options',MODEL,DEF,EltSel) returns the stresses evaluated at elements of Model selected by EltSel.

The CritFcn part of the command string is used to select a criterion. Currently supported criteria are

sI, sII, sIIIprincipal stresses from max to min. sI is the default.
misesReturns the von Mises stress (note that the plane strain case is not currently handled consistently).
-comp iReturns the stress components of index i. This component index is giving in the engineering rather than tensor notation (before applying the TensorTopology transformation).

Supported command Options (to select a restitution method, ...) are

The fecom ColorDataStress directly calls fe_stress and displays the result. For example, run the basic element test q4p testsurstress, then display various stresses using

% Using stress display commands
 fecom('ColorDataStress atcenter')
 fecom('ColorDataStress mises')
 fecom('ColorDataStress sII atcenter')

To obtain strain computations, use the strain material as shown below.

% Accessing stress computation data (older calls)
 [model,def]=hexa8('testload stress');'dbval 100 strain','dbval 112 strain');'dbval 111 d3 -3');
 data=fe_stress('stress atcenter',model,def)


For stress processing, one must often distinguish the raw stress components associated with the element formulation and the desired output. CritFcn are callback functions that take a local variable r1 of dimensions (stress components × nodes × deformations) and to replace this variable with the desired stress quantity(ies). For example

% Sample declaration of a user defined stress criterium computation
 function out=first_comp(r1)

would be a function taking the first component of a computed stress. sdtweb fe_stress(''Principal''' provides stress evaluations classical for mechanics.

For example, a list of predefined CritFcn callback :

Redefining the CritFcn callback is in particular used in the StressCut functionality, see section 4.7.

See also

fe_mk, feplot, fecom

©1991-2019 by SDTools
Previous Up Next