findStructuralIC
Find initial displacement and velocity assigned to geometric region
Domain-specific structural workflow is not recommended. New features might not be compatible with this workflow. For help migrating your existing code to the unified finite element workflow, see Migration from Domain-Specific to Unified Workflow.
Description
returns the initial displacement and velocity assigned to the specified
region.sica
= findStructuralIC(initialConditions
,RegionType
,RegionID
)
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: [3x1 double] InitialVelocity: [3x1 double]
SICACell1.InitialDisplacement
ans = 3×1
0
0
0
SICACell1.InitialVelocity
ans = 3×1
0
0
0
Find Initial Displacement Set as Previously Obtained Static Solution
Use a static solution as an initial condition for a dynamic structural model. Check and plot the initial displacement.
Create a static model.
staticmodel = createpde("structural","static-solid");
Create the geometry and include it in the model. Plot the geometry.
gm = multicuboid(0.06,0.005,0.01); staticmodel.Geometry = gm; pdegplot(staticmodel,"FaceLabels","on","FaceAlpha",0.5) view(50,20)
Specify Young's modulus, Poisson's ratio, and the mass density.
structuralProperties(staticmodel,"YoungsModulus",210E9, ... "PoissonsRatio",0.3,... "MassDensity",7800);
Apply the boundary condition and static load.
structuralBC(staticmodel,"Face",5,"Constraint","fixed"); structuralBoundaryLoad(staticmodel,"Face",3, ... "SurfaceTraction",[0;1E6;0]); generateMesh(staticmodel,"Hmax",0.02); Rstatic = solve(staticmodel);
Create a dynamic model and assign geometry.
dynamicmodel = createpde("structural","transient-solid"); gm = multicuboid(0.06,0.005,0.01); dynamicmodel.Geometry = gm;
Apply the boundary condition.
structuralBC(dynamicmodel,"Face",5,"Constraint","fixed");
Specify the initial condition using the static solution.
generateMesh(dynamicmodel,"Hmax",0.02);
structuralIC(dynamicmodel,Rstatic)
ans = NodalStructuralICs with properties: InitialDisplacement: [189x3 double] InitialVelocity: [189x3 double]
Check the initial condition specification for dynamicmodel
.
sica = findStructuralIC(dynamicmodel.InitialConditions,"Cell",1)
sica = NodalStructuralICs with properties: InitialDisplacement: [189x3 double] InitialVelocity: [189x3 double]
Plot the z
-component of the initial displacement.
pdeplot3D(dynamicmodel,"ColorMapData",sica.InitialDisplacement(:,3)) title("Initial Displacement in the Z-direction")
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
MATLAB Command
You clicked a link that corresponds to this MATLAB command:
Run the command by entering it in the MATLAB Command Window. Web browsers do not support MATLAB commands.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
You can also select a web site from the following list:
How to Get Best Site Performance
Select the China site (in Chinese or English) for best site performance. Other MathWorks country sites are not optimized for visits from your location.
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)