h2lib-tests 13.1.304__py3-none-any.whl → 13.1.502__py3-none-any.whl

h2lib_tests/test_files/DTU_10_MW/htc/DTU_10MW_RWT_yaw_acturator.htc

@@ -0,0 +1,732 @@
+;DTU_10MW_RWT, version 9, 25-09-2017, mhha
+; 
+begin simulation;
+  time_stop    1000.0;
+  solvertype   1 ;    (newmark)
+  on_no_convergence continue ;
+  convergence_limits 1E3 1.0 1E-7 ;
+  logfile ./log/DTU_10MW_RWT.log ;
+  log_deltat .1;
+  begin newmark;
+    deltat    0.02;  
+  end newmark;
+  ; visualization visu.hdf5 0 30;
+end simulation;
+;
+;---------------------------------------------------------------------------------------------------------------------------------------------------------------- 
+begin new_htc_structure;
+;   beam_output_file_name  ./log/DTU_10MW_RWT_beam.dat;                    Optional - Calculated beam properties of the bodies are written to file
+;   body_output_file_name  ./log/DTU_10MW_RWT_body.dat;                    Optional - Body initial position and orientation are written to file
+;   body_eigenanalysis_file_name ./eig/DTU_10MW_RWT_body_eigen.dat;
+;   structure_eigenanalysis_file_name ./eig/DTU_10MW_RWT_strc_eigen.dat ;
+;-------------------------------------------------------------------------------------------------------------------------------
+;-------------------------------------------------------------------------------------------------------------------------------
+  begin main_body;         tower 115m
+    name        tower ;            
+    type        timoschenko ;
+    nbodies     1 ;
+    node_distribution     c2_def ;
+    damping_posdef   0.0 0.0 0.0 4.12E-03 4.12E-03 4.5E-04  ; Mx My Mz Kx Ky Kz , M´s raises overall level, K´s raises high freguency level "tuned by Larh"
+     begin timoschenko_input;
+      filename ./data/DTU_10MW_RWT_Tower_st.dat;
+      set 1 1 ; 
+    end timoschenko_input;
+    begin c2_def;              Definition of centerline (main_body coordinates)
+      nsec 11;
+      sec	1	0	0	0.00	0	;  x,y,z,twist
+      sec	2	0	0	-11.50	0	;
+      sec	3	0	0	-23.00	0	;
+      sec	4	0	0	-34.50	0	;
+      sec	5	0	0	-46.00	0	;
+      sec	6	0	0	-57.50	0	;
+      sec	7	0	0	-69.00	0	;
+      sec	8	0	0	-80.50	0	;
+      sec	9	0	0	-92.00	0	 ;
+      sec	10	0	0	-103.50	0	;
+      sec	11	0	0	-115.63	0	;	  
+     end c2_def ;
+    end main_body;
+;
+  begin main_body;
+    name        towertop ;              
+    type        timoschenko ;
+    nbodies     1 ;
+    node_distribution     c2_def ;
+    damping_posdef  0.0  0.0  0.0  7.00E-03  7.00E-03  7.00E-03  ;   "changed by Larh"	
+	concentrated_mass	2.0	0.0	2.6870E+00	3.0061E-01	4.4604E+05	4.1060E+06	4.1060E+05	4.1060E+06	;	Nacelle mass and inertia "corrected by Anyd 25/4/13"
+	begin timoschenko_input;
+      filename ./data/DTU_10MW_RWT_Towertop_st.dat ;
+      set 1 2 ;                
+    end timoschenko_input;
+    begin c2_def;              Definition of centerline (main_body coordinates)
+      nsec 2;
+      sec 1 0.0 0.0  0.0    0.0 ; x,y,z,twist
+      sec 2 0.0 0.0 -2.75   0.0 ; 
+    end c2_def ;
+  end main_body;
+;
+  begin main_body;
+    name        shaft ;              
+    type        timoschenko ;
+    nbodies     1 ;
+    node_distribution     c2_def ;
+	damping_posdef  0.0 0.0 0.0 4.65E-04  4.65E-04  3.983E-03 ; "tuned by Anyd 23/5/13 to 31.45 log decr. damping for free free with stiff rotor and tower"
+    concentrated_mass	1.0	0.0	0.0	0.0	0.0	0.0	0.0	3.751E+06	;	generator equivalent slow shaft "re_tuned by Anyd 20/2/13"  
+    concentrated_mass	5.0	0.0	0.0	0.0	1.0552E+05	0.0	0.0	3.257E+05	;	hub mass and inertia;	"re_tuned by Anyd 20/2/13"  
+	begin timoschenko_input;
+      filename ./data/DTU_10MW_RWT_Shaft_st.dat ;
+      set 1 1 ;                
+    end timoschenko_input;
+    begin c2_def;              Definition of centerline (main_body coordinates)
+      nsec 5;
+      sec 1 0.0 0.0 0.0 0.0 ; Tower top x,y,z,twist
+      sec 2 0.0 0.0 1.5 0.0 ; 
+      sec 3 0.0 0.0 3.0 0.0 ; 
+      sec 4 0.0 0.0 4.4	0.0 ; Main bearing
+      sec 5 0.0 0.0 7.1 0.0 ; Rotor centre
+    end c2_def ;
+  end main_body;	
+;
+  begin main_body;
+    name        hub1 ;              
+    type        timoschenko ;
+    nbodies     1 ;
+    node_distribution     c2_def ;
+    damping_posdef  0.0  0.0  0.0  3.00E-06  3.00E-06  2.00E-05;   "changed by Larh"	
+	begin timoschenko_input;
+      filename ./data/DTU_10MW_RWT_Hub_st.dat ;
+      set 1 2 ;                
+    end timoschenko_input;
+    begin c2_def;              Definition of centerline (main_body coordinates)
+      nsec 2;
+      sec 1 0.0 0.0 0.0 0.0 ; x,y,z,twist
+      sec 2 0.0 0.0 2.8 0.0 ; 
+    end c2_def ;
+  end main_body;
+;
+  begin main_body;
+    name           hub2 ;
+    copy_main_body hub1;
+  end main_body;
+;
+  begin main_body;
+    name           hub3 ;
+    copy_main_body hub1 ;
+  end main_body;
+;
+  begin main_body;
+    name        blade1 ;        
+    type        timoschenko ;
+    nbodies     10 ;
+    node_distribution    c2_def;
+	damping_posdef   0.0 0.0 0.0 1.53e-3 2.55e-3 3.3e-4 ; " 3% damping tuned by tkim 23/03/13 unable to fit 3rd and higher mode"
+    begin timoschenko_input ;
+      filename ./data/DTU_10MW_RWT_Blade_st.dat;
+      set 1 1 ;                set subset
+    end timoschenko_input;
+    begin c2_def;              Definition of centerline (main_body coordinates)
+		nsec 27 ;
+			sec	1	0.00000E+00	7.00600E-05	4.44089E-16	-1.45000E+01	;
+			sec	2	-2.06477E-05	-1.22119E-02	3.00000E+00	-1.45000E+01	;
+			sec	3	-7.28810E-03	-2.49251E-02	6.00000E+00	-1.44851E+01	;
+			sec	4	-1.89235E-02	-2.73351E-02	7.00004E+00	-1.44610E+01	;
+			sec	5	-5.41282E-02	-2.82163E-02	8.70051E+00	-1.43388E+01	;
+			sec	6	-1.26633E-01	-2.13210E-02	1.04020E+01	-1.40201E+01	;
+			sec	7	-2.25666E-01	-1.28378E-02	1.22046E+01	-1.33904E+01	;
+			sec	8	-2.88563E-01	-7.70659E-03	1.32065E+01	-1.29371E+01	;
+			sec	9	-3.99194E-01	-4.88317E-03	1.50100E+01	-1.19445E+01	;
+			sec	10	-5.76634E-01	-1.80296E-02	1.82151E+01	-9.98243E+00	;
+			sec	11	-7.07136E-01	-5.01772E-02	2.14178E+01	-8.45147E+00	;
+			sec	12	-7.91081E-01	-9.41228E-02	2.46189E+01	-7.46417E+00	;
+			sec	13	-8.37195E-01	-1.48880E-01	2.78193E+01	-6.72916E+00	;
+			sec	14	-8.53948E-01	-2.14514E-01	3.10194E+01	-6.08842E+00	;
+			sec	15	-8.49367E-01	-2.90618E-01	3.42197E+01	-5.49322E+00	;
+			sec	16	-7.93920E-01	-4.62574E-01	4.02204E+01	-4.39222E+00	;
+			sec	17	-7.16284E-01	-6.88437E-01	4.66217E+01	-3.09315E+00	;
+			sec	18	-6.34358E-01	-9.60017E-01	5.30232E+01	-1.75629E+00	;
+			sec	19	-5.53179E-01	-1.28424E+00	5.94245E+01	-5.00650E-01	;
+			sec	20	-4.75422E-01	-1.66402E+00	6.58255E+01	6.01964E-01	;
+			sec	21	-4.03180E-01	-2.10743E+00	7.22261E+01	1.55560E+00	;
+			sec	22	-3.30085E-01	-2.65630E+00	7.90266E+01	2.51935E+00	;
+			sec	23	-3.10140E-01	-2.78882E+00	8.05267E+01	2.72950E+00	;
+			sec	24	-2.86719E-01	-2.92517E+00	8.20271E+01	2.93201E+00	;
+			sec	25	-2.55823E-01	-3.06577E+00	8.35274E+01	3.11874E+00	;
+			sec	26	-2.07891E-01	-3.20952E+00	8.50277E+01	3.28847E+00	;
+			sec	27	-8.98940E-02	-3.33685E+00	8.63655E+01	3.42796E+00	;
+     end c2_def ;                                 
+   end main_body;
+;
+  begin main_body;
+    name           blade2 ;
+    copy_main_body blade1;
+  end main_body;
+;
+  begin main_body;
+    name           blade3 ;
+    copy_main_body blade1 ;
+  end main_body;
+;-------------------------------------------------------------------------------------------------------------------------------                   
+;
+  begin orientation;
+    begin base;
+      body   tower;
+      inipos        0.0 0.0 0.0 ;         initial position of node 1
+      body_eulerang 0.0 0.0 0.0;
+    end base;
+; 
+    begin relative;
+      body1  tower last;
+      body2  towertop 1;
+      body2_eulerang 0.0 0.0 0.0; 
+    end relative;
+;
+    begin relative;
+      body1  towertop last;
+      body2  shaft 1;
+      body2_eulerang 90.0 0.0 0.0; 
+      body2_eulerang 5.0 0.0 0.0;    5 deg tilt angle
+      mbdy2_ini_rotvec_d1 0.0 0.0 -1.0 0.2 ; mbdy2_ini_rotvec_d1 0.0 0.0 -1.0 [init_wr]; 
+    end relative;
+;
+    begin relative;
+      body1  shaft last;         
+      body2  hub1 1;
+      body2_eulerang -90.0 0.0 0.0;    
+      body2_eulerang 0.0 180.0 0.0;    
+      body2_eulerang 2.5 0.0 0.0;      2.5deg cone angle
+    end relative;
+;
+    begin relative;
+      body1  shaft last;         
+      body2  hub2 1;
+      body2_eulerang -90.0 0.0 0.0;    
+      body2_eulerang 0.0 60.0 0.0;   
+      body2_eulerang 2.5 0.0 0.0;      2.5deg cone angle
+    end relative;
+;
+    begin relative;
+      body1  shaft last;         
+      body2  hub3 1;
+      body2_eulerang -90.0 0.0 0.0;    
+      body2_eulerang 0.0 -60.0 0.0;    
+      body2_eulerang 2.5 0.0 0.0;      2.5deg cone angle
+    end relative;
+;
+    begin relative;
+      body1  hub1 last;         
+      body2  blade1 1;
+      body2_eulerang 0.0 0.0 0;    
+    end relative;
+;
+    begin relative;
+      body1  hub2 last;         
+      body2  blade2 1;
+      body2_eulerang 0.0 0.0 0.0;    
+    end relative;
+;
+    begin relative;
+      body1  hub3 last;         
+      body2  blade3 1;
+      body2_eulerang 0.0 0.0 0.0;    
+    end relative;
+;
+ 	end orientation;
+;-------------------------------------------------------------------------------------------------------------------------------
+begin constraint;   
+;
+    begin fix0;  fixed to ground in translation and rotation of node 1
+      body tower;
+    end fix0;
+;
+     ; begin fix1;
+		   ; body1 tower last ;
+		   ; body2 towertop 1;
+		 ; end fix1;
+;
+	begin bearing2;
+		name	yaw_rot;
+		mbdy1	tower last;
+		mbdy2	towertop 1;
+		bearing_vector	2 0 0 -1;
+	end bearing2;
+	begin bearing1;                       free bearing
+		name  shaft_rot;
+		body1 towertop last;
+		body2 shaft 1;
+		bearing_vector 2 0.0 0.0 -1.0;        x=coo (0=global.1=body1.2=body2) vector in body2 coordinates where the free rotation is present
+	end bearing1; 
+;
+     begin fix1;
+		   body1 shaft last ;
+		   body2 hub1 1;
+		 end fix1;
+;
+     begin fix1;
+		   body1 shaft last ;
+		   body2 hub2 1;
+		 end fix1;
+;
+     begin fix1;
+		   body1 shaft last ;
+		   body2 hub3 1;
+		 end fix1; 
+;	
+    begin bearing2;
+      name pitch1;		
+      body1 hub1 last;
+     body2 blade1 1;
+			bearing_vector 2 0.0 0.0 -1.0;
+   end bearing2;
+;
+    begin bearing2;
+      name pitch2;		
+      body1 hub2 last;
+      body2 blade2 1;
+			bearing_vector 2 0.0 0.0 -1.0;
+    end bearing2;
+;
+    begin bearing2;
+      name pitch3;		
+      body1 hub3 last;
+      body2 blade3 1;
+			bearing_vector 2 0.0 0.0 -1.0;
+    end bearing2;
+end constraint;
+;
+end new_htc_structure;
+;---------------------------------------------------------------------------------------------------------------------------------------------------------------- 
+begin wind ;
+  density                 1.225 ;
+  wsp                     6.0   ;
+  tint                    0.1   ;
+  horizontal_input        1     ;
+  windfield_rotations     0.0 0.0 0.0 ;    yaw, tilt, rotation
+  center_pos0             0.0 0.0 -119 ; hub heigth
+  shear_format            1 0.2 ; 
+  turb_format             0     ;  0=none, 1=mann,2=flex
+  tower_shadow_method     0     ;  0=none, 1=potential flow, 2=jet
+  scale_time_start       0.0 ; 
+  ;wind_ramp_factor   0.0 40.0 0.6 1.0 ;
+;  
+begin mann;
+  create_turb_parameters 33.6 .1 3.9 1 0;
+  filename_u /turb/turb_u.bin;
+  filename_v /turb/turb_v.bin;
+  filename_w /turb/turb_w.bin;
+  box_dim_u  4096 2;
+  box_dim_v  32 3;
+  box_dim_w  16 4;
+  dont_scale 1;
+end mann;
+; 
+  begin tower_shadow_potential_2;
+    tower_mbdy_link tower;
+    nsec  2;
+    radius      0.0  4.15 ;
+    radius     115.63 2.75 ;
+  end tower_shadow_potential_2;
+end wind;
+;
+begin aerodrag ;
+  begin aerodrag_element ;
+    mbdy_name tower;
+    aerodrag_sections uniform 10 ;
+    nsec 2 ;
+    sec 0.0 0.6 8.3 ;  tower bottom
+    sec 115.63 0.6 5.5 ;  tower top
+  end aerodrag_element;
+;
+  begin aerodrag_element ;        Nacelle drag side
+    mbdy_name shaft;
+    aerodrag_sections uniform 2 ;
+    nsec 2 ;
+    sec 0.0   0.8 10.0 ;  
+    sec 7.01  0.8 10.0 ;  
+  end aerodrag_element;
+end aerodrag;
+;
+begin aero ;
+  nblades  3;
+  hub_vec shaft -3 ;         rotor rotation vector (normally shaft composant directed from pressure to sustion side)
+  link 1 mbdy_c2_def blade1;
+  link 2 mbdy_c2_def blade2;
+  link 3 mbdy_c2_def blade3;
+  ae_filename        ./data/DTU_10MW_RWT_ae.dat ;
+  pc_filename        ./data/DTU_10MW_RWT_pc.dat ;
+  induction_method   1 ;     0=none, 1=normal
+  aerocalc_method    1 ;     0=ingen aerodynamic, 1=med aerodynamic
+  aerosections       50 ; def. 50
+  ae_sets            1 1 1;
+  tiploss_method     1 ;     0=none, 1=prandtl
+  dynstall_method    2;     0=none, 1=stig øye method,2=mhh method
+;  
+end aero ;
+;-------------------------------------------------------------------------------------------------  
+begin dll;
+;
+  begin type2_dll;
+    name dtu_we_controller ;
+    filename  ./control/dtu_we_controller.dll ;
+    dll_subroutine_init init_regulation_advanced ;
+    dll_subroutine_update update_regulation ;
+    arraysizes_init  100 1 ;
+    arraysizes_update  100 100 ;
+	begin init ;
+       ; Overall parameters
+      constant   1 10000.0    	; Rated power [kW]                         
+      constant   2   0.628    	; Minimum rotor (LSS) speed [rad/s]
+      constant   3   1.005    	; Rated rotor (LSS) speed [rad/s]
+      constant   4  15.6E+06  	; Maximum allowable generator torque [Nm]
+      constant   5  100.0    	; Minimum pitch angle, theta_min [deg], 
+								; if |theta_min|>90, then a table of <wsp,theta_min> is read ;
+								; from a file named 'wptable.n', where n=int(theta_min)
+      constant   6  82.0    	; Maximum pitch angle [deg]
+      constant   7  10.0    	; Maximum pitch velocity operation [deg/s]
+      constant   8   0.4    	; Frequency of generator speed filter [Hz]
+      constant   9   0.7    	; Damping ratio of speed filter [-]
+      constant  10   1.80   	; Frequency of free-free DT torsion mode [Hz], if zero no notch filter used
+      ; Partial load control parameters
+      constant  11   13013100.0 ; Optimal Cp tracking K factor [Nm/(rad/s)^2], ;
+                                ; Qg=K*Omega^2, K=eta*0.5*rho*A*Cp_opt*R^3/lambda_opt^3                     
+      constant  12   0.683456e8 ; Proportional gain of torque controller [Nm/(rad/s)]
+      constant  13   0.153367e8 ; Integral gain of torque controller [Nm/rad]
+      constant  14   0.0    	; Differential gain of torque controller [Nm/(rad/s^2)]
+;     Full load control parameters
+      constant  15   1      	; Generator control switch [1=constant power, 0=constant torque]
+      constant  16   1.06713    ; Proportional gain of pitch controller [rad/(rad/s)]
+      constant  17   0.242445   ; Integral gain of pitch controller [rad/rad]
+      constant  18   0.0    	; Differential gain of pitch controller [rad/(rad/s^2)]
+      constant  19   0.4e-8 	; Proportional power error gain [rad/W]
+      constant  20   0.4e-8 	; Integral power error gain [rad/(Ws)]
+      constant  21 	 11.4       ; Coefficient of linear term in aerodynamic gain scheduling, KK1 [deg]
+      constant  22 	 402.9      ; Coefficient of quadratic term in aerodynamic gain scheduling, KK2 [deg^2] &
+								; (if zero, KK1 = pitch angle at double gain)
+      constant  23   1.3    	; Relative speed for double nonlinear gain [-]
+;     Cut-in simulation parameters
+      constant  24  -1    ; Cut-in time [s], no cut-in is simulated if zero or negative
+      constant  25  1.0   ; Time delay for soft start of torque [1/1P]
+;     Cut-out simulation parameters
+      constant  26  -1    ; Shut-down time [s], no shut-down is simulated if zero or negative
+      constant  27   5.0  ; Time of linear torque cut-out during a generator assisted stop [s]
+      constant  28  1     ; Stop type [1=normal, 2=emergency]
+      constant  29  1.0   ; Time delay for pitch stop after shut-down signal [s]
+      constant  30  3     ; Maximum pitch velocity during initial period of stop [deg/s]
+      constant  31  3.0   ; Time period of initial pitch stop phase [s] (maintains pitch speed specified in constant 30)
+      constant  32  4     ; Maximum pitch velocity during final phase of stop [deg/s]
+;     Expert parameters (keep default values unless otherwise given)
+      constant  33   2.0  	; Time for the maximum torque rate = Maximum allowable generator torque/(constant 33 + 0.01s) [s]
+      constant  34   2.0  	; Upper angle above lowest minimum pitch angle for switch [deg], if equal then hard switch
+      constant  35  95.0  	; Percentage of the rated speed when the torque limits are fully opened [%]
+      constant  36   2.0  	; Time constant of 1st order filter on wind speed used for minimum pitch [1/1P]
+      constant  37   1.0  	; Time constant of 1st order filter on pitch angle used for gain scheduling [1/1P]
+;     Drivetrain damper
+      constant  38   0.0  	; Proportional gain of active DT damper [Nm/(rad/s)], requires frequency in input 10
+;	  Over speed
+	  constant  39  25.0  	; Overspeed percentage before initiating turbine controller alarm (shut-down) [%]
+;     Additional non-linear pitch control term (not used when all zero)
+	  constant  40   0.0  	; Rotor speed error scaling factor [rad/s]
+	  constant  41   0.0  	; Rotor acceleration error scaling factor [rad/s^2]
+	  constant  42   0.0  	; Pitch rate gain [rad/s]
+;     Storm control command
+	  constant 43   28.0  	; Wind speed 'Vstorm' above which derating of rotor speed is used [m/s]
+	  constant 44   28.0  	; Cut-out wind speed (only used for derating of rotor speed in storm) [m/s]	  
+;     Safety system parameters
+	  constant 45   30.0  ; Overspeed percentage before initiating safety system alarm (shut-down) [%]
+	  constant 46    1.5  ; Max low-pass filtered tower top acceleration level [m/s^2]
+;     Turbine parameter
+	  constant 47  178.0  ; Nominal rotor diameter [m]
+;     Parameters for rotor inertia reduction in variable speed region
+      constant 48    0.0  ; Proportional gain on rotor acceleration in variable speed region [Nm/(rad/s^2)] (not used when zero)
+;     Parameters for alternative partial load controller with PI regulated TSR tracking
+      constant 49    7.8  ; Optimal tip speed ratio [-] (only used when K=constant 11 = 0 otherwise  Qg=K*Omega^2 is used)
+;     Parameters for adding aerodynamic drivetrain damping on gain scheduling
+      constant 50    0.0  ; Aerodynamic DT damping coefficient at the operational point of zero pitch angle [Nm/(rad/s)] (not used when zero)
+      constant 51    0.0  ; Coefficient of linear term in aerodynamic DT damping scheduling, KK1 [deg]
+      constant 52    0.0  ; Coefficient of quadratic term in aerodynamic DT damping scheduling, KK2 [deg^2]
+;     Torque exclusion zone
+      constant 53     0.0 ; Exclusion zone: Lower speed limit [rad/s] (Default 0 used if zero)	  
+      constant 54     0.0 ; Exclusion zone: Generator torque at lower limit [Nm] (Default 0 used if zero)	  
+      constant 55     0.0 ; Exclusion zone: Upper speed limit [rad/s] (if =< 0 then exclusion zone functionality is inactive)               	  
+      constant 56     0.0 ; Exclusion zone: Generator torque at upper limit [Nm] (Default 0 used if zero) 	  
+      constant 57     0.0 ; Time constant of reference switching at exclusion zone [s] (Default 0 used if zero)	  
+;     DT torsion mode damper	  
+      constant 58     0.0 ; Frequency of notch filter [Hz] (Default 10 x input 10 used if zero)	  
+      constant 59     0.0 ; Damping of BP filter [-] (Default 0.02 used if zero) 	  
+      constant 60     0.0 ; Damping of notch filter [-] (Default 0.01 used if zero) 	  
+      constant 61     0.0 ; Phase lag of damper [s] =>  max 40*dt (Default 0 used if zero) 	  
+;     Fore-aft Tower mode damper	  
+      constant 62     0.0 ; Frequency of BP filter [Hz] (Default 10 used if zero)\\ 	  
+      constant 63     0.0 ; Frequency of notch fiter [Hz] (Default 10 used if zero)\\ 	  
+      constant 64     0.0 ; Damping of BP filter [-] (Default 0.02 used if zero)\\	  
+      constant 65     0.0 ; Damping of notch filter [-] (Default 0.01 used if zero)\\	  
+      constant 66     0.0 ; Gain of damper [-] (Default 0 used if zero)\\ 	  
+      constant 67     0.0 ; Phase lag of damper [s] =>  max 40*dt (Default 0 used if zero)\\ 	  
+      constant 68     0.0 ; Time constant of 1st order filter on PWR used for fore-aft Tower mode damper GS [Hz] (Default 10 used if zero)	  
+      constant 69     0.0 ; Lower PWR limit used for fore-aft Tower mode damper GS [-] (Default 0 used if zero)	  
+      constant 70     0.0 ; Upper PWR limit used for fore-aft Tower mode damper GS [-] (Default 0 used if zero) 	  
+;     Side-to-side Tower mode filter	  
+      constant 71     0.0 ; Frequency of Tower side-to-sede notch filter [Hz] (Default 100 used if zero)	  
+      constant 72     0.0 ; Damping of notch filter [-] (Default 0.01 used if zero)	  
+      constant 73     0.0 ; Max low-pass filtered tower top acceleration level before initiating safety system alarm (shut-down) [m/s^2] (Default 1.1 x input 46 used if zero)	  
+      constant 74     0.0 ; Time constant of 1st order filter on tower top acceleration [1/1P] (Default 1 used if zero)	  
+;     Pitch deviation monitor parameters	  
+      constant 75 1005020 ; Parameters for pitch deviation monitoring. The format is 1,nnn,mmm 	  
+                          ; where 'nnn' [s] is the period of the moving average and 'mmm' is threshold of the deviation [0.1 deg] (functionality is inactive if value $<$ 1,000,000)	  
+;     Gear ratio	  
+      constant 76     0.0 ; Gear ratio used for the calculation of the LSS rotational speeds and the HSS generator torque reference [-] (Default 1 if zero)	  
+	end init ;
+;
+    begin output ;
+      general time ; [s]     
+      constraint bearing1 shaft_rot 1 only 2 ; Drivetrain speed [rad/s]
+      constraint bearing2 pitch1 1 only 1; [rad]         
+      constraint bearing2 pitch2 1 only 1; [rad]                               
+      constraint bearing2 pitch3 1 only 1; [rad]                               
+      wind free_wind 1 0.0 0.0 -119      ; Global coordinates at hub height
+      dll inpvec 2 2                     ; Elec. power from generator servo .dll
+      dll inpvec 2 8                     ; Grid state flag from generator servo .dll
+      mbdy state acc towertop   1 1.0 global only 1 ; Tower top x-acceleration [m/s^2]
+      mbdy state acc towertop   1 1.0 global only 2 ; Tower top y-acceleration [m/s^2]
+    end output;    
+  end type2_dll;
+;
+   begin type2_dll;
+     name generator_servo ;
+     filename  ./control/generator_servo.dll ;
+     dll_subroutine_init init_generator_servo ;
+     dll_subroutine_update update_generator_servo ;
+     arraysizes_init  100 1 ;
+     arraysizes_update 100 100 ;
+ 	begin init ;
+       constant 1  20.0    ; Frequency of 2nd order servo model of generator-converter system [Hz]   
+       constant 2  0.9     ; Damping ratio 2nd order servo model of generator-converter system [-]
+       constant 3 15.6E+06 ; Maximum allowable LSS torque (pull-out torque) [Nm]
+       constant 4 0.94     ; Generator efficiency [-]
+       constant 5 1.0      ; Gearratio [-]
+       constant 6 0.0      ; Time for half value in softstart of torque [s] 
+       constant 7 -1       ; Time for grid loss [s] (never if lower than zero)
+     end init ;
+;
+     begin output;
+       general time                          ;   Time [s]    
+       dll inpvec 1 1                        ;   Electrical torque reference [Nm]  
+       constraint bearing1 shaft_rot 1 only 2;   Generator LSS speed [rad/s]   
+       mbdy momentvec shaft 1 1 shaft only 3 ;   Shaft moment [kNm] (Qshaft)
+     end output;
+;
+     begin actions;    
+        mbdy moment_int shaft 1 -3 shaft towertop 2 ;   Generator LSS torque [Nm]
+     end actions;
+   end type2_dll;
+;
+   begin type2_dll;
+     name mech_brake ;
+     filename  ./control/mech_brake.dll ;
+     dll_subroutine_init init_mech_brake ;
+     dll_subroutine_update update_mech_brake ;
+     arraysizes_init    100 1 ;
+     arraysizes_update  100 100 ;
+ 	begin init ;
+      constant 1  9.36E+06   ; Fully deployed maximum brake torque [Nm] (0.6*max torque)
+      constant 2     100.0   ; Parameter alpha used in Q = tanh(omega*alpha), typically 1e2/Omega_nom
+      constant 3       0.5   ; Delay time for before brake starts to deploy [s]
+      constant 4       0.6   ; Time for brake to become fully deployed [s]
+     end init ;
+;
+     begin output;
+	   general time                           ; Time [s]
+	   constraint bearing1 shaft_rot 1 only 2 ; Generator LSS speed [rad/s]
+	   dll inpvec 1 25                        ; Command to deploy mechanical disc brake [0,1]
+     end output;
+;
+     begin actions;    
+        mbdy moment_int shaft 1 -3 shaft towertop 2 ;   Brake LSS torque [Nm]
+     end actions;
+   end type2_dll;
+;
+  begin type2_dll;
+    name servo_with_limits ;
+    filename  ./control/servo_with_limits.dll ;
+    dll_subroutine_init init_servo_with_limits ;
+    dll_subroutine_update update_servo_with_limits ;
+    arraysizes_init  100 1 ;
+    arraysizes_update  100 100 ;
+	begin init ;
+      constant 1   3    ; Number of blades [-]
+      constant 2   1.0  ; Frequency of 2nd order servo model of pitch system [Hz]
+      constant 3   0.7  ; Damping ratio 2nd order servo model of pitch system [-]
+      constant 4  10.0  ; Max. pitch speed [deg/s]
+      constant 5  15.0  ; Max. pitch acceleration [deg/s^2]
+      constant 6  -5.0  ; Min. pitch angle [deg] 
+      constant  7 90.0  ; Max. pitch angle [deg] 	  
+	  constant  8 -1    ; Time for pitch runaway [s]
+	  constant  9 -1    ; Time for stuck blade 1 [s]
+	  constant 10 0.0   ; Angle of stuck blade 1 [deg] (if > 90 deg then blade is stuck at instantaneous angle)
+	end init ;
+    begin output;
+      general time        ;  Time                         [s]     
+       dll inpvec 1 2     ;  Pitch1 demand angle          [rad]
+       dll inpvec 1 3     ;  Pitch2 demand angle          [rad]
+       dll inpvec 1 4     ;  Pitch3 demand angle          [rad]
+       dll inpvec 1 26    ;  Flag for emergency pitch stop         [0=off/1=on]
+    end output;           
+;
+    begin actions;    
+      constraint bearing2 angle pitch1 ; Angle pitch1 bearing    [rad]
+      constraint bearing2 angle pitch2 ; Angle pitch2 bearing    [rad]
+      constraint bearing2 angle pitch3 ; Angle pitch3 bearing    [rad]
+    end actions;                      
+  end type2_dll;
+  ;
+  begin type2_dll;
+      name	yaw_servo;
+      deltat	0.01;
+      filename	./control/servo_with_limits2.dll;
+      dll_subroutine_init	init_servo_with_limits;
+      dll_subroutine_update	update_servo_with_limits;
+      arraysizes_init	100 1;
+      arraysizes_update	100 100;
+      begin init;
+        constant	1 1;	Number of bearings [-]
+        constant	2 0.5;	Frequency of 2nd order servo model of yaw system [Hz]
+        constant	3 0.7;	Damping ratio 2nd order servo model of yaw system [-]
+        constant	4 1;	Max. yaw speed [deg/s]
+        constant	5 1;	Max. yaw acceleration [deg/s^2]
+        constant	6 -100000000000000000000;	Min. yaw angle [deg]
+        constant	7 100000000000000000000;	Max. yaw angle [deg]
+        constant	8 100000000000000000000;	Time for yaw runaway [s]
+        constant	9 -1;	Time for stuck yaw [s]
+        constant	10 0;	Angle of stuck yaw [deg]
+      end init;
+      begin output;
+        general time;	Time [s]
+        general variable 1 0 yaw_setpoint rad;   Demanded nacelle yaw angle [rad]
+      end output;
+      begin actions;
+        constraint	bearing2 angle yaw_rot;	Yaw bearing angle [rad];
+      end actions;
+    end type2_dll;
+;
+;	--- DLL for tower-blade tip distance -- ;
+;-------------------------------------------
+begin type2_dll;
+  name towerclearance_mblade ;
+  filename  ./control/towerclearance_mblade.dll ;
+  dll_subroutine_init initialize ;
+  dll_subroutine_update update ;
+  arraysizes_init  3 1 ;
+  arraysizes_update  15 6 ;
+  begin init ;    Variables passed into initialization function
+    constant  1 4.15  ; Tower radius at tower bottom [m] 
+    constant  2 2.75  ; Tower radius at tower top [m]
+    constant  3    3  ; Number of points to check [-]
+  end init ;
+  begin output;   Variables passed into update function
+    mbdy state pos tower    1 0.0 global ; [1,2,3] global coordinates of tower base
+    mbdy state pos tower   10 1.0 global ; [4,5,6] global coordinates of tower top
+    mbdy state pos blade1  26 1.0 global  ; [7,8,9] global coordinates of point 1 (blade 1 tip)
+    mbdy state pos blade2  26 1.0 global  ; [10,11,12] global coordinates of point 2 (blade 2 tip)
+	mbdy state pos blade3  26 1.0 global  ; [13,14,15] global coordinates of point 3 (blade 3 tip)
+  end output;           
+end type2_dll;
+;
+end dll;
+;---------------------------------------------------------------------------------------------------------------------------------------------------------------- 
+;
+begin output;
+  filename ./res/DTU_10MW_RWT_ver09 ;
+  data_format  gtsdf;
+  buffer 250 ;
+;
+  general time;  
+  constraint bearing1 shaft_rot 2; angle and angle velocity 
+  constraint bearing2 pitch1 5;    angle and angle velocity 
+  constraint bearing2 pitch2 5;    angle and angle velocity 
+  constraint bearing2 pitch3 5;    angle and angle velocity 
+  aero omega ;
+  aero torque;
+  aero power;
+  aero thrust;
+  wind free_wind 1 0.0 0.0 -119; local wind at fixed position: coo (1=global,2=non-rotation rotor coo.), pos x, pos y, pos z
+  mbdy state pos shaft 4 1 global # rotor position;
+  ; Moments:
+  mbdy momentvec tower  1 1  tower # tower base ;
+  mbdy momentvec tower  10 2 tower # tower yaw bearing ;
+  mbdy momentvec shaft  4 1  shaft # main bearing ;
+  mbdy momentvec blade1 2  2 blade1 # blade 1 root ;
+  mbdy momentvec blade2 2  2 blade2 # blade 2 root ;
+  mbdy momentvec blade3 2  2 blade3 # blade 3 root ;
+  mbdy momentvec blade1 13 1 local # blade 1 50% local e coo ;
+  mbdy momentvec blade2 13 1 local # blade 2 50% local e coo ;
+  mbdy momentvec blade3 13 1 local # blade 3 50% local e coo ;
+  ; Displacements and accellerations
+  mbdy state pos tower 10 1.0 global only 1 # Tower top FA displ; 
+  mbdy state pos tower 10 1.0 global only 2 # Tower top SS displ; 
+  mbdy state acc tower 10 1.0 global only 1 # Tower top FA acc; 
+  mbdy state acc tower 10 1.0 global only 2 # Tower top SS acc; 
+;  
+  mbdy state pos blade1  26 1.0 blade1 # blade 1 tip pos ;
+  mbdy state pos blade2  26 1.0 blade2 # blade 2 tip pos ;
+  mbdy state pos blade3  26 1.0 blade3 # blade 3 tip pos ;
+  mbdy state pos blade1  26 1.0 global # gl blade 1 tip pos ;  
+; - Monitor Aerodynamics - ;
+  aero windspeed 3 1 1 72.5;
+  aero alfa 1 72.5; 
+  aero alfa 2 72.5; 
+  aero alfa 3 72.5; 
+  aero cl 1 72.5;
+  aero cl 2 72.5;
+  aero cl 3 72.5;
+  aero cd 1 72.5;
+  aero cd 2 72.5;
+  aero cd 3 72.5;    
+; DLL outputs and into HAWC2 
+  dll inpvec 1  1 # Generator torque reference            [Nm]   ;
+  dll inpvec 1  2 # Pitch angle reference of blade 1      [rad]  ;
+  dll inpvec 1  3 # Pitch angle reference of blade 2      [rad]  ;
+  dll inpvec 1  4 # Pitch angle reference of blade 3      [rad]  ;
+  dll inpvec 1  5 # Power reference                       [W]    ;
+  dll inpvec 1  6 # Filtered wind speed                   [m/s]  ;
+  dll inpvec 1  7 # Filtered rotor speed                  [rad/s];
+  dll inpvec 1  8 # Filtered rotor speed error for torque [rad/s];
+  dll inpvec 1  9 # Bandpass filtered rotor speed         [rad/s];
+  dll inpvec 1 10 # Proportional term of torque contr.    [Nm]   ;
+  dll inpvec 1 11 # Integral term of torque controller    [Nm]   ;
+  dll inpvec 1 12 # Minimum limit of torque               [Nm]   ;
+  dll inpvec 1 13 # Maximum limit of torque               [Nm]   ;
+  dll inpvec 1 14 # Torque limit switch based on pitch    [-]    ;
+  dll inpvec 1 15 # Filtered rotor speed error for pitch  [rad/s];
+  dll inpvec 1 16 # Power error for pitch                 [W]    ;
+  dll inpvec 1 17 # Proportional term of pitch controller [rad]  ;
+  dll inpvec 1 18 # Integral term of pitch controller     [rad]  ;
+  dll inpvec 1 19 # Minimum limit of pitch                [rad]  ;
+  dll inpvec 1 20 # Maximum limit of pitch                [rad]  ;
+  dll inpvec 1 21 # Torque reference from DT dammper      [Nm]  ;
+  dll inpvec 1 22 # Status signal                         [-]  ;
+  dll inpvec 1 23 # Total added pitch rate                [rad/s]  ;  
+  dll inpvec 1 24 # Filtered Mean pitch for gain sch      [rad]  ;  
+  dll inpvec 1 25 # Flag for mechnical brake              [0=off/1=on] ;
+  dll inpvec 1 26 # Flag for emergency pitch stop         [0=off/1=on] ;
+  dll inpvec 1 27 # LP filtered acceleration level        [m/s^2] ;
+  dll inpvec 1 31 #  Monitored average of reference pitch     [rad] ;
+  dll inpvec 1 32 #  Monitored ave. of actual pitch (blade 1) [rad] ;
+; Input from generator model
+   dll inpvec 2 1  # Mgen LSS [Nm] ;
+   dll inpvec 2 2  # Pelec    [W]  ;
+   dll inpvec 2 3  # Mframe   [Nm] ;
+   dll inpvec 2 4  # Mgen HSS [Nm] ;
+   dll inpvec 2 8  # Grid flag [0=run/1=stop];
+; Input from mechanical brake
+   dll inpvec 3 1 # Brake torque [Nm] ;
+; Input from pitch servo
+   dll inpvec 4 1 # pitch 1 [rad];
+   dll inpvec 4 2 # pitch 2 [rad];
+   dll inpvec 4 3 # pitch 3 [rad];
+; yaw set point
+   general variable 1 0 yaw_setpoint rad;
+   dll inpvec 5 1 # yaw angle [rad];
+; Check tower clearence
+   dll inpvec 6 1 # Bltip tow min d [m];
+end output;
+;
+exit;

h2lib_tests/test_files/minimal/htc/minimal_nooutput.htc

@@ -3,6 +3,7 @@ begin simulation;
   solvertype   1 ;    (newmark)
   on_no_convergence continue ;
   convergence_limits 1E3 1.0 1E-7 ;
+  logfile log.log;
   log_deltat 1;
   begin newmark;
     deltat    0.1;  

h2lib_tests/test_static_solver_wrapper.py

@@ -0,0 +1,20 @@
+import pytest
+
+from h2lib._h2lib import H2Lib
+from h2lib_tests.test_files import tfp
+
+from .test_write_htc import write_dtu10mw_only_tower
+
+
+def test_init_static_solver(write_dtu10mw_only_tower):
+    with H2Lib() as h2:
+        # Load a basic model.
+        model_path = tfp + "DTU_10_MW/"
+        h2.init("htc/DTU_10MW_RWT_only_tower.htc", model_path)
+
+        # First execution is fine.
+        h2.init_static_solver()
+
+        # The next should give an error.
+        with pytest.raises(RuntimeError, match="STATIC_SOLVER_ALREADY_INITIALIZED"):
+            h2.init_static_solver()

h2lib_tests/test_topology_h2lib.py

@@ -0,0 +1,143 @@
+from h2lib._h2lib import H2Lib
+
+from numpy import testing as npt
+from h2lib_tests.test_files import tfp
+import numpy as np
+import pytest
+from .test_write_htc import (
+    write_dtu10mw_only_tower,
+    write_dtu10mw_only_tower_rotated,
+    write_dtu10mw_only_tower_encrypted
+)
+
+
+@pytest.fixture(scope="module")
+def h2_dtu_10mw_only_tower(write_dtu10mw_only_tower):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_tower.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="module")
+def h2_dtu_10mw_only_tower_rotated(write_dtu10mw_only_tower_rotated):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_tower_rotated.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="module")
+def h2_dtu_10mw_only_tower_encrypted(write_dtu10mw_only_tower_encrypted):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_tower_encrypted.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    yield h2
+    h2.close()
+
+
+def test_number_of_bodies_and_constraints(
+    h2_dtu_10mw_only_tower,
+):
+    nbdy, ncst = h2_dtu_10mw_only_tower.get_number_of_bodies_and_constraints()
+    assert nbdy == 3
+    assert ncst == 9
+
+
+def test_number_of_bodies_and_constraints_encrypted(
+    h2_dtu_10mw_only_tower_encrypted,
+):
+    with pytest.raises(RuntimeError, match="STRUCTURE_IS_CONFIDENTIAL"):
+        h2_dtu_10mw_only_tower_encrypted.get_number_of_bodies_and_constraints()
+
+
+def test_get_number_of_elements(h2_dtu_10mw_only_tower):
+    nelem = h2_dtu_10mw_only_tower.get_number_of_elements()
+    npt.assert_array_equal(nelem, np.array([3, 3, 4]))
+
+
+def test_get_number_of_elements_encrypted(
+    h2_dtu_10mw_only_tower_encrypted,
+):
+    # This test is not really needed, since the check for confidential structure
+    # is already done by test_number_of_bodies_and_constraints_encrypted().
+    with pytest.raises(RuntimeError, match="STRUCTURE_IS_CONFIDENTIAL"):
+        h2_dtu_10mw_only_tower_encrypted.get_number_of_elements()
+
+
+def test_get_timoshenko_location(
+    h2_dtu_10mw_only_tower,
+):
+    # Test first element.
+    l, r1, r12, tes = h2_dtu_10mw_only_tower.get_timoshenko_location(ibdy=0, ielem=0)
+    assert l - 11.5 < 1e-14
+    npt.assert_array_equal(r1, np.array([0.0, 0.0, 0]))
+    npt.assert_array_almost_equal_nulp(r12, np.array([0.0, 0.0, -11.5]))
+    npt.assert_array_equal(
+        tes,
+        np.array([[-1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, -1.0]]),
+    )
+
+    # Test last element.
+    l, r1, r12, tes = h2_dtu_10mw_only_tower.get_timoshenko_location(ibdy=2, ielem=3)
+    assert l - 12.13 < 1e-14
+    npt.assert_array_almost_equal_nulp(r1, np.array([0.0, 0.0, -34.5]))
+    npt.assert_array_almost_equal_nulp(r12, np.array([0.0, 0.0, -12.13]), nulp=3)
+    npt.assert_array_equal(
+        tes,
+        np.array([[-1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, -1.0]]),
+    )
+
+
+def test_get_timoshenko_location_body_does_not_exist(
+    h2_dtu_10mw_only_tower,
+):
+    with pytest.raises(IndexError, match="BODY_DOES_NOT_EXIST"):
+        h2_dtu_10mw_only_tower.get_timoshenko_location(ibdy=1000, ielem=0)
+
+
+def test_get_timoshenko_location_element_does_not_exist(
+    h2_dtu_10mw_only_tower,
+):
+    with pytest.raises(IndexError, match="ELEMENT_DOES_NOT_EXIST"):
+        h2_dtu_10mw_only_tower.get_timoshenko_location(ibdy=0, ielem=1000)
+
+
+def test_get_timoshenko_location_encrypted(
+    h2_dtu_10mw_only_tower_encrypted,
+):
+    with pytest.raises(RuntimeError, match="STRUCTURE_IS_CONFIDENTIAL"):
+        h2_dtu_10mw_only_tower_encrypted.get_timoshenko_location(ibdy=0, ielem=0)
+
+
+def test_get_body_rotation_tensor_1(h2_dtu_10mw_only_tower):
+    amat = h2_dtu_10mw_only_tower.get_body_rotation_tensor(ibdy=0)
+    npt.assert_array_equal(amat, np.eye(3))
+
+
+def test_get_body_rotation_tensor_2(
+    h2_dtu_10mw_only_tower_rotated, write_dtu10mw_only_tower_rotated
+):
+    amat = h2_dtu_10mw_only_tower_rotated.get_body_rotation_tensor(ibdy=0)
+    _, alpha = write_dtu10mw_only_tower_rotated
+    alpha_rad = np.deg2rad(alpha)
+    sa = np.sin(alpha_rad)
+    ca = np.cos(alpha_rad)
+    npt.assert_array_almost_equal_nulp(
+        amat, np.array([[1.0, 0.0, 0.0], [0.0, ca, -sa], [0.0, sa, ca]])
+    )
+
+
+def test_get_body_rotation_tensor_body_does_not_exist(h2_dtu_10mw_only_tower):
+    with pytest.raises(IndexError, match="BODY_DOES_NOT_EXIST"):
+        h2_dtu_10mw_only_tower.get_body_rotation_tensor(ibdy=1000)
+
+
+def test_get_body_rotation_tensor_encrypted(h2_dtu_10mw_only_tower_encrypted):
+    with pytest.raises(RuntimeError, match="STRUCTURE_IS_CONFIDENTIAL"):
+        h2_dtu_10mw_only_tower_encrypted.get_body_rotation_tensor(ibdy=0)

h2lib_tests/test_write_htc.py

@@ -0,0 +1,59 @@
+import pytest
+from wetb.hawc2.htc_file import HTCFile
+from h2lib_tests.test_files import tfp
+
+
+@pytest.fixture(scope="module")
+def write_dtu10mw_only_tower():
+    # Start from DTU_10MW_RWT and delete everything except the tower.
+    htc = HTCFile(tfp + "DTU_10_MW/htc/DTU_10MW_RWT.htc")
+    htc.set_name("DTU_10MW_RWT_only_tower")
+    for key1 in htc["new_htc_structure"].keys():
+        if key1.startswith("main_body"):
+            if "tower" not in htc["new_htc_structure"][key1]["name"].values:
+                htc["new_htc_structure"][key1].delete()
+        if key1 == "orientation":
+            for key2 in htc["new_htc_structure"]["orientation"].keys():
+                if key2.startswith("relative"):
+                    htc["new_htc_structure"]["orientation"][key2].delete()
+        if key1 == "constraint":
+            for key2 in htc["new_htc_structure"]["constraint"].keys():
+                if key2 != "fix0":
+                    htc["new_htc_structure"]["constraint"][key2].delete()
+    htc["wind"].delete()
+    htc["aerodrag"].delete()
+    htc["aero"].delete()
+    htc["dll"].delete()
+    htc["output"].delete()
+    # Reduce simulation time.
+    htc.simulation.time_stop = 10.0
+    # Change number of bodies in the tower.
+    htc.new_htc_structure.main_body.nbodies = 3
+    # Save the new file.
+    htc.save()
+    return htc
+
+
+@pytest.fixture(scope="module")
+def write_dtu10mw_only_tower_rotated(write_dtu10mw_only_tower):
+    # Start from the DTU_10MW_RWT_only_tower and rotate the tower.
+    htc = write_dtu10mw_only_tower.copy()
+    htc.set_name("DTU_10MW_RWT_only_tower_rotated")
+    alpha = 30.0
+    htc.new_htc_structure.orientation.base.body_eulerang = [
+        alpha,
+        0.0,
+        0.0,
+    ]
+    htc.save()
+    return (htc, alpha)
+
+
+@pytest.fixture(scope="module")
+def write_dtu10mw_only_tower_encrypted(write_dtu10mw_only_tower):
+    # Start from the DTU_10MW_RWT_only_tower and then encrypt the tower.
+    htc = write_dtu10mw_only_tower.copy()
+    htc.set_name("DTU_10MW_RWT_only_tower_encrypted")
+    # Only the tower is left.
+    htc.new_htc_structure.main_body.timoschenko_input.filename = "./data/DTU_10MW_RWT_Tower_st.dat.v3.enc"
+    htc.save()

{h2lib_tests-13.1.304.dist-info → h2lib_tests-13.1.502.dist-info}/METADATA

@@ -1,7 +1,7 @@
 Metadata-Version: 2.1
 Name: h2lib-tests
-Version: 13.1.304
-Summary: Tests and test_files for test h2lib (13.1.3)
+Version: 13.1.502
+Summary: Tests and test_files for test h2lib (13.1.5+4-g20eca1a)
 Download-URL: 
 Author: Mads M Pedersen
 Author-email: 

{h2lib_tests-13.1.304.dist-info → h2lib_tests-13.1.502.dist-info}/RECORD

@@ -6,6 +6,9 @@ h2lib_tests/test_h2lib.py,sha256=ZKp8pfQQpCUQ7IBXypTmvewbJwTQowSiCEt6t4R9TgY,133
 h2lib_tests/test_h2rotor.py,sha256=kaE9AFDYSPaMFSgBAbeegd7-l6Z_4_BsQaWie0k32q4,6210
 h2lib_tests/test_mpi.py,sha256=emBgXRAvvFFOsVrAziiQCUZvEF9HS5Wc-x4KqqltQP0,6835
 h2lib_tests/test_multiprocessinterface.py,sha256=h2o4havtK6IuMXsplNjGUa3VxOnbpEYGxdrrAKQilj0,1470
+h2lib_tests/test_static_solver_wrapper.py,sha256=xCWDkxtW8K2kD6l6MRgyPS1kz91IpbRA0Cc3BG7o3Es,604
+h2lib_tests/test_topology_h2lib.py,sha256=cClNtAIc1wR_gZa7eDO4ZySvhScqCEIHVMKZEhk29QE,4980
+h2lib_tests/test_write_htc.py,sha256=t78STsVoQ2qR5WB-jkYhaMOy0jWHHsDfuO2nfyPj4ag,2242
 h2lib_tests/test_files/__init__.py,sha256=9e6ZUPb42e0wf2E1rutdcTM8hROcWFRVPXtZriU3ySw,50
 h2lib_tests/test_files/my_test_cls.py,sha256=7ZDsFkxrLfOY6q00U5Y-daxfuhATK-K5H04RP-VmQdE,850
 h2lib_tests/test_files/DTU_10_MW/control/dtu_we_controller.dll,sha256=C5T_CuAFtIuDgCXSYAoNu24yKPwj2nWOeORacJbLN9s,1134592
@@ -19,6 +22,9 @@ h2lib_tests/test_files/DTU_10_MW/control/mech_brake.so,sha256=YG-cTEp83x2J2jpRMg
 h2lib_tests/test_files/DTU_10_MW/control/mech_brake_64.dll,sha256=kG5SEM0A6M1eCSjSZb2LBIbI-Gi4gWfqLGZZ2PxUfjY,893952
 h2lib_tests/test_files/DTU_10_MW/control/servo_with_limits.dll,sha256=RsYVurATKF1-1mVKiMDoORSUGEU3Th8qBAl-a0BNrZg,839680
 h2lib_tests/test_files/DTU_10_MW/control/servo_with_limits.so,sha256=6nERT3HsDZPrxfisn8T6dOFGf0U4PF4t8pJlrOPg8Hk,174824
