interpolateVonMisesStress
在任意空间位置插值冯·米塞斯应力
语法
描述
intrpVMStress= interpolateVonMisesStress (structuralresults，xq，yq）xq而且yq．对于瞬态和频响结构模型，interpolateVonMisesStress分别为所有时间步长或频率步长插值冯·米塞斯应力。
intrpVMStress= interpolateVonMisesStress (structuralresults，xq，yq，的zq）xq，yq,的zq．
intrpVMStress= interpolateVonMisesStress (structuralresults，querypoints）querypoints．
例子
平面应变问题的内插von Mises应力
建立平面应变问题的结构分析模型。
Structuralmodel = 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(结构模型);
创建一个网格，并将冯·米塞斯的应力插入网格。
V = linspace(-1,1,151);[X,Y] = meshgrid(v);intrpVMStress = interpolateVonMisesStress(strucalresults,X,Y);
将冯·米塞斯的应力重塑成网格的形状并绘制出来。
VMStress =重塑(intrpVMStress,size(X));p = pcolor(X,Y,VMStress);p.EdgeColor =“没有”；colorbar
三维静力结构分析问题的Von Mises应力插值
求解一个表示受拉双金属索的静态结构模型，并在索的横截面上插值冯米塞斯应力。
为解决一个实体(3-D)问题创建一个静态结构模型。
Structuralmodel = 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);结构结果=求解(结构模型)
strucalresults = staticstrucalresults with properties:位移:[1x1 FEStruct]应变:[1x1 FEStruct]应力:[1x1 FEStruct] VonMisesStress: [22402x1 double]网格:[1x1 FEMesh]
定义电缆跨中截面的坐标。
[X,Y] = meshgrid(linspace(-0.015,0.015,50));Z = ones(size(X))*0.025;
插入冯·米塞斯应力并绘制结果。
IntrpVMStress = interpolateVonMisesStress(strucalresults,X,Y,Z);冲浪(X, Y,重塑(IntrpVMStress、大小(X)))
或者，您可以使用查询点矩阵来指定网格。
querypoints = [X(:)，Y(:)，Z(:)]';IntrpVMStress =.．.interpolateVonMisesStress (structuralresults querypoints);冲浪(X, Y,重塑(IntrpVMStress、大小(X)))
三维结构动力学问题的内插von Mises应力
在简谐激励下，在梁的几何中心插值冯·米塞斯应力。
为一个三维问题创建一个瞬态动态模型。
Structuralmodel = 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;strucalresults = solve(strucalmodel,tlist);
在梁的几何中心插入冯·米塞斯应力。
coordsMidSpan = [0;0;0.005];intrpStress = interpolateStress(strucalresults,coordsMidSpan);
画出梁的几何中心的冯·米塞斯应力。
图(strucalresults . solutiontimes,intrpStress.sxx)“梁中心的冯·米塞斯应力”）
输入参数
输出参数
版本历史
在R2017b中引入
您也可以从以下列表中选择一个网站: