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

pyfemtet/opt/samples/femprj_sample/cad_ex01_NX.py

@@ -1,135 +0,0 @@
-"""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.
-        export_curves=False,
-        export_surfaces=False,
-        export_solids=True,
-        export_flattened_assembly=False,
-    )
-
-    # 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)

pyfemtet/opt/samples/femprj_sample/cad_ex01_NX_test_result.reccsv

@@ -1,23 +0,0 @@
-"{""femprj_path"": null, ""model_name"": ""Analysis Model""}",prm,prm_lb,prm_ub,prm,prm_lb,prm_ub,prm,prm_lb,prm_ub,obj,obj_direction,obj,obj_direction,,cns,cns_lb,cns_ub,,,,
-,,,,,,,,,,,,,,,,,,,,,
-trial,A,A_lower_bound,A_upper_bound,B,B_lower_bound,B_upper_bound,C,C_lower_bound,C_upper_bound,von Mises (Pa),von Mises (Pa)_direction,mass (kg),mass (kg)_direction,non_domi,C>B,C>B_lower_bound,C>B_upper_bound,feasible,hypervolume,message,time
-1,10.0,1.0,59.0,10.0,1.0,40.0,20.0,5.0,59.0,12187480208.414427,minimize,5.615241407210086,minimize,True,10.0,1,,True,0.0,initial,2024-07-24 09:50:32.097923
-2,35.722192083428126,1.0,59.0,7.084726977255024,1.0,40.0,13.423704098154943,5.0,59.0,7519676665.53638,minimize,8.728900507395284,minimize,False,6.338977120899918,1,,True,0.0,,2024-07-24 09:50:46.512275
-3,42.06820951217064,1.0,59.0,1.8027952775362954,1.0,40.0,57.37513201674769,5.0,59.0,22480675646.096836,minimize,9.011481291190496,minimize,False,55.572336739211394,1,,True,0.0,,2024-07-24 09:51:01.360991
-4,49.28167316642446,1.0,59.0,9.28122531645277,1.0,40.0,14.818548229183433,5.0,59.0,4473829483.306867,minimize,11.38617808276437,minimize,True,5.537322912730662,1,,True,4203391875.452818,,2024-07-24 09:51:10.100850
-5,11.637461571499161,1.0,59.0,12.865447475421972,1.0,40.0,33.33684730814085,5.0,59.0,9250865273.829187,minimize,7.638439032091176,minimize,True,20.471399832718873,1,,True,7878103684.844186,,2024-07-24 09:51:18.443756
-6,26.052811081242716,1.0,59.0,12.357936467723635,1.0,40.0,38.04005631500849,5.0,59.0,5608846597.7734165,minimize,10.015718884236131,minimize,True,25.68211984728486,1,,True,18883770151.979397,,2024-07-24 09:51:26.735550
-7,9.090643917818426,1.0,59.0,12.393641292873507,1.0,40.0,24.78353953785935,5.0,59.0,9677909990.12775,minimize,6.4477508178425245,minimize,True,12.389898244985844,1,,True,23875008147.791668,,2024-07-24 09:51:34.794313
-8,21.441337720814797,1.0,59.0,18.072617642299363,1.0,40.0,45.21168857033972,5.0,59.0,4042233022.9095497,minimize,11.546379519585152,minimize,True,27.139070928040354,1,,True,26863123829.53503,,2024-07-24 09:51:50.555086
-9,37.73538456111445,1.0,59.0,2.051023324497365,1.0,40.0,5.434182199522157,5.0,59.0,49467429429.04644,minimize,7.690690579272188,minimize,False,3.383158875024792,1,,True,26863123829.53503,,2024-07-24 09:52:00.042133
-10,34.35890124528135,1.0,59.0,6.2478110218036935,1.0,40.0,23.49586366183886,5.0,59.0,7426521617.115658,minimize,8.713041181770238,minimize,True,17.24805264003517,1,,True,26999590015.41306,,2024-07-24 09:52:08.391160
-11,33.07122876558135,1.0,59.0,16.015501119793537,1.0,40.0,24.772488687883417,5.0,59.0,4967659892.43169,minimize,10.571545999177301,minimize,True,8.75698756808988,1,,True,29376121818.648865,,2024-07-24 09:52:17.118345
-12,44.49232391461052,1.0,59.0,6.594534954244819,1.0,40.0,31.754020417753264,5.0,59.0,5722567823.730824,minimize,10.860949688986132,minimize,False,25.159485463508446,1,,True,29376121818.648865,,2024-07-24 09:52:25.653325
-13,58.768958841351264,1.0,59.0,5.331491062142676,1.0,40.0,29.653341259152,5.0,59.0,6504123574.132211,minimize,12.846059384800473,minimize,False,24.321850197009326,1,,True,29376121818.648865,,2024-07-24 09:52:34.300756
-14,41.931667740603416,1.0,59.0,15.432432848846828,1.0,40.0,25.249518595235475,5.0,59.0,3944370145.7178416,minimize,11.794198824024951,minimize,True,9.817085746388647,1,,True,29898453080.38972,,2024-07-24 09:52:49.958108
-15,36.816310624574804,1.0,59.0,10.138554917176371,1.0,40.0,18.890867424341323,5.0,59.0,5890926471.473402,minimize,9.793323153548432,minimize,True,8.752312507164952,1,,True,29905143965.590454,,2024-07-24 09:52:58.152601
-16,50.4663695486681,1.0,59.0,5.90091414254359,1.0,40.0,29.425620895390857,5.0,59.0,6756288747.227509,minimize,11.575200835069767,minimize,False,23.524706752847266,1,,True,29905143965.590454,,2024-07-24 09:53:06.781909
-17,33.00960738761119,1.0,59.0,3.647820207249005,1.0,40.0,57.34057201063078,5.0,59.0,11701146103.834673,minimize,8.422141022544057,minimize,False,53.692751803381775,1,,True,29905143965.590454,,2024-07-24 09:53:15.773832
-18,29.89239146401912,1.0,59.0,21.093464343210364,1.0,40.0,38.1038646366275,5.0,59.0,4562525040.83442,minimize,12.253735414187783,minimize,False,17.010400293417135,1,,True,29905143965.590454,,2024-07-24 09:53:24.116460
-19,39.65771866413639,1.0,59.0,11.726584835682985,1.0,40.0,44.63063781981377,5.0,59.0,4738124374.9653225,minimize,12.384862721051968,minimize,False,32.904052984130786,1,,True,29905143965.590454,,2024-07-24 09:53:32.782403
-20,10.443255153670096,1.0,59.0,7.644383519738286,1.0,40.0,19.30571184104179,5.0,59.0,12772033981.056494,minimize,4.957198530620907,minimize,True,11.661328321303506,1,,True,30246653770.046047,,2024-07-24 09:53:41.944991

pyfemtet/opt/samples/femprj_sample/cad_ex01_SW.py

@@ -1,131 +0,0 @@
-"""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)