+h2lib_tests/test_files/DTU_10_MW/control/servo_with_limits2.dll,sha256=RsYVurATKF1-1mVKiMDoORSUGEU3Th8qBAl-a0BNrZg,839680
+h2lib_tests/test_files/DTU_10_MW/control/servo_with_limits2.so,sha256=6nERT3HsDZPrxfisn8T6dOFGf0U4PF4t8pJlrOPg8Hk,174824
+h2lib_tests/test_files/DTU_10_MW/control/servo_with_limits2_64.dll,sha256=tz0G4nfYQjXdskVU0Cv2nMN2eNPtZlBglgBZr2xrRcU,1055232
 h2lib_tests/test_files/DTU_10_MW/control/servo_with_limits_64.dll,sha256=tz0G4nfYQjXdskVU0Cv2nMN2eNPtZlBglgBZr2xrRcU,1055232
 h2lib_tests/test_files/DTU_10_MW/control/towerclearance_mblade.dll,sha256=T7LoFiQ9AYIaO-P9OMmM5NC1XmHEWj1lThDLlJHO1ys,141312
 h2lib_tests/test_files/DTU_10_MW/control/towerclearance_mblade.so,sha256=bfgcVPyw7YegMW2g5nOGpseNjspaZHk8qEKNeJ-tPL8,138032
