解决gydF4y2Ba

解决传热或结构分析的问题gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

structuralStaticResultsgydF4y2Ba=解决(gydF4y2BastructuralStaticgydF4y2Ba)gydF4y2Ba返回解决静态结构分析模型中表示gydF4y2BastructuralStaticgydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

structuralModalResultsgydF4y2Ba=解决(gydF4y2BastructuralModalgydF4y2Ba“FrequencyRange”,gydF4y2Ba(ω,₂)gydF4y2Ba)gydF4y2Ba返回所有的模态分析模型的解决方案模式的频率范围gydF4y2Ba(ω,₂)gydF4y2Ba。定义gydF4y2BaωgydF4y2Ba比预期的最低频率和略小gydF4y2Ba₂gydF4y2Ba比最高的期望频率。例如,如果最低的期望频率为零,然后用一个小的负面价值gydF4y2BaωgydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

structuralTransientResultsgydF4y2Ba=解决(gydF4y2BastructuralTransientgydF4y2Ba,gydF4y2BatlistgydF4y2Ba)gydF4y2Ba返回解决瞬态结构动力学模型中表示gydF4y2BastructuralTransientgydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

structuralFrequencyResponseResultsgydF4y2Ba=解决(gydF4y2BastructuralFrequencyResponsegydF4y2Ba,gydF4y2Ba弗利斯特雨gydF4y2Ba)gydF4y2Ba返回表示的频率响应模型的解决方案gydF4y2BastructuralFrequencyResponsegydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

structuralTransientResultsgydF4y2Ba=解决(gydF4y2BastructuralTransientgydF4y2Ba,gydF4y2BatlistgydF4y2Ba“ModalResults”,gydF4y2BamodalresultsgydF4y2Ba)gydF4y2Ba和gydF4y2BastructuralFrequencyResponseResultsgydF4y2Ba=解决(gydF4y2BastructuralFrequencyResponsegydF4y2Ba,gydF4y2Ba弗利斯特雨gydF4y2Ba“ModalResults”,gydF4y2BamodalresultsgydF4y2Ba)gydF4y2Ba解决瞬态和频率响应结构模型,分别采用模态叠加法加快计算速度。首先,执行计算固有频率和振型模态分析在一个特定的频率范围。然后,使用这个语法来调用模态叠加法。结果的准确性依赖于模式的模态分析结果。gydF4y2Ba

例子gydF4y2Ba

thermalSteadyStateResultsgydF4y2Ba=解决(gydF4y2BathermalSteadyStategydF4y2Ba)gydF4y2Ba返回解决稳态热模型中表示gydF4y2BathermalSteadyStategydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

thermalTransientResultsgydF4y2Ba=解决(gydF4y2BathermalTransientgydF4y2Ba,gydF4y2BatlistgydF4y2Ba)gydF4y2Ba返回解决瞬态热模型中表示gydF4y2BathermalTransientgydF4y2Ba在《纽约时报》gydF4y2BatlistgydF4y2Ba。gydF4y2Ba

例子gydF4y2Ba

全部折叠gydF4y2Ba

解决一个三维稳态热的问题。gydF4y2Ba

关于这个问题创建一个热模型。gydF4y2Ba

thermalmodel = createpde (gydF4y2Ba“热”gydF4y2Ba);gydF4y2Ba

导入和图块几何。gydF4y2Ba

importGeometry (thermalmodelgydF4y2Ba“Block.stl”gydF4y2Ba);pdegplot (thermalmodelgydF4y2Ba“FaceLabel”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba,gydF4y2Ba“FaceAlpha”gydF4y2Ba,0.5)轴gydF4y2Ba平等的gydF4y2Ba

指定材料属性。gydF4y2Ba

thermalProperties (thermalmodelgydF4y2Ba“ThermalConductivity”gydF4y2Ba,80);gydF4y2Ba

施加一个恒定的温度gydF4y2Ba 1gydF4y2Ba 0gydF4y2Ba 0gydF4y2Ba ∘gydF4y2Ba CgydF4y2Ba 左边的块(脸1)和一个恒定的温度gydF4y2Ba 3gydF4y2Ba 0gydF4y2Ba 0gydF4y2Ba ∘gydF4y2Ba CgydF4y2Ba 右边的块(3)脸。所有其他面临默认是绝缘的。gydF4y2Ba

thermalBC (thermalmodelgydF4y2Ba“脸”gydF4y2Ba,1gydF4y2Ba“温度”gydF4y2Ba,100);thermalBC (thermalmodelgydF4y2Ba“脸”gydF4y2Ba3,gydF4y2Ba“温度”gydF4y2Ba,300);gydF4y2Ba

网格的几何和解决这个问题。gydF4y2Ba

generateMesh (thermalmodel);thermalresults =解决(thermalmodel)gydF4y2Ba
thermalresults = SteadyStateThermalResults属性:温度:[12691 x1双]XGradients: [12691 x1双]YGradients: [12691 x1双]ZGradients: [12691 x1双]网:[1 x1 FEMesh]gydF4y2Ba

解算器发现节点位置的温度和温度梯度。要访问这些值,使用gydF4y2Bathermalresults.TemperaturegydF4y2Ba,gydF4y2Bathermalresults.XGradientsgydF4y2Ba,等等。例如,绘制温度节点的位置。gydF4y2Ba

pdeplot3D (thermalmodelgydF4y2Ba“ColorMapData”gydF4y2Bathermalresults.Temperature)gydF4y2Ba

解决一个二维瞬态热的问题。gydF4y2Ba

关于这个问题创建一个瞬态热模型。gydF4y2Ba

thermalmodel = createpde (gydF4y2Ba“热”gydF4y2Ba,gydF4y2Ba瞬态的gydF4y2Ba);gydF4y2Ba

创建几何和包括在模型中。gydF4y2Ba

于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”gydF4y2Ba;ns = char (gydF4y2Ba“于SQ1”gydF4y2Ba,gydF4y2Ba“D1”gydF4y2Ba);ns = ns ';dl = decsg (gd、科幻、ns);geometryFromEdges (thermalmodel dl);pdegplot (thermalmodelgydF4y2Ba“EdgeLabels”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba,gydF4y2Ba“FaceLabels”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba)xlim ([-1.5 - 4.5]) ylim([-0.5 - 3.5])轴gydF4y2Ba平等的gydF4y2Ba

正方形,分配这些热性能:gydF4y2Ba

  • 导热系数是gydF4y2Ba 1gydF4y2Ba 0gydF4y2Ba WgydF4y2Ba /gydF4y2Ba (gydF4y2Ba 米gydF4y2Ba ∘gydF4y2Ba CgydF4y2Ba )gydF4y2Ba

  • 质量密度是gydF4y2Ba 2gydF4y2Ba kgydF4y2Ba ggydF4y2Ba /gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba

  • 比热是gydF4y2Ba 0gydF4y2Ba 。gydF4y2Ba 1gydF4y2Ba JgydF4y2Ba /gydF4y2Ba (gydF4y2Ba kgydF4y2Ba ggydF4y2Ba ∘gydF4y2Ba CgydF4y2Ba )gydF4y2Ba

thermalProperties (thermalmodelgydF4y2Ba“ThermalConductivity”gydF4y2Ba10gydF4y2Ba…gydF4y2Ba“MassDensity”gydF4y2Ba2,gydF4y2Ba…gydF4y2Ba“SpecificHeat”gydF4y2Ba,0.1,gydF4y2Ba…gydF4y2Ba“脸”gydF4y2Ba1);gydF4y2Ba

钻石的地区分配这些热性能:gydF4y2Ba

  • 导热系数是gydF4y2Ba 2gydF4y2Ba WgydF4y2Ba /gydF4y2Ba (gydF4y2Ba 米gydF4y2Ba ∘gydF4y2Ba CgydF4y2Ba )gydF4y2Ba

  • 质量密度是gydF4y2Ba 1gydF4y2Ba kgydF4y2Ba ggydF4y2Ba /gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba

  • 比热是gydF4y2Ba 0gydF4y2Ba 。gydF4y2Ba 1gydF4y2Ba JgydF4y2Ba /gydF4y2Ba (gydF4y2Ba kgydF4y2Ba ggydF4y2Ba ∘gydF4y2Ba CgydF4y2Ba )gydF4y2Ba

thermalProperties (thermalmodelgydF4y2Ba“ThermalConductivity”gydF4y2Ba2,gydF4y2Ba…gydF4y2Ba“MassDensity”gydF4y2Ba,1gydF4y2Ba…gydF4y2Ba“SpecificHeat”gydF4y2Ba,0.1,gydF4y2Ba…gydF4y2Ba“脸”gydF4y2Ba2);gydF4y2Ba

假设菱形区域密度的热源gydF4y2Ba 4gydF4y2Ba WgydF4y2Ba /gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba 。gydF4y2Ba

internalHeatSource (thermalmodel 4gydF4y2Ba“脸”gydF4y2Ba2);gydF4y2Ba

施加一个恒定的温度gydF4y2Ba 0gydF4y2Ba ∘gydF4y2Ba CgydF4y2Ba 广场的两侧板。gydF4y2Ba

thermalBC (thermalmodelgydF4y2Ba“温度”gydF4y2Ba0,gydF4y2Ba“边缘”gydF4y2Ba[1 2 7 8]);gydF4y2Ba

设置初始温度gydF4y2Ba 0gydF4y2Ba ∘gydF4y2Ba CgydF4y2Ba 。gydF4y2Ba

thermalIC (thermalmodel 0);gydF4y2Ba

网格几何。gydF4y2Ba

generateMesh (thermalmodel);gydF4y2Ba

这个问题的动力学是非常快。温度达到稳定状态在大约0.1秒。捕捉有趣的部分动态,设置时间的解决方案gydF4y2Balogspace (2, 1, 10)gydF4y2Ba。这个命令返回10对数间隔的解决方案*在0.01和0.1之间。gydF4y2Ba

tlist = logspace (2, 1, 10);gydF4y2Ba

解方程。gydF4y2Ba

thermalresults =解决(thermalmodel tlist)gydF4y2Ba
thermalresults = TransientThermalResults属性:温度:[1481 x10双]SolutionTimes: [1 x10双]XGradients: [1481 x10双]YGradients: [1481 x10双]ZGradients:[]网:[1 x1 FEMesh]gydF4y2Ba

情节的解决方案与等温线用等高线图。gydF4y2Ba

T = thermalresults.Temperature;pdeplot (thermalmodelgydF4y2Ba“XYData”gydF4y2BaT (10),gydF4y2Ba“轮廓”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba,gydF4y2Ba“ColorMap”gydF4y2Ba,gydF4y2Ba“热”gydF4y2Ba)gydF4y2Ba

解决一个静态结构模型代表一个双金属电缆张力。gydF4y2Ba

创建一个静态结构模型为解决固体(3 d)问题。gydF4y2Ba

structuralmodel = createpde (gydF4y2Ba“结构”gydF4y2Ba,gydF4y2Ba“static-solid”gydF4y2Ba);gydF4y2Ba

创建几何和包括在模型中。绘制几何图形。gydF4y2Ba

通用=多缸(0.05 [0.01,0.015]);structuralmodel。几何=通用;pdegplot (structuralmodelgydF4y2Ba“FaceLabels”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba,gydF4y2Ba“CellLabels”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba,gydF4y2Ba“FaceAlpha”gydF4y2Ba,0.5)gydF4y2Ba

指定杨氏模量和泊松比为每个金属。gydF4y2Ba

structuralProperties (structuralmodelgydF4y2Ba“细胞”gydF4y2Ba,1gydF4y2Ba“YoungsModulus”gydF4y2Ba110 e9,gydF4y2Ba…gydF4y2Ba“PoissonsRatio”gydF4y2Ba,0.28);structuralProperties (structuralmodelgydF4y2Ba“细胞”gydF4y2Ba2,gydF4y2Ba“YoungsModulus”gydF4y2Ba210 e9,gydF4y2Ba…gydF4y2Ba“PoissonsRatio”gydF4y2Ba,0.3);gydF4y2Ba

指定面临1和4是固定的边界。gydF4y2Ba

structuralBC (structuralmodelgydF4y2Ba“脸”gydF4y2Ba(1、4),gydF4y2Ba“约束”gydF4y2Ba,gydF4y2Ba“固定”gydF4y2Ba);gydF4y2Ba

指定表面引力面临2和5。gydF4y2Ba

structuralBoundaryLoad (structuralmodelgydF4y2Ba“脸”gydF4y2Ba(2、5),gydF4y2Ba“SurfaceTraction”gydF4y2Ba,(0,0,100));gydF4y2Ba

生成一个网格,并解决问题。gydF4y2Ba

generateMesh (structuralmodel);structuralresults =解决(structuralmodel)gydF4y2Ba
structuralresults = StaticStructuralResults属性:位移:[1 x1 struct]应变:[1 x1 struct]压力:[1 x1 struct] VonMisesStress: [22281 x1双]网:[1 x1 FEMesh]gydF4y2Ba

解决发现的价值位移、应力、应变,·冯·米塞斯应力节点位置。要访问这些值,使用gydF4y2Bastructuralresults.DisplacementgydF4y2Ba,gydF4y2Bastructuralresults.StressgydF4y2Ba,等等。位移、应力和应变值的节点的位置作为结构数组返回字段代表他们的组件。gydF4y2Ba

structuralresults.DisplacementgydF4y2Ba
ans =gydF4y2Ba结构体字段:gydF4y2Ba用户体验(22281 x1双):uy: [22281 x1双]是乌斯:[22281 x1双]级:[22281 x1双]gydF4y2Ba
structuralresults.StressgydF4y2Ba
ans =gydF4y2Ba结构体字段:gydF4y2Basxx (22281 x1双):syy: [22281 x1双]szz: [22281 x1双]syz: [22281 x1双]sxz: [22281 x1双]sxy: [22281 x1双]gydF4y2Ba
structuralresults.StraingydF4y2Ba
ans =gydF4y2Ba结构体字段:gydF4y2Baexx (22281 x1双):eyy: [22281 x1双]ezz: [22281 x1双]eyz: [22281 x1双]exz: [22281 x1双]exy: [22281 x1双]gydF4y2Ba

情节的变形形状gydF4y2BazgydF4y2Ba分的正常压力。gydF4y2Ba

pdeplot3D (structuralmodelgydF4y2Ba“ColorMapData”gydF4y2Bastructuralresults.Stress.szz,gydF4y2Ba…gydF4y2Ba“变形”gydF4y2Bastructuralresults.Displacement)gydF4y2Ba

解决瞬态响应的薄3 d板下谐波负荷中心。gydF4y2Ba

创建一个瞬态动力学模型的三维问题。gydF4y2Ba

structuralmodel = createpde (gydF4y2Ba“结构”gydF4y2Ba,gydF4y2Ba“transient-solid”gydF4y2Ba);gydF4y2Ba

创建几何和包括在模型中。绘制几何图形。gydF4y2Ba

通用= multicuboid ([5, 0.05), 5, 0.05, 0.01);structuralmodel.Geometry =通用;pdegplot (structuralmodelgydF4y2Ba“FaceLabels”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba,gydF4y2Ba“FaceAlpha”gydF4y2Ba,0.5)gydF4y2Ba

放大的脸小板上的标签的中心。gydF4y2Ba

