主要内容

evaluateReaction

Evaluate reaction forces on boundary

Description

example

F= evaluateReaction(structuralresults,RegionType,区域)评估在指定的边界上的反作用力RegionType区域. The function uses the global Cartesian coordinate system. For transient and frequency response structural models,evaluateReactionevaluates reaction forces for all time- and frequency-steps, respectively.

Examples

collapse all

Create a static structural model.

structuralmodel = createpde('structural','static-solid');

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

structuralmodel.Geometry = multicuboid(0.01,0.01,0.05); pdegplot(structuralmodel,'FaceLabels','on','FaceAlpha',0.5);

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

指定Young的模量和Poisson的比例。

structuralProperties(structuralmodel,'Youngsmodulus',210E9,...'PoissonsRatio',0.3);

固定条的一端,并向另一端施加压力。

结构BC(结构模型,'Face',1,'Constraint','fixed')
ANS =具有属性的结构BC:regionType:'face'regionId:1矢量化:'OFF'边界约束和强制置换位移位移:[] xdisplacement:[] ydisplacement:[] ydisplacement:[] zdisplacement:[] zdisplacement:[] zdisplacement:[]约束:“固定” radius:[] radius:[] radius:[]参考:[]参考:[]参考:[]参考:[]参考:[]:[]标签:[]边界载荷力:[]表面效果:[]压力:[]转换速度:[]标签:[]
structuralBoundaryLoad(structuralmodel,'Face',2,'Pressure', 100)
ans = StructuralBC with properties: RegionType: 'Face' RegionID: 2 Vectorized: 'off' Boundary Constraints and Enforced Displacements Displacement: [] XDisplacement: [] YDisplacement: [] ZDisplacement: [] Constraint: [] Radius: [] Reference: [] Label: [] Boundary Loads Force: [] SurfaceTraction: [] Pressure: 100 TranslationalStiffness: [] Label: []

Generate a mesh and solve the problem.

generateMesh(structuralmodel,'hmax',0.003);structuralResults = solve(structuralModel);

Compute the reaction forces on the fixed end.

reaction = evaluateReaction(structuralresults,'Face',1)
reaction =struct with fields:Fx: -1.3620e-06 Fy: 2.2303e-06 Fz: 0.0103

Evaluate the reaction forces at the fixed end of a beam subject to harmonic excitation.

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

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

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

gm =多蛋布(0.06,0.005,0.01);bunstoralModel.DEMETRY = GM;PDEGPLOT(结构模型,'FaceLabels','on','FaceAlpha',0.5) view(50,20)

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

指定杨氏模量,泊松比和材料的质量密度。

structuralProperties(structuralmodel,'Youngsmodulus',210E9,...'PoissonsRatio',0.3,...'MassDensity',7800);

Fix one end of the beam.

结构BC(结构模型,'Face',5,'Constraint','fixed');

沿着沿着的正弦位移y-direction on the end opposite the fixed end of the beam.

结构BC(结构模型,'Face',3,...'YDisplacement',1E-4,...'频率',50);

生成网格。

generateMesh(structuralmodel,'hmax',0.01);

指定零初始位移和速度。

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

Solve the model.

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

Compute the reaction forces on the fixed end.

reaction = evaluateReaction(structuralresults,'Face',5)
reaction =struct with fields:FX:[101x1 double] FY:[101x1 double] fz:[101x1 double]

Input Arguments

collapse all

Solution of the structural analysis problem, specified as aStaticStructuralResults,构造术, orFrequencyStructuralResultsobject. Createstructuralresultsby using thesolvefunction.

Example:structuralResults = solve(structuralModel)

几何区域类型,指定为'Edge'for a 2-D model or'Face'对于3-D模型.

Example:evaluateReaction(structuralresults,'Face',2)

Data Types:char|string

几何区域ID,指定为正整数的向量。通过使用PDEGPLOT.

Example:evaluateReaction(structuralresults,'Face',2)

Data Types:double

Output Arguments

collapse all

反作用力,作为结构阵列返回。阵列场代表集成的反作用力和表面牵引矢量,通过使用边界和外部正常的应力状态计算出来。

版本嗨story

Introduced in R2017b