加载中
加载中
表情图片
评为精选
鼓励
加载中...
分享
加载中...
文件下载
加载中...
修改排序
加载中...
文号932716
2183
1
收藏
Python脚本计算内弹道（似乎有问题）
ZiRui2024/06/06喷气推进 IP:广东 
Python
from tkinter import Tk, Label, Entry, Button, StringVar
from tkinter.messagebox import showerror
from tkinter.filedialog import askdirectory
from math import sqrt

root = Tk()
root.title("推力室内弹道")


def calculate():
    try:
        gamma = float(gamma_entry.get())
        mass_flow = float(mass_flow_entry.get())
        pressure = float(pressure_entry.get())
        temperature = float(temperature_entry.get())
        gas_constant = float(gas_constant_entry.get())
        pressure_export = float(pressure_export_entry.get())
        characteristic_length = float(characteristic_length_entry.get())
        diameter_convergence_ratio_length = float(diameter_convergence_ratio_length_entry.get())
        convergence_half_angle = float(convergence_half_angle_entry.get())

        G = sqrt(gamma) * (2 / (gamma + 1)) ** ((gamma + 1) / (2 * (gamma - 1)))
        throat_area = mass_flow / (G * (pressure / sqrt(gas_constant * temperature)))
        compute_throat_diameter = sqrt(throat_area / 3.14159) * 2

        expansion_area_ratio = ((2 / (gamma + 1)) ** (1 / (gamma - 1)) * (
                (pressure / pressure_export) ** (1 / gamma))) / sqrt(
            ((gamma + 1) / (gamma - 1)) * (1 - ((pressure_export / pressure) ** ((gamma - 1) / gamma))))
        compute_export_diameter = sqrt(throat_area * expansion_area_ratio / 3.14159) * 2

        chamber_volume = throat_area * characteristic_length
        compute_chamber_radius = round(compute_throat_diameter * diameter_convergence_ratio_length + 0.4, 0) / 2
        compute_convergence_length = (90 / convergence_half_angle - 1) * \
                                     (compute_chamber_radius - (compute_throat_diameter / 2))
        chamber_volume -= float(compute_convergence_length) * 3.14159 * float(compute_chamber_radius ** 2 +
                                                                  (compute_throat_diameter / 2) ** 2 +
                                                                  compute_throat_diameter / 2 * compute_chamber_radius) / 3000
        compute_cylinder_length = chamber_volume * 1000 / (compute_chamber_radius ** 2 * 3.14159)

        chamber_diameter.set(f"{compute_chamber_radius * 2}")
        cylinder_length.set(f"{round(compute_cylinder_length + 0.4, 0)}")
        convergence_length.set(f"{round(compute_convergence_length + 0.4, 0)}")
        throat_diameter.set(f"{round(compute_throat_diameter, 3)}")
        expansion_diameter.set(f"{round(compute_export_diameter, 3)}")
        expansion_ratio.set(f"{round(expansion_area_ratio, 3)}")
    except Exception as e:
        showerror("错误", f"计算过程中发生错误：\n{str(e)}")
        print(e)


mass_flow_label = Label(root, text="质量流量（Kg/s）：")
mass_flow_label.grid(row=0, column=0, padx=5, pady=5)
mass_flow_entry = Entry(root)
mass_flow_entry.grid(row=0, column=1, padx=5, pady=5)

gamma_label = Label(root, text="比热比（γ）：")
gamma_label.grid(row=1, column=0, padx=5, pady=5)
gamma_entry = Entry(root)
gamma_entry.grid(row=1, column=1, padx=5, pady=5)

pressure_label = Label(root, text="压力（MPa）：")
pressure_label.grid(row=2, column=0, padx=5, pady=5)
pressure_entry = Entry(root)
pressure_entry.grid(row=2, column=1, padx=5, pady=5)

temperature_label = Label(root, text="温度（K）：")
temperature_label.grid(row=3, column=0, padx=5, pady=5)
temperature_entry = Entry(root)
temperature_entry.grid(row=3, column=1, padx=5, pady=5)