图pdegplot (structuralmodel,gydF4y2Ba“FaceLabels”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba,gydF4y2Ba“FaceAlpha”gydF4y2Ba,0.25)轴([-0.2 0.2 -0.2 0.2 -0.1 0.1])gydF4y2Ba

指定杨氏模量,泊松比,质量密度的材料。gydF4y2Ba

structuralProperties (structuralmodelgydF4y2Ba“YoungsModulus”gydF4y2Ba210 e9,gydF4y2Ba…gydF4y2Ba“PoissonsRatio”gydF4y2Ba,0.3,gydF4y2Ba…gydF4y2Ba“MassDensity”gydF4y2Ba,7800);gydF4y2Ba

指定所有的脸瘦的外围3 d板是固定的边界。gydF4y2Ba

structuralBC (structuralmodelgydF4y2Ba“约束”gydF4y2Ba,gydF4y2Ba“固定”gydF4y2Ba,gydF4y2Ba“脸”gydF4y2Ba,8);gydF4y2Ba

应用一个正弦压力载荷的小脸上在板的中心。gydF4y2Ba

structuralBoundaryLoad (structuralmodelgydF4y2Ba“脸”gydF4y2Ba12gydF4y2Ba“压力”gydF4y2Ba5 e7,gydF4y2Ba“频率”gydF4y2Ba25);gydF4y2Ba

生成一个网格和线性元素。gydF4y2Ba

generateMesh (structuralmodelgydF4y2Ba“GeometricOrder”gydF4y2Ba,gydF4y2Ba“线性”gydF4y2Ba,gydF4y2Ba“Hmax”gydF4y2Ba,0.2);gydF4y2Ba

指定零初始位移和速度。gydF4y2Ba

structuralIC (structuralmodelgydF4y2Ba“位移”gydF4y2Ba(0,0,0),gydF4y2Ba“速度”gydF4y2Ba,(0,0,0));gydF4y2Ba

解决模型。gydF4y2Ba

tlist = linspace (0, 1300);structuralresults =解决(structuralmodel tlist)gydF4y2Ba
structuralresults = TransientStructuralResults属性:位移:[1 x1 struct]速度:[1 x1 struct]加速度:[1 x1 struct] SolutionTimes:[1×双]网:[1 x1 FEMesh]gydF4y2Ba

解算器发现的价值位移,速度,加速度的节点位置。要访问这些值,使用gydF4y2Bastructuralresults.DisplacementgydF4y2Ba,gydF4y2Bastructuralresults.VelocitygydF4y2Ba,等等。位移、速度和加速度值作为结构数组返回字段代表他们的组件。gydF4y2Ba

structuralresults.DisplacementgydF4y2Ba
ans =gydF4y2Ba结构体字段:gydF4y2Ba用户体验:[1873 x300双]uy: [1873 x300双]是乌斯:[1873 x300双]级:[1873 x300双]gydF4y2Ba
structuralresults.VelocitygydF4y2Ba
ans =gydF4y2Ba结构体字段:gydF4y2Bavx(1873×双):v: [1873 x300双]vz: [1873 x300双]级:[1873 x300双]gydF4y2Ba
structuralresults.AccelerationgydF4y2Ba
ans =gydF4y2Ba结构体字段:gydF4y2Baax: [1873 x300双]ay: [1873 x300双]阿兹:[1873 x300双]级:[1873 x300双]gydF4y2Ba

找到最基本的二维悬臂梁(最低)模式,假设患病率的平面应力状态。gydF4y2Ba

指定以下梁的几何和结构属性,以及一个平面应力单元厚度。gydF4y2Ba

长度= 5;身高= 0.1;E = 3 e7;ν= 0.3;ρ= 0.3/386;gydF4y2Ba

创建一个模型平面应力模型,指定一个几何,并生成网格。gydF4y2Ba

structuralmodel = createpde (gydF4y2Ba“结构”gydF4y2Ba,gydF4y2Ba“modal-planestress”gydF4y2Ba);gdm =[3、4 0;长度;长度;0,0,0;高度;高度);g = decsg (gdm,gydF4y2Ba“S1 'gydF4y2Ba,(gydF4y2Ba“S1 'gydF4y2Ba)');geometryFromEdges (structuralmodel g);gydF4y2Ba

定义一个最大元素大小(五个元素通过梁厚度)。gydF4y2Ba

hmax =身高/ 5;msh = generateMesh (structuralmodel,gydF4y2Ba“Hmax”gydF4y2Ba,hmax);gydF4y2Ba

指定结构的属性和边界约束。gydF4y2Ba

structuralProperties (structuralmodelgydF4y2Ba“YoungsModulus”gydF4y2Ba,E,gydF4y2Ba…gydF4y2Ba“MassDensity”gydF4y2Baρ,gydF4y2Ba…gydF4y2Ba“PoissonsRatio”gydF4y2Baν);structuralBC (structuralmodelgydF4y2Ba“边缘”gydF4y2Ba4gydF4y2Ba“约束”gydF4y2Ba,gydF4y2Ba“固定”gydF4y2Ba);gydF4y2Ba

计算分析基本频率(赫兹)用梁理论。gydF4y2Ba

我=身高^ 3/12;analyticalOmega1 = 3.516 * sqrt (E *我/(长度^ 4 *(ρ*高)))/(2 *π)gydF4y2Ba
analyticalOmega1 = 126.9498gydF4y2Ba

指定一个频率范围,包括一个分析计算频率和解决模型。gydF4y2Ba

modalresults =解决(structuralmodel,gydF4y2Ba“FrequencyRange”gydF4y2Ba(0,1 e6))gydF4y2Ba
modalresults = ModalStructuralResults属性:NaturalFrequencies: [32 x1双]ModeShapes: [1 x1 struct]网:[1 x1 FEMesh]gydF4y2Ba

解算器发现固有频率和模态位移值在节点的位置。要访问这些值,使用gydF4y2Bamodalresults.NaturalFrequenciesgydF4y2Ba和gydF4y2Bamodalresults.ModeShapesgydF4y2Ba。gydF4y2Ba

modalresults.NaturalFrequencies /(2 *π)gydF4y2Ba
ans =gydF4y2Ba32×1gydF4y2Ba10gydF4y2Ba5gydF4y2Ba×0.0013 0.0079 0.0222 0.0433 0.0711 0.0983 0.1055 0.1462 0.1930 0.2455⋮gydF4y2Ba
modalresults.ModeShapesgydF4y2Ba
ans =gydF4y2Ba结构体字段:gydF4y2Ba用户体验:[6511 x32双]uy: [6511 x32双]gydF4y2Ba

画出gydF4y2BaygydF4y2Ba基频分量的解决方案。gydF4y2Ba

pdeplot (structuralmodelgydF4y2Ba“XYData”gydF4y2Bamodalresults.ModeShapes.uy(: 1)标题([gydF4y2Ba“第一模式频率”gydF4y2Ba,gydF4y2Ba…gydF4y2Banum2str (modalresults.NaturalFrequencies(1) /(2 *π)),gydF4y2Ba“赫兹”gydF4y2Ba])轴gydF4y2Ba平等的gydF4y2Ba

执行音叉的频率响应分析。gydF4y2Ba

首先,创建一个结构模型的模态分析固体音叉。gydF4y2Ba

模型= createpde (gydF4y2Ba“结构”gydF4y2Ba,gydF4y2Ba“frequency-solid”gydF4y2Ba);gydF4y2Ba

导入音叉几何。gydF4y2Ba

importGeometry(模型,gydF4y2Ba“TuningFork.stl”gydF4y2Ba);gydF4y2Ba

指定杨氏模量,泊松比,质量密度模型线性弹性材料的行为。指定所有物理性质的单位一致。gydF4y2Ba

structuralProperties(模型,gydF4y2Ba“YoungsModulus”gydF4y2Ba210 e9,gydF4y2Ba…gydF4y2Ba“PoissonsRatio”gydF4y2Ba,0.3,gydF4y2Ba…gydF4y2Ba“MassDensity”gydF4y2Ba,8000);gydF4y2Ba

识别面孔应用边界约束和载荷通过绘制的几何面标签。gydF4y2Ba

图(gydF4y2Ba“单位”gydF4y2Ba,gydF4y2Ba“归一化”gydF4y2Ba,gydF4y2Ba“outerposition”gydF4y2Ba,(0 0 1 1))pdegplot(模型,gydF4y2Ba“FaceLabels”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba15)视图(-50)标题gydF4y2Ba“几何与脸标签”gydF4y2Ba

施加足够的边界限制,以防止刚体运动应用载荷作用下。通常,您举行一个音叉用手或挂载到桌子上。创建一个简单的近似边界条件,解决一个地区交界处附近的表面和处理(面临21和22)。gydF4y2Ba

structuralBC(模型,gydF4y2Ba“脸”gydF4y2Ba(21、22),gydF4y2Ba“约束”gydF4y2Ba,gydF4y2Ba“固定”gydF4y2Ba);gydF4y2Ba

指定的压力加载在一个齿(面对11)作为一个短矩形压力脉冲。在频域,这种压力脉冲是一个单位荷载均匀分布在所有频率。gydF4y2Ba

structuralBoundaryLoad(模型,gydF4y2Ba“脸”gydF4y2Ba11gydF4y2Ba“压力”gydF4y2Ba1);弗利斯特雨= linspace (0、4000、150);网= generateMesh(模型,gydF4y2Ba“Hmax”gydF4y2Ba,0.005);R =解决(模型2 *π*弗利斯特雨)gydF4y2Ba
R = FrequencyStructuralResults属性:位移:[1 x1 struct]速度:[1 x1 struct]加速度:[1 x1 struct] SolutionFrequencies: [1 x150双]网:[1 x1 FEMesh]gydF4y2Ba

齿尖的振动频率,这是面对12。发现节点的脸上和情节gydF4y2BaygydF4y2Ba位移对频率的分量,利用其中的一个节点。gydF4y2Ba

excitedTineTipNodes = findNodes(网、gydF4y2Ba“地区”gydF4y2Ba,gydF4y2Ba“脸”gydF4y2Ba12);tipDisp = R.Displacement.uy (excitedTineTipNodes (1):);图绘制(弗利斯特雨、abs (tipDisp))包含(gydF4y2Ba“频率”gydF4y2Ba);ylabel (gydF4y2Ba“| Y-Displacement |”gydF4y2Ba);gydF4y2Ba

找到一个三维悬臂梁的挠度下非均匀热负荷。指定的热负荷结构模型上使用瞬态热分析的解决方案相同的几何和网格。gydF4y2Ba

瞬态热模型分析gydF4y2Ba

创建一个瞬态热模型。gydF4y2Ba

thermalmodel = createpde (gydF4y2Ba“热”gydF4y2Ba,gydF4y2Ba瞬态的gydF4y2Ba);gydF4y2Ba

创建和绘制几何。gydF4y2Ba

通用= multicuboid (0.5, 0.1, 0.05);thermalmodel。几何=通用;pdegplot (thermalmodelgydF4y2Ba“FaceLabels”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba,gydF4y2Ba“FaceAlpha”gydF4y2Ba,0.5)gydF4y2Ba

生成一个网格。gydF4y2Ba

网= generateMesh (thermalmodel);gydF4y2Ba

指定材料的热性能。gydF4y2Ba

thermalProperties (thermalmodelgydF4y2Ba“ThermalConductivity”gydF4y2Ba,5 e - 3,gydF4y2Ba…gydF4y2Ba“MassDensity”gydF4y2Ba2.7 * 10 ^ (6),gydF4y2Ba…gydF4y2Ba“SpecificHeat”gydF4y2Ba10);gydF4y2Ba

指定恒定温度应用到梁左右结束。gydF4y2Ba

thermalBC (thermalmodelgydF4y2Ba“脸”gydF4y2Ba3,gydF4y2Ba“温度”gydF4y2Ba,100);thermalBC (thermalmodelgydF4y2Ba“脸”gydF4y2Ba5,gydF4y2Ba“温度”gydF4y2Ba,0);gydF4y2Ba

指定热源在整个几何。gydF4y2Ba

internalHeatSource (thermalmodel 10);gydF4y2Ba

设置初始温度。gydF4y2Ba

thermalIC (thermalmodel 0);gydF4y2Ba

解决模型。gydF4y2Ba

tlist = [0:1e-4:2e-4];thermalresults =解决(thermalmodel tlist)gydF4y2Ba
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]gydF4y2Ba