pyfemtet/opt/samples/femprj_sample/cad_ex01_SW_test_result.reccsv

@@ -1,23 +0,0 @@
-"{""femprj_path"": null, ""model_name"": ""Analysis Model""}",prm,prm_lb,prm_ub,prm,prm_lb,prm_ub,prm,prm_lb,prm_ub,obj,obj_direction,obj,obj_direction,,cns,cns_lb,cns_ub,,,,
-,,,,,,,,,,,,,,,,,,,,,
-trial,A,A_lower_bound,A_upper_bound,B,B_lower_bound,B_upper_bound,C,C_lower_bound,C_upper_bound,von Mises (Pa),von Mises (Pa)_direction,mass (kg),mass (kg)_direction,non_domi,C>B,C>B_lower_bound,C>B_upper_bound,feasible,hypervolume,message,time
-1,10.0,1.0,59.0,10.0,1.0,40.0,20.0,5.0,59.0,10713610908.24202,minimize,5.615241407210081,minimize,True,10.0,1,,True,0.0,initial,2024-07-24 09:54:50.307882
-2,35.722192083428126,1.0,59.0,7.084726977255024,1.0,40.0,13.423704098154943,5.0,59.0,6451312936.468065,minimize,8.728900507478476,minimize,True,6.338977120899918,1,,True,145876137626.11978,,2024-07-24 09:54:53.688085
-3,42.06820951217064,1.0,59.0,1.8027952775362954,1.0,40.0,57.37513201674769,5.0,59.0,24065344414.778633,minimize,9.011481292235954,minimize,False,55.572336739211394,1,,True,145876137626.11978,,2024-07-24 09:54:57.979731
-4,49.28167316642446,1.0,59.0,9.28122531645277,1.0,40.0,14.818548229183433,5.0,59.0,5503645118.524387,minimize,11.386178083846733,minimize,False,5.537322912730662,1,,True,145876137626.11978,,2024-07-24 09:55:01.056168
-5,11.637461571499161,1.0,59.0,12.865447475421972,1.0,40.0,33.33684730814085,5.0,59.0,11438869243.373928,minimize,7.638439033366426,minimize,False,20.471399832718873,1,,True,145876137626.11978,,2024-07-24 09:55:03.859850
-6,26.052811081242716,1.0,59.0,12.357936467723635,1.0,40.0,38.04005631500849,5.0,59.0,5495681889.607063,minimize,10.01571888509523,minimize,True,25.68211984728486,1,,True,153516179853.76138,,2024-07-24 09:55:06.837354
-7,9.090643917818426,1.0,59.0,12.393641292873507,1.0,40.0,24.78353953785935,5.0,59.0,9483855035.111774,minimize,6.447750817595856,minimize,True,12.389898244985844,1,,True,153999396463.2158,,2024-07-24 09:55:09.733585
-8,37.99457827883769,1.0,59.0,1.5209697385859142,1.0,40.0,23.897679644327795,5.0,59.0,30406250170.856964,minimize,7.778758622447119,minimize,False,22.37670990574188,1,,True,153999396463.2158,,2024-07-24 09:55:13.314997
-9,34.35890124528135,1.0,59.0,6.2478110218036935,1.0,40.0,12.678079497056089,5.0,59.0,8761031684.825426,minimize,8.28729418409825,minimize,True,6.4302684752523955,1,,True,156804653691.83817,,2024-07-24 09:55:17.128276
-10,16.48548223083144,1.0,59.0,15.540089806204016,1.0,40.0,47.73934394149744,5.0,59.0,5443514483.782774,minimize,10.442763870961809,minimize,True,32.19925413529342,1,,True,157123857054.0115,,2024-07-24 09:55:20.740157
-11,31.61117579574186,1.0,59.0,6.349996259293115,1.0,40.0,19.822188557253973,5.0,59.0,9074520374.76584,minimize,8.13162247977363,minimize,True,13.47219229796086,1,,True,157127957770.34027,,2024-07-24 09:55:23.753393
-12,45.77159271804414,1.0,59.0,6.286691200972086,1.0,40.0,11.516115144447737,5.0,59.0,7146383394.662491,minimize,10.180750385294035,minimize,False,5.22942394347565,1,,True,157127957770.34027,,2024-07-24 09:55:27.065773
-13,25.856967098172454,1.0,59.0,13.177162500962542,1.0,40.0,39.81237900447662,5.0,59.0,7865665149.924187,minimize,10.363696138818305,minimize,False,26.635216503514076,1,,True,157127957770.34027,,2024-07-24 09:55:30.765081
-14,21.401950687804057,1.0,59.0,11.148262705918398,1.0,40.0,32.95590234333009,5.0,59.0,7221556874.706432,minimize,8.579324905603093,minimize,True,21.80763963741169,1,,True,157191679594.55545,,2024-07-24 09:55:33.521637
-15,12.670878704766015,1.0,59.0,4.379466144113291,1.0,40.0,10.277664459492886,5.0,59.0,28528452444.874542,minimize,3.9668848228126214,minimize,True,5.898198315379595,1,,True,157421947465.85098,,2024-07-24 09:55:36.562223
-16,7.797785002291775,1.0,59.0,8.535545896262494,1.0,40.0,57.34057201063078,5.0,59.0,13421903355.914164,minimize,6.845257989137097,minimize,False,48.80502611436829,1,,True,157421947465.85098,,2024-07-24 09:55:39.711790
-17,6.384676158645501,1.0,59.0,2.9595622079140185,1.0,40.0,21.66854065738089,5.0,59.0,40100488412.84058,minimize,2.582077175952352,minimize,True,18.708978449466873,1,,True,177068059660.24457,,2024-07-24 09:55:43.519329
-18,1.2749269398791796,1.0,59.0,4.355873741157909,1.0,40.0,50.335120552067295,5.0,59.0,66671888609.68277,minimize,3.061340758008378,minimize,False,45.97924681090939,1,,True,177068059660.24457,,2024-07-24 09:55:46.497664
-19,6.676161230791632,1.0,59.0,3.799493060906398,1.0,40.0,21.166583460985745,5.0,59.0,29298549516.19457,minimize,2.980158111433043,minimize,True,17.36709040007935,1,,True,177547992085.38452,,2024-07-24 09:55:50.730197
-20,8.206418867626418,1.0,59.0,9.303411630645828,1.0,40.0,25.37658417274892,5.0,59.0,13800596362.955494,minimize,5.449166488749558,minimize,True,16.073172542103094,1,,True,188206553729.3952,,2024-07-24 09:55:53.707442

