pyfemtet 0.4.1__py3-none-any.whl → 0.4.3__py3-none-any.whl

pyfemtet/FemtetPJTSample/her_ex40_parametric.py

@@ -1,3 +1,9 @@
+"""Single-objective optimization: Resonant frequency of a circular patch antenna
+
+Using Femtet’s electromagnetic wave analysis solver,
+we explain an example of setting the resonant frequency
+of a circular patch antenna to a specific value.
+"""
 from time import sleep
 
 import numpy as np
@@ -8,9 +14,8 @@ from optuna.integration.botorch import BoTorchSampler
 from pyfemtet.opt import OptunaOptimizer, FEMOpt
 
 
-
 class SParameterCalculator:
-    """S パラメータ・共振周波数計算用クラスです。"""
+    """This class is for calculating S-parameters and resonance frequencies."""
     
     def __init__(self):
         self.freq = []
@@ -20,51 +25,55 @@ class SParameterCalculator:
         self.minimum_S = None
 
     def get_result_from_Femtet(self, Femtet):
-        """Femtet の解析結果から周波数と S パラメータの関係を取得します。"""
+        """Obtain the relationship between frequency and S-parameter from the Femtet analysis results."""
 
-        # 前準備
+        # Preparation
         Femtet.OpenCurrentResult(True)
         Gogh = Femtet.Gogh
         
-        # 各モードに対して周波数と S(1,1) を取得する
+        # Obtain the frequency and S(1,1) for each mode
         mode = 0
         freq_list = []
         dB_S_list = []
-        for mode in tqdm(range(Gogh.Hertz.nMode), '結果取得中'):
-            # Femtet 結果画面のモード設定
+        for mode in tqdm(range(Gogh.Hertz.nMode), 'Obtaining frequency and S-parameter'):
+            # Femtet result screen mode settings
             Gogh.Hertz.Mode = mode
             sleep(0.01)
-            # 周波数の取得
+            # Get frequency
             freq = Gogh.Hertz.GetFreq().Real
-            # S パラメータの取得
+            # Get S-parameters
             comp_S = Gogh.Hertz.GetSMatrix(0, 0)
             norm = np.linalg.norm((comp_S.Real, comp_S.Imag))
             dB_S = 20 * np.log10(norm)
-            # 結果の取得
+            # Get results
             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 に対して最初のピークを与える freq を計算します。"""
+        """Compute the frequency that gives the first peak for S-parameter."""
         x = -np.array(self.S)
         peaks, _ = find_peaks(x, height=None, threshold=None, distance=None, prominence=0.5, width=None, wlen=None, rel_height=0.5, plateau_size=None)
         from pyfemtet.core import SolveError
         if len(peaks) == 0:
-            raise SolveError('ピークが検出されませんでした。')
+            raise SolveError('No peaks detected.')
         self.resonance_frequency = self.freq[peaks[0]]
         self.minimum_S = self.S[peaks[0]]
 
     def get_resonance_frequency(self, Femtet):
-        """共振周波数を計算します。
-
-        Femtet : マクロを使用するためのインスタンスです。詳しくは "Femtet マクロヘルプ / CFemtet クラス" をご覧ください。
-            目的関数は第一引数に Femtet インスタンスを取る必要があります。
+        """Calculate the resonant frequency.
 
-        f : 計算された共振周波数です。
-            目的関数は単一の float を返す必要があります。
+        Note:
+            The objective or constraint function
+            must take a Femtet as its first argument
+            and must return a single float.
 
+        Params:
+            Femtet: An instance for using Femtet macros. For more information, see "Femtet Macro Help / CFemtet Class".
+        
+        Returns:
+            float: A resonance frequency.
         """
         self.get_result_from_Femtet(Femtet)
         self.calc_resonance_frequency()
@@ -73,14 +82,16 @@ class SParameterCalculator:
         
 
 def antenna_is_smaller_than_substrate(Femtet):
-    """アンテナのサイズと基板のサイズの関係を計算します。
-
-    Femtet : マクロを使用するためのインスタンスです。詳しくは "Femtet マクロヘルプ / CFemtet クラス" をご覧ください。
-        拘束関数は第一引数に Femtet インスタンスを取る必要があります。
+    """Calculate the relationship between antenna size and board size.
 
