ヒータプレートの伝熱解析
CAEは、SolidworksSimulation、CATIA、Abaqusなどのソフトがあります。
一方で、フリーウェアCAE（別名、OpenCAE）と呼ばれる、フリーウェアで開発したソフトがあります。
今回は、SalomeMecaというフリーウェアを使用して解析した事例を紹介します。
Monoistに掲載されていた構造物の熱伝導解析の記事をSalomeMecaで非定常熱伝導解析を行いました。
シーズヒータを内蔵したプレートををシーズヒータで１８００秒(３０分)加熱したときのプレートの表面温度を計算しています。記事から読み取れない点は推定でモデル化を行っています。
・記事名：実例で学ぶステップアップ設計者CAE（3）
設計者はどんな視点で設計者CAEを進めていくべきか【ケース2：構造物の熱伝導解析】
・記事の筆者：土橋美博氏
・掲載メディア：MONOist
・解析ソフトウェア：SolidWorks（バージョンは不明）
URL: 設計者はどんな視点で設計者CAEを進めていくべきか【ケース2：構造物の熱伝導解析】：実例で学ぶステップアップ設計者CAE（3）（1/3 ページ） - MONOist (itmedia.co.jp)
＜動作確認環境＞
Salome-Meca 2021(Windows11)
１．１ モデル寸法・物性値
ヒータプレートは、プレートとシーズヒータで構成されたアセンブリ品です。
シーズヒータは、中空のパイプに対してパイプ中央部に電熱線（ニクロム線）を配置して中空部にアルミナを充填した発熱体で、両端部は電熱線がなく、端面はガイシなどの絶縁体で封止をしています。
シーズヒータは、シーズヒータのモデルは記事に詳細がないことから、本ブログでは概略計算の位置づけでシーズヒータはステンレス製の棒状発熱体とみなしました。
シーズヒータとプレートのすきまは、記事と同様にないものとしました。
・プレート：150ｍｍｘ150ｍｍｘ30ｍｍ（ヒータ穴あり）
・シーズヒータ：φ10ｘ120ｍｍ(200W) 4本
・物性値は下表を参照ください。熱伝達率は記事と同一数値を使用しました。
ステンレスの物性値は、http://www.jssa.gr.jp/contents/faq-article/q6の数値から参照しました。
１．２ 解析モデル
メッシュ：四面体１次要素
拘束条件：プレート上面のみ自然対流による熱伝達（周囲は断熱）
発熱源はシーズヒータ全体で発熱させる発熱量：シーズヒータ１本あたり200Wにて30分加熱
SalomeMeca上では単位体積の発熱量0.021W/m3として入力初期温度：室温20℃、シーズヒータ１本あたり200Wにて30分加熱
解析：非定常熱伝導解析
出力：プレート上面の温度
２．モデルの作成
２．１ プレートのモデル作成
plateの形状をFreeCADで作成しエクスポートからstep方式で出力します。
出力は、メニューバーからのファイル(F)をクリックし、エクスポートをクリックします。
ダイアログボックスがでますので、名前を付けて保存します。
SalomeMecaからモデルをインポートします。
インポートは、メニューバーからFileをクリックし、Import→STEPをクリックします。
クリックすると、ダイアログボックスがでるので対象のモデルを選択します。
２．２ 形状作成
(1)Geometry上でヒータの形状を作成します。
(2)Plate上に熱伝達面やヒータとの接触面をExplodeやPartationなどを用いて、Face（面）として作成します。
この操作は、SalomeMecaでモデルを作成した場合と同じです。
２．３ メッシュの作成
Mesh上にてメッシュを作成します。
CreateMeshを用いて、Plateとheaterのメッシュを作成します。Plateは、Create submeshを用いてヒータと接触する面をヒータのメッシュサイズと同一にします。
メッシュのアセンブリコマンド（Create a Compound)を用いて、プレートとヒータを一つのメッシュにします。このときに、プレートとヒータの接合面のメッシュ（接点）を接合する設定にします。
初期値１e-05に対して許容差を緩くしています。トライアンドエラーで行いました。
３．解析結果
ヒータプレート上面の任意１ケ所赤〇部の温度を確認しました。温度結果は以下です。
469秒：179℃
900秒：307℃
1800秒: 558℃
面温度は、記事の解析では表面温度が200℃を基準にシーズヒータをON-OFF制御するのに対して、本解析はヒータが常にONであることから表面温度が上がり続ける結果になりました。
表面の温度分布は、記事と同様にヒータの発熱体が入っている位置の表面温度が上昇していることから意図した解析になっていると考えます。
469℃の時の表面温度は179℃であり、プレートとシーズヒータが均一に469秒間800W加熱された場合141.8℃上昇し温度161.8℃に対して近いことから、単位や入力する発熱量のオーダミスはないものと考えます。この表面温度が均一加熱した温度に対して高いのはシーズヒータに近いためと推察します。
４．まとめ、できたこと、できなかったこと
実務設計を想定したアセンブリ部品での非定常熱伝導解析を行うことができました。
本解析のモデルは、FreeCADで作成したモデルをstep形式へエクスポートしてからSalomeMecaへインポートして作成しました。step形式へエクスポートできるCADソフトがあれば、設計したモデルからSalomeMecaへインポートを行い解析を行うことが可能です。Fusion360やCATIA、SolidWorksなどたいていの3D-CADソフトにはこのエクスポート機能がありますので使用可能です。
今後やっていきたいことは、発熱量の制御で、記事にあった任意の位置の温度で発熱量の関数を作れるかを折をみて調べたいと考えています。
ヒータとプレートの接触面の接点は、接触面のメッシュ長さを同じにして作成し、Toreranceで接点をコンピュータの自動計算でつなぎました。シーズヒータの場合、パイプと充填剤のアルミナと発熱体と合成することから、ほかの方法でやれないか調べたいと考えています。
５．補足
SalomeMecaは、フランスのeDFが開発したオープンソースのソフトウェアです。このソフトは、ソルバー（解析）をおこなうと、CodeAsterとメッシュやモデルを作成するSalomeが一体化したソフトです。