pyfemtet/opt/samples/femprj_sample/constrained_pipe.py

@@ -1,96 +0,0 @@
-"""A sample to implement constrained optimization.
-
-This section describes the types of constraints and
-the steps to run optimization on models that require them.
-
-"""
-
-from pyfemtet.opt import FEMOpt, OptunaOptimizer
-from pyfemtet.opt.optimizer import PoFBoTorchSampler
-
-
-def mises_stress(Femtet):
-    """Calculate the von Mises stress as the objective function.
-
-    This function is called automatically by the FEMOpt
-    object while the optimization is running.
-
-    Args:
-        Femtet: When defining an objective or constraint
-            function using PyFemtet, the first argument
-            must take a Femtet instance.
-
-    Returns:
-        float: A single float representing the expression value you want to constrain.
-    """
-    return Femtet.Gogh.Galileo.GetMaxStress_py()[2]
-
-
-def radius_diff(Femtet, opt):
-    """Calculate the difference between the outer and inner radii of the pipe.
-
-    This constraint is called to ensure that the
-    inner radius of the pipe does not exceed the
-    outer radius while the optimization is running.
-
-    Note:
-        If you are using BoTorchSampler of OptunaOptimizer
-        and use strict constraints, be aware that accessing
-        the Femtet can be very slow, as it requires repeated
-        calculations to propose parameters.
-        We recommend that you do not access the Femtet,
-        but rather get the parameters and perform the
-        calculations via the Optimizer object, as in this
-        function example.
-
-        NOT recommended::
-
-            p = Femtet.GetVariableValue('p')
-
-        instead, use optimizer::
-
-            params = opt.get_parameter()
-            p = params['p']
-
-    Args:
-        Femtet: When defining an objective or constraint
-            function using PyFemtet, the first argument
-            must take a Femtet instance.
-        opt: This object allows you to obtain the outer
-            radius and inner radius values without going
-            through Femtet.
-    """
-    params = opt.get_parameter()
-    internal_r = params['internal_r']
-    external_r = params['external_r']
-    return external_r - internal_r
-
-if __name__ == '__main__':
-    # Setup optimization method
-    opt = OptunaOptimizer(
-        sampler_class=PoFBoTorchSampler,
-        sampler_kwargs=dict(
-            n_startup_trials=3,  # The first three samples are randomly sampled.
-        )
-    )
-    femopt = FEMOpt(opt=opt)
-
-    # Add parameters
-    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)
-
-    # Add the strict constraint not to exceed the
-    # outer radius while the optimization is running.
-    femopt.add_constraint(
-        radius_diff,  # Constraint function (returns external radius - internal radius).
-        name='wall thickness',  # You can name the function anything you want.
-        lower_bound=1,  # Lower bound of constraint function (set minimum wall thickness is 1).
-        args=(femopt.opt,)  # Additional arguments passed to the function.
-    )
-
-    # Add the objective
-    femopt.add_objective(mises_stress, name='Mises Stress')
-
-    # Run optimization.
-    femopt.set_random_seed(42)
-    femopt.optimize(n_trials=10)