-    Sx/2 - ant_r : 基板サイズとアンテナサイズの差です。
-        拘束関数は単一の float を返す必要があります。
+    This function is used to constrain the model
+    from breaking down while changing parameters.
 
+    Params:
+        Femtet: An instance for using Femtet macros.
+    
+    Returns:
+        float: Difference between the board size and antenna size. Must be equal to or grater than 1 mm.
     """
     ant_r = Femtet.GetVariableValue('ant_r')
     Sx = Femtet.GetVariableValue('sx')
@@ -88,14 +99,16 @@ def antenna_is_smaller_than_substrate(Femtet):
 
 
 def port_is_inside_antenna(Femtet):
-    """給電ポートの位置とアンテナのサイズの関係を計算します。
-
-    Femtet : マクロを使用するためのインスタンスです。詳しくは "Femtet マクロヘルプ / CFemtet クラス" をご覧ください。
-        拘束関数は第一引数に Femtet インスタンスを取る必要があります。
+    """Calculate the relationship between the feed port location and antenna size.
 
-    ant_r - xf : アンテナ辺縁と給電ポートの位置の差です。
-        拘束関数は単一の float を返す必要があります。
+    This function is used to constrain the model
+    from breaking down while changing parameters.
 
+    Params:
+        Femtet: An instance for using Femtet macros.
+    
+    Returns:
+        float: Difference between the antenna edge and the position of the feed port. Must be equal to or grater than 1 mm.
     """
     ant_r = Femtet.GetVariableValue('ant_r')
     xf = Femtet.GetVariableValue('xf')
@@ -103,10 +116,10 @@ def port_is_inside_antenna(Femtet):
 
 
 if __name__ == '__main__':
-    # S パラメータ計算用クラス
+    # Define the object for calculating S-parameters and resonance frequencies.
     s = SParameterCalculator()
 
-    # 最適化手法を定義するオブジェクトを用意
+    # Define mathematical optimization object.
     opt = OptunaOptimizer(
         sampler_class=BoTorchSampler,
         sampler_kwargs=dict(
@@ -114,22 +127,21 @@ if __name__ == '__main__':
         )
     )
     
-    # 最適化処理を行うオブジェクトを用意
-    femopt = FEMOpt(
-        opt=opt,
-    )  # ここで起動している Femtet が紐づけされます
+    # Define FEMOpt object (This process integrates mathematical optimization and FEM.).
+    femopt = FEMOpt(opt=opt)
     
-    # 設計変数の設定
-    femopt.add_parameter('ant_r', 10, 5, 20, memo='円形アンテナの半径')
-    femopt.add_parameter('sx', 50, 40, 60, memo='基板のサイズ')
-    femopt.add_parameter('xf', 5, 1, 20, memo='給電ポートの偏心量')
-
-    # 拘束の設定
-    femopt.add_constraint(antenna_is_smaller_than_substrate, 'アンテナサイズ', lower_bound=1)
-    femopt.add_constraint(port_is_inside_antenna, 'ポート位置', lower_bound=1)
-
-    # 目的の設定
-    femopt.add_objective(s.get_resonance_frequency, '第一共振周波数（GHz）', direction=3.0)
+    # Add design variables (Use variable names set in Femtet) to the optimization problem.
+    femopt.add_parameter('ant_r', 10, 5, 20)
+    femopt.add_parameter('sx', 50, 40, 60)
+    femopt.add_parameter('xf', 5, 1, 20)
+
+    # Add constraint to the optimization problem.
+    femopt.add_constraint(antenna_is_smaller_than_substrate, 'board_antenna_clearance', lower_bound=1)
+    femopt.add_constraint(port_is_inside_antenna, 'antenna_port_clearance', lower_bound=1)
+
+    # Add objective to the optimization problem.
+    # The target frequency is 3 GHz.
+    femopt.add_objective(s.get_resonance_frequency, 'First_resonant_frequency(GHz)', direction=3.0)
 
     femopt.set_random_seed(42)
     femopt.optimize(n_trials=20)

pyfemtet/FemtetPJTSample/wat_ex14_parallel_parametric.py

