pyfemtet 0.9.5__py3-none-any.whl → 1.0.0b0__py3-none-any.whl

pyfemtet/opt/samples/femprj_sample_jp/cad_ex01_NX_jp.py

@@ -1,129 +0,0 @@
-"""外部 CAD (NX) 連携
-
-Femtet の応力解析ソルバ、および
-Siemens 社製 CAD ソフト NX を用いて
-軽量かつ高剛性な H 型鋼の設計を行います。
-
-事前準備として、下記の手順を実行してください。
-- NX のインストール
-- C:\temp フォルダを作成する
-    - Note: NX が .x_t ファイルをこのフォルダに保存します。
-- 以下のファイルを同じフォルダに配置
-    - cad_ex01_NX_jp.py (このファイル)
-    - cad_ex01_NX.prt
-    - cad_ex01_NX_jp.femprj
-"""
-
-import os
-
-from win32com.client import constants
-
-from pyfemtet.opt import FEMOpt
-from pyfemtet.opt.interface import FemtetWithNXInterface
-from pyfemtet.core import ModelError
-
-
-here, me = os.path.split(__file__)
-os.chdir(here)
-
-
-def von_mises(Femtet):
-    """モデルの最大フォン・ミーゼス応力を取得します。
-
-    Note:
-        目的関数または制約関数は、
-        第一引数としてFemtetを受け取り、
-        戻り値としてfloat型を返す必要があります。
-
-    Warning:
-        CAD 連携機能では、意図しない位置に境界条件が設定される可能性があります。
-
-        この例では、境界条件が意図したとおりに割り当てられている場合、
-        最大変位は常に負になります。最大変位が正の場合、境界条件の割り当てが
-        失敗したとみなし、ModelError を送出します。
-
-        最適化中に ModelError、MeshError、または SolveError が発生した場合、
-        最適化プロセスは試行を失敗とみなし、次のトライアルにスキップします。
-    """
-
-    # 簡易的な境界条件の正しさチェック
-    dx, dy, dz = Femtet.Gogh.Galileo.GetMaxDisplacement_py()
-    if dz >= 0:
-        raise ModelError('境界条件の設定が間違っています。')
-
-    # ミーゼス応力計算
-    Gogh = Femtet.Gogh
-    Gogh.Galileo.Potential = constants.GALILEO_VON_MISES_C
-    succeed, (x, y, z), mises = Gogh.Galileo.GetMAXPotentialPoint_py(constants.CMPX_REAL_C)
-
-    return mises
-
-
-def mass(Femtet):
-    """モデルの質量を取得します。"""
-    return Femtet.Gogh.Galileo.GetMass('H_beam')
-
-
-def C_minus_B(Femtet, opt):
-    """C 寸法と B 寸法の差を計算します。
-
-    別の例では、次のスニペットを使用して設計変数にアクセスします。
-
-        A = Femtet.GetVariableValue('A')
-    
-    ただし、CAD 連携機能を使用する場合、設計変数が .femprj ファイルに
-    設定されていないため、この方法は機能しません。
-
-    CAD 連携機能を使用する場合、以下の方法で設計変数にアクセスすることができます。
-
-        # add_parameter() で追加したパラメータの変数名をキーとする辞書を得る方法
-        params: dict = opt.get_parameter()
-        A = params['A']
-
-    又は
-
-        # add_parameter() で追加した順のパラメータの値の配列を得る方法
-        values: np.ndarray = opt.get_parameter('values')
-        A, B, C = values
-
-    目的関数と拘束関数は、最初の引数の後に任意の変数を取ることができます。
-    FEMOpt のメンバ変数 opt には get_parameter() というメソッドがあります。
-    このメソッドによって add_parameter() で追加された設計変数を取得できます。
-    opt を第 2 引数として取ることにより、目的関数または拘束関数内で
-    get_parameter() を実行して設計変数を取得できます。
-    """
-    A, B, C = opt.get_parameter('values')
-    return C - B
-
-
-if __name__ == '__main__':
-
-    # NX-Femtet 連携オブジェクトの初期化
-    # この処理により、Python プロセスは Femtet に接続を試みます。
-    fem = FemtetWithNXInterface(
-        prt_path='cad_ex01_NX.prt',
-        open_result_with_gui=False,
-        export_curves=False,
-        export_surfaces=False,
-        export_solids=True,
-        export_flattened_assembly=False,
-    )
-
-    # FEMOpt オブジェクトの初期化 (最適化問題とFemtetとの接続を行います)
-    femopt = FEMOpt(fem=fem)
-
-    # 設計変数を最適化問題に追加 (femprj ファイルに登録されている変数を指定してください)
-    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(von_mises, name='von Mises (Pa)')
-    femopt.add_objective(mass, name='mass (kg)')
-
-    # 最適化を実行
-    femopt.set_random_seed(42)
-    femopt.optimize(n_trials=20)

pyfemtet/opt/samples/femprj_sample_jp/cad_ex01_SW_jp.py

@@ -1,125 +0,0 @@
-"""外部 CAD (SOLIDWORKS) 連携
-
-Femtet の応力解析ソルバ、および
-Dassault Systemes 社製 CAD ソフト SOLIDWORKS を用いて
-軽量かつ高剛性な H 型鋼の設計を行います。
-
-事前準備として、下記の手順を実行してください。
-- SOLIDWORKS のインストール
-- C:\temp フォルダを作成する
-    - Note: SOLIDWORKS が .x_t ファイルをこのフォルダに保存します。
-- 以下のファイルを同じフォルダに配置
-    - cad_ex01_SW_jp.py (このファイル)
-    - cad_ex01_SW.SLDPRT
-    - cad_ex01_SW_jp.femprj
-"""
-
-import os
-
-from win32com.client import constants
-
-from pyfemtet.opt import FEMOpt
-from pyfemtet.opt.interface import FemtetWithSolidworksInterface
-from pyfemtet.core import ModelError
-
-
-here, me = os.path.split(__file__)
-os.chdir(here)
-
-
-def von_mises(Femtet):
-    """モデルの最大フォン・ミーゼス応力を取得します。
-
-    Note:
-        目的関数または制約関数は、
-        第一引数としてFemtetを受け取り、
-        戻り値としてfloat型を返す必要があります。
-
-    Warning:
-        CAD 連携機能では、意図しない位置に境界条件が設定される可能性があります。
-
-        この例では、境界条件が意図したとおりに割り当てられている場合、
-        最大変位は常に負になります。最大変位が正の場合、境界条件の割り当てが
-        失敗したとみなし、ModelError を送出します。
-
-        最適化中に ModelError、MeshError、または SolveError が発生した場合、
-        最適化プロセスは試行を失敗とみなし、次のトライアルにスキップします。
-    """
-
-    # 簡易的な境界条件の正しさチェック
-    dx, dy, dz = Femtet.Gogh.Galileo.GetMaxDisplacement_py()
-    if dz >= 0:
-        raise ModelError('境界条件の設定が間違っています。')
-
-    # ミーゼス応力計算
-    Gogh = Femtet.Gogh
-    Gogh.Galileo.Potential = constants.GALILEO_VON_MISES_C
-    succeed, (x, y, z), mises = Gogh.Galileo.GetMAXPotentialPoint_py(constants.CMPX_REAL_C)
-
-    return mises
-
-
-def mass(Femtet):
-    """モデルの質量を取得します。"""
-    return Femtet.Gogh.Galileo.GetMass('H_beam')
-
-
-def C_minus_B(Femtet, opt):
-    """C 寸法と B 寸法の差を計算します。
-
-    別の例では、次のスニペットを使用して設計変数にアクセスします。
-
-        A = Femtet.GetVariableValue('A')
-    
-    ただし、CAD 連携機能を使用する場合、設計変数が .femprj ファイルに
-    設定されていないため、この方法は機能しません。
-
-    CAD 連携機能を使用する場合、以下の方法で設計変数にアクセスすることができます。
-
-        # add_parameter() で追加したパラメータの変数名をキーとする辞書を得る方法
-        params: dict = opt.get_parameter()
-        A = params['A']
-
-    又は
-
-        # add_parameter() で追加した順のパラメータの値の配列を得る方法
-        values: np.ndarray = opt.get_parameter('values')
-        A, B, C = values
-
-    目的関数と拘束関数は、最初の引数の後に任意の変数を取ることができます。
-    FEMOpt のメンバ変数 opt には get_parameter() というメソッドがあります。
-    このメソッドによって add_parameter() で追加された設計変数を取得できます。
-    opt を第 2 引数として取ることにより、目的関数または拘束関数内で
-    get_parameter() を実行して設計変数を取得できます。
-    """
-    A, B, C = opt.get_parameter('values')
-    return C - B
-
-
-if __name__ == '__main__':
-
-    # NX-Femtet 連携オブジェクトの初期化
-    # この処理により、Python プロセスは Femtet に接続を試みます。
-    fem = FemtetWithSolidworksInterface(
-        sldprt_path='cad_ex01_SW.SLDPRT',
-        open_result_with_gui=False,
-    )
-
-    # FEMOpt オブジェクトの初期化 (最適化問題とFemtetとの接続を行います)
-    femopt = FEMOpt(fem=fem)
-
-    # 設計変数を最適化問題に追加 (femprj ファイルに登録されている変数を指定してください)
-    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(von_mises, name='von Mises (Pa)')
-    femopt.add_objective(mass, name='mass (kg)')
-
-    # 最適化を実行
-    femopt.set_random_seed(42)
-    femopt.optimize(n_trials=20)

pyfemtet/opt/samples/femprj_sample_jp/constrained_pipe_jp.py

@@ -1,93 +0,0 @@
-"""拘束付き最適化を実装するサンプル。
-
-このセクションでは、拘束の種類と、拘束を必要とするモデルで
-最適化を実行する手順について説明します。
-
-"""
-
-from optuna_integration import BoTorchSampler
-from pyfemtet.opt import FEMOpt, OptunaOptimizer
-
-
-def mises_stress(Femtet):
-    """フォンミーゼス応力を目的関数として計算します。
-
-    この関数は、最適化の実行中に FEMOpt オブジェクトによって
-    自動的に呼び出されます。
-
-    引数:
-        Femtet: PyFemtet を使用して目的関数または拘束関数を
-            定義する場合、最初の引数は Femtet インスタンスを
-            取る必要があります。
-
-    戻り値:
-        float: 目的または拘束関数は単一の float を返すよう定義してください。
-    """
-    return Femtet.Gogh.Galileo.GetMaxStress_py()[2]
-
-
-def radius_diff(Femtet, opt):
-    """パイプの外側の半径と内側の半径の差を計算します。
-
-    この拘束は、最適化の実行中にパイプの内側の半径が
-    外側の半径を超えないようにするために呼び出されます。
-
-    注意:
-        OptunaOptimizer の BoTorchSampler を使用していて、
-        strict な拘束を使用する場合、パラメータを提案するため
-        に繰り返し計算が必要になるため、Femtet へのアクセスが
-        非常に遅くなる可能性があることに注意してください。
-        この関数の例のように、Femtet にアクセスするのではなく、
-        Optimizer オブジェクトを介してパラメータを取得して計算
-        を実行することをお勧めします。
-
-        非推奨::
-
-            p = Femtet.GetVariableValue('p')
-
-        代わりに::
-
-            params = opt.get_parameter()
-            p = params['p']
-
-    引数:
-        Femtet: PyFemtet を使用して目的関数または拘束関数を
-            定義する場合、最初の引数は Femtet インスタンスを
-            取る必要があります。
-        opt: このオブジェクトを使用すると、Femtet を経由せず
-            に外側の半径と内側の半径の値を取得できます。
-    """
-    params = opt.get_parameter()
-    internal_r = params['internal_r']
-    external_r = params['external_r']
-    return external_r - internal_r
-
-
-if __name__ == '__main__':
-    # 最適化手法のセットアップ
-    opt = OptunaOptimizer(
-        sampler_class=BoTorchSampler,
-        sampler_kwargs=dict(
-            n_startup_trials=3,  # 最初の 3 回はランダムサンプリングを行います。
-        )
-    )
-    femopt = FEMOpt(opt=opt)
-
-    # 変数の追加
-    femopt.add_parameter("external_r", 10, lower_bound=0.1, upper_bound=10)
-    femopt.add_parameter("internal_r", 5, lower_bound=0.1, upper_bound=10)
-
-    # 最適化の実行中に外側の半径を超えないように strict 拘束を追加します。
-    femopt.add_constraint(
-        radius_diff,  # 拘束関数 (ここでは 外半径 - 内半径).
-        name='管厚さ',  # 拘束関数にはプログラム上の名前とは別に自由な名前を付与できます.
-        lower_bound=1,  # 拘束関数の下限 (ここでは管の厚みを最低 1 とする).
-        args=(femopt.opt,)  # 拘束関数に渡される、Femtet 以外の追加の引数.
-    )
-
-    # 目的関数の追加
-    femopt.add_objective(mises_stress, name='ミーゼス応力')
-
-    # 最適化の実行
-    femopt.set_random_seed(42)
-    femopt.optimize(n_trials=10)

pyfemtet/opt/samples/femprj_sample_jp/gal_ex58_parametric_jp.py

@@ -1,70 +0,0 @@
-"""単目的最適化: スプリングバックを加味した曲げ
-
-Femtet の応力解析ソルバを利用して、スプリングバックを考慮した
-目標の材料曲げ角度を達成するために必要な曲げ角度を決定します。
-※ 弾塑性解析は特別オプション機能です。
-
-対応プロジェクト：gal_ex58_parametric_jp.femprj
-"""
-import numpy as np
-from win32com.client import constants
-from optuna.integration.botorch import BoTorchSampler
-
-from pyfemtet.opt import FEMOpt, OptunaOptimizer
-
-
-def bending(Femtet):
-    """材料の曲げ角度を取得します。
-
-    Note:
-        目的関数または制約関数は、
-        第一引数としてFemtetを受け取り、
-        戻り値としてfloat型を返す必要があります。
-
-    Params:
-        Femtet: Femtet をマクロで操作するためのインスタンスです。詳細な情報については、「Femtet マクロヘルプ」をご覧ください。
-
-    Returns:
-        float: 曲げ。
-    """
-    Gogh = Femtet.Gogh
-
-    # モードを除荷後に設定
-    Gogh.Galileo.Mode = Gogh.Galileo.nMode - 1
-
-    # 計測対象点の変位を取得
-    Gogh.Galileo.Vector = constants.GALILEO_DISPLACEMENT_C
-    succeed, (x, y, z) = Gogh.Galileo.GetVectorAtPoint_py(200, 0, 0)
-
-    # 曲げ起点 (100, 0) と変形後の点を結ぶ線分が X 軸となす角度を計算
-    bending_point = np.array((100, 0))
-    bended_point = np.array((200 + 1000 * x.Real, 1000 * z.Real))
-    dx, dz = bended_point - bending_point
-    degree = np.arctan2(-dz, dx)
-
-    return degree * 360 / (2*np.pi)  # 単位: 度
-
-
-if __name__ == '__main__':
-
-    # 数値最適化問題の初期化 (最適化手法を決定します)
-    opt = OptunaOptimizer(
-        sampler_class=BoTorchSampler,
-        sampler_kwargs=dict(
-            n_startup_trials=3,
-        )
-    )
-
-    # FEMOpt オブジェクトの初期化 (最適化問題とFemtetとの接続を行います)
-    femopt = FEMOpt(opt=opt)
-
-    # 設計変数を最適化問題に追加 (femprj ファイルに登録されている変数を指定してください)
-    femopt.add_parameter("rot", 90, lower_bound=80, upper_bound=100)
-
-    # 目的関数を最適化問題に追加
-    # 曲げ角度の目標は 90° です。
-    femopt.add_objective(bending, name='曲げ角度（度）', direction=90)
-
-    # 最適化を実行
-    femopt.set_random_seed(42)
-    femopt.optimize(n_trials=10)

pyfemtet/opt/samples/femprj_sample_jp/gau_ex08_parametric_jp.py

@@ -1,57 +0,0 @@
-"""単目的最適化: 有限長ヘリカルコイルの自己インダクタンス
-
-Femtet の磁場解析ソルバを利用して、
-有限長ヘリカルコイルの自己インダクタンスを
-目標の値にする設計を行います。
-
-対応プロジェクト: gau_ex08_parametric_jp.femprj
-"""
-from optuna.integration.botorch import BoTorchSampler
-from pyfemtet.opt import FEMOpt, OptunaOptimizer
-
-
-def inductance(Femtet):
-    """自己インダクタンスを取得します。
-
-    Note:
-        目的関数または制約関数は、
-        第一引数としてFemtetを受け取り、
-        戻り値としてfloat型を返す必要があります。
-
-    Params:
-        Femtet: Femtet をマクロで操作するためのインスタンスです。詳細な情報については、「Femtet マクロヘルプ」をご覧ください。
-
-    Returns:
-        float: 自己インダクタンスです。
-    """
-    Gogh = Femtet.Gogh
-
-    coil_name = Gogh.Gauss.GetCoilList()[0]
-    return Gogh.Gauss.GetL(coil_name, coil_name)  # 単位: F
-
-
-if __name__ == '__main__':
-
-    # 数値最適化問題の初期化 (最適化手法を決定します)
-    opt = OptunaOptimizer(
-        sampler_class=BoTorchSampler,
-        sampler_kwargs=dict(
-            n_startup_trials=5,
-        )
-    )
-
-    # FEMOpt オブジェクトの初期化 (最適化問題とFemtetとの接続を行います)。
-    femopt = FEMOpt(opt=opt)
-
-    # 設計変数を最適化問題に追加 (femprj ファイルに登録されている変数を指定してください)。
-    femopt.add_parameter("helical_pitch", 6, lower_bound=4.2, upper_bound=8)
-    femopt.add_parameter("coil_radius", 10, lower_bound=1, upper_bound=10)
-    femopt.add_parameter("n_turns", 5, lower_bound=1, upper_bound=5)
-
-    # 目的関数を最適化問題に追加
-    # 目標の自己インダクタンスは 0.1 μF です。
-    femopt.add_objective(inductance, name='自己インダクタンス (F)', direction=1e-7)
-
-    # 最適化を実行
-    femopt.set_random_seed(42)
-    femopt.optimize(n_trials=20)

pyfemtet/opt/samples/femprj_sample_jp/gau_ex12_parametric_jp.py

@@ -1,52 +0,0 @@
-"""パラメトリック解析出力設定を目的関数とする最適化
-
-Femtet のパラメトリック解析の結果出力設定機能で出力される値を
-最適化の目的関数として使用する方法をデモします。
-この機能により、目的関数をコーディングすることなく
-最適化を実施できます。
-
-
-注意:
-
-    この機能を使う際は、以下のことに注意してください。
-
-    - パラメトリック解析のスイープテーブルが削除されます。
-    - 複素数やベクトルを出力する出力設定は、第一の値のみが
-      目的関数として使用されます。（複素数の場合は実数、
-      ベクトル値の場合は X 成分など）
-
-
-対応するプロジェクト: gau_ex12_parametric.femprj
-"""
-
-from pyfemtet.opt import FEMOpt, FemtetInterface
-
-
-if __name__ == '__main__':
-
-    # Femtet の設定を参照するため、Femtet と接続を
-    # 行うためのオブジェクトを初期化します。
-    fem = FemtetInterface()
-
-    # パラメトリック解析の結果出力設定を目的関数にします。
-    # number は Femtet パラメトリック解析ダイアログの
-    # 結果出力設定タブのテーブルの番号で、direction は
-    # その目的関数の目標です(FEMOpt.add_objective と同様)。
-
-    # 相互インダクタンス
-    fem.use_parametric_output_as_objective(number=1, direction=1.5e-7)
-
-    # コイル中央の磁界の強さ
-    fem.use_parametric_output_as_objective(number=2, direction='minimize')
-
-    # 最適化用オブジェクトを初期化します。
-    # さきほど初期化した fem を渡します。
-    femopt = FEMOpt(fem=fem)
-
-    # パラメータを設定します。
-    femopt.add_parameter('in_radius', 10, 5, 10)
-    femopt.add_parameter('out_radius', 20, 20, 25)
-
-    # 最適化を実行します。
-    femopt.set_random_seed(42)  # 乱数シードの固定
-    femopt.optimize(n_trials=20)

pyfemtet/opt/samples/femprj_sample_jp/her_ex40_parametric_jp.py