pyfemtet/opt/samples/femprj_sample/constrained_pipe_test_result.reccsv

@@ -1,13 +0,0 @@
-"",prm,prm_lb,prm_ub,prm,prm_lb,prm_ub,obj,obj_direction,,cns,cns_lb,cns_ub,,,,
-,,,,,,,,,,,,,,,,
-trial,external_r,external_r_lower_bound,external_r_upper_bound,internal_r,internal_r_lower_bound,internal_r_upper_bound,Mises Stress,Mises Stress_direction,non_domi,wall thickness,wall thickness_lower_bound,wall thickness_upper_bound,feasible,hypervolume,message,time
-1,10.0,0.1,10.0,5.0,0.1,10.0,2369132.9299968677,minimize,False,5.0,1,,True,-1.0,initial,2024-09-17 14:52:32.436730
-2,7.34674002393291,0.1,10.0,6.026718993550662,0.1,10.0,5797829.750527166,minimize,False,1.3200210303822484,1,,True,-1.0,,2024-09-17 14:52:42.927138
-3,8.341182143924176,0.1,10.0,2.202157195714934,0.1,10.0,1977631.590419345,minimize,False,6.139024948209242,1,,True,-1.0,,2024-09-17 14:52:54.382893
-4,9.999999999999995,0.1,10.0,0.6764719819158711,0.1,10.0,1514026.80705058,minimize,False,9.323528018084124,1,,True,-1.0,,2024-09-17 14:53:32.128331
-5,7.61174007452904,0.1,10.0,0.10000000000001566,0.1,10.0,398224.5447358758,minimize,True,7.511740074529024,1,,True,-1.0,,2024-09-17 14:54:08.541865
-6,5.026961232239628,0.1,10.0,0.10000000000001101,0.1,10.0,485322.32959281537,minimize,False,4.926961232239617,1,,True,-1.0,,2024-09-17 14:54:39.695340
-7,1.1000000916874046,0.1,10.0,0.1,0.1,10.0,613872.6231008903,minimize,False,1.0000000916874046,1,,True,-1.0,,2024-09-17 14:55:07.711879
-8,3.122874622702614,0.1,10.0,0.1,0.1,10.0,497544.0108231131,minimize,False,3.0228746227026138,1,,True,-1.0,,2024-09-17 14:55:41.036763
-9,9.999999999999998,0.1,10.0,3.7062659141309533,0.1,10.0,2092780.7380693094,minimize,False,6.293734085869045,1,,True,-1.0,,2024-09-17 14:56:28.269065
-10,6.601571760977393,0.1,10.0,0.10000000004171199,0.1,10.0,456263.04356135987,minimize,False,6.501571760935681,1,,True,-1.0,,2024-09-17 14:57:04.427906

pyfemtet/opt/samples/femprj_sample/gal_ex58_parametric.py

@@ -1,74 +0,0 @@
-"""Single-objective optimization: bending with consideration for springback.
-
-Using Femtet's stress analysis solver, we will determine the bending angle
-to achieve the desired material bend angle with consideration for springback.
-Elasto-plastic analysis is available in an optional package.
-
-Corresponding project: gal_ex58_parametric.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):
-    """Get the material bend angle.
-
-    Note:
-        The objective or constraint function should take Femtet
-        as its first argument and return a float as the output.
-
-    Params:
-        Femtet: This is an instance for manipulating Femtet with macros. For detailed information, please refer to "Femtet Macro Help".
-
-    Returns:
-        float: material bend angle.
-    """
-    Gogh = Femtet.Gogh
-
-    # Set the mode after unloading.
-    Gogh.Galileo.Mode = Gogh.Galileo.nMode - 1
-
-    # Obtain the displacement of the measurement target point.
-    Gogh.Galileo.Vector = constants.GALILEO_DISPLACEMENT_C
-    succeed, (x, y, z) = Gogh.Galileo.GetVectorAtPoint_py(200, 0, 0)
-
-    # Calculate the angle formed by the line segment
-    # connecting the bending origin (100, 0) and the
-    # deformed point with the X-axis.
-    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)  # unit: degree
-
-
-if __name__ == '__main__':
-
-    # Initialize the numerical optimization problem.
-    # (determine the optimization method)
-    opt = OptunaOptimizer(
-        sampler_class=BoTorchSampler,
-        sampler_kwargs=dict(
-            n_startup_trials=3,
-        )
-    )
-
-    # Initialize the FEMOpt object.
-    # (establish connection between the optimization problem and Femtet)
-    femopt = FEMOpt(opt=opt)
-
-    # Add design variables to the optimization problem.
-    # (Specify the variables registered in the femprj file.)
-    femopt.add_parameter("rot", 90, lower_bound=80, upper_bound=100)
-
-    # Add the objective function to the optimization problem.
-    # The target bending angle is 90 degrees.
-    femopt.add_objective(bending, name='final angle (degree)', direction=90)
-
-    # Run optimization.
-    femopt.set_random_seed(42)
-    femopt.optimize(n_trials=10)

pyfemtet/opt/samples/femprj_sample/gal_ex58_parametric_test_result.reccsv

@@ -1,13 +0,0 @@
-"",prm,prm_lb,prm_ub,obj,obj_direction,,,,,
-,,,,,,,,,,
-trial,rot,rot_lower_bound,rot_upper_bound,final angle (degree),final angle (degree)_direction,non_domi,feasible,hypervolume,message,time
-1,90.0,80.0,100.0,81.62540666214329,90,False,True,-1.0,initial,2024-07-24 10:57:24.248429
-2,87.49080237694724,80.0,100.0,79.02571753406903,90,False,True,-1.0,,2024-07-24 10:57:27.596882
-3,99.01428612819832,80.0,100.0,90.96145916607375,90,False,True,-1.0,,2024-07-24 10:57:31.010370
-4,99.99999999999999,80.0,100.0,91.89033160850201,90,False,True,-1.0,,2024-07-24 10:57:34.878708
-5,96.87624029068405,80.0,100.0,88.88281351640448,90,False,True,-1.0,,2024-07-24 10:57:38.443873
-6,97.86932823598721,80.0,100.0,89.88148176887373,90,False,True,-1.0,,2024-07-24 10:57:42.063686
-7,80.0,80.0,100.0,71.5309188592805,90,False,True,-1.0,,2024-07-24 10:57:45.658223
-8,97.6356229310297,80.0,100.0,89.64880326597084,90,False,True,-1.0,,2024-07-24 10:57:49.220276
-9,98.08539472199536,80.0,100.0,90.07549754912081,90,True,True,-1.0,,2024-07-24 10:57:52.882134
-10,93.84354938661555,80.0,100.0,85.71599925662619,90,False,True,-1.0,,2024-07-24 10:57:56.448858

pyfemtet/opt/samples/femprj_sample/gau_ex08_parametric.py

@@ -1,58 +0,0 @@
-"""Single-objective optimization: Self-inductance of a finite-length helical coil.
-
-Using Femtet's magnetic field analysis solver, design to achieve
-the target value for the self-inductance of a finite-length helical coil.
-
-Corresponding project: gau_ex08_parametric.femprj
-"""
-from optuna.integration.botorch import BoTorchSampler
-from pyfemtet.opt import FEMOpt, OptunaOptimizer
-
-
-def inductance(Femtet):
-    """Obtain the self-inductance.
-
-    Note:
-        The objective or constraint function should take Femtet
-        as its first argument and return a float as the output.
-
-    Params:
-        Femtet: This is an instance for manipulating Femtet with macros. For detailed information, please refer to "Femtet Macro Help".
-
-    Returns:
-        float: Self-inductance.
-    """
-    Gogh = Femtet.Gogh
-
-    coil_name = Gogh.Gauss.GetCoilList()[0]
-    return Gogh.Gauss.GetL(coil_name, coil_name)  # unit: F
-
-
-if __name__ == '__main__':
-
-    # Initialize the numerical optimization problem.
-    # (determine the optimization method)
-    opt = OptunaOptimizer(
-        sampler_class=BoTorchSampler,
-        sampler_kwargs=dict(
-            n_startup_trials=5,
-        )
-    )
-
-    # Initialize the FEMOpt object.
-    # (establish connection between the optimization problem and Femtet)
-    femopt = FEMOpt(opt=opt)
-
-    # Add design variables to the optimization problem.
-    # (Specify the variables registered in the femprj file.)
-    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)
-
-    # Add the objective function to the optimization problem.
-    # The target inductance is 0.1 uF.
-    femopt.add_objective(inductance, name='self-inductance (F)', direction=1e-7)
-
-    # Run optimization.
-    femopt.set_random_seed(42)
-    femopt.optimize(n_trials=20)

pyfemtet/opt/samples/femprj_sample/gau_ex08_parametric_test_result.reccsv

@@ -1,23 +0,0 @@
-"{""femprj_path"": null, ""model_name"": ""ex8""}",prm,prm_lb,prm_ub,prm,prm_lb,prm_ub,prm,prm_lb,prm_ub,obj,obj_direction,,,,,
-,,,,,,,,,,,,,,,,
-trial,helical_pitch,helical_pitch_lower_bound,helical_pitch_upper_bound,coil_radius,coil_radius_lower_bound,coil_radius_upper_bound,n_turns,n_turns_lower_bound,n_turns_upper_bound,self-inductance (F),self-inductance (F)_direction,non_domi,feasible,hypervolume,message,time
-1,6.0,4.2,8.0,10.0,1.0,10.0,5.0,1.0,5.0,2.2519131001615743e-07,1e-07,False,True,-1.0,initial,2024-07-24 09:58:15.358912
-2,5.623252451619978,4.2,8.0,9.556428757689245,1.0,10.0,3.9279757672456204,1.0,5.0,1.566960058810769e-07,1e-07,False,True,-1.0,,2024-07-24 09:58:19.258262
-3,4.420717726239158,4.2,8.0,8.795585311974417,1.0,10.0,3.404460046972835,1.0,5.0,1.248686556706198e-07,1e-07,False,True,-1.0,,2024-07-24 09:58:23.877333
-4,6.890675795624973,4.2,8.0,1.185260448662222,1.0,10.0,4.879639408647977,1.0,5.0,1.1068853339246818e-08,1e-07,False,True,-1.0,,2024-07-24 09:58:26.798210
-5,7.363282035041603,4.2,8.0,2.9110519961044856,1.0,10.0,1.7272998688284025,1.0,5.0,5.928345320411823e-09,1e-07,False,True,-1.0,,2024-07-24 09:58:29.046402
-6,4.506720164610959,4.2,8.0,8.847748295293755,1.0,10.0,3.369256043237217,1.0,5.0,1.2297792684115126e-07,1e-07,False,True,-1.0,,2024-07-24 09:58:33.300913
-7,4.726319192706824,4.2,8.0,8.983609183728657,1.0,10.0,3.2249528933303577,1.0,5.0,1.1498139019761039e-07,1e-07,False,True,-1.0,,2024-07-24 09:58:37.272948
-8,4.962549980320558,4.2,8.0,9.171627422480721,1.0,10.0,2.851110895485008,1.0,5.0,9.983020941290325e-08,1e-07,True,True,-1.0,,2024-07-24 09:58:41.746230
-9,5.050776505064675,4.2,8.0,9.300811078672245,1.0,10.0,2.4947121172622975,1.0,5.0,8.3155685279894e-08,1e-07,False,True,-1.0,,2024-07-24 09:58:45.329267
-10,5.000826194581878,4.2,8.0,9.199169607996211,1.0,10.0,2.8627720197139612,1.0,5.0,9.972903605466685e-08,1e-07,False,True,-1.0,,2024-07-24 09:58:49.405862
-11,4.95860020375077,4.2,8.0,9.178798062739858,1.0,10.0,2.8575303624554946,1.0,5.0,1.0072301643095888e-07,1e-07,False,True,-1.0,,2024-07-24 09:58:53.695530
-12,5.122546460811877,4.2,8.0,9.083931701720807,1.0,10.0,2.8617894824562873,1.0,5.0,9.647304150013922e-08,1e-07,False,True,-1.0,,2024-07-24 09:58:58.340271
-13,4.764485683472439,4.2,8.0,9.862539991891884,1.0,10.0,2.8711470086976076,1.0,5.0,1.1533811097380779e-07,1e-07,False,True,-1.0,,2024-07-24 09:59:02.339150
-14,4.968945595366566,4.2,8.0,9.161892011657752,1.0,10.0,2.8734870374362838,1.0,5.0,1.0056359149401165e-07,1e-07,False,True,-1.0,,2024-07-24 09:59:06.390307
-15,4.859845797581561,4.2,8.0,9.26151325278322,1.0,10.0,2.8664730982646485,1.0,5.0,1.0218819608040895e-07,1e-07,False,True,-1.0,,2024-07-24 09:59:10.251527
-16,4.973161997562328,4.2,8.0,9.184057459793724,1.0,10.0,2.8650256176749336,1.0,5.0,1.0070520417340285e-07,1e-07,False,True,-1.0,,2024-07-24 09:59:14.654137
-17,4.969144305852682,4.2,8.0,9.169512917581903,1.0,10.0,2.8604997178931355,1.0,5.0,9.943667954213088e-08,1e-07,False,True,-1.0,,2024-07-24 09:59:18.304112
-18,4.971098140849773,4.2,8.0,9.178412506598546,1.0,10.0,2.861253162496574,1.0,5.0,9.956172119870154e-08,1e-07,False,True,-1.0,,2024-07-24 09:59:22.125248
-19,4.971607724219878,4.2,8.0,9.179995137812535,1.0,10.0,2.861038617389135,1.0,5.0,1.0069749441418149e-07,1e-07,False,True,-1.0,,2024-07-24 09:59:26.682442
-20,4.9713095460170535,4.2,8.0,9.179603131437648,1.0,10.0,2.8608051498628724,1.0,5.0,9.93079291332147e-08,1e-07,False,True,-1.0,,2024-07-24 09:59:30.361099