模型= createpde;importGeometry(模型,“Block.stl”）;applyBoundaryCondition(模型,“边界条件”面对= 1:4,u = 0);specifyCoefficients(模型、m = 0 d = 0, c = 1 = 0, f = 2);generateMesh(模型);结果= solvepde(模型)

u = results. nodesolution;msh = results.Mesh;pdeplot3D(模型、ColorMapData = u)

Thermalmodel = createpde(thermal=“稳态”）;

importGeometry (thermalmodel“Block.stl”）;pdegplot (thermalmodel FaceLabels =“上”, FaceAlpha = 0.5)轴平等的

thermalProperties (thermalmodel ThermalConductivity = 80);

thermalBC (thermalmodel,面对= 1,温度= 100);thermalBC (thermalmodel,面对= 3,温度= 300);

generateMesh (thermalmodel);热结果=解决(热模型)

pdeplot3D (thermalresults.Mesh ColorMapData = thermalresults.Temperature)

Thermalmodel = createpde(thermal=“稳态”）;

importGeometry (thermalmodel“Block.stl”）;pdegplot (thermalmodel FaceLabels =“上”, FaceAlpha = 0.5)标题(“铜块，厘米”)轴平等的

thermalProperties (thermalmodel ThermalConductivity = 4);

thermalBC (thermalmodel,面对= 1,温度= 373);thermalBC (thermalmodel,面对= 3,温度= 573);

thermalBC (thermalmodel面对= 4,HeatFlux = -20);

generateMesh (thermalmodel);热结果=解决(热模型)

[qx,qy,qz] = evaluateHeatFlux(热流结果);图pdeplot3D (thermalresults。网格，FlowData=[qx qy qz])

[X,Y,Z] =网格(1:26:100,1:6:20,1:11:50);[qx,qy,qz] = evaluateHeatFlux(热流结果，X,Y,Z);

qx =重塑(qx，大小(X));qy =重塑(qy，大小(Y));qz =重塑(qz，大小(Z));图quiver3 (X, Y, Z, qx, qy,求出)

querypoints = [X(:) Y(:) Z(:)]';[qx,qy,qz] = evaluateHeatFlux(热流结果，查询点);qx =重塑(qx，大小(X));qy =重塑(qy，大小(Y));qz =重塑(qz，大小(Z));图quiver3 (X, Y, Z, qx, qy,求出)

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

importGeometry (structuralmodel“SquareBeam.stl”）;pdegplot (structuralmodel FaceLabels =“上”, FaceAlpha = 0.5)

structuralProperties (structuralmodel PoissonsRatio = 0.3,.．.YoungsModulus = 210 e3);

structuralBC (structuralmodel面对= 6,约束=“固定”）;

structuralBoundaryLoad (structuralmodel面对= 5,.．.SurfaceTraction = (0, 0; 2));

generateMesh (structuralmodel);结构结果= solve(结构模型);

图pdeplot3D (structuralresults。网,.．.ColorMapData = structuralresults。VonMisesStress,.．.变形= structuralresults.Displacement)

图pdeplot3D (structuralresults。网,.．.ColorMapData = structuralresults。VonMisesStress,.．.变形= structuralresults。位移,.．.DeformationScaleFactor = 500)

图pdeplot3D (structuralresults。网,.．.ColorMapData = structuralresults.VonMisesStress)

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 e9 = 210,.．.PoissonsRatio = 0.3,.．.MassDensity = 7800);

structuralBC (structuralmodel面对= 5,约束=“固定”）;

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);

vmStress = evaluateVonMisesStress(strucalresults);

图pdeplot3D (structuralresults。网,ColorMapData= vmStress(:,end)) title(《冯·米塞斯在光束上的最后一阶应力》）

Emagmodel = createpde(电磁=“静电”）;

gm = importGeometry(emagmodel，“PlateHoleSolid.stl”）;pdegplot(通用、FaceLabels =“上”, FaceAlpha = 0.3)

emagmodel。真空介电常数= 8.8541878128E-12;

electromagneticProperties (emagmodel RelativePermittivity = 1);

electromagneticSource (emagmodel ChargeDensity = 5 e-9);

electromagneticBC (emagmodel、电压= 0,面对=三6);electromagneticBC (emagmodel,电压= 1000,面对= 7);

generateMesh (emagmodel);

R = solve(emagmodel)

图pdeplot3D (R.Mesh ColorMapData = R.ElectricPotential)

pdeplot3D (R.Mesh FlowData = [R.ElectricField.Ex.．.R.ElectricField.Ey.．.R.ElectricField.Ez])

模型= createpde;importGeometry(模型,“Tetrahedron.stl”）;网格= generateMesh(模型，Hmax=20，几何顺序=“线性”）;

pdeplot3D(网)

pdeplot3D(模型)

pdeplot3D (mesh.Nodes mesh.Elements)

pdeplot3D(网格,NodeLabels =“上”12)视图(101)

pdeplot3D(网格,ElementLabels =“上”12)视图(101)