@@ -28,6 +34,7 @@ h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Blade_st.dat,sha256=7stkwSf7R
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Hub_st.dat,sha256=UPIVkaThyxRsEeMIyjhXIr_n9kr1SAGQYudgj_ofJ04,916
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Shaft_st.dat,sha256=Jz_QA1hJUY6ZG4re1oFj0_7MkXu8az9jn2-WLGb0jY0,905
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Tower_st.dat,sha256=0asFvlyhL1CZltWve-PNGgjfpfdMGRCCwMf3U0K_x_c,8367
+h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Tower_st.dat.v3.enc,sha256=2MIV4QV0Y5XZ6hHsUdO95pG4b1HVlNWeuIzSM6n-88A,8422
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_Towertop_st.dat,sha256=iW8RzjgQe7D3y9DEQMD74e0PWqJ4d7meYGTT6Fsl4D0,955
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_ae.dat,sha256=LkAlElRW_Fl1V4GD8hQZIgLD4IW0o3uxXL2olnEU3dI,1531
 h2lib_tests/test_files/DTU_10_MW/data/DTU_10MW_RWT_pc.dat,sha256=2Nc_ofcqfCB3iu8dSFRj8Ub6AhaQADdW0O8BXBWj5Uk,26888
@@ -35,6 +42,7 @@ h2lib_tests/test_files/DTU_10_MW/data/operation.dat,sha256=QH3ULApoKiZJfT4Z65aGN
 h2lib_tests/test_files/DTU_10_MW/htc/DTU_10MW_RWT.htc,sha256=80y88seN4BL_gAiM-RYaTlZxKCYJMr27IBdCX38O6DE,32313
 h2lib_tests/test_files/DTU_10_MW/htc/DTU_10MW_RWT_mann_turb.htc,sha256=-jFK-ra64GeY3fQrP8FxeZ25DkQMdS8yQ9sO78SS9G8,32006
 h2lib_tests/test_files/DTU_10_MW/htc/DTU_10MW_RWT_no_output.htc,sha256=PxkaacH8iVp0ewk6x8fIjkITPUkiVar5bxZCWqZSoP4,25887
+h2lib_tests/test_files/DTU_10_MW/htc/DTU_10MW_RWT_yaw_acturator.htc,sha256=FET165Hd3w6r0hwOeMU1WKARiBToIWL_tXPjYsjIatQ,33776
 h2lib_tests/test_files/IEA-15-240-RWT/IEA_15MW_RWT_Blade_st_FPM.st,sha256=ytwFSx1cFSF25vZ8WZpYFiV2aL3DTvetPrR7oEqFYpc,16605
 h2lib_tests/test_files/IEA-15-240-RWT/IEA_15MW_RWT_Blade_st_noFPM.st,sha256=Pdb8rfJoV7zOAmbPn_xUyGgQk-xIjQ5fCZPRzsu9Q-M,20945
 h2lib_tests/test_files/IEA-15-240-RWT/IEA_15MW_RWT_Dummy_st.dat,sha256=fJGn9bFdMwCV1ceoxT2SUbuDMuNJGCaeqHWUqFlyfKE,532
@@ -74,12 +82,12 @@ h2lib_tests/test_files/minimal/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5
 h2lib_tests/test_files/minimal/htc/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 h2lib_tests/test_files/minimal/htc/minimal.htc,sha256=qdtxlFNczVxmH3h7RDefUHUSuDESeWgj1uSf_rQjjT4,985
 h2lib_tests/test_files/minimal/htc/minimal_mann_turb.htc,sha256=hT5e9AQ7orN1dBwdNlJIwHaj9wuj_DRALTOb4EfhQw4,1432
-h2lib_tests/test_files/minimal/htc/minimal_nooutput.htc,sha256=qdtxlFNczVxmH3h7RDefUHUSuDESeWgj1uSf_rQjjT4,985
+h2lib_tests/test_files/minimal/htc/minimal_nooutput.htc,sha256=z5KPGWjjryBA2uWe0WQ17W6g6Q21tWi-QjEAQtERhPo,1005
 h2lib_tests/test_files/minimal/res/minimal_mann_turb.hdf5,sha256=Q3cs3bZyplZjBponctO2EARIaP6TpSWP7KI8SfXWRAY,41400
 h2lib_tests/test_files/minimal/turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_u,sha256=byiorJmXDL6uKFbyfXthHTjJdm6ELvLR2lS202KrhRI,1048576
 h2lib_tests/test_files/minimal/turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_v,sha256=cxK5Rfgfm3gyJsEYi_KlmYY8DIIl_G0aizN2jt18Glc,1048576
 h2lib_tests/test_files/minimal/turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_w,sha256=xs61jAwhP3fIR1P5Oa8ovEt2baLoF8uCNs6pKIT8L4o,1048576
-h2lib_tests-13.1.304.dist-info/METADATA,sha256=On7mnogtZ5MNXvUszTaV95PSOZdV1i7YS0RPbHzSRAs,349
-h2lib_tests-13.1.304.dist-info/WHEEL,sha256=R0nc6qTxuoLk7ShA2_Y-UWkN8ZdfDBG2B6Eqpz2WXbs,91
-h2lib_tests-13.1.304.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
-h2lib_tests-13.1.304.dist-info/RECORD,,
+h2lib_tests-13.1.502.dist-info/METADATA,sha256=VlhajO_jSvh6Ff37KcgGpcJmCCRapAnmN8NGeCq17RY,360
+h2lib_tests-13.1.502.dist-info/WHEEL,sha256=GV9aMThwP_4oNCtvEC2ec3qUYutgWeAzklro_0m4WJQ,91
+h2lib_tests-13.1.502.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
+h2lib_tests-13.1.502.dist-info/RECORD,,

{h2lib_tests-13.1.304.dist-info → h2lib_tests-13.1.502.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (72.1.0)
+Generator: setuptools (75.1.0)
 Root-Is-Purelib: true
 Tag: py3-none-any