interpolateStrain
插值应变在任意空间位置
语法
描述
intrpStrain= interpolateStrain (structuralresults,xq,yq）xq而且yq．对于瞬态和频响结构模型，interpolateStrain插值应变的所有时间或频率步长，分别。
intrpStrain= interpolateStrain (structuralresults,xq,yq,的zq）xq,yq,的zq．
intrpStrain= interpolateStrain (structuralresults,querypoints）querypoints．
例子
平面应变问题的插值应变
为平面应变问题建立结构分析模型。
结构模型= createpde(“结构性”,“static-planestrain”）;
在模型中包含正方形几何。画出几何图形。
geometryFromEdges (structuralmodel @squareg);pdegplot (structuralmodel“EdgeLabels”,“上”)轴平等的
指定杨氏模量和泊松比。
structuralProperties (structuralmodel“PoissonsRatio”, 0.3,.．.“YoungsModulus”210年e3);
指定x边1的强制位移的-分量。
structuralBC (structuralmodel“XDisplacement”, 0.001,“边缘”1);
指定边3为固定边界。
structuralBC (structuralmodel“约束”,“固定”,“边缘”3);
生成一个网格并解决问题。
generateMesh (structuralmodel);结构结果= solve(结构模型);
创建一个网格并插值x- - -y-网格的法向应变分量。
V = linspace(-1,1,101);[X,Y] = meshgrid(v);intrpStrain = interpolateStrain(strucalresults,X,Y);
重塑x法向应变的-分量到网格的形状，并绘制它。
exx =重塑(intrpStrain.exx,size(X));px = pcolor(X,Y,exx);px。EdgeColor =“没有”；colorbar
重塑y法向应变的-分量到网格的形状，并绘制它。
eyy =重塑(intrpStrain.eyy,size(Y));图py = pcolor(X,Y, Y);py。EdgeColor =“没有”；colorbar
三维静力结构分析问题的插值应变法
求解受拉双金属拉索的静态结构模型，插值拉索截面上的应变。
为解决实体(3d)问题创建一个静态结构模型。
结构模型= createpde(“结构性”,“static-solid”）;
创建几何图形并将其包含在模型中。画出几何图形。
Gm =多缸([0.01,0.015]，0.05);structuralmodel。几何= gm;pdegplot (structuralmodel“FaceLabels”,“上”,.．.“CellLabels”,“上”,.．.“FaceAlpha”, 0.5)
指定每种金属的杨氏模量和泊松比。
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);结构结果= solve(结构模型)
strucalresults = staticstrucalresults with properties:位移:[1x1 FEStruct]应变:[1x1 FEStruct]应力:[1x1 FEStruct] VonMisesStress: [22402x1 double] Mesh: [1x1 FEMesh]
确定缆绳跨中截面的坐标。
[X,Y] = meshgrid(linspace(-0.015,0.015,50));Z = ones(size(X))*0.025;
插值应变并画出结果。
intrpStrain = interpolateStrain(strucalresults,X,Y,Z);冲浪(X, Y,重塑(intrpStrain.ezz、大小(X)))
或者，您可以通过使用查询点矩阵来指定网格。
querypoints = [X(:)，Y(:)，Z(:)]';intrpStrain = interpolateStrain(结构结果，查询点);冲浪(X, Y,重塑(intrpStrain.ezz、大小(X)))
三维结构动力学问题的插值应变法
在谐波激励下插值梁几何中心处的应变。
为一个三维问题创建一个瞬态动态模型。
结构模型= createpde(“结构性”,“transient-solid”）;
创建一个几何图形，并将其包含在模型中。画出几何图形。
Gm = multicuboid(0.06,0.005,0.01);structuralmodel。几何= gm;pdegplot (structuralmodel“FaceLabels”,“上”,“FaceAlpha”, 0.5)视图(50、20)
指定材料的杨氏模量、泊松比和质量密度。
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;Structuralresults = solve(structuralmodel,tlist);
在梁的几何中心插值应变。
coordsMidSpan = [0;0;0.005];intrpStrain = interpolateStrain(strucalresults,coordsMidSpan);
画出梁几何中心的法向应变。
图(structuralresults.SolutionTimes,intrpStrain.exx)标题(“梁中心x方向法向应变”）
输入参数
输出参数
版本历史
在R2017b中引入
您也可以从以下列表中选择网站: