Orthotropic Conductivity on Face

1 次查看（过去 30 天）
Paul Safier
Paul Safier 2025-5-2
评论： Paul Safier 2025-5-3
Hi. My domain is as shown in the images attached.
I am trying to place orthotropic conductivity on face 3 and uniform conductivity on faces 1&2. The code/solution fails with the following syntax:
kappa = 1;
orthoK = [0.75*kappa,0 ; 0, 0.85*kappa];
model.MaterialProperties([1 2]) = materialProperties(ThermalConductivity=kappa, ...
MassDensity=1,SpecificHeat=1);
model.MaterialProperties(3) = materialProperties(ThermalConductivity=orthoK, ...
MassDensity=1,SpecificHeat=1);
The code also fails with this test:
orthoK = [kappa,0 ; 0, kappa];
but would run fine with this (not what I need though):
orthoK = [kappa];
It also fails with this syntax that I have seen in the documentation online:
ortho = [0.75*kappa 0.85*kappa];
Can anyone tell me the proper syntax?
Thanks!

请先登录，再回答此问题。

采纳的回答

Torsten
Torsten 2025-5-2
编辑：Torsten 2025-5-2
I found
orthoK = [0.75*kappa; 0.85*kappa];
in the documentation - thus a column instead of a row vector.
  7 个评论
显示 5更早的评论
Paul Safier
Paul Safier 2025-5-3
I'm using:
model = femodel(AnalysisType="thermalSteady",Geometry=gm);
Paul Safier
Paul Safier 2025-5-3
I made a small test problem and was able to get results using this format that I gleaned from the full pdf documentation.
orthoK = [0.75*kappa ; 0 ; 0 ; 0.8*kappa];
However, my problem has a nonlinear boundary condition and is more involved so there is something peculiar about my problem that does not permit a proper solution (at the nonlinear BC) when using an orthotropic conductivity. For the small test problem, I also tried a nonlinear BC and it worked as well, so the issue is specific to my problem. Oh well, more troubleshooting to be had. Thanks for the help, @Torsten.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Heat Transfer 的更多信息

标签

产品


版本

R2023b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by