@@ -1,23 +1,26 @@
-"""基板上の発熱体
+"""Parallel computing / Multi-objective optimization: heating element on board
 
-wat_ex14_parametric.femprj に対し熱伝導解析を行い、
-チップの発熱による温度上昇を最小にしつつ
-基板寸法を最小にする
-基板寸法・チップ配置寸法を探索します。
+Perform thermal conduction analysis on wat_ex14_parametric.femprj
+and search for board dimensions and chip placement dimensions that
+minimize the board dimensions while minimizing the temperature rise
+due to chip heat generation.
 """
-
 from pyfemtet.opt import FEMOpt
 
 
 def max_temperature(Femtet, body_name):
-    """Femtet の解析結果からチップの最高温度を取得します。
+    """Get the maximum temperature of the chip.
 
-    Femtet : マクロを使用するためのインスタンスです。詳しくは "Femtet マクロヘルプ / CFemtet クラス" をご覧ください。
-        目的関数は第一引数に Femtet インスタンスを取る必要があります。
+    Note:
+        The objective or constraint function
+        must take a Femtet as its first argument
+        and must return a single float.
 
-    max_temp : 計算された最高温度です。
-        目的関数は単一の float を返す必要があります。
+    Params:
+        Femtet: An instance for using Femtet macros. For more information, see "Femtet Macro Help / CFemtet Class".
 
+    Returns:
+        float: Max-temperature.
     """
     Gogh = Femtet.Gogh
 
@@ -27,14 +30,13 @@ def max_temperature(Femtet, body_name):
 
 
 def substrate_size(Femtet):
-    """Femtet の設計変数から基板サイズを取得します。
-
-    Femtet : マクロを使用するためのインスタンスです。詳しくは "Femtet マクロヘルプ / CFemtet クラス" をご覧ください。
-        目的関数は第一引数に Femtet インスタンスを取る必要があります。
-
-    subs_w * subs_d : XY 平面における基板の占有面積です。
-        目的関数は単一の float を返す必要があります。
+    """Calculate the substrate size.
 
+    Params:
+        Femtet: An instance for using Femtet macros. For more information, see "Femtet Macro Help / CFemtet Class".
+    
+    Returns:
+        float: The area occupied by the board in the XY plane.
     """
     subs_w = Femtet.GetVariableValue('substrate_w')
     subs_d = Femtet.GetVariableValue('substrate_d')
@@ -44,20 +46,20 @@ def substrate_size(Femtet):
 
 if __name__ == '__main__':
 
-    # 最適化処理を行うオブジェクトを用意
+    # Define FEMOpt object (This process integrates mathematical optimization and FEM.).
     femopt = FEMOpt()
 
-    # 設計変数の設定
-    femopt.add_parameter("substrate_w", 40, lower_bound=22, upper_bound=40, memo='基板サイズ X')
-    femopt.add_parameter("substrate_d", 60, lower_bound=33, upper_bound=60, memo='基板サイズ Y')
+    # Add design variables (Use variable names set in Femtet) to the optimization problem.
+    femopt.add_parameter("substrate_w", 40, lower_bound=22, upper_bound=40)
+    femopt.add_parameter("substrate_d", 60, lower_bound=33, upper_bound=60)
 
-    # 目的関数の設定
-    femopt.add_objective(max_temperature, name='メインチップ温度', args=('MAINCHIP',))
-    femopt.add_objective(max_temperature, name='サブチップ温度', args=('SUBCHIP',))
-    femopt.add_objective(substrate_size, name='基板サイズ')
+    # Add objective to the optimization problem.
+    femopt.add_objective(max_temperature, name='main chip temp.', args=('MAINCHIP',))
+    femopt.add_objective(max_temperature, name='sub chip temp.', args=('SUBCHIP',))
+    femopt.add_objective(substrate_size, name='substrate size')
 
-    # 最適化の実行
+    # Run optimization.
     femopt.set_random_seed(42)
-    # femopt.main(n_trials=20)
-    femopt.optimize(n_trials=20, n_parallel=3)  # ここのみ wat_ex14_parametric.py から変更しました。
+    # femopt.optimize(n_trials=20)
+    femopt.optimize(n_trials=20, n_parallel=3)  # Change only this line.
     femopt.terminate_all()

pyfemtet/FemtetPJTSample/wat_ex14_parametric.py

@@ -1,23 +1,26 @@
-"""基板上の発熱体
+"""Multi-objective optimization: heating element on board
 
-wat_ex14_parametric.femprj に対し熱伝導解析を行い、
-チップの発熱による温度上昇を最小にしつつ
-基板寸法を最小にする
-基板寸法・チップ配置寸法を探索します。
+Perform thermal conduction analysis on wat_ex14_parametric.femprj
+and search for board dimensions and chip placement dimensions that
+minimize the board dimensions while minimizing the temperature rise
+due to chip heat generation.
 """
-
 from pyfemtet.opt import FEMOpt
 
 
 def max_temperature(Femtet, body_name):
-    """Femtet の解析結果からチップの最高温度を取得します。
+    """Get the maximum temperature of the chip.
 
-    Femtet : マクロを使用するためのインスタンスです。詳しくは "Femtet マクロヘルプ / CFemtet クラス" をご覧ください。
-        目的関数は第一引数に Femtet インスタンスを取る必要があります。
+    Note:
+        The objective or constraint function
+        must take a Femtet as its first argument
+        and must return a single float.
 
-    max_temp : 計算された最高温度です。
-        目的関数は単一の float を返す必要があります。
+    Params:
+        Femtet: An instance for using Femtet macros. For more information, see "Femtet Macro Help / CFemtet Class".
 
+    Returns:
+        float: Max-temperature.
     """
     Gogh = Femtet.Gogh
 
@@ -27,14 +30,13 @@ def max_temperature(Femtet, body_name):
 
 
 def substrate_size(Femtet):
-    """Femtet の設計変数から基板サイズを取得します。
-
-    Femtet : マクロを使用するためのインスタンスです。詳しくは "Femtet マクロヘルプ / CFemtet クラス" をご覧ください。
-        目的関数は第一引数に Femtet インスタンスを取る必要があります。
-
-    subs_w * subs_d : XY 平面における基板の占有面積です。
-        目的関数は単一の float を返す必要があります。
+    """Calculate the substrate size.
 
+    Params:
+        Femtet: An instance for using Femtet macros. For more information, see "Femtet Macro Help / CFemtet Class".
+    
+    Returns:
+        float: The area occupied by the board in the XY plane.
     """
     subs_w = Femtet.GetVariableValue('substrate_w')
     subs_d = Femtet.GetVariableValue('substrate_d')
@@ -44,19 +46,19 @@ def substrate_size(Femtet):
 
 if __name__ == '__main__':
 
-    # 最適化処理を行うオブジェクトを用意
+    # Define FEMOpt object (This process integrates mathematical optimization and FEM.).
     femopt = FEMOpt()
 
-    # 設計変数の設定
-    femopt.add_parameter("substrate_w", 40, lower_bound=22, upper_bound=40, memo='基板サイズ X')
-    femopt.add_parameter("substrate_d", 60, lower_bound=33, upper_bound=60, memo='基板サイズ Y')
+    # Add design variables (Use variable names set in Femtet) to the optimization problem.
+    femopt.add_parameter("substrate_w", 40, lower_bound=22, upper_bound=40)
+    femopt.add_parameter("substrate_d", 60, lower_bound=33, upper_bound=60)
 
-    # 目的関数の設定
-    femopt.add_objective(max_temperature, name='メインチップ温度', args=('MAINCHIP',))
-    femopt.add_objective(max_temperature, name='サブチップ温度', args=('SUBCHIP',))
-    femopt.add_objective(substrate_size, name='基板サイズ')
+    # Add objective to the optimization problem.
+    femopt.add_objective(max_temperature, name='main chip temp.', args=('MAINCHIP',))
+    femopt.add_objective(max_temperature, name='sub chip temp.', args=('SUBCHIP',))
+    femopt.add_objective(substrate_size, name='substrate size')
 
-    # 最適化の実行
+    # Run optimization.
     femopt.set_random_seed(42)
     femopt.optimize(n_trials=20)
     femopt.terminate_all()

pyfemtet/__init__.py

@@ -1 +1 @@
-__version__ = "0.4.1"
+__version__ = "0.4.3"

pyfemtet/opt/_femopt_core.py

