Luận án tiến sĩ trường nhiệt độ và trường ẩm trong khuôn tươi

using System; using System.Collections.Generic; using System.ComponentModel; using System.Data; using System.Drawing; using System.Linq; using System.Text; using System.Windows.Forms; using Excel = Microsoft.Office.Interop.Excel; using System.Windows.Forms.DataVisualization.Charting; namespace Tien_nhiet { public partial class Form1 : Form { public Form1() { InitializeComponent(); } public Excel.ApplicationClass app; public Excel.Workbook workbook; public Excel.Worksheet worksheet; public double T_rot; // chieu day vat duc public double s_duc; //chieu day khuon public double s_khuon; //do am ban dau public double w_0_am; //nhiet ban dau public double T_0; //buoc sai khong gian public double denta_x,denta_y; //buoc sai thoi gian public double denta_tg; //độ chênh áp suất biên public double denta_ap; //Các thông số nhiệt lý public double c_l = 0.229; public double c_kho = 0.188; public double c_w = 1.0; public double H = 92.61; public double K_l = 0.115; public double L = 539.0; public double T_l = 605; public double T_s = 515.0; public double T_c = 100.0; public double ro_al = 2796000.0; public double ro_k = 1490000.0; public double anpha_k = 3.4 * Math.Pow(10.0, -7.0); public double[,] T; public double[] epsilon; public double[] T_vd; public double w; public double F_0; public double tong_tg; public int buoc_kg; int N_buoc; int so_buoc_kg; int so_buoc_tg; int save; double P_0,D_cat, Epxilon_xop, Phi, Ro_cat, C_cat, Lamda_cat, Miu_khi, Lamda_khi, C_a, C_w, H_bemat, H_ngoai, R_kk, R_w, K_d, S_b; double[,] A, A_p, T_cat_save, T_g_save,P_a_save,P_w_save,W_save; double[] B, nghiemT, B_p, nghiemP, T_cat_0, T_cat_1, T_g_0, T_g_1, P_a_0, P_a_1, P_w_0, P_w_1,W_0,W_1; double[] truongam,truongL,U, Ro_g, V, Re, S_c, D_w, K_c, P_r, Anpha_g, Ro_a, Ro_w, Cg_Ro_g, h_p, L_w ; double dap1, dap2,dot1,dot2,dot3,dot4; int buoc_save; int kogd; double tg_save; public class AutoClosingMessageBox { System.Threading.Timer _timeoutTimer; string _caption; AutoClosingMessageBox(string text, string caption, int timeout) { _caption = caption; _timeoutTimer = new System.Threading.Timer(OnTimerElapsed, null, timeout, System.Threading.Timeout.Infinite); MessageBox.Show(text, caption); } public static void Show(string text, string caption, int timeout) { new AutoClosingMessageBox(text, caption, timeout); } void OnTimerElapsed(object state) { IntPtr mbWnd = FindWindow(null, _caption); if (mbWnd != IntPtr.Zero) SendMessage(mbWnd, WM_CLOSE, IntPtr.Zero, IntPtr.Zero); _timeoutTimer.Dispose(); } const int WM_CLOSE = 0x0010; [System.Runtime.InteropServices.DllImport("user32.dll", SetLastError = true)] static extern IntPtr FindWindow(string lpClassName, string lpWindowName); [System.Runtime.InteropServices.DllImport("user32.dll", CharSet = System.Runtime.InteropServices.CharSet.Auto)] static extern IntPtr SendMessage(IntPtr hWnd, UInt32 Msg, IntPtr wParam, IntPtr lParam); } public void khoitaothamso(double s1, double s2, double s3, double s4, double s5, double s6, int s7, double s8, double s9, double s10) { w_0_am = s1; T_0 = s2; T_rot = s3; denta_tg = s4; denta_x = s5; tong_tg = s6; buoc_kg = s7; tg_save = s8; denta_y = s9; denta_ap = s10; } public double PW0(double T) { if ((T > 273.0) & (T < 373.0)) { return 0.005 * Math.Pow((Math.E), 0.045 * T); } else { return P_0; } } public void ktraapsuat(double s1, double s2) { dap1 = s1; dap2 = s2; } public void ktranhietdo(double s1, double s2, double s3, double s4) { dot1 = s1; dot2 = s2; dot3 = s3; dot4 = s4; } public void nhapheso() { P_0 = 101325.0; D_cat = 0.0002; Epxilon_xop = 0.436; Phi = 1.44; Ro_cat = 1465; C_cat = 1134.0; Lamda_cat = 0.574; Miu_khi = 2.19 * Math.Pow(10.0, -5.0); Lamda_khi = 0.0317; C_a = 1005.0; C_w = 1846.0; H_bemat = 8370.0; H_ngoai = 20.9; R_kk = 287.1; R_w = 461.5; K_d = (Math.Pow(Epxilon_xop, 3.0)) * D_cat*D_cat / (180.0 * (1 - Epxilon_xop) * (1 - Epxilon_xop) * Phi * Phi); S_b = 6.0 * Phi * (1 - Epxilon_xop) / D_cat; } public void giaihept1(double[,] C, double[] D, int n, double[] nghiem,int erros) { double[,] a; double[] b; erros = 0; a = new double[n, n + 1]; b = new double[n + 1]; for (int i1 = 0; i1 < n; i1++) for (int j1 = 0; j1 < n; j1++) { a[i1, j1] = C[i1, j1]; } for (int i1 = 0; i1 < n; i1++) { a[i1, n] = D[i1]; } double t1; int n1, i, j, k; for (i = 0; i < n; i++) { t1 = Math.Abs(a[i, i]); n1 = i; for (k = i + 1; k < n; k++) { if (Math.Abs(a[k, i]) > t1) { t1 = Math.Abs(a[k, i]); n1 = k; } } int i1, j1, k1; if (Math.Abs(t1) == 0) { //CWnd::MessageBox(_T("pt vo nghiem")); // MessageBox.Show("pt vo nghiem"); erros = 1; } for (i1 = 0; i1 < n + 1; i1++) { b[i1] = a[i, i1]; a[i, i1] = a[n1, i1]; a[n1, i1] = b[i1]; } double tt1 = a[i, i]; for (i1 = i; i1 < n + 1; i1++) { a[i, i1] = a[i, i1] / tt1; } for (j1 = i + 1; j1 < n; j1++) { tt1 = a[j1, i]; for (k1 = i; k1 < n + 1; k1++) { a[j1, k1] = a[j1, k1] - tt1 * a[i, k1]; } } } for (i = n - 1; i > 0; i ) { for (j = 0; j < i; j++) { a[j, n] = a[j, n] - a[j, i] * a[i, n]; a[j, i] = 0; } } for (j = 0; j < n; j++) { nghiem[j] = a[j, n]; } } public void openexcel() { string fileexcel = Application.StartupPath + "\\database.xls"; app = new Excel.ApplicationClass(); app.Visible = false; app.DisplayAlerts = false; workbook = app.Workbooks.Open(fileexcel, Type.Missing, Type.Missing, Type.Missing, Type.Missing, Type.Missing, Type.Missing, Type.Missing, Type.Missing, Type.Missing, Type.Missing, Type.Missing, Type.Missing, Type.Missing, Type.Missing); worksheet = (Excel.Worksheet)workbook.Sheets["Sheet1"]; } public int sokq; public string[,] mangkq; public string luufile; void solve() { nhapheso(); N_buoc=(int)(Math.Floor(tong_tg/denta_tg)); buoc_save = (int)(Math.Floor(tg_save / denta_tg)); save=0; T_cat_0 = new double[ buoc_kg + 1]; T_cat_1 = new double[ buoc_kg + 1]; T_cat_save = new double[200000, buoc_kg + 1]; T_g_0 = new double[buoc_kg + 1]; T_g_1 = new double[buoc_kg + 1]; T_g_save = new double[200000, buoc_kg + 1]; P_a_0 = new double[buoc_kg + 1]; P_a_1 = new double[buoc_kg + 1]; P_a_save = new double[200000, buoc_kg + 1]; P_w_0 = new double[buoc_kg + 1]; P_w_1 = new double[buoc_kg + 1]; P_w_save = new double[200000, buoc_kg + 1]; W_0 = new double[buoc_kg + 1]; W_1 = new double[buoc_kg + 1]; W_save = new double[200000, buoc_kg + 1]; U = new double[ buoc_kg + 1]; Ro_g = new double[ buoc_kg + 1]; V = new double[ buoc_kg + 1]; Re = new double[ buoc_kg + 1]; S_c = new double[ buoc_kg + 1]; D_w = new double[ buoc_kg + 1]; K_c = new double[ buoc_kg + 1]; P_r = new double[ buoc_kg + 1]; Anpha_g = new double[ buoc_kg + 1]; Ro_a=new double[ buoc_kg+1]; Ro_w=new double[ buoc_kg+1]; Cg_Ro_g = new double[buoc_kg + 1]; h_p = new double[buoc_kg + 1]; L_w = new double[ buoc_kg + 1]; truongL = new double[ buoc_kg + 1]; for (int i = 0; i < buoc_kg + 1; i++) { if (i == 0) { T_cat_0[ i] = T_rot; T_g_0[ i] = T_rot; P_w_0[ i] = P_0; P_a_0[i] = 250.0; W_0[ i] = 0.0; } else { T_cat_0[ i] = T_0; T_g_0[ i] = T_0; P_w_0[ i] = PW0(T_cat_0[ i]); P_a_0[ i] = P_0 - P_w_0[ i]; W_0[ i] = w_0_am; } } A= new double[2*buoc_kg+2,2*buoc_kg+2]; B= new double[2*buoc_kg+2]; nghiemT=new double[2*buoc_kg+2]; A_p = new double[2 * buoc_kg + 2, * buoc_kg + 2]; B_p = new double[2 * buoc_kg + 2]; nghiemP = new double[2 * buoc_kg + 2]; for (int j = 0; j < N_buoc-1; j++) { for (int i = 1; i < buoc_kg + 1; i++) { U[ i] = (-K_d / (Miu_khi * denta_x)) * (P_w_0[ i] + P_a_0[ i] - P_w_0[ i - 1] - P_a_0[ i - 1]); if (U[i] < 0) { U[i] = 0; } Ro_w[i] = (P_w_0[i] / (R_w * T_g_0[i])); Ro_a[i] = (P_a_0[i] / (R_kk * T_g_0[i])); Ro_g[i] = Ro_w[i] + Ro_a[i]; V[i] = Miu_khi / Ro_g[i]; Re[i] = (Math.Abs(U[i])* D_cat) / V[ i]; D_w[i] = 2.65 * (Math.Pow(10, -5.0)) * (Math.Pow((T_g_0[i] / 273), 1.78)); S_c[i] = V[i] / D_w[ i]; Cg_Ro_g[i] = C_a * Ro_a[i] + C_w * Ro_w[i]; Anpha_g[i]=Lamda_khi/(1.0*Cg_Ro_g[i]); P_r[i] = V[i] / Anpha_g[i]; if (W_0[i]