每个时间步的温度分布。gydF4y2Ba

为gydF4y2Ban = 1:元素个数(thermalresults.SolutionTimes)图pdeplot3D (thermalmodel,gydF4y2Ba“ColorMapData”gydF4y2Bathermalresults.Temperature (:, n))标题([gydF4y2Ba温度时间= 'gydF4y2Banum2str (tlist (n)))) caxis (100 [0])gydF4y2Ba结束gydF4y2Ba

结构分析与热负荷gydF4y2Ba

创建一个静态结构模型。gydF4y2Ba

structuralmodel = createpde (gydF4y2Ba“结构”gydF4y2Ba,gydF4y2Ba“static-solid”gydF4y2Ba);gydF4y2Ba

包括热模型的几何一样。gydF4y2Ba

structuralmodel。几何=通用;gydF4y2Ba

使用相同的网格,用于获得热解决方案。gydF4y2Ba

structuralmodel。网=网;gydF4y2Ba

指定杨氏模量、泊松比、热膨胀系数。gydF4y2Ba

structuralProperties (structuralmodelgydF4y2Ba“YoungsModulus”gydF4y2Ba1 e10,gydF4y2Ba…gydF4y2Ba“PoissonsRatio”gydF4y2Ba,0.3,gydF4y2Ba…gydF4y2Ba”gydF4y2BaCTE的gydF4y2Ba11.7 e-6);gydF4y2Ba

应用一个固定边界条件面临5。gydF4y2Ba

structuralBC (structuralmodelgydF4y2Ba“脸”gydF4y2Ba5,gydF4y2Ba“约束”gydF4y2Ba,gydF4y2Ba“固定”gydF4y2Ba);gydF4y2Ba

应用一个身体负荷使用瞬态热模型的解决方案。默认情况下,gydF4y2BastructuralBodyLoadgydF4y2Ba最后一次使用解决方案的步骤。gydF4y2Ba

structuralBodyLoad (structuralmodelgydF4y2Ba“温度”gydF4y2Ba,thermalresults);gydF4y2Ba

指定参考温度。gydF4y2Ba

structuralmodel。ReferenceTemperature = 10;gydF4y2Ba

解决结构模型。gydF4y2Ba

thermalstressresults =解决(structuralmodel);gydF4y2Ba

情节相对应的梁的变形形状的最后一步瞬态热模型的解决方案。gydF4y2Ba

pdeplot3D (structuralmodelgydF4y2Ba“ColorMapData”gydF4y2Bathermalstressresults.Displacement.Magnitude,gydF4y2Ba…gydF4y2Ba“变形”gydF4y2Bathermalstressresults.Displacement)标题([gydF4y2Ba的热膨胀解决方案时间= 'gydF4y2Banum2str (tlist(结束)))caxis ([0 3 e - 3])gydF4y2Ba

现在指定身体负荷作为热模型解决方案的步骤。金宝搏官方网站对于每一个身体负荷,解决结构模型和策划相应的梁的变形形状。gydF4y2Ba

为gydF4y2Ban = 1:元素个数(thermalresults.SolutionTimes) structuralBodyLoad (structuralmodel,gydF4y2Ba“温度”gydF4y2Bathermalresults,gydF4y2Ba“步伐”gydF4y2Ban);thermalstressresults =解决(structuralmodel);图pdeplot3D (structuralmodel,gydF4y2Ba“ColorMapData”gydF4y2Ba,gydF4y2Ba…gydF4y2Bathermalstressresults.Displacement.Magnitude,gydF4y2Ba…gydF4y2Ba“变形”gydF4y2Ba,gydF4y2Ba…gydF4y2Bathermalstressresults.Displacement)标题(gydF4y2Ba的热解决方案时间=结果gydF4y2Banum2str (tlist (n)))) caxis ([0 3 e - 3])gydF4y2Ba结束gydF4y2Ba

解决了瞬态响应的中心三维光束谐波负载下它的一个角落。gydF4y2Ba

模态分析gydF4y2Ba

创建一个三维的模态分析模型问题。gydF4y2Ba

modelM = createpde (gydF4y2Ba“结构”gydF4y2Ba,gydF4y2Ba“modal-solid”gydF4y2Ba);gydF4y2Ba

创建几何和包括在模型中。绘制几何图形和显示边缘和顶点标签。gydF4y2Ba

通用= multicuboid (0.05, 0.003, 0.003);modelM。几何=通用;pdegplot (modelMgydF4y2Ba“EdgeLabels”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba,gydF4y2Ba“VertexLabels”gydF4y2Ba,gydF4y2Ba“上”gydF4y2Ba);视图(95 [5])gydF4y2Ba

生成一个网格。gydF4y2Ba

msh = generateMesh (modelM);gydF4y2Ba

指定杨氏模量,泊松比,质量密度的材料。gydF4y2Ba

structuralProperties (modelMgydF4y2Ba“YoungsModulus”gydF4y2Ba210 e9,gydF4y2Ba…gydF4y2Ba“PoissonsRatio”gydF4y2Ba,0.3,gydF4y2Ba…gydF4y2Ba“MassDensity”gydF4y2Ba,7800);gydF4y2Ba

指定最小约束梁的一端,以防止刚体模式。例如,指定4边和顶点7是固定的边界。gydF4y2Ba

structuralBC (modelMgydF4y2Ba“边缘”gydF4y2Ba4gydF4y2Ba“约束”gydF4y2Ba,gydF4y2Ba“固定”gydF4y2Ba);structuralBC (modelMgydF4y2Ba“顶点”gydF4y2Ba7gydF4y2Ba“约束”gydF4y2Ba,gydF4y2Ba“固定”gydF4y2Ba);gydF4y2Ba

解决这个问题的频率范围从0到500000。推荐的方法是使用一个值略小于预期的最低频率。因此,使用gydF4y2Ba-0.1gydF4y2Ba而不是gydF4y2Ba0gydF4y2Ba。gydF4y2Ba

Rm =解决(modelM,gydF4y2Ba“FrequencyRange”gydF4y2Ba(-0.1,500000));gydF4y2Ba

瞬态分析gydF4y2Ba

创建一个三维瞬态分析模型问题。gydF4y2Ba

细胞系= createpde (gydF4y2Ba“结构”gydF4y2Ba,gydF4y2Ba“transient-solid”gydF4y2Ba);gydF4y2Ba

使用相同的几何和网格的模态分析。gydF4y2Ba

细胞系。几何=通用;细胞系。要看更多有关憩苑网=;gydF4y2Ba

指定相同的值的杨氏模量,泊松比,质量密度的材料。gydF4y2Ba

structuralProperties(细胞系,gydF4y2Ba“YoungsModulus”gydF4y2Ba210 e9,gydF4y2Ba…gydF4y2Ba“PoissonsRatio”gydF4y2Ba,0.3,gydF4y2Ba…gydF4y2Ba“MassDensity”gydF4y2Ba,7800);gydF4y2Ba

指定相同的最小约束梁的一端,以防止刚体模式。gydF4y2Ba

structuralBC(细胞系,gydF4y2Ba“边缘”gydF4y2Ba4gydF4y2Ba“约束”gydF4y2Ba,gydF4y2Ba“固定”gydF4y2Ba);structuralBC(细胞系,gydF4y2Ba“顶点”gydF4y2Ba7gydF4y2Ba“约束”gydF4y2Ba,gydF4y2Ba“固定”gydF4y2Ba);gydF4y2Ba

角落上施加一个正弦力相反的约束边和顶点。gydF4y2Ba

structuralBoundaryLoad(细胞系,gydF4y2Ba“顶点”gydF4y2Ba5,gydF4y2Ba“力”gydF4y2Ba(0,0,10),gydF4y2Ba“频率”gydF4y2Ba,7600);gydF4y2Ba

指定零初始位移和速度。gydF4y2Ba

structuralIC(细胞系,gydF4y2Ba“速度”gydF4y2Ba(0,0,0),gydF4y2Ba“位移”gydF4y2Ba,(0,0,0));gydF4y2Ba

指定的相对和绝对误差解算器。gydF4y2Ba

modelD.SolverOptions。RelativeTolerance = 1 e-5;modelD.SolverOptions。AbsoluteTolerance = 1 e-9;gydF4y2Ba

解决使用模态模型的结果。gydF4y2Ba

tlist = linspace (0, 0.004,120);Rdm =解决(tlist细胞系,gydF4y2Ba“ModalResults”gydF4y2BaRm)gydF4y2Ba
Rdm = TransientStructuralResults属性:位移:[1 x1 struct]速度:[1 x1 struct]加速度:[1 x1 struct] SolutionTimes: [1 x120双]网:[1 x1 FEMesh]gydF4y2Ba

插入和情节在梁的中心位移。gydF4y2Ba

intrpUdm = interpolateDisplacement (Rdm, 0, 0, 0.0015);情节(Rdm.SolutionTimes intrpUdm.uz)网格gydF4y2Ba在gydF4y2Ba包含(gydF4y2Ba“时间”gydF4y2Ba);ylabel (gydF4y2Ba“梁位移的中心”gydF4y2Ba)gydF4y2Ba

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

静态结构分析模型,指定为一个gydF4y2BaStructuralModelgydF4y2Ba对象。包含几何模型、网格、结构性能的材料,身体负荷,边界载荷和边界条件。gydF4y2Ba

例子:gydF4y2Bastructuralmodel = createpde(“结构”、“static-solid”)gydF4y2Ba

模态分析的结构模型,指定为一个gydF4y2BaStructuralModelgydF4y2Ba对象。包含几何模型、网格、结构性能的材料,身体负荷,边界载荷和边界条件。gydF4y2Ba

例子:gydF4y2Bastructuralmodel = createpde(“结构”、“modal-solid”)gydF4y2Ba

瞬态结构分析模型,指定为一个gydF4y2BaStructuralModelgydF4y2Ba对象。包含几何模型、网格、结构性能的材料,身体负荷,边界载荷和边界条件。gydF4y2Ba

例子:gydF4y2Bastructuralmodel = createpde(“结构”、“transient-solid”)gydF4y2Ba

频率响应分析的结构模型,指定为一个gydF4y2BaStructuralModelgydF4y2Ba对象。包含几何模型、网格、结构性能的材料,身体负荷,边界载荷和边界条件。gydF4y2Ba

例子:gydF4y2Bastructuralmodel = createpde(“结构”、“frequency-solid”)gydF4y2Ba

解决方案的时候,指定为一个真正的单调递增或递减的矢量值。gydF4y2Ba

例子:gydF4y2Ba0:20gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba

解决方案的频率,指定为一个真正的单调递增或递减的矢量值。gydF4y2Ba

例子:gydF4y2Balinspace (0, 4000, 150)gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba

频率范围,指定为一个向量的两个元素。定义gydF4y2BaωgydF4y2Ba比预期的最低频率和略小gydF4y2Ba₂gydF4y2Ba比最高的期望频率。例如,如果最低的期望频率为零,然后用一个小的负面价值gydF4y2BaωgydF4y2Ba。gydF4y2Ba

例子:gydF4y2Ba(-0.1,1000)gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba

模态分析结果,指定为一个gydF4y2BaModalStructuralResultsgydF4y2Ba对象。gydF4y2Ba

例子:gydF4y2Bamodalresults =解决(structuralmodel FrequencyRange, [0,1 e6])gydF4y2Ba

稳态热分析模型,指定为一个gydF4y2BaThermalModelgydF4y2Ba对象。包含几何模型、网格、热性能的材料,内部热源,边界条件和初始条件。gydF4y2Ba

例子:gydF4y2Bathermalmodel = createpde(“热”、“稳态”)gydF4y2Ba

瞬态热分析模型,指定为一个gydF4y2BaThermalModelgydF4y2Ba对象。包含几何模型、网格、热性能的材料,内部热源,边界条件和初始条件。gydF4y2Ba

例子:gydF4y2Bathermalmodel = createpde(“热”、“瞬时”)gydF4y2Ba

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

静态结构分析结果,作为一个返回gydF4y2BaStaticStructuralResultsgydF4y2Ba对象。gydF4y2Ba

模态结构分析结果,作为一个返回gydF4y2BaModalStructuralResultsgydF4y2Ba对象。gydF4y2Ba

瞬态结构分析结果,作为一个返回gydF4y2BaTransientStructuralResultsgydF4y2Ba对象。gydF4y2Ba

频率响应结构分析结果,作为一个返回gydF4y2BaFrequencyStructuralResultsgydF4y2Ba对象。gydF4y2Ba

稳态热分析的结果,作为一个返回gydF4y2BaSteadyStateThermalResultsgydF4y2Ba对象。gydF4y2Ba

瞬态热分析的结果,作为一个返回gydF4y2BaTransientThermalResultsgydF4y2Ba对象。gydF4y2Ba

提示gydF4y2Ba

  • 当你使用模态分析结果来解决一个瞬态结构动力学模型,gydF4y2BamodalresultsgydF4y2Ba在偏微分方程参数必须创建工具箱™R2019a或更新的版本。gydF4y2Ba

  • 频率响应模型的阻尼,结果是复杂的。使用等功能gydF4y2Ba腹肌gydF4y2Ba和gydF4y2Ba角gydF4y2Ba获得实值结果,如大小和相位。gydF4y2Ba

介绍了R2017agydF4y2Ba