@@ -374,7 +374,6 @@ class History:
     prm_names = []
     obj_names = []
     cns_names = []
-    local_data = pd.DataFrame()
     is_restart = False
     is_processing = False
     _future = None
@@ -397,6 +396,9 @@ class History:
         self.cns_names = cns_names
         self.additional_metadata = additional_metadata or ''
 
+        # 初期化
+        self.local_data = pd.DataFrame()
+
         # 最適化実行中かどうか
         self.is_processing = client is not None
 

pyfemtet/opt/femprj_sample/cad_ex01_NX.py

@@ -0,0 +1,132 @@
+"""External CAD (NX) Integration
+
+Using Femtet's stress analysis solver and Siemens' CAD software NX,
+design a lightweight and high-strength H-shaped beam.
+
+As a preliminary step, please perform the following procedures:
+- Install NX
+- Create a C:\temp folder
+    - Note: NX will save a .x_t file in this folder.
+- Place the following files in the same folder:
+    - cad_ex01_NX.py (this file)
+    - cad_ex01_NX.prt
+    - cad_ex01_NX.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):
+    """Obtain the maximum von Mises stress of the model.
+
+    Note:
+        The objective or constraint function should take Femtet
+        as its first argument and return a float as the output.
+
+    Warning:
+        CAD integration may assign boundary conditions to unintended locations.
+
+        In this example, if the boundary conditions are assigned as intended,
+        the maximum z displacement is always negative.
+        If the maximum displacement is not negative, it is assumed that
+        boundary condition assignment has failed.
+        Then this function raises a ModelError.
+
+        If a ModelError, MeshError, or SolveError occurs during optimization,
+        the optimization process considers the attempt a failure and skips to
+        the next trial.
+    """
+
+    # Simple check for the correctness of boundary conditions.
+    dx, dy, dz = Femtet.Gogh.Galileo.GetMaxDisplacement_py()
+    if dz >= 0:
+        raise ModelError('Assigning unintended boundary conditions.')
+
+    # Von Mises stress calculation.
+    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):
+    """Obtain model mass."""
+    return Femtet.Gogh.Galileo.GetMass('H_beam')
+
+
+def C_minus_B(Femtet, opt):
+    """Calculate the difference between C and B dimensions.
+
+    Another example uses the following snippet to access design variables:
+
+        A = Femtet.GetVariableValue('A')
+    
+    However, when performing CAD integration, this method does not work
+    because the variables are not set in the .femprj file.
+
+    In CAD integration, design variables are obtained in the following way.
+
+        # How to obtain a dictionary with the variable names of parameters
+        # added by add_parameter() as keys.
+        params: dict = opt.get_parameter()
+        A = params['A']
+
+    Or
+
+        # How to obtain an array of values of parameters added in the order
+        # by add_parameter().
+        values: np.ndarray = opt.get_parameter('values')
+        A, B, C = values
+
+    Objective functions and constraint functions can take arbitrary variables
+    after the first argument.
+    The FEMOpt member variable `opt` has a method called get_parameter().
+    This method allows you to retrieve design variables added by add_parameter().
+    By taking `opt` as the second argument, you can execute get_parameter()
+    within the objective or constraint function to retrieve design variables.
+    """
+    A, B, C = opt.get_parameter('values')
+    return C - B
+
+
+if __name__ == '__main__':
+
+    # Initialize NX-Femtet integration object.
+    # At this point, Python is connected to the Femtet.
+    fem = FemtetWithNXInterface(
+        prt_path='cad_ex01_NX.prt',
+        open_result_with_gui=False,  # To calculate von Mises stress, set this argument to False. See Femtet Macro Help.
+    )
+
+    # Initialize the FEMOpt object.
+    # (establish connection between the optimization problem and Femtet)
+    femopt = FEMOpt(fem=fem)
+
+    # Add design variables to the optimization problem.
+    # (Specify the variables registered in the femprj file.)
+    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)
+
+    # Add the constraint function to the optimization problem.
+    femopt.add_constraint(C_minus_B, 'C>B', lower_bound=1, args=femopt.opt)
+
+    # Add the objective function to the optimization problem.
+    femopt.add_objective(von_mises, name='von Mises (Pa)')
+    femopt.add_objective(mass, name='mass (kg)')
+
+    # Run optimization.
+    femopt.set_random_seed(42)
+    femopt.optimize(n_trials=20)
+    femopt.terminate_all()

pyfemtet/opt/femprj_sample/cad_ex01_SW.py

@@ -0,0 +1,132 @@
+"""External CAD (SOLIDWORKS) Integration
+
+Using Femtet's stress analysis solver and Dassault Systemes' CAD software SOLIDWORKS,
+design a lightweight and high-strength H-shaped beam.
+
+As a preliminary step, please perform the following procedures:
+- Install SOLIDWORKS
+- Create a C:\temp folder
+    - Note: SOLIDWORKS will save a .x_t file in this folder.
+- Place the following files in the same folder:
+    - cad_ex01_SW.py (this file)
+    - cad_ex01_SW.SLDPRT
+    - cad_ex01_SW.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):
+    """Obtain the maximum von Mises stress of the model.
+
+    Note:
+        The objective or constraint function should take Femtet
+        as its first argument and return a float as the output.
+
+    Warning:
+        CAD integration may assign boundary conditions to unintended locations.
+
+        In this example, if the boundary conditions are assigned as intended,
+        the maximum z displacement is always negative.
+        If the maximum displacement is not negative, it is assumed that
+        boundary condition assignment has failed.
+        Then this function raises a ModelError.
+
+        If a ModelError, MeshError, or SolveError occurs during optimization,
+        the optimization process considers the attempt a failure and skips to
+        the next trial.
+    """
+
+    # Simple check for the correctness of boundary conditions.
+    dx, dy, dz = Femtet.Gogh.Galileo.GetMaxDisplacement_py()
+    if dz >= 0:
+        raise ModelError('Assigning unintended boundary conditions.')
+
+    # Von Mises stress calculation.
+    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):
+    """Obtain model mass."""
+    return Femtet.Gogh.Galileo.GetMass('H_beam')
+
+
+def C_minus_B(Femtet, opt):
+    """Calculate the difference between C and B dimensions.
+
+    Another example uses the following snippet to access design variables:
+
+        A = Femtet.GetVariableValue('A')
+    
+    However, when performing CAD integration, this method does not work
+    because the variables are not set in the .femprj file.
+
+    In CAD integration, design variables are obtained in the following way.
+
+        # How to obtain a dictionary with the variable names of parameters
+        # added by add_parameter() as keys.
+        params: dict = opt.get_parameter()
+        A = params['A']
+
+    Or
+
+        # How to obtain an array of values of parameters added in the order
+        # by add_parameter().
+        values: np.ndarray = opt.get_parameter('values')
+        A, B, C = values
+
+    Objective functions and constraint functions can take arbitrary variables
+    after the first argument.
+    The FEMOpt member variable `opt` has a method called get_parameter().
+    This method allows you to retrieve design variables added by add_parameter().
+    By taking `opt` as the second argument, you can execute get_parameter()
+    within the objective or constraint function to retrieve design variables.
+    """
+    A, B, C = opt.get_parameter('values')
+    return C - B
+
+
+if __name__ == '__main__':
+
+    # Initialize NX-Femtet integration object.
+    # At this point, Python is connected to the Femtet.
+    fem = FemtetWithSolidworksInterface(
+        sldprt_path='cad_ex01_SW.SLDPRT',
+        open_result_with_gui=False,  # To calculate von Mises stress, set this argument to False. See Femtet Macro Help.
+    )
+
+    # Initialize the FEMOpt object.
+    # (establish connection between the optimization problem and Femtet)
+    femopt = FEMOpt(fem=fem)
+
+    # Add design variables to the optimization problem.
+    # (Specify the variables registered in the femprj file.)
+    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)
+
+    # Add the constraint function to the optimization problem.
+    femopt.add_constraint(C_minus_B, 'C>B', lower_bound=1, args=femopt.opt)
+
+    # Add the objective function to the optimization problem.
+    femopt.add_objective(von_mises, name='von Mises (Pa)')
+    femopt.add_objective(mass, name='mass (kg)')
+
+    # Run optimization.
+    femopt.set_random_seed(42)
+    femopt.optimize(n_trials=20)
+    femopt.terminate_all()