gas_constant_label = Label(root, text="气体常数（J/（Kg·K））：")
gas_constant_label.grid(row=4, column=0, padx=5, pady=5)
gas_constant_entry = Entry(root)
gas_constant_entry.grid(row=4, column=1, padx=5, pady=5)

pressure_export_label = Label(root, text="出口压力（MPa）：")
pressure_export_label.grid(row=5, column=0, padx=5, pady=5)
pressure_export_entry = Entry(root)
pressure_export_entry.grid(row=5, column=1, padx=5, pady=5)

characteristic_length_label = Label(root, text="特征长度（m）：")
characteristic_length_label.grid(row=6, column=0, padx=5, pady=5)
characteristic_length_entry = Entry(root)
characteristic_length_entry.grid(row=6, column=1, padx=5, pady=5)

diameter_convergence_ratio_length_label = Label(root, text="直径收敛比（Dc/Dt）：")
diameter_convergence_ratio_length_label.grid(row=7, column=0, padx=5, pady=5)
diameter_convergence_ratio_length_entry = Entry(root)
diameter_convergence_ratio_length_entry.grid(row=7, column=1, padx=5, pady=5)

convergence_half_angle_label = Label(root, text="敛段半角（°）：")
convergence_half_angle_label.grid(row=8, column=0, padx=5, pady=5)
convergence_half_angle_entry = Entry(root)
convergence_half_angle_entry.grid(row=8, column=1, padx=5, pady=5)


def save():
    with open(f"{askdirectory()}/液机推力室内弹道计算.txt", 'w+', encoding="utf-8") as file:
        file.write(f"推力室形状：\n")
        file.write(f"燃烧室直径Dc（mm）：{chamber_diameter.get()}\n")
        file.write(f"圆筒段长度Lcyl（mm）：{cylinder_length.get()}\n")
        file.write(f"喷喉直径Dt（mm）：{throat_diameter.get()}\n")
        file.write(f"收敛段长度（mm）：{convergence_length.get()}\n")
        file.write(f"出口直径De（mm）：{expansion_diameter.get()}\n")
        file.write(f"膨胀比Ae/At：{expansion_ratio.get()}\n")


calculate_button = Button(root, text="计算", command=calculate)
calculate_button.grid(row=9, column=0, columnspan=2, padx=5, pady=5, sticky="ew")

save_button = Button(root, text="保存", command=save)
save_button.grid(row=10, column=0, columnspan=2, padx=5, pady=5, sticky="ew")

chamber_diameter = StringVar()
chamber_diameter_label = Label(root, text='燃烧室直径Dc（mm）：')
chamber_diameter_label.grid(row=11, column=0)
chamber_diameter_value = Label(root, textvariable=chamber_diameter)
chamber_diameter_value.grid(row=11, column=1)

cylinder_length = StringVar()
cylinder_length_label = Label(root, text='圆筒段长度Lcyl（mm）：')
cylinder_length_label.grid(row=12, column=0)
cylinder_length_value = Label(root, textvariable=cylinder_length)
cylinder_length_value.grid(row=12, column=1)

throat_diameter = StringVar()
throat_diameter_label = Label(root, text='喷喉直径Dt（mm）：')
throat_diameter_label.grid(row=13, column=0)
throat_diameter_value = Label(root, textvariable=throat_diameter)
throat_diameter_value.grid(row=13, column=1)

convergence_length = StringVar()
convergence_length_label = Label(root, text='收敛段长度（mm）：')
convergence_length_label.grid(row=14, column=0)
convergence_length_value = Label(root, textvariable=convergence_length)
convergence_length_value.grid(row=14, column=1)

expansion_diameter = StringVar()
expansion_diameter_label = Label(root, text='出口直径De（mm）：')
expansion_diameter_label.grid(row=15, column=0)
expansion_diameter_value = Label(root, textvariable=expansion_diameter)
expansion_diameter_value.grid(row=15, column=1)

expansion_ratio = StringVar()
expansion_ratio_label = Label(root, text='膨胀比Ae/At：')
expansion_ratio_label.grid(row=16, column=0)
expansion_ratio_value = Label(root, textvariable=expansion_ratio)
expansion_ratio_value.grid(row=16, column=1)

