findStructuralIC

(To be removed) Find initial displacement and velocity assigned to geometric region

findStructuralIC will be removed. Use the CellIC, FaceIC, EdgeIC, and VertexIC properties of femodel instead. (since R2023a) For more information on updating your code, see Version History.

Syntax

sica = findStructuralIC(initialConditions,RegionType,RegionID)

Description

sica = findStructuralIC(initialConditions,RegionType,RegionID) returns the initial displacement and velocity assigned to the specified region.

Examples

Find Initial Conditions for Cells of 3-D Geometry

Find the initial displacement and velocity assigned to the cells of a 3-D geometry.

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

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

Create the geometry consisting of the three nested cylinders and include it in the model. Plot the geometry.

gm = multicylinder([5 10 15],2);
structuralmodel = createpde("structural","transient-solid");
structuralmodel.Geometry = gm;
pdegplot(structuralmodel,"CellLabels","on","FaceAlpha",0.5)

Set the initial conditions for each cell. When you specify only the initial velocity or initial displacement, structuralIC assumes that the omitted parameter is zero.

structuralIC(structuralmodel,"Displacement",[0;0;0],...
                             "Velocity",[0;0;0], ...
                             "Cell",1);
structuralIC(structuralmodel,"Displacement",[0;0.1;0], ...
                             "Cell",2);
structuralIC(structuralmodel,"Velocity",[0;0.2;0], ...
                             "Cell",3);

Check the initial condition specification for cell 1.

SICACell1 = findStructuralIC(structuralmodel.InitialConditions,"Cell",1)

SICACell1 = 
  GeometricStructuralICs with properties:

             RegionType: 'Cell'
               RegionID: 1
    InitialDisplacement: [3×1 double]
        InitialVelocity: [3×1 double]

SICACell1.InitialDisplacement

ans = 3×1

     0
     0
     0

SICACell1.InitialVelocity

ans = 3×1

     0
     0
     0

Input Arguments

initialConditions — Initial conditions

InitialConditions property of a StructuralModel object

Initial conditions of a transient structural model, specified as the InitialConditions property of a StructuralModel object.

Example: structuralmodel.InitialConditions

RegionType — Geometric region type

"Face" | "Edge" | "Vertex" | "Cell" for a 3-D model

Geometric region type, specified as "Face", "Edge", or "Vertex" for a 2-D model or 3-D model, or "Cell" for a 3-D model.

Data Types: char

RegionID — Geometric region ID

vector of positive integers

Geometric region ID, specified as a vector of positive integers. Find the region IDs by using pdegplot.

Data Types: double

Output Arguments

sica — Structural initial condition assignment

GeometricStructuralICs object | NodalStructuralICs object

Structural initial condition for a particular region, returned as a GeometricStructuralICs or NodalStructuralICs object. For details, see GeometricStructuralICs Properties and NodalStructuralICs Properties.

Version History

Introduced in R2018a

R2025a: To be removed

findStructuralIC will be removed. Use the CellIC, FaceIC, EdgeIC, and VertexIC properties of femodel instead.

For example, you can check the initial conditions for all faces as follows.

model = femodel(AnalysisType="structuralTransient", ...
                Geometry="Block.stl");
model.CellIC = cellIC(Displacement=[0;0;0],Velocity=[0;0;0]);
model.CellIC

The unified finite element model workflow defines the type of a problem and all of its parameters as the properties of an femodel object. This object enables you to specify physical parameters for structural, thermal, and electromagnetic types of analyses. The solver in the unified workflow uses only the parameters (properties) appropriate for the current analysis type while ignoring all other properties. If you switch the analysis type by setting the AnalysisType property of the model, the solver uses the appropriate set of properties corresponding to the new analysis type.

For more help migrating your existing code that uses StructuralModel to the unified finite element workflow, see Migration from Domain-Specific to Unified Workflow.