pdeplot
二维问题的图解或网格
语法
描述
pdeplot (为二维结构分析模型绘制von Mises应力，并显示变形形状。模型, results.VonMisesStress XYData的“变形”,results.Displacement)
例子
二维网格图
创建PDE模型。包括内置函数的几何形状lshapeg．网格几何和绘图。
模型= createpde;geometryFromEdges(模型、@lshapeg);网格= generateMesh(模型);pdeplot(模型)
或者，您可以使用网作为输入参数。
pdeplot(网)
另一种方法是使用网格的节点和元素作为输入参数pdeplot．
pdeplot (mesh.Nodes mesh.Elements)
显示节点标签。使用xlim而且ylim放大特定节点。
pdeplot(模型,“NodeLabels”，“上”xlim([-0.2,0.2])
显示元素标签。
pdeplot(模型,“ElementLabels”，“上”xlim([-0.2,0.2])
解决方案图
创建PDE模型解的彩色2-D和3-D图。
创建PDE模型。包括内置函数的几何形状lshapeg．网格几何。
模型= createpde;geometryFromEdges(模型、@lshapeg);generateMesh(模型);
在所有边上设置零狄利克雷边界条件。
applyBoundaryCondition(模型,“边界条件”，.．.“边缘”1: model.Geometry.NumEdges,.．.“u”, 0);
指定系数并求解偏微分方程。
specifyCoefficients(模型,“米”0,.．.' d '0,.．.“c”, 1.．.“一个”0,.．.“f”1);结果= solvepde(模型)
results = StationaryResults with properties: nodesolution: [1177x1 double] XGradients: [1177x1 double] YGradients: [1177x1 double] ZGradients: [] Mesh: [1x1 FEMesh]
在节点位置访问解决方案。
u = results. nodesolution;
画出二维解。
pdeplot(模型,“XYData”, u)
画出三维解。
pdeplot(模型,“XYData”u“ZData”, u)
溶液颤振图
将PDE溶液的梯度绘制为颤振图。
创建PDE模型。包括内置函数的几何形状lshapeg．网格几何。
模型= createpde;geometryFromEdges(模型、@lshapeg);generateMesh(模型);
在所有边上设置零狄利克雷边界条件。
applyBoundaryCondition(模型,“边界条件”，.．.“边缘”1: model.Geometry.NumEdges,.．.“u”, 0);
指定系数并求解偏微分方程。
specifyCoefficients(模型,“米”0,.．.' d '0,.．.“c”, 1.．.“一个”0,.．.“f”1);结果= solvepde(模型)
results = StationaryResults with properties: nodesolution: [1177x1 double] XGradients: [1177x1 double] YGradients: [1177x1 double] ZGradients: [] Mesh: [1x1 FEMesh]
在节点位置访问解决方案的梯度。
ux = results.XGradients;uy = results.YGradients;
把梯度画成箭袋图。
pdeplot(模型,“FlowData”(用户体验,uy))
综合图
将二维偏微分方程的解画成三维“喷气机”着色和网格，并包括一个箭袋情节。获取轴对象的句柄。
创建PDE模型。包括内置函数的几何形状lshapeg．网格几何。
模型= createpde;geometryFromEdges(模型、@lshapeg);generateMesh(模型);
在所有边上设置零狄利克雷边界条件。
applyBoundaryCondition(模型,“边界条件”，.．.“边缘”1: model.Geometry.NumEdges,.．.“u”, 0);
指定系数并求解偏微分方程。
specifyCoefficients(模型,“米”0,.．.' d '0,.．.“c”, 1.．.“一个”0,.．.“f”1);结果= solvepde(模型)
results = StationaryResults with properties: nodesolution: [1177x1 double] XGradients: [1177x1 double] YGradients: [1177x1 double] ZGradients: [] Mesh: [1x1 FEMesh]
在节点位置访问解决方案及其梯度。
u = results. nodesolution;ux = results.XGradients;uy = results.YGradients;
用三维图形绘制解“喷气机”上色和一个网格，并包括渐变作为一个箭筒图。
H = pdeploy(模型，“XYData”u“ZData”u.．.“FaceAlpha”, 0.5,.．.“FlowData”(用户体验,uy),.．.“ColorMap”，“喷气机”，.．.“网”，“上”）
h = 3x1图形数组:Patch Quiver ColorBar
瞬态热模型解
求解一个二维瞬态热问题。
为这个问题创建一个瞬态热模型。
Thermalmodel = createpde(“热”，瞬态的）;
创建几何图形并将其包含在模型中。
Sq1 = [3;4;0;3;3;0;0;0;3;3); D1 = [2; 4; 0.5; 1.5; 2.5; 1.5; 1.5; 0.5; 1.5; 2.5]; gd = [SQ1 D1]; sf =“于SQ1 + D1”；Ns = char(“于SQ1”，“D1”）;Ns = Ns ';Dl = decsg(gd,sf,ns);geometryFromEdges (thermalmodel dl);pdegplot (thermalmodel“EdgeLabels”，“上”，“FaceLabels”，“上”) xlim([-1.5 4.5]) ylim([-0.5 3.5]平等的
对于正方形区域，分配这些热特性:
导热系数为
质量密度为
比热为
thermalProperties (thermalmodel“ThermalConductivity”10.．.“MassDensity”，2，.．.“SpecificHeat”, 0.1,.．.“脸”1);
对于金刚石区域，分配这些热特性:
导热系数为
质量密度为
比热为
thermalProperties (thermalmodel“ThermalConductivity”，2，.．.“MassDensity”, 1.．.“SpecificHeat”, 0.1,.．.“脸”2);
假设菱形区域是一个密度为的热源 ．
internalHeatSource (thermalmodel 4“脸”2);
涂上恒温的 在方板的两侧。
thermalBC (thermalmodel“温度”0,“边缘”，[1 2 7 8]);
设置初始温度为0°C。
thermalIC (thermalmodel 0);
生成网格。
generateMesh (thermalmodel);
这个问题的动力学非常快。温度在0.1秒左右达到稳定状态。要捕获动态中有趣的部分，请将解决方案时间设置为logspace (2, 1, 10)．该命令返回10个在0.01到0.1之间的对数间隔解时间。
Tlist = logspace(-2，-1,10);
解这个方程。
Thermalresults = solve(thermalmodel,tlist);
用等高线画出解的等温线。
T =热结果温度;pdeplot (thermalmodel“XYData”T (10),“轮廓”，“上”，“ColorMap”，“热”）
平面应变问题的平面变形图
为静态平面应变问题建立结构分析模型。
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(结构模型);
使用默认比例因子绘制变形形状。默认情况下,pdeplot内部根据几何尺寸和变形大小确定比例因子。
pdeplot (structuralmodel.．.“XYData”, structuralresults。VonMisesStress,.．.“变形”, structuralresults。位移,.．.“ColorMap”，“喷气机”）
用比例因子500绘制变形的形状。
pdeplot (structuralmodel.．.“XYData”, structuralresults。VonMisesStress,.．.“变形”, structuralresults。位移,.．.“DeformationScaleFactor”, 500,.．.“ColorMap”，“喷气机”）
在不缩放的情况下绘制变形形状。
pdeplot (structuralmodel“XYData”, structuralresults。VonMisesStress,.．.“ColorMap”，“喷气机”）
结构模型模态分析解
假定平面应力条件下，求二维悬臂梁的基本(最低)模态。
指定梁的以下几何和结构特性，以及单位平面应力厚度。
长度= 5;高度= 0.1;E = 3e7;Nu = 0.3;Rho = 0.3/386;
创建一个平面应力模型，分配一个几何，并生成一个网格。
Structuralmodel = createpde(“结构”，“modal-planestress”）;GDM = [3;4;0;length;length;0;0; height;height];G = decsg(gdm，“S1 '，（“S1 ') ');geometryFromEdges (structuralmodel g);
定义一个最大元素大小(五个元素通过光束厚度)。
Hmax =身高/5;msh = generateMesh (structuralmodel,“Hmax”, hmax);
指定结构属性和边界约束。
structuralProperties (structuralmodel“YoungsModulus”, E,.．.“MassDensity”ρ,.．.“PoissonsRatio”ν);structuralBC (structuralmodel“边缘”4“约束”，“固定”）;
利用光束理论计算分析基频(Hz)。
I =身高^3/12;analyticalOmega1 = 3.516*根号(E*I/(长度^4*(rho*高度)))/(2*pi)
analyticalOmega1 = 126.9498
指定一个频率范围，其中包括分析计算的频率并求解模型。
Modalresults = solve(结构模型，“FrequencyRange”(0,1 e6))
modalresults = ModalStructuralResults with properties: NaturalFrequencies: [32x1 double] ModeShapes: [1x1 FEStruct] Mesh: [1x1 FEMesh]
求解器在节点位置找到固有频率和模态位移值。要访问这些值，请使用modalresults。NaturalFrequencies而且modalresults。ModeShapes．
modalresults.NaturalFrequencies /(2 *π)
ans =32×1105× 0.0013 0.0079 0.0222 0.0433 0.0711 0.0983 0.1055 0.1462 0.1930 0.2455
modalresults。ModeShapes
ans = FEStruct with properties: ux: [6511x32 double] uy: [6511x32 double] Magnitude: [6511x32 double]
画出y-基频解的分量。
pdeplot (structuralmodel“XYData”modalresults.ModeShapes.uy(: 1)标题([“有频率的第一模式”，.．.num2str (modalresults.NaturalFrequencies(1) /(2 *π)),“赫兹”])轴平等的
二维静电分析模型的解
求解一个电磁问题，求出代表带孔板的二维几何图形的电势和场分布。
创建用于静电分析的电磁模型。
Emagmodel = createpde(“电磁”，“静电”）;
导入并绘制表示带孔板的几何图形。
importGeometry (emagmodel“PlateHolePlanar.stl”）;pdegplot (emagmodel“EdgeLabels”，“上”）
指定国际单位制中的真空介电常数。
emagmodel。真空介电常数= 8.8541878128E-12;
指定材料的相对介电常数。
electromagneticProperties (emagmodel“RelativePermittivity”1);
在矩形和圆的边框上应用电压边界条件。
electromagneticBC (emagmodel“电压”0,“边缘”1:4);electromagneticBC (emagmodel“电压”, 1000,“边缘”5);
指定整个几何的电荷密度。
electromagneticSource (emagmodel“ChargeDensity”5 e-9);
生成网格。
generateMesh (emagmodel);
求解模型。
R = solve(emagmodel)
R = ElectrostaticResults with properties: ElectricPotential: [1218x1 double] ElectricField: [1x1 FEStruct] ElectricFluxDensity: [1x1 FEStruct] Mesh: [1x1 FEMesh]
画出电势和电场。
pdeplot (emagmodel“XYData”, R。电位,.．.“FlowData”, (R.ElectricField.Ex.．.R.ElectricField.Ey])轴平等的
[p,e,t]网格和解图
画出p, e, t网。使用2-D和3-D彩色图显示解决方案。
创建几何，网格，边界条件，PDE系数和解决方案。
[p,e,t] = initmesh(“lshapeg”）;U =汇编(“lshapeb”, p, e t - 1, 0, 1);
绘制网格。
pdeplot (p, e, t)
将解绘制为二维彩色图。
pdeplot (p, e t“XYData”, u)
将解绘制为三维彩色图。
pdeplot (p, e t“XYData”u“ZData”, u)
输入参数
模型- - - - - -模型对象
PDEModel对象|ThermalModel对象|StructuralModel对象|ElectromagneticModel对象
模型对象，指定为PDEModel对象,ThermalModel对象,StructuralModel对象,或ElectromagneticModel对象。
例子:模型= createpde(1)
例子:Thermalmodel = createpde('thermal'，'steadystate')
例子:Structuralmodel = createpde('structural'，'static-solid')
例子:Emagmodel = createpde('电磁'，'静磁')
网- - - - - -网格物体
网的属性PDEModel对象|的输出generateMesh
对象，指定为网的属性PDEModel对象或作为的输出generateMesh．
例子:模型。网
节点- - - - - -节点坐标
2 -NumNodes矩阵
节点坐标，指定为2 *NumNodes矩阵。NumNodes为节点数。
元素- - - - - -以节点id表示的元素连通性矩阵
3 * -NumElements矩阵|6日-NumElements矩阵
以节点id表示的元素连通性矩阵，指定为3-by-NumElements或6日-NumElements矩阵。线性网格只包含角节点。对于线性网格，连通性矩阵每个二维元素有三个节点。二次网格包含角节点和元素每条边中间的节点。对于二次网格，连通性矩阵每个二维元素有六个节点。
名称-值参数
的可选逗号分隔对名称,值参数。的名字参数名称和价值对应的值。的名字必须出现在引号内。您可以以任意顺序指定多个名称和值对参数Name1, Value1,…,的家．
pdeplot(模型、“XYData”u“ZData”,u)
当你使用PDEModel对象,pdeplot(模型、“XYData”u“ZData”,u)将表面绘图着色设置为解决方案u，并将3d图的高度设置为u．在这里u是一个NodalSolution返回的PDE结果的属性solvepde或solvepdeeig．
当你使用[p, e t]表示,pdeplot (p, e t XYData, u, ZData, u)将表面绘图着色设置为解决方案u并为解决方案设置3d图的高度u．在这里u遗留求解器是否返回解决方案，例如assempde．
提示
项中指定至少一个FlowData(向量场图)，XYData(彩色曲面图)，或ZData(三维高度图)名-值对。否则,pdeplot绘制没有数据的网格。
输出参数
更多关于
您也可以从以下列表中选择一个网站:
