turning off altitude

Author: SteveDoyle2

pyNastran/f06/dev/flutter/utils.py

@@ -15,7 +15,8 @@
 def load_f06_op2(f06_filename: str, log: SimpleLogger,
                  in_units: str,
                  out_units: str,
-                 use_rhoref: bool) -> tuple[OP2, dict[int, FlutterResponse]]:
+                 use_rhoref: bool,
+                 make_alt: bool=False) -> tuple[OP2, dict[int, FlutterResponse]]:
     """
     load a Vg-Vf plot from:
      - OP2 / F06
@@ -43,6 +44,7 @@ def load_f06_op2(f06_filename: str, log: SimpleLogger,
                 f06_units=in_units_dict,
                 out_units=out_units_dict,
                 use_rhoref=use_rhoref,
+                make_alt=make_alt,
                 log=log)
         except Exception as e:
             log.error(str(e))

pyNastran/f06/flutter_response.py

@@ -575,10 +575,6 @@ def _set_pknl_results(self,
 
         vel_units_in = in_units['velocity']
         q_units_in = in_units['dynamic_pressure']
-        if _is_q_units_consistent(density_units_in, vel_units_in, q_units_in):
-            q = 0.5 * rho * vel**2
-        else:  # pragma: no cover
-            raise NotImplementedError((density_units_in, vel_units_in, q_units_in))
 
         # eas  = (2 * q / rho_ref)**0.5
         # eas = V * sqrt(rho / rhoSL)
@@ -590,27 +586,35 @@ def _set_pknl_results(self,
         altitude_units = in_units['altitude']
 
         # print('density_units_in=%r density_units2=%r' % (density_units_in, density_units2))
-        kdensityi = convert_density(1., density_units_in, 'slug/ft^3')
-
         resultsi = results[:, :, :9]
         assert resultsi.shape[2] == 9, resultsi.shape
 
-        rho_in_slug_ft3 = rho * kdensityi
+        rho_in_slug_ft3 = rho_ref * rho
         ft_to_alt_unit = convert_altitude(1., 'ft', altitude_units)
         if self.make_alt:
+            nrho = len(rho)
             alt_ft = []
             for idensity, densityi in enumerate(rho_in_slug_ft3.ravel()):
                 try:
                     alt_fti = get_alt_for_density(densityi, density_units='slug/ft^3',
-                                                  alt_units='ft', nmax=20)
+                                                  alt_units='ft', nmax=50)
                 except Exception:
-                    raise RuntimeError(f'failed to find altitude for density[{idensity}]='
-                                       f'{rho.ravel()[idensity]:g}; density_units_in={density_units_in!r}')
+                    raise RuntimeError(f'Case {idensity+1}/{nrho}: failed to find altitude for density='
+                                       f'{rho.ravel()[idensity]:g}; density_units_in={density_units_in!r} ({densityi:g} slug/ft^3)\n'
+                                       '  output_rho_range = rho_ref * input_rho_range\n'
+                                       f'  rho_ref:             {rho_ref:g} {density_units_in}\n'
+                                       f'  input_rho_range:     [{rho.min():g}, {rho.max()}]\n'
+                                       f'  => output_rho_range: [{rho_in_slug_ft3.min():g}, {rho_in_slug_ft3.max()}] slug/ft^3')
                 alt_ft.append(alt_fti)
             alt = np.array(alt_ft, dtype=rho.dtype).reshape(vel.shape) * ft_to_alt_unit
         else:
             alt = np.full(vel.shape, np.nan, dtype=vel.dtype)
 
+        # get dynamic pressure
+        if _is_q_units_consistent(density_units_in, vel_units_in, q_units_in):
+            q = 0.5 * rho * vel**2
+        else:  # pragma: no cover
+            raise NotImplementedError((density_units_in, vel_units_in, q_units_in))
         results2 = np.dstack([resultsi, eas, q, alt])
         # results2[:, :, self.idensity] = rho
         return results2
@@ -809,6 +813,9 @@ def get_flutter_crossings(self,
             easi, dampi, freqi = remove_eas_range(
                 easi, dampi, freqi, eas_range)
 
+            if len(freqi) == 0:
+                continue
+
             for damping_targeti, damping_requiredi in damping_crossings_dict.items():
                 if dampi.max() < damping_requiredi:
                     continue
@@ -1796,6 +1803,8 @@ def plot_vg_vf(self, fig=None, damp_axes=None, freq_axes=None,
             vel_calc, damping_calc, freq_calc = remove_eas_range(
                 vel, damping, freq, eas_range)
 
+            if len(freq) == 0:
+                continue
             jcolor, color, linestyle2, symbol2, texti, is_removedi = _increment_jcolor(
                 mode, jcolor, color, linestyle, symbol,
                 freq, damping, freq_tol=freq_tol, freq_tol_remove=freq_tol_remove,
@@ -3653,14 +3662,22 @@ def _get_divergence(self,
         dampi = self.results[imode, :, self.idamping].flatten()
         easi = self.results[imode, :, self.ieas].flatten()
         freqi = self.results[imode, :, self.ifreq].flatten()
-
         easi, dampi, freqi = remove_excluded_points(
             easi, dampi, freqi, point_removal)
         easi, dampi, freqi = remove_eas_range(
             easi, dampi, freqi, eas_range)
 
+        if len(freqi) == 0:
+            continue
+
         for freq_targeti, freq_requiredi in freq_crossings_dict.items():
-            if freqi.min() > freq_requiredi:
+            try:
+                freqi_min = freqi.min()
+            except:
+                warnings.warn(f'mode={mode}; freqi={freqi}')
+                raise
+
+            if freqi_min > freq_requiredi:
                 continue
 
             # dfreq sign is flipped because get_zero_crossings

pyNastran/f06/parse_flutter.py

@@ -39,6 +39,7 @@ def make_flutter_response(f06_filename: PathLike,
                           f06_units=None, out_units=None,
                           use_rhoref: bool=False,
                           read_flutter: bool=True,
+                          make_alt: bool=True,
                           log: Optional[SimpleLogger]=None) -> tuple[dict[int, FlutterResponse],
                                                                      dict[str, Any]]:
     """
@@ -57,7 +58,10 @@ def make_flutter_response(f06_filename: PathLike,
         False: assume the density in the table is absolute density
         True: assume the density should be defined by sea level density,
               so density is a density ratio
-
+    read_flutter : bool; default=True
+        don't parse the flutter to extract eigenvalues/mass
+    make_alt : bool; default=True
+        make altitude an output variable
 
     Returns
     -------
@@ -225,7 +229,8 @@ def make_flutter_response(f06_filename: PathLike,
                         in_units=f06_units,
                         use_rhoref=use_rhoref,
                         eigenvector=eigenvectors_array,
-                        eigr_eigi_velocity=eigr_eigi_velocity)
+                        eigr_eigi_velocity=eigr_eigi_velocity,
+                        make_alt=make_alt)
                     response.set_out_units(out_units)
                     #_remove_neutrinos(response, log)
                     flutters[subcase] = response
@@ -426,7 +431,8 @@ def make_flutter_response(f06_filename: PathLike,
                 in_units=f06_units,
                 use_rhoref=use_rhoref,
                 eigenvector=eigenvectors_array,
-                eigr_eigi_velocity=eigr_eigi_velocity)
+                eigr_eigi_velocity=eigr_eigi_velocity,
+                make_alt=make_alt)
             response.set_out_units(out_units)
             flutters[subcase] = response