@@ -1,138 +0,0 @@
-"""単目的最適化: 円形パッチアンテナの共振周波数
-
-Femtet の電磁波解析ソルバを利用して、円形パッチアンテナの
-電磁波調和解析を行い、共振特性を目標の値にする設計を行います。
-
-対応プロジェクト: her_ex40_parametric_jp.femprj
-"""
-from time import sleep
-
-import numpy as np
-from scipy.signal import find_peaks
-from tqdm import tqdm
-from optuna.integration.botorch import BoTorchSampler
-
-from pyfemtet.core import SolveError
-from pyfemtet.opt import OptunaOptimizer, FEMOpt
-
-
-class SParameterCalculator:
-    """Sパラメータ計算用クラス"""
-    
-    def __init__(self):
-        self.freq = []
-        self.S = []
-        self.interpolated_function = None
-        self.resonance_frequency = None
-        self.minimum_S = None
-
-    def _get_freq_and_S_parameter(self, Femtet):
-        """周波数とSパラメータの関係を取得します。"""
-
-        Gogh = Femtet.Gogh
-        
-        freq_list = []
-        dB_S_list = []
-        for mode in tqdm(range(Gogh.Hertz.nMode), '周波数と S(1, 1) の関係を取得'):
-            # 周波数モード設定
-            Gogh.Hertz.Mode = mode
-            sleep(0.01)
-
-            # 周波数を取得
-            freq = Gogh.Hertz.GetFreq().Real
-
-            # S(1, 1) を取得
-            comp_S = Gogh.Hertz.GetSMatrix(0, 0)
-            norm = np.linalg.norm((comp_S.Real, comp_S.Imag))
-            dB_S = 20 * np.log10(norm)
-
-            # 結果を保存
-            freq_list.append(freq)
-            dB_S_list.append(dB_S)
-
-        self.freq = freq_list
-        self.S = dB_S_list
-
-    def _calc_resonance_frequency(self):
-        """Sパラメータの第一ピークを与える周波数を取得します。"""
-        peaks, _ = find_peaks(-np.array(self.S), height=None, threshold=None, distance=None, prominence=0.5, width=None, wlen=None, rel_height=0.5, plateau_size=None)
-        if len(peaks) == 0:
-            raise SolveError('S(1,1) のピークを取得できませんでした。')
-        self.resonance_frequency = self.freq[peaks[0]]
-        self.minimum_S = self.S[peaks[0]]
-
-    def get_resonance_frequency(self, Femtet):
-        """パッチアンテナの共振周波数を計算します。
-
-        Note:
-            目的関数または制約関数は、
-            第一引数としてFemtetを受け取り、
-            戻り値としてfloat型を返す必要があります。
-
-        Params:
-            Femtet: Femtet をマクロで操作するためのインスタンスです。詳細な情報については、「Femtet マクロヘルプ」をご覧ください。
-        
-        Returns:
-            float: パッチアンテナの共振周波数。
-        """
-        self._get_freq_and_S_parameter(Femtet)
-        self._calc_resonance_frequency()
-        return self.resonance_frequency  # 単位: Hz
-        
-
-def antenna_is_smaller_than_substrate(Femtet, opt):
-    """アンテナの大きさと基板の大きさの関係を計算します。
-
-    この関数は、変数の更新によってモデル形状が破綻しないように
-    変数の組み合わせを拘束するために使われます。
-
-    Params:
-        Femtet: Femtet をマクロで操作するためのインスタンスです。詳細な情報については、「Femtet マクロヘルプ」をご覧ください。
-    
-    Returns:
-        float: 基板エッジとアンテナエッジの間隙。1 mm 以上が必要です。
-    """
-    params = opt.get_parameter()
-    r = params['antenna_radius']
-    w = params['substrate_w']
-    return w / 2 - r  # 単位: mm
-
-
-def port_is_inside_antenna(Femtet, opt):
-    """給電ポートの位置とアンテナの大きさの関係を計算します。"""
-    params = opt.get_parameter()
-    r = params['antenna_radius']
-    x = params['port_x']
-    return r - x  # 単位: mm。1 mm 以上が必要です。
-
-
-if __name__ == '__main__':
-    # 周波数特性を計算するためのオブジェクトを初期化
-    s = SParameterCalculator()
-
-    # 数値最適化問題の初期化 (最適化手法を決定します)
-    opt = OptunaOptimizer(
-        sampler_class=BoTorchSampler,
-        sampler_kwargs=dict(
-            n_startup_trials=10,
-        )
-    )
-    
-    # FEMOpt オブジェクトの初期化 (最適化問題とFemtetとの接続を行います)
-    femopt = FEMOpt(opt=opt)
-    
-    # 設計変数を最適化問題に追加 (femprj ファイルに登録されている変数を指定してください)
-    femopt.add_parameter('antenna_radius', 10, 5, 20)
-    femopt.add_parameter('substrate_w', 50, 40, 60)
-    femopt.add_parameter('port_x', 5, 1, 20)
-
-    # 拘束関数を最適化問題に追加
-    femopt.add_constraint(antenna_is_smaller_than_substrate, 'アンテナと基板エッジの間隙', lower_bound=1, args=(opt,))
-    femopt.add_constraint(port_is_inside_antenna, 'アンテナエッジと給電ポートの間隙', lower_bound=1, args=(opt,))
-
-    # 目的関数を最適化問題に追加
-    # 共振周波数の目標は 3.0 GHz です。
-    femopt.add_objective(s.get_resonance_frequency, '第一共振周波数(Hz)', direction=3.0 * 1e9)
-
-    femopt.set_random_seed(42)
-    femopt.optimize(n_trials=15)