解决
解决结构分析、传热或电磁分析问题
语法
描述
结果=解决(有限元法“DecayRange”,(λ₁,lambda2))有限元法所有模式衰减范围(λ₁,lambda2)。由此产生的模式使您可以:
用模态叠加法来加速瞬态热分析。
提取减少模态系统使用,例如,在仿真软件金宝app®。
结果=解决(有限元法,tlist“ModalResults”,thermalModalR)结果=解决(有限元法,tlist“ModalResults”,structuralModalR)结果=解决(有限元法,弗利斯特雨“ModalResults”,structuralModalR)
structuralStaticResults=解决(structuralStatic)structuralStatic。
structuralTransientResults=解决(structuralTransient,tlist)structuralTransient中指定的时间tlist。
structuralFrequencyResponseResults=解决(structuralFrequencyResponse,弗利斯特雨)structuralFrequencyResponse在指定的频率弗利斯特雨。
structuralModalResults=解决(structuralModal“FrequencyRange”,(ω,₂))(ω,₂)。定义ω略低于最低频率和预期₂是略高于期望频率最高。例如,如果最低的期望频率为零,然后用一个小的负面价值ω。
structuralTransientResults=解决(structuralTransient,tlist“ModalResults”,structuralModalR)structuralFrequencyResponseResults=解决(structuralFrequencyResponse,弗利斯特雨“ModalResults”,structuralModalR)
thermalSteadyStateResults=解决(thermalSteadyState)thermalSteadyState。
thermalTransientResults=解决(thermalTransient,tlist)thermalTransient中指定的时间tlist。
thermalModalResults=解决(thermalModal“DecayRange”,(λ₁,lambda2))thermalModal所有模式衰减范围(λ₁,lambda2)。由此产生的模式使您可以:
用模态叠加法来加速瞬态热分析。
提取减少模态系统使用,例如,在仿真软件。金宝app
thermalModalResults=解决(thermalModal“快照”,Tmatrix)
thermalTransientResults=解决(thermalTransient,tlist“ModalResults”,thermalModalR)
emagStaticResults=解决(emagmodel)emagmodel。
emagHarmonicResults=解决(emagmodel“频率”,ω)emagmodel在指定的频率ω。
例子
解决静态结构模型
解决一个静态结构模型代表一个双金属电缆张力。
为一个坚实的创建一个静态结构模型(三维)的问题。
structuralmodel = createpde (“结构性”,“static-solid”);
创建几何和包括在模型中。绘制几何图形。
通用=多缸(0.05 [0.01 - 0.015]);structuralmodel。几何=通用;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);structuralresults =解决(structuralmodel)
structuralresults = StaticStructuralResults属性:位移:[1 x1 FEStruct]应变:[1 x1 FEStruct]压力:[1 x1 FEStruct] VonMisesStress: [22402 x1双]网:[1 x1 FEMesh]
解决发现的价值位移、应力、应变和•冯•米塞斯应力节点位置。要访问这些值,使用structuralresults.Displacement,structuralresults.Stress,等等。位移、应力和应变值在返回节点的位置FEStruct代表他们的组件对象的属性。注意,属性的FEStruct对象是只读的。
structuralresults.Displacement
ans = FEStruct属性:用户体验:[22402 x1双]uy: [22402 x1双]是乌斯:[22402 x1双]级:[22402 x1双]
structuralresults.Stress
ans = FEStruct属性:sxx: [22402 x1双]syy: [22402 x1双]szz: [22402 x1双]syz: [22402 x1双]sxz: [22402 x1双]sxy: [22402 x1双]
structuralresults.Strain
ans = FEStruct属性:exx: [22402 x1双]eyy: [22402 x1双]ezz: [22402 x1双]eyz: [22402 x1双]exz: [22402 x1双]exy: [22402 x1双]
情节的变形形状z分的正常压力。
pdeplot3D (structuralmodel…“ColorMapData”structuralresults.Stress.szz,…“变形”structuralresults.Displacement)
解决瞬态结构模型
解决瞬态响应的薄3 d板下谐波负荷中心。
创建一个瞬态动力学模型的三维问题。
structuralmodel = createpde (“结构性”,“transient-solid”);
创建几何和包括在模型中。绘制几何图形。
通用= multicuboid ([5, 0.05), 5, 0.05, 0.01);structuralmodel。几何=通用;pdegplot (structuralmodel“FaceLabels”,“上”,“FaceAlpha”,0.5)
放大的脸小板上的标签的中心。
图pdegplot (structuralmodel,“FaceLabels”,“上”,“FaceAlpha”,0.25)轴([-0.2 0.2 -0.2 0.2 -0.1 0.1])
指定杨氏模量、泊松比和材料的质量密度。
structuralProperties (structuralmodel“YoungsModulus”210 e9,…“PoissonsRatio”,0.3,…“MassDensity”,7800);
指定所有的脸瘦的外围3 d板是固定的边界。
structuralBC (structuralmodel“约束”,“固定”,“面子”,8);
应用一个正弦压力载荷的小脸上在板的中心。
structuralBoundaryLoad (structuralmodel“面子”12…“压力”5 e7,…“频率”25);
生成一个网格和线性元素。
generateMesh (structuralmodel“GeometricOrder”,“线性”,“Hmax”,0.2);
指定零初始位移和速度。
structuralIC (structuralmodel“位移”(0,0,0),“速度”,(0,0,0));
解决模型。
tlist = linspace (0, 1300);structuralresults =解决(structuralmodel tlist);
解算器发现的价值位移,速度,加速度的节点位置。要访问这些值,使用structuralresults.Displacement,structuralresults.Velocity,等等。位移、速度和加速度值返回FEStruct代表他们的组件对象的属性。注意,属性的FEStruct对象是只读的。
structuralresults.Displacement
ans = FEStruct属性:用户体验:[1873 x300双]uy: [1873 x300双]是乌斯:[1873 x300双]级:[1873 x300双]
structuralresults.Velocity
ans = FEStruct属性:vx:(1873 x300双v: [1873 x300双]vz: [1873 x300双]级:[1873 x300双]
structuralresults.Acceleration
ans = FEStruct属性:ax: [1873 x300双]ay: [1873 x300双]阿兹:[1873 x300双]级:[1873 x300双]
频率响应分析
执行音叉的频率响应分析。
首先,创建一个结构模型的模态分析固体音叉。
模型= createpde (“结构性”,“frequency-solid”);
导入音叉几何。
importGeometry(模型,“TuningFork.stl”);
指定杨氏模量、泊松比和质量密度模型线性弹性材料的行为。指定所有物理性质的单位一致。
structuralProperties(模型,“YoungsModulus”210 e9,…“PoissonsRatio”,0.3,…“MassDensity”,8000);
识别面孔应用边界约束和载荷通过绘制的几何面标签。
图(“单位”,“归一化”,“outerposition”,(0 0 1 1))pdegplot(模型,“FaceLabels”,“上”15)视图(-50)标题(“几何与标签”)
施加足够的边界限制,以防止刚体运动应用载荷作用下。通常,您举行一个音叉用手或挂载到桌子上。创建一个简单的近似边界条件,解决一个地区交界处附近的表面和处理(面临21和22)。
structuralBC(模型,“面子”(21、22),“约束”,“固定”);
指定的压力加载在一个齿(面对11)作为一个短矩形压力脉冲。在频域,这种压力脉冲是一个单位荷载均匀分布在所有频率。
structuralBoundaryLoad(模型,“面子”11“压力”1);弗利斯特雨= linspace (0、4000、150);网= generateMesh(模型,“Hmax”,0.005);R =解决(模型2 *π*弗利斯特雨);
齿尖的振动频率,这是面对12。发现节点的脸上和情节y位移对频率的分量,利用其中的一个节点。
excitedTineTipNodes = findNodes(网、“地区”,“面子”12);tipDisp = R.Displacement.uy (excitedTineTipNodes (1):);图绘制(弗利斯特雨、abs (tipDisp))包含(“频率”);ylabel (“| Y-Displacement |”);
模态分析结构模型解决方案
找到最基本的二维悬臂梁(最低)模式,假设患病率的平面应力状态。
指定几何和梁的结构属性,以及一个平面应力单元厚度。
长度= 5;身高= 0.1;E = 3 e7;ν= 0.3;ρ= 0.3/386;
创建一个模态平面应力模型,指定一个几何,并生成网格。
structuralmodel = createpde(结构=“modal-planestress”);gdm =[3、4 0;长度;长度;0,0,0;高度;高度);g = decsg (gdm,“S1 ',(“S1 ')');geometryFromEdges (structuralmodel g);
定义一个最大元素大小(五个元素通过梁厚度)。
hmax =身高/ 5;msh = generateMesh (structuralmodel Hmax = Hmax);
指定结构的属性和边界约束。
structuralProperties (structuralmodel YoungsModulus = E,…MassDensity =ρ,…PoissonsRatio =ν);structuralBC (structuralmodel = 4,约束=“固定”);
计算分析基本频率(赫兹)用梁理论。
我=身高^ 3/12;analyticalOmega1 = 3.516 * sqrt (E *我/(长度^ 4 *(ρ*高)))/(2 *π)
analyticalOmega1 = 126.9498
指定一个频率范围,包括一个分析计算频率和解决模型。
modalresults =解决(structuralmodel FrequencyRange = [0,1 e6])
modalresults = ModalStructuralResults属性:NaturalFrequencies: [32 x1双]ModeShapes: [1 x1 FEStruct]网:[1 x1 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属性:用户体验:[6511 x32双]uy: [6511 x32双]级:[6511 x32双)
画出y基频分量的解决方案。
pdeplot (modalresults.Mesh XYData = modalresults.ModeShapes.uy(: 1))标题([“第一模式频率”,…num2str (modalresults.NaturalFrequencies(1) /(2 *π)),“赫兹”])轴平等的
采用模态叠加法解决瞬态结构模型
解出瞬态响应的中心三维光束谐波负载下它的一个角落。
模态分析
创建一个三维的模态分析模型问题。
structuralmodel = createpde (“结构性”,“modal-solid”);
创建几何和包括在模型中。绘制几何图形和显示边缘和顶点标签。
通用= multicuboid (0.05, 0.003, 0.003);structuralmodel。几何=通用;pdegplot (structuralmodel“EdgeLabels”,“上”,…“VertexLabels”,“上”);视图(95 [5])
生成一个网格。
msh = generateMesh (structuralmodel);
指定杨氏模量、泊松比和材料的质量密度。
structuralProperties (structuralmodel“YoungsModulus”210 e9,…“PoissonsRatio”,0.3,…“MassDensity”,7800);
指定最小约束梁的一端,以防止刚体模式。例如,指定4边和顶点7是固定的边界。
structuralBC (structuralmodel“边缘”4“约束”,“固定”);structuralBC (structuralmodel“顶点”7“约束”,“固定”);
解决这个问题的频率范围从0到500000。推荐的方法是使用一个值略低于预期的最低频率。因此,使用-0.1而不是0。
Rm =解决(structuralmodel,“FrequencyRange”(-0.1,500000));
瞬态分析
开关瞬态模型的分析类型。
structuralmodel。AnalysisType =“transient-solid”;
角落上施加一个正弦力相反的约束边和顶点。
structuralBoundaryLoad (structuralmodel“顶点”5,…“力”(0,0,10),…“频率”,7600);
指定零初始位移和速度。
structuralIC (structuralmodel“速度”(0,0,0),…“位移”,(0,0,0));
指定的相对和绝对误差解算器。
structuralmodel.SolverOptions。RelativeTolerance = 1 e-5;structuralmodel.SolverOptions。AbsoluteTolerance = 1 e-9;
解决使用模态模型的结果。
tlist = linspace (0, 0.004,120);Rdm =解决(structuralmodel tlist,“ModalResults”Rm);
插入和情节在梁的中心位移。
intrpUdm = interpolateDisplacement (Rdm, 0, 0, 0.0015);情节(Rdm.SolutionTimes intrpUdm.uz)网格在包含(“时间”);ylabel (“梁位移的中心”)
扩张下悬臂梁的热负荷
找到一个三维悬臂梁的挠度下非均匀热负荷。指定的热负荷结构模型上使用瞬态热分析的解决方案相同的几何和网格。
瞬态热模型分析
创建一个瞬态热模型。
thermalmodel = createpde (“热”,“瞬态”);
创建和绘制几何。
通用= multicuboid (0.5, 0.1, 0.05);thermalmodel。几何=通用;pdegplot (thermalmodel“FaceLabels”,“上”,“FaceAlpha”,0.5)
生成一个网格。
网= generateMesh (thermalmodel);
指定材料的热性能。
thermalProperties (thermalmodel“ThermalConductivity”,5 e - 3,…“MassDensity”2.7 * 10 ^ (6),…“SpecificHeat”10);
指定恒定温度应用于左、右两端的梁。
thermalBC (thermalmodel“面子”3,“温度”,100);thermalBC (thermalmodel“面子”5,“温度”,0);
指定热源在整个几何。
internalHeatSource (thermalmodel 10);
设置初始温度。
thermalIC (thermalmodel 0);
解决模型。
tlist = [0:1e-4:2e-4];thermalresults =解决(thermalmodel tlist)
thermalresults = TransientThermalResults属性:温度:[3870 x3双]SolutionTimes: [0 1.0000 2.0000 e-04 e-04] XGradients: [3870 x3双]YGradients: [3870 x3双]ZGradients: [3870 x3双]网:[1 x1 FEMesh]
每个时间步的温度分布。
为n = 1:元素个数(thermalresults.SolutionTimes)图pdeplot3D (thermalmodel,…“ColorMapData”,…thermalresults.Temperature (:, n))标题([“温度时间= "…num2str (tlist (n)))) caxis (100 [0])结束
结构分析与热负荷
创建一个静态结构模型。
structuralmodel = createpde (“结构性”,“static-solid”);
包括热模型的几何一样。
structuralmodel。几何=通用;
使用相同的网格,用于获得热解决方案。
structuralmodel。网=网;
指定杨氏模量、泊松比和热膨胀系数。
structuralProperties (structuralmodel“YoungsModulus”1 e10,…“PoissonsRatio”,0.3,…“CTE”11.7 e-6);
应用一个固定边界条件面临5。
structuralBC (structuralmodel“面子”5,“约束”,“固定”);
应用一个身体负荷使用瞬态热模型的解决方案。默认情况下,structuralBodyLoad最后一次使用解决方案的步骤。
structuralBodyLoad (structuralmodel“温度”,thermalresults);
指定参考温度。
structuralmodel。ReferenceTemperature = 10;
解决结构模型。
thermalstressresults =解决(structuralmodel);
情节相对应的梁的变形形状的最后一步瞬态热模型的解决方案。
pdeplot3D (structuralmodel…“ColorMapData”,…thermalstressresults.Displacement.Magnitude,…“变形”,…thermalstressresults.Displacement)标题(“在解决方案时间=热膨胀”…num2str (tlist(结束)))caxis ([0 3 e - 3])
现在指定身体负荷作为热模型解决方案的步骤。金宝搏官方网站对于每一个身体负荷,解决结构模型和策划相应的梁的变形形状。
为n = 1:元素个数(thermalresults.SolutionTimes) structuralBodyLoad (structuralmodel,…“温度”,…thermalresults,…“步伐”n);thermalstressresults =解决(structuralmodel);图pdeplot3D (structuralmodel,…“ColorMapData”,…thermalstressresults.Displacement.Magnitude,…“变形”,…thermalstressresults.Displacement)标题(“热结果在解决时间= "…num2str (tlist (n)))) caxis ([0 3 e - 3])结束
解决稳态热模型
解决一个三维稳态热的问题。
关于这个问题创建一个热模型。
thermalmodel = createpde(热=“稳态”);
导入和图块几何。
importGeometry (thermalmodel“Block.stl”);pdegplot (thermalmodel FaceLabels =“上”,FaceAlpha = 0.5)轴平等的
指定材料属性。
thermalProperties (thermalmodel ThermalConductivity = 80);
100年施加一个恒定的温度°C的左边块(脸1)和300年的一个恒定的温度°C的右侧块(3)脸。所有其他面临默认是绝缘的。
thermalBC (thermalmodel,面对= 1,温度= 100);thermalBC (thermalmodel,面对= 3,温度= 300);
网格的几何和解决这个问题。
generateMesh (thermalmodel);thermalresults =解决(thermalmodel)
thermalresults = SteadyStateThermalResults属性:温度:[12691 x1双]XGradients: [12691 x1双]YGradients: [12691 x1双]ZGradients: [12691 x1双]网:[1 x1 FEMesh]
解算器发现节点位置的温度和温度梯度。要访问这些值,使用thermalresults.Temperature,thermalresults.XGradients,等等。例如,绘制温度节点的位置。
pdeplot3D (thermalresults.Mesh ColorMapData = thermalresults.Temperature)
解决瞬态热模型
解决一个二维瞬态热的问题。
关于这个问题创建一个瞬态热模型。
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、科幻、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);
施加一个恒定的温度为0°C的平方板。
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 =解决(thermalmodel tlist);
情节的解决方案与等温线用等高线图。
T = thermalresults.Temperature;msh = thermalresults.Mesh;要看更多有关憩苑pdeplot (, XYData = T(:, 10),轮廓=“上”ColorMap =“热”)
采用模态叠加法解决瞬态热模型
解决一个瞬态热的问题,首先获取模式形状为特定的衰减范围,然后利用模态叠加法。
模态分解
首先,创建一个模态热模型。
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、科幻、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);
施加一个恒定的温度为0°C的平方板。
thermalBC (thermalmodel“温度”0,“边缘”[1 2 7 8]);
设置初始温度为0°C。
thermalIC (thermalmodel 0);
生成网格。
generateMesh (thermalmodel);
计算热模型的固有模式衰减范围(100、10000) 。
RModal =解决(thermalmodel,“DecayRange”(100、10000))
RModal = ModalThermalResults属性:DecayRates: x1双[164]ModeShapes: [1481 x164双]ModeType:“固有模式”网:[1 x1 FEMesh]
瞬态分析
了解模式形状,您现在可以用模态叠加法解决瞬态热的问题。首先,开关瞬态的热模型分析类型。
thermalmodel。AnalysisType =“瞬态”;
这个问题的动力学是非常快。温度大概在0.1秒到达稳定状态。捕捉最活跃的动力的一部分,设置时间的解决方案logspace (-1100)。这个命令返回100对数间隔的解决方案*在0.01和0.1之间。
tlist = logspace (2, 1, 10);
解方程。
Rtransient =解决(thermalmodel tlist,“ModalResults”,RModal);
情节的解决方案与等温线用等高线图。
T = Rtransient.Temperature;pdeplot (thermalmodel“XYData”T(:,结束),…“轮廓”,“上”,…“ColorMap”,“热”)
适当的正交分解的快照
获得POD模式的线性瞬态的热模型使用几个实例解决方案(快照)。
创建一个瞬态热模型。
thermalmodel = createpde (“热”,“瞬态”);
创建一个单位平方几何,包括在模型中。
geometryFromEdges (thermalmodel @squareg);
画出几何图形,显示标签边缘。
pdegplot (thermalmodel“EdgeLabels”,“上”)xlim ([-1.1 - 1.1]) ylim ([-1.1 - 1.1])
指定导热系数、质量密度和比热的材料。
thermalProperties (thermalmodel“ThermalConductivity”,400,…“MassDensity”,1300,…“SpecificHeat”,600);
设置温度的右边缘One hundred.。
thermalBC (thermalmodel“边缘”2,“温度”,100);
设置一个初始值0的温度。
thermalIC (thermalmodel 0);
生成一个网格。
generateMesh (thermalmodel);
解决三种不同值的热源模型和收集快照。
tlist = 0:10:600;snapShotIDs = [1:10 59 60 61];Tmatrix = [];heatVariation = (10000 15000 20000);为q = heatVariation internalHeatSource (thermalmodel q);结果=解决(thermalmodel tlist);Tmatrix = [Tmatrix results.Temperature (:, snapShotIDs)];结束
开关模态的热模型分析类型。
thermalmodel。AnalysisType =“模态”;
计算圆荚体模式。
RModal =解决(thermalmodel,“快照”Tmatrix)
RModal = ModalThermalResults属性:DecayRates: x1双[6]ModeShapes: [1541 x6双]SnapshotsAverage: x1双[1541]ModeType:“PODModes”网:[1 x1 FEMesh]
解决二维静电分析模型
解决电磁问题和找到电势场分布的二维几何代表一盘有一个洞。
创建一个静电电磁模型分析。
emagmodel = createpde(电磁=“静电”);
导入和绘制几何图形代表一盘有一个洞。
importGeometry (emagmodel“PlateHolePlanar.stl”);pdegplot (emagmodel EdgeLabels =“上”)
指定在SI单位制真空介电常数的值。
emagmodel。VacuumPermittivity = 8.8541878128 e-12;
指定材料的相对介电常数。
electromagneticProperties (emagmodel RelativePermittivity = 1);
应用框架边缘上的电压边界条件矩形和圆。
electromagneticBC (emagmodel电压= 0,边缘= 1:4);electromagneticBC (emagmodel、电压= 1000,边缘= 5);
指定整个几何的电荷密度。
electromagneticSource (emagmodel ChargeDensity = 5 e-9);
生成网格。
generateMesh (emagmodel);
解决模型。
R =解决(emagmodel)
R = ElectrostaticResults属性:电位:x1双[1218]电场:[1 x1 FEStruct] ElectricFluxDensity: [1 x1 FEStruct]网:[1 x1 FEMesh]
情节的电势和字段。
pdeplot (R.Mesh XYData = R.ElectricPotential,…FlowData = [R.ElectricField.Ex…R.ElectricField.Ey])轴平等的
解决三维静磁分析模型
解决三维电磁问题在几何代表一个盘子在它的中心有一个洞。情节产生的潜力和磁场分布。
创建一个电磁模型静磁分析。
emagmodel = createpde (“电磁”,“静磁”);
导入和绘制几何图形代表一盘有一个洞。
importGeometry (emagmodel“PlateHoleSolid.stl”);pdegplot (emagmodel“FaceLabels”,“上”,“FaceAlpha”,0.3)
指定真空磁导率值在SI单位制。
emagmodel。VacuumPermeability = 1.2566370614 e-6;
指定材料的相对磁导率。
electromagneticProperties (emagmodel“RelativePermeability”,5000);
应用磁势接壤边界条件的侧脸,脸上的洞。
electromagneticBC (emagmodel“磁势”(0,0,0),“面子”、三6);electromagneticBC (emagmodel“磁势”(0,0,0.01),“面子”7);
指定整个几何的电流密度。
electromagneticSource (emagmodel“CurrentDensity”[0,0,0.5]);
生成网格。
generateMesh (emagmodel);
解决模型。
R =解决(emagmodel)
R = MagnetostaticResults属性:磁势:[1 x1 FEStruct]磁场:[1 x1 FEStruct] MagneticFluxDensity: [1 x1 FEStruct]网:[1 x1 FEMesh]
画出z磁势的分量。
pdeplot3D (emagmodel“ColormapData”R.MagneticPotential.Az)
情节的磁场。
pdeplot3D (emagmodel“FlowData”,(R.MagneticField.Hx…R.MagneticField.Hy…R.MagneticField.Hz])
解决三维直流电传导模型
解决直流导电问题在几何代表3 d板在它的中心有一个洞。情节的电势和电流密度的组件。
创建一个直流传导电磁模型分析。
emagmodel = createpde (“电磁”,“传导”);
导入和阴谋一个几何图形代表一盘有一个洞。
通用= importGeometry (emagmodel,“PlateHoleSolid.stl”);pdegplot(通用、“FaceLabels”,“上”,“FaceAlpha”,0.3)
指定材料的电导率。
electromagneticProperties (emagmodel“导”6 e4);
应用电压边界条件在左,右,前,和底部的面孔。
electromagneticBC (emagmodel“电压”0,“面子”、三6);
指定表面电流密度在脸上接壤的洞。
electromagneticBC (emagmodel“SurfaceCurrentDensity”,100,“面子”7);
生成网格。
generateMesh (emagmodel);
解决模型。
R =解决(emagmodel)
R = ConductionResults属性:电位:x1双[4359]电场:[1 x1 FEStruct] CurrentDensity: [1 x1 FEStruct]网:[1 x1 FEMesh]
情节的电势。
图pdeplot3D (emagmodel,“ColorMapData”R.ElectricPotential)
画出x电流密度的分量。
图pdeplot3D (emagmodel,“ColorMapData”R.CurrentDensity.Jx)标题(“x分量电流密度”)
画出y电流密度的分量。
图pdeplot3D (emagmodel,“ColorMapData”R.CurrentDensity.Jy)标题(“电流密度分量”)
画出z电流密度的分量。
图pdeplot3D (emagmodel,“ColorMapData”R.CurrentDensity.Jz)标题(“z分量的电流密度”)
直流导电溶液作为静磁电流密度模型
使用解决方案获得通过执行指定直流导电分析电流密度对静磁模型。
创建一个直流传导电磁模型分析。
emagmodel = createpde (“电磁”,“传导”);
导入和阴谋一个几何图形代表一盘有一个洞。
通用= importGeometry (emagmodel,“PlateHoleSolid.stl”);pdegplot(通用、“FaceLabels”,“上”,“FaceAlpha”,0.3)
指定材料的电导率。
electromagneticProperties (emagmodel“导”6 e4);
应用电压边界条件在左,右,前,和底部的面孔。
electromagneticBC (emagmodel“电压”0,“面子”、三6);
指定表面电流密度在脸上接壤的洞。
electromagneticBC (emagmodel“SurfaceCurrentDensity”,100,“面子”7);
生成网格。
generateMesh (emagmodel);
解决模型。
R =解决(emagmodel);
模型的分析类型改为静磁。
emagmodel。AnalysisType =“静磁”;
这个模型已经生成的二次网你直流导电分析。三维静磁模型,网格必须是线性的。生成一个新的线性网格。的generateMesh函数创建一个线性网格默认情况下如果模型是三维静磁。
generateMesh (emagmodel);
指定真空磁导率值在SI单位制。
emagmodel。VacuumPermeability = 1.2566370614 e-6;
指定材料的相对磁导率。
electromagneticProperties (emagmodel“RelativePermeability”,5000);
应用磁势接壤边界条件的侧脸,脸上的洞。
electromagneticBC (emagmodel“磁势”(0,0,0),“面子”、三6);electromagneticBC (emagmodel“磁势”(0,0,0.01),“面子”7);
指定使用直流电流密度对整个几何传导的解决方案。
electromagneticSource (emagmodel“CurrentDensity”,R);
解决模型。
Rmagnetostatic =解决(emagmodel);
画出x- - -z制造磁性的潜力。
pdeplot3D (emagmodel“ColormapData”Rmagnetostatic.MagneticPotential.Ax)
pdeplot3D (emagmodel“ColormapData”Rmagnetostatic.MagneticPotential.Az)
解决二维静磁模型与永久磁铁
解决一个铜的静磁模型与一个永久的钕磁铁的中心广场。
创建单元正方形和一个圆的几何中心。
L = 0.8;r = 0.25;平方= [3 4 - L L L - L - L - L L L) ';中国保监会= r (1 0 0 0 0 0 0 0 0) ';gd =(平方,保监会);科幻小说=“平方+保监会”;ns = char (“平方”,“保监会”);ns = ns ';g = decsg (gd、科幻、ns);
情节的脸和边缘的几何标签。
pdegplot (g,“FaceLabels”,“上”,“EdgeLabels”,“上”)
创建一个静磁模型,包括几何模型。
emagmodel = createpde (“电磁”,“静磁”);geometryFromEdges (emagmodel g);
指定真空磁导率值在SI单位制。
emagmodel。VacuumPermeability = 1.2566370614 e-6;
指定的铜的相对渗透率广场。
electromagneticProperties (emagmodel“面子”,1…“RelativePermeability”1);
指定的钕圆的相对渗透率。
electromagneticProperties (emagmodel“面子”2,…“RelativePermeability”,1.05);
指定钕磁铁的磁化强度。
M = 1 e6;
指定磁化圆形脸的正面x方向。为二维磁化模型是一个列向量的两个元素。
dir = (1, 0);electromagneticSource (emagmodel“面子”2,“磁化”,M * dir);
应用磁势在广场边缘构造边界条件。
electromagneticBC (emagmodel“边缘”1:4,“磁势”,0);
生成的网格细网格圆的边缘附近。
generateMesh (emagmodel“对冲”{8,0.007});图pdemesh (emagmodel)
解决模型,找到产生的磁场B和H。在这里, ,在那里 绝对磁导率的材料, 真空磁导率, 磁化。
R =解决(emagmodel);Bmag =√R.MagneticFluxDensity.Bx。^ 2 + R.MagneticFluxDensity.By。^ 2);Hmag =√R.MagneticField.Hx。^ 2 + R.MagneticField.Hy。^ 2);
情节的磁场B。
图pdeplot (emagmodel,“XYData”Bmag,…“FlowData”,(R.MagneticFluxDensity.Bx…R.MagneticFluxDensity.By])
情节的磁场H。
图pdeplot (emagmodel,“XYData”Hmag,…“FlowData”,(R.MagneticField.HxR.MagneticField.Hy])
解决二维谐波电磁模型
对于一个电磁谐波分析问题,找到了x- - -y制造的电场。解决这个问题在一个域组成的正方形和一个圆孔。
创建一个电磁谐波分析模型。
emagmodel = createpde (“电磁”,“谐”);
定义一个圆广场,把它们在一个矩阵,并创建一套公式,减去的圆广场。
平方= [3、4、5、5、5、5、5、5、5、5]';C = [1, 0, 0, 1] ';C = [C; 0(长度(平方)——(C), 1)];通用=[平方,C];科幻小说=“SQ-C”;
创建几何。
ns = char (“平方”,“C”);ns = ns ';g = decsg(通用汽车、科幻、ns);
包括几何模型和绘制几何标签。
geometryFromEdges (emagmodel g);pdegplot (emagmodel“EdgeLabels”,“上”)xlim ([-5.5 - 5.5]) ylim ([-5.5 - 5.5])
真空介电常数和磁导率的值指定为1。
emagmodel。VacuumPermittivity = 1;emagmodel。VacuumPermeability = 1;
指定相对介电常数、相对渗透率和导电性的材料。
electromagneticProperties (emagmodel“RelativePermittivity”,1…“RelativePermeability”,1…“导”,0);
应用吸收边界条件的厚度2边缘的广场上。使用默认的衰减率吸收区域。
electromagneticBC (emagmodel“边缘”1:4,…“FarField”,“吸收”,…“厚度”2);
指定一个洞的边缘电场。
E = @(位置、状态)(1,0)* exp(1 * 2 *π* location.y);electromagneticBC (emagmodel“边缘”8,“电场”,E);
生成一个网格。
generateMesh (emagmodel“Hmax”、1/2 ^ 3);
解决模型的频率 。
结果=解决(emagmodel,“频率”,2 *π);
情节真实的一部分x分量产生的电场。
图pdeplot (emagmodel,“XYData”真正的(result.ElectricField.Ex));标题(“真正的电场分量的一部分”)
情节真实的一部分y分量产生的电场。
图pdeplot (emagmodel,“XYData”真正的(result.ElectricField.Ey));标题(“真正的电场分量的一部分”)
输入参数
输出参数
提示
当你使用模态分析结果来解决一个瞬态结构动力学模型,
modalresults在偏微分方程参数必须创建工具箱™R2019a或更新。频率响应模型的阻尼,结果是复杂的。使用等功能
腹肌和角获得实值结果,如大小和相位。