root.mainloop()

算出来的数据可能有问题，给大家提供一个可能正确的方法

来自：航空航天 / 喷气推进
1
全部 只看作者
新版本公告
~~空空如也
ZiRui 作者
8个月4天前 IP:广东
934967
1楼

收敛段长度计算有误，已改正并改为PyQt5：

Python
from PyQt5.QtWidgets import QApplication, QMainWindow, QLabel, QLineEdit, QPushButton, QVBoxLayout, QHBoxLayout, \
    QWidget, QFileDialog, QMessageBox
from PyQt5.QtCore import Qt
from math import sqrt, tan, radians


class RocketChamberDesign(QMainWindow):
    def __init__(self):
        super().__init__()
        self.setWindowTitle("推力室内弹道")
        self.setGeometry(100, 100, 600, 400)

        self.setStyleSheet("""
            QMainWindow {
                background-color: #f0f0f0;
            }

            QLabel {
                font-size: 14px;
                color: #333333;
            }

            QLineEdit {
                font-size: 14px;
                border: 1px solid #cccccc;
                padding: 5px;
                background-color: white;
                border-radius: 5px;
            }

            QPushButton {
                font-size: 16px;
                padding: 10px 20px;
                border: none;
                background-color: #007acc;
                color: white;
                border-radius: 5px;
            }

            QPushButton:hover {
                background-color: #005fa3;
            }
        """)

        self.initUI()

    def initUI(self):
        central_widget = QWidget(self)
        main_layout = QVBoxLayout()
        central_widget.setLayout(main_layout)
        self.setCentralWidget(central_widget)

        # Labels and LineEdits
        self.create_input("质量流量（Kg/s）：", main_layout)
        self.create_input("气体比热比（γ）：", main_layout)
        self.create_input("燃室压力（MPa）：", main_layout)
        self.create_input("燃室温度（K）：", main_layout)
        self.create_input("气体常数（J/(Kg·K)）：", main_layout)
        self.create_input("出口压力（MPa）：", main_layout)
        self.create_input("特征长度（m）：", main_layout)
        self.create_input("径收敛比（Dc/Dt）：", main_layout)
        self.create_input("敛段半角（°）：", main_layout)

        # Button layout
        button_layout = QHBoxLayout()
        calculate_button = QPushButton("计算", self)
        calculate_button.clicked.connect(self.calculate)
        button_layout.addWidget(calculate_button, alignment=Qt.AlignHCenter)

        save_button = QPushButton("保存", self)
        save_button.clicked.connect(self.save)
        button_layout.addWidget(save_button, alignment=Qt.AlignHCenter)

        # Add the button layout to the main layout
        main_layout.addLayout(button_layout)

        # Output labels and line edits
        output_layout = QVBoxLayout()
        self.create_output("燃烧室直径Dc（mm）：", output_layout)
        self.create_output("圆筒段长度Lcyl（mm）：", output_layout)
        self.create_output("喷喉直径Dt（mm）：", output_layout)
        self.create_output("收敛段长度（mm）：", output_layout)
        self.create_output("出口直径De（mm）：", output_layout)
        self.create_output("膨胀比Ae/At：", output_layout)

        # Add the output layout to the main layout
        main_layout.addLayout(output_layout)

    def create_input(self, label_text, layout):
        input_layout = QHBoxLayout()

        label = QLabel(label_text, self)
        label.setAlignment(Qt.AlignRight | Qt.AlignVCenter)
        input_layout.addWidget(label)

        line_edit = QLineEdit(self)
        line_edit.setAlignment(Qt.AlignLeft)
        input_layout.addWidget(line_edit)

        layout.addLayout(input_layout)

        attr_name = self._format_attribute_name(label_text)
        setattr(self, attr_name, line_edit)
        print(f"Created attribute: {attr_name}")

    def create_output(self, label_text, layout):
        output_layout = QHBoxLayout()

        label = QLabel(label_text, self)
        label.setAlignment(Qt.AlignRight | Qt.AlignVCenter)
        output_layout.addWidget(label)

        line_edit = QLineEdit(self)
        line_edit.setReadOnly(True)
        output_layout.addWidget(line_edit)

        layout.addLayout(output_layout)

        attr_name = self._format_attribute_name(label_text) + '_output'
        setattr(self, attr_name, line_edit)
        print(f"Created attribute: {attr_name}")

    @staticmethod
    def _format_attribute_name(text):
        # Replace non-alphanumeric characters with underscores and remove leading numbers
        return ''.join([char if char.isalnum() else '_' for char in text]).lstrip('_')

    def calculate(self):
        try:
            gamma = float(getattr(self, self._format_attribute_name("比热比（γ）：")).text())
            mass_flow = float(getattr(self, self._format_attribute_name("质量流量（Kg/s）：")).text())
            pressure = float(getattr(self, self._format_attribute_name("压力（MPa）：")).text())
            temperature = float(getattr(self, self._format_attribute_name("温度（K）：")).text())
            gas_constant = float(getattr(self, self._format_attribute_name("气体常数（J/(Kg·K)）：")).text())
            pressure_export = float(getattr(self, self._format_attribute_name("出口压力（MPa）：")).text())
            characteristic_length = float(getattr(self, self._format_attribute_name("特征长度（m）：")).text())
            diameter_convergence_ratio_length = float(
                getattr(self, self._format_attribute_name("直径收敛比（Dc/Dt）：")).text())
            convergence_half_angle = float(getattr(self, self._format_attribute_name("敛段半角（°）：")).text())

            G = sqrt(gamma) * (2 / (gamma + 1)) ** ((gamma + 1) / (2 * (gamma - 1)))
            throat_area = mass_flow / (G * (pressure / sqrt(gas_constant * temperature)))
            compute_throat_diameter = sqrt(throat_area / 3.14159) * 2

            expansion_area_ratio = ((2 / (gamma + 1)) ** (1 / (gamma - 1)) * (
                    (pressure / pressure_export) ** (1 / gamma))) / sqrt(
                ((gamma + 1) / (gamma - 1)) * (1 - ((pressure_export / pressure) ** ((gamma - 1) / gamma))))
            compute_export_diameter = sqrt(throat_area * expansion_area_ratio / 3.14159) * 2

            chamber_volume = throat_area * characteristic_length
            compute_chamber_radius = round(compute_throat_diameter * diameter_convergence_ratio_length + 0.4, 0) / 2
            compute_convergence_length = (compute_chamber_radius * 2 - compute_throat_diameter) / (
                    2 * tan(radians(convergence_half_angle)))
            chamber_volume -= float(compute_convergence_length) * 3.14159 * float(compute_chamber_radius ** 2 +
                                                                                  (compute_throat_diameter / 2) ** 2 +
                                                                                  compute_throat_diameter / 2 * compute_chamber_radius) / 3000
            compute_cylinder_length = chamber_volume * 1000 / (compute_chamber_radius ** 2 * 3.14159)

            getattr(self, self._format_attribute_name("燃烧室直径Dc（mm）：") + '_output').setText(
                f'{compute_chamber_radius * 2}')
            getattr(self, self._format_attribute_name("圆筒段长度Lcyl（mm）：") + '_output').setText(
                f'{round(compute_cylinder_length + 0.4, 0)}')
            getattr(self, self._format_attribute_name("喷喉直径Dt（mm）：") + '_output').setText(
                f'{round(compute_throat_diameter, 3)}')
            getattr(self, self._format_attribute_name("收敛段长度（mm）：") + '_output').setText(
                f'{round(compute_convergence_length + 0.4, 0)}')
            getattr(self, self._format_attribute_name("出口直径De（mm）：") + '_output').setText(
                f'{round(compute_export_diameter, 3)}')
            getattr(self, self._format_attribute_name("膨胀比Ae/At：") + '_output').setText(
                f'{round(expansion_area_ratio, 3)}')
        except Exception as e:
            QMessageBox.critical(self, "错误", f"计算过程中发生错误：\n{str(e)}")
            print(e)

    def save(self):
        filename, _ = QFileDialog.getSaveFileName(self, "保存文件", "", "Text Files (*.txt)")
        if filename:
            with open(filename, 'w', encoding="utf-8") as file:
                file.write("推力室设计参数：\n")
                file.write(f"燃烧室直径Dc（mm）：{getattr(self, '燃烧室直径Dc_mm___output').text()}\n")
                file.write(f"圆筒段长度Lcyl（mm）：{getattr(self, '圆筒段长度Lcyl_mm___output').text()}\n")
                file.write(f"喷喉直径Dt（mm）：{getattr(self, '喷喉直径Dt_mm___output').text()}\n")
                file.write(f"收敛段长度（mm）：{getattr(self, '收敛段长度_mm___output').text()}\n")
                file.write(f"出口直径De（mm）：{getattr(self, '出口直径De_mm___output').text()}\n")
                file.write(f"膨胀比Ae/At：{getattr(self, '膨胀比Ae_At__output').text()}")


if __name__ == '__main__':
    import sys

    app = QApplication(sys.argv)
    main_window = RocketChamberDesign()
    main_window.show()
    app.exec_()
    app.quit()


引用
评论
加载评论中，请稍候...
200字以内，仅用于支线交流，主线讨论请采用回复功能。
折叠评论

想参与大家的讨论？现在就 登录 或者 注册。

所属专业
喷气推进
所属分类
测试仿真
上级专业
航空航天
同级专业
航空技术航天技术火箭燃料
ZiRui
进士 机友 笔友
文章
6
回复
38
学术分
0
2023/12/12注册，6时18分前活动

<div class="btn btn-default btn-xs" onclick="changeUsername()">修改</div> <a href="XXXXXXXXXXXXXXXXXXXXXXX">Web</a>

主体类型：个人
所属领域：无
认证方式：手机号
IP归属地：广东
名片私信
作者最新文章
Python脚本+PyQt5计算内弹道（更正收敛半角转收敛段长度，GUI由AI编写）
喷气推进ZiRui 8个月4天前
[新人向]拉瓦尔管简单计算
喷气推进ZiRui 9个月6天前
Python脚本计算内弹道（似乎有问题）
喷气推进ZiRui 10个月1天前
高氯酸铵分解事故
化学ZiRui 10个月4天前
化学火箭发动机燃料计算软件“Rocket Propulsion Analysis”教程，及发动机设计
喷气推进ZiRui 1年1个月前
相似文章推荐
[贴图]火箭上搭载摄像MP4进行航拍实验！
喷气推进AsyncMaster 17年2个月前
环形气动塞式喷管塞锥型面设计+实例计算
喷气推进LS 5年4个月前
精心打造的液体发动机
喷气推进拉风 13年7个月前
单组元（过氧化氢）液机火箭发射成功！
喷气推进乖雪狼 12年9个月前
过氧化氢发动机试车
喷气推进乖雪狼 14年0个月前
两股直流互击喷注器计算软件
喷气推进就是氧化铁喽 6个月1天前
实际喷管基于理论的各指标修正系数评估和损耗分析
喷气推进lyhhjahj 3年3个月前
含有交叉耦合系数的陀螺传感器补偿公式
喷气推进xuetianwuran 13年9个月前
液体火箭发动机推力室设计
喷气推进LSEQY 3个月8天前
年终献礼---RNX静推力及发射试验
喷气推进乖老汉 13年5个月前
插入公式
评论控制
加载中...
文号：{{pid}}
投诉或举报
加载中...
{{tip}}
请选择违规类型：
{{reason.type}}

空空如也

提交成功！
笔记
{{note.content}}
{{n.user.username}}
{{fromNow(n.toc)}} {{n.status === noteStatus.disabled ? "已屏蔽" : ""}} {{n.status === noteStatus.unknown ? "正在审核" : ""}} {{n.status === noteStatus.deleted ? '已删除' : ''}}
  • 编辑
  • 删除
  • {{n.status === 'disabled' ? "解除屏蔽" : "屏蔽" }}
我也是有底线的