PyFemtet 0.3.8__tar.gz → 0.3.10__tar.gz

{pyfemtet-0.3.8 → pyfemtet-0.3.10}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: pyfemtet
-Version: 0.3.8
+Version: 0.3.10
 Summary: Design parameter optimization using Femtet.
 Home-page: https://github.com/pyfemtet/pyfemtet
 License: BSD-3-Clause

pyfemtet-0.3.10/pyfemtet/FemtetPJTSample/Sldworks_ex01/Sldworks_ex01.py

@@ -0,0 +1,99 @@
+"""外部 CAD ソフト Solidworks との連携
+
+外部 CAD ソフト Solidworks で作成したパラメトリックモデルを
+Femtet にインポートして最適化を行います。
+
+事前準備として、以下の作業を行ってください。
+    - Sldworks_ex01.SLDPRT と Sldworks_ex01.femprj を同じフォルダに配置する
+    - C:\temp フォルダを作成する
+        ※ 最適化の実行中に Solidworks から .x_t ファイルがエクスポートされます。
+"""
+
+import os
+from pyfemtet.opt import FEMOpt, OptunaOptimizer
+from pyfemtet.opt.interface import FemtetWithSolidworksInterface
+from pyfemtet.core import ModelError
+
+
+here, me = os.path.split(__file__)
+os.chdir(here)
+
+
+def disp(Femtet):
+    """解析結果から Z 方向最大変位を取得します。
+    目的関数の第一引数は Femtet インスタンスである必要があります。
+    """
+    _, _, ret = Femtet.Gogh.Galileo.GetMaxDisplacement_py()
+
+    # CAD 連携の場合、境界条件をトポロジーIDに基づいて割り当てているので、
+    # 意図しない箇所に境界条件が割り当てられることがあります。
+    # この問題では意図通りの割り当てができていれば最大変位は必ず負なので
+    # 答えが正の場合は境界条件割り当てに失敗しているとみなして
+    # ModelError を送出するようにします。
+    # 最適化ルーチン中に ModelError, MeshError, SolveError が送出された場合
+    # Optimizer はその試行を失敗とみなし、スキップして次の試行を行います。
+    if ret >= 0:
+        raise ModelError('境界条件の割り当てが間違えています。')
+    return ret
+
+
+def volume(Femtet):
+    """解析結果からモデル体積を取得します。
+    目的関数の第一引数は Femtet インスタンスである必要があります。
+    """
+    _, ret = Femtet.Gogh.CalcVolume_py([0])
+    return ret
+
+
+def C_minus_B(_, opt):
+    """C 寸法と B 寸法の差を計算します。
+
+    拘束関数の第一引数は Femtet インスタンスである必要がありますが、
+    この例では使用していません。
+
+    ほかの例では設計変数にアクセスする際以下のスニペットを用いますが、
+    CAD 連携を行う場合、Femtet に設計変数が設定されていないためです。
+        A = Femtet.GetVariableValue('A')
+
+    pyfemtet.opt では目的関数・拘束関数の第二引数以降に任意の変数を設定できます。
+    この例では、第二引数に FEMOpt のメンバー変数 opt を取り、
+    そのメソッド get_parameter() を用いて設計変数を取得しています。
+
+    """
+    A, B, C = opt.get_parameter('values')
+    return C - B
+
+
+if __name__ == '__main__':
+
+    # Solidworks-Femtet 連携オブジェクトを用意
+    # ここで起動している Femtet と紐づけがされます。
+    fem = FemtetWithSolidworksInterface(
+        sldprt_path='Sldworks_ex01.SLDPRT',
+        femprj_path='Sldworks_ex01.femprj',
+    )
+
+    opt = OptunaOptimizer(
+        sampler_kwargs=dict(
+            n_startup_trials=5,
+        )
+    )
+
+    # 最適化処理を行うオブジェクトに連携オブジェクトを紐づけ
+    femopt = FEMOpt(fem=fem, opt=opt)
+
+    # 設計変数の設定
+    femopt.add_parameter('A', 10, lower_bound=1, upper_bound=59)
+    femopt.add_parameter('B', 10, lower_bound=1, upper_bound=40)
+    femopt.add_parameter('C', 20, lower_bound=5, upper_bound=59)
+
+    # 拘束関数の設定
+    femopt.add_constraint(C_minus_B, 'C>B', lower_bound=1, args=femopt.opt)
+
+    # 目的関数の設定
+    femopt.add_objective(disp, name='変位', direction=0)
+    femopt.add_objective(volume, name='体積', direction='minimize')
+
+    # 最適化の実行
+    femopt.set_random_seed(42)
+    femopt.main(n_trials=20)

pyfemtet-0.3.10/pyfemtet/__init__.py

@@ -0,0 +1 @@
+__version__ = "0.3.10"

{pyfemtet-0.3.8 → pyfemtet-0.3.10}/pyfemtet/opt/base.py

@@ -2,6 +2,7 @@ from abc import ABC, abstractmethod
 from typing import List, Iterable
 
 import os
+import sys
 import datetime
 import inspect
 import ast
@@ -842,7 +843,7 @@ class OptunaOptimizer(AbstractOptimizer):
         try:
             _, _y, c = self.f(x)
         except (ModelError, MeshError, SolveError) as e:
-            logger.info(str(type(e)))
+            logger.info(e)
             logger.info('以下の変数で FEM 解析に失敗しました。')
             print(self.get_parameter('dict'))
 
@@ -1330,7 +1331,7 @@ class FEMOpt:
         # monitor 用 worker を起動
         logger.info('Launching monitor server. This may take a few seconds.')
         self.monitor_process_worker_name = datetime.datetime.now().strftime("Monitor-%Y%m%d-%H%M%S")
-        cmd = f'dask worker {self.client.scheduler.address} --name {self.monitor_process_worker_name} --no-nanny'
+        cmd = f'{sys.executable} -m dask worker {self.client.scheduler.address} --name {self.monitor_process_worker_name} --no-nanny'
         current_n_workers = len(self.client.nthreads().keys())
         Popen(cmd, shell=True)  # , stdout=PIPE) --> cause stream error
 

{pyfemtet-0.3.8 → pyfemtet-0.3.10}/pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "pyfemtet"
-version = "0.3.8"  # ignored by versioning plugin
+version = "0.3.10"  # ignored by versioning plugin
 description = "Design parameter optimization using Femtet."
 authors = ["kazuma.naito <kazuma.naito@murata.com>"]
 readme = "README.md"

pyfemtet-0.3.8/pyfemtet/__init__.py

@@ -1 +0,0 @@
-__version__ = "0.3.8"