主要内容

evaluatePrincipalStrain

评估节点位置的主应变

描述

例子

pStrain= evaluatePrincipalStrain (structuralresults用的应变值计算节点位置的主应变structuralresults.对于瞬态和频响结构模型,evaluatePrincipalStrain分别计算所有时间步长或频率步长的主应变。

例子

全部折叠

求解受拉双金属索的静力结构模型,计算八面体剪切应变。

创建一个结构模型。

Structuralmodel = createpde(“结构性”“static-solid”);

创建几何图形并将其包含在模型中。画出几何图形。

Gm =多缸([0.01,0.015],0.05);structuralmodel。几何= gm;pdegplot (structuralmodel“FaceLabels”“上”...“CellLabels”“上”...“FaceAlpha”, 0.5)

图中包含一个轴对象。坐标轴对象包含颤动、贴片、线条类型的3个对象。

指定每种金属的杨氏模量和泊松比。

structuralProperties (structuralmodel“细胞”, 1“YoungsModulus”110 e9,...“PoissonsRatio”, 0.28);structuralProperties (structuralmodel“细胞”2,“YoungsModulus”210 e9,...“PoissonsRatio”, 0.3);

指定面1和面4为固定边界。

structuralBC (structuralmodel“面子”(1、4),“约束”“固定”);

指定面2和面5的表面牵引。

structuralBoundaryLoad (structuralmodel“面子”(2、5),...“SurfaceTraction”, (0, 0, 100));

生成一个网格并解决问题。

generateMesh (structuralmodel);结构结果=求解(结构模型)
strucalresults = staticstrucalresults with properties:位移:[1x1 FEStruct]应变:[1x1 FEStruct]应力:[1x1 FEStruct] VonMisesStress: [22402x1 double]网格:[1x1 FEMesh]

计算节点位置的主应变。

pStrain = evaluatePrincipalStrain(结构结果);

用主应变求应变的一、二不变量。

I1 = pStrain。e1 + pStrain。e2 +pStrain.e3; I2 = pStrain.e1.*pStrain.e2 +...pStrain.e2。* pStrain。e3 +...pStrain.e3。* pStrain.e1;tauOct =√(2*(I1;^ 2 * I2) / 3;pdeplot3D (structuralmodel“ColorMapData”tauOct)

计算谐波激励下梁的主应变和八面体剪切应变。

为一个三维问题创建一个瞬态动态模型。

Structuralmodel = createpde(“结构性”“transient-solid”);

创建一个几何图形并将其包含在模型中。画出几何图形。

Gm = multicuboid(0.06,0.005,0.01);structuralmodel。几何= gm;pdegplot (structuralmodel“FaceLabels”“上”“FaceAlpha”, 0.5)视图(50、20)

图中包含一个轴对象。坐标轴对象包含颤动、贴片、线条类型的3个对象。

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

structuralProperties (structuralmodel“YoungsModulus”210 e9,...“PoissonsRatio”, 0.3,...“MassDensity”, 7800);

固定梁的一端。

structuralBC (structuralmodel“面子”5,“约束”“固定”);

施加正弦位移沿y-方向上与横梁固定端相对的一端。

structuralBC (structuralmodel“面子”3,...“YDisplacement”1的军医,...“频率”, 50);

生成一个网格。

generateMesh (structuralmodel“Hmax”, 0.01);

指定零初始位移和速度。

structuralIC (structuralmodel“位移”(0, 0, 0),“速度”, (0, 0, 0));

求解模型。

Tlist = 0:0.002:0.2;strucalresults = solve(strucalmodel,tlist);

计算梁的主应变。

pStrain = evaluatePrincipalStrain(结构结果);

用主应变求第一不变量和第二不变量。

I1 = pStrain。e1 + pStrain。e2 +pStrain.e3; I2 = pStrain.e1.*pStrain.e2 +...pStrain.e2。* pStrain。e3 +...pStrain.e3。* pStrain.e1;

利用应力不变量计算八面体剪切应变。

tauOct =√(2*(I1;^ 2 * I2) / 3;

画出结果。

图pdeplot3D (structuralmodel,“ColorMapData”结束,tauOct (:,)

输入参数

全部折叠

结构分析问题的解决方法,指定为StaticStructuralResultsTransientStructuralResults,或FrequencyStructuralResults对象。创建structuralresults通过使用解决函数。

例子:结构结果=求解(结构模型)

输出参数

全部折叠

节点位置的主应变,作为结构数组返回。

版本历史

在R2017b中引入

全部展开