# Lattice Result Result lres = lsys.results['Test Alpha'] display_markdown(lres) #%% # Plot Lift Distribution axl = None axl = pres.plot_strip_lift_force_distribution(ax=axl, label='pyapm Strip') axl = lres.plot_strip_lift_force_distribution(ax=axl, label='pyvlm Strip') axl = pres.plot_trefftz_lift_force_distribution(ax=axl, label='pyapm Trefftz') axl = lres.plot_trefftz_lift_force_distribution(ax=axl, label='pyvlm Trefftz') #%% # Plot Strip Drag Distribution axd = None axd = pres.plot_strip_drag_force_distribution(ax=axd, label='pyapm Strip') axd = lres.plot_strip_drag_force_distribution(ax=axd, label='pyvlm Strip') axd = pres.plot_trefftz_drag_force_distribution(ax=axd, label='pyapm Trefftz') axd = lres.plot_trefftz_drag_force_distribution(ax=axd, label='pyvlm Trefftz') #%% # Plot Wash Distribution axw = None axw = pres.plot_trefftz_wash_distribution(ax=axw, label='pyapm Trefftz') axw = lres.plot_trefftz_wash_distribution(ax=axw, label='pyvlm Trefftz') #%% # MSH File Output mshfilepath = '../results/' + psys.name + '.msh' panelresult_to_msh(pres, mshfilepath)
psys = panelsystem_from_json(jsonfilepath)
#%%
# Display System
display_markdown(psys)
#%%
# Display Results
for case in psys.results:
pres = psys.results[case]
display_markdown(pres)
#%%
# Mesh File Output
ptrm = psys.results['Positive 1g Cruise']
panelresult_to_msh(ptrm, '../results/Aircraft.msh')
display_markdown(ptrm)
#%%
# Deflect Elevator
pres = ptrm.to_result('Positive 1g Cruise + Controls 25 deg')
pres.set_controls(elevator=25.0, aileron=25.0, rudder=25.0)
panelresult_to_msh(pres, '../results/Aircraft_Trim.msh')
display_markdown(pres)
#%%
# Plot Lift Distribution
axl = ptrm.plot_trefftz_lift_force_distribution()
axl = ptrm.plot_strip_lift_force_distribution(ax=axl)
axl = pres.plot_trefftz_lift_force_distribution(ax=axl)
axl = pres.plot_strip_lift_force_distribution(ax=axl)