Программа по расчету количество проходов для уплотнения грунта катком (неправильно выдает результат)

@ДмитрийБатраков · Регистрация: 21.03.2018

Студворк — интернет-сервис помощи студентам

Пишется программа для определения количества проходов катком для уплотнения грунта, но почему-то выдает всегда только "Уплотнение не возможно". Я новичек в этом. Помогите найти ошибки в коде. И как потом еще реализовать вывод промежуточных результатов в datagridview?

C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
private void MainButton_Click(object sender, EventArgs e)
        {
            {//расчеты
                double q, s, c, dt, dz, ql, k1, f0, K, ay, tu, a, b,
              v, r, eps, t, t0, td, sig, sig0, sigz, tt, kp, mDef, imp, mdef0, mdefz, pm0, pmz,
             zkd, zkd1, zkd2, usl, usl1, usl2, KyP, hor,
             ky0n, exo, exz, pm, et0, SigR0, SigRz, kyzn, te0, ampl, p,
             sig1, df, dp, dd, ee, SigR, Sig0m, SigZm, g, pi, tau, d2, zkdt, Ksig0, nvozd, 
             Kez, Kedef, khitr0, khitrz, ket,
             KuSt, aa, Sig0Pr, mm, Npr, EkWl, A1, A2, sigP, sig0P, Kzvibr, Kzst, Kpq, n,
             qlv, tn, tr, signag, sigraz, t1, t2, mdef;
 
                int  ksy, dgl, Nv, Nsh, Uq0, UqZ, Nst;
                //где
                //    Nst - число проходов в статическом режиме;
                //    Kvibr - коэф - т пересчета глубины распр - я напр - й в вибр-м режиме уплотнения;
                //    Kzst - коэф - т пересчета глубины распр - я напр - й в статическом режиме упл-я}
                //{ Ksig0 = 0.727 - к - т пересчета sig0}
                // Исходные данные
                double Wl = Convert.ToDouble(tb_vlGr.Text); //влажность  
                double Hu = Convert.ToDouble(tb_GlUp.Text); //Глубина уплотнения
                double KyZ = Convert.ToDouble(tb_KuVgI.Text); //Коэффициент уплотнения в глубине
                double Ky0 = Convert.ToDouble(tb_KuNpI.Text); //Коэффициент уплотнения на поверхности
                double Ky0t = Convert.ToDouble(tb_KuNpT.Text); //Коэффициент уплотнения на поверхности тр.
                double KyZt = Convert.ToDouble(tb_KuVgT.Text); //Коэффициент уплотнения в глубине тр.
 
                double m = Convert.ToDouble(tb_MasPerM.Text); //Нагрузка на передний мост, кг
                double m2 = Convert.ToDouble(tb_MasZadM.Text); //Нагрузка на задний мост, кг
                double bb = Convert.ToDouble(tb_Shir.Text) / 1000; //Ширина вибровальца, м
                double d = Convert.ToDouble(tb_Diam.Text) / 1000; //Диаметр переднего вибровальца, м
 
                double F = Convert.ToDouble(tb_ChMin.Text); //Частота Гц
                double P = Convert.ToDouble(tb_VinS1.Text); //Вынуждающая сила кН
                double Ampl = Convert.ToDouble(tb_Am1.Text); //Амплитуда кН
 
                double dGl = Convert.ToDouble(cb_SodGlCh.Text); //Содержание глинистых частиц, %
                double Vk = Convert.ToDouble(tb_Skor.Text) * 1000 / 3600; //Собственная скорость катка, км / ч
 
                K= 0.0602; ay= 0.0017; tu= 0.016; a= 0.5076; b= -16.328; ksy= 2;
                Ksig0 = 0.727;
                g = 9.81;
                pi= 3.14;
                tau= 0.016;
                dz= 0.01; dt= 0.00005;
                bool fi= false; n= 1; t= dt; zkd= 0; ee= 0.005;
                Nst= 1; Kzvibr= 1; Kzst= 5;
                Kedef= 1.00;
                signag = 0;
                td = 0;
                kp = 0;
                pm0 = 0;
                sigz = 0;
                pmz = 0;
                n = 1;
                EkWl = 0;
                Npr = 0;
                switch (cb_TipGr.Text.ToString())
                {
                    case "Песок":
                        Ksig0 = 0.143;
                       
                        break;
                    case "Супесь":
                        Ksig0 = 0.53715;
                        
                        break;
                  
                }
                q= m * g / 1000;
               
                ql= q / bb;
              
               double R = d / 2;
                Kpq= 5.4 - 0.2 * P / q;
                qlv= Kpq * ((P + q) / bb);
                
                while (true)
                {
                    if (n > 50)
                    { //Упл не возможно
                        MessageBox.Show("Уплотнение не возможно", "Внимание!", MessageBoxButtons.OK, MessageBoxIcon.Error);
                        goto l7;
                    }
                    KyP = 0.7 * Ky0 + 0.3 * KyZ;
                    switch (cb_TipGr.Text.ToString())
                    {
                        case "Песок":
                         
                            mdef = 12 * Math.Pow(KyP, 15.506) * 1.626 * Math.Pow((Hu / (d / 2)), -0.606);
                            break;
                        case "Супесь":
                           
                            mdef = 16.97 * Math.Pow(KyP, 7.506) * 0.726 * Math.Pow((Hu / (d / 2)), -1.179);
                            break;
                        default:
                            mdef = Kedef * (74.28 * Math.Pow(KyP, 2) - 34.51 * KyP - 9.96);
                            break;
                    }
                   
                    hor = 4 * Math.Sqrt((qlv * R) / (1000 * mdef));
                    Kez = 1;
                   
                    if ((P == 0) || (F == 0))
                    {
                        sig0 = (1 / 1000) * Math.Sqrt(((ql) * 1000 * mdef) / R);
                        t0 = Vk / (2 * hor);
                    }
                    else
                    {
                        sig0 = Ksig0 * 0.5 * Math.Sqrt(qlv * mdef / R / 1000);
                        t0 = 1 / (4 * F);
                    }//Мпа
                    t1 = t0 - dt;
                    te0 = t0 / 2;
                    tt = 2 * t0;
                    nvozd = hor * F / Vk; // Число воздействий за 1 проход катка
                    c = (Ky0 - KyZ) / (Hu / hor); // Коэффициент изменения плотности по глубине
                    SigR = 0.1639 + 0.1815 * Math.Pow(KyP, 10.206) * Math.Pow(dGl, 0.8344); //предел точности Мпа
                    while (t1 >= dt)
                    {
                        signag = sig0 * t1 / t0;
                        td = 2 * (t0 - t1);
                    }
                    l9:
                    zkd = 0;
                    zkd1 = zkd;
                    if (Ky0 == KyZ)
                    {
                        zkd = Math.Sqrt(td / (a * K * signag * Math.Exp(b * Math.Log10(Ky0))));
                        goto l3;
                    }
                    else
                    while(true)
                    {
                        l5:
                        zkd2 = zkd1 + dz;
                        usl1 = Math.Pow((Ky0 - c * zkd1), b) * (ay * zkd1 + K * signag * Math.Pow(zkd1, 2) + (ay * (Ky0 - c * zkd1)) / (c * (b + 1))) - (ay * Math.Pow(Ky0, (b + 1))) / (c * (b + 1)) - td / a;
                        usl2 = Math.Pow((Ky0 - c * zkd2), b) * (ay * zkd2 + K * signag * Math.Pow(zkd2, 2) + (ay * (Ky0 - c * zkd2)) / (c * (b + 1))) - (ay * Math.Pow(Ky0, (b + 1))) / (c * (b + 1)) - td / a;
                            if ((Math.Abs(usl1) < ee) || (Math.Abs(usl2) < ee))
                        {
                            goto l1;
                        }
                        else
                        {
                            zkd1 = zkd2;
                            goto l5;
                        }
                    }
                    l1:
                    if (Math.Abs(usl1) < Math.Abs(usl2))
                    { zkd = zkd1; }
                    else
                    { zkd = zkd2; }
 
                    l3:
                        if ((zkd) > (Hu / hor))
                        {
                            sigz = signag;
                            goto l2;
                        }
                        else
                        {
                            t1 = t1 - dt;
                        }
                    l2:
                    if (n == 1)
                    {
                        kp = 3.6;
                    }
                    if (n == 2)
                    {
                        kp = 1.8;
                    }
                    if (n == 3)
                    {
                        kp = 1.2;
                    }
 
                    if (n > 3)
                    {
                        kp = 1;
                    }
                    et0 = 6.447 + 191.9 * (Math.Exp(14.046 * Math.Log10(Ky0)));
                    switch (cb_TipGr.Text.ToString())
                    {
                        case "Песок":
                            pm0 = (1.298 + 42.141 * Math.Pow(Ky0, 25.183));
                            pmz = (1.298 + 42.141 * Math.Pow(KyZ, 25.183));
                            break;
                        case "Супесь":
                            Ksig0 = 0.53715;
                            pm0 = (3.298 + 81.141 * Math.Pow(Ky0, 15.183));
                            pmz = (3.298 + 81.141 * Math.Pow(KyZ, 15.183));
                            break;
                    }
                    khitr0 = 0.588;
                    exo = kp * khitr0 * sig0 * (1 / pm0 + (1 / (et0 * ksy)) * Math.Log10(1 + ksy * te0 * nvozd));
                    switch (cb_TipGr.Text.ToString())
                    {
                        case "Песок":
                           
                        pmz = (1.298 + 42.141 * Math.Pow(KyZ, 25.183));
                        break;
                        case "Супесь":
                            Ksig0 = 0.53715;
                       
                        pmz = (3.298 + 81.141 * Math.Pow(KyZ, 15.183));
                        break;
                    }
 
                    khitrz = 0.588;
                    et0 = 6.447 + 191.9 * (Math.Exp(14.046 * Math.Log10(KyZ)));
                    exz = Kez * kp * khitrz * sigz * (1 / pmz + (1 / (et0 * ksy)) * Math.Log10(1 + ksy * te0 * nvozd));
                    Ky0 = Ky0 / (1 - exo);
                    if (Ky0 >= 1.1)
                    { Ky0 = 1.1; }
                    KyZ = KyZ / (1 - exz);
                    if (KyZ >= 1.1)
                    { KyZ = 1.1; }
                    c = (Ky0 - KyZ) / (Hu / hor);
                    KyP = 0.7 * Ky0 + 0.3 * KyZ;
                    SigR0 = 0.1639 + 0.1815 * Math.Pow(Ky0t, 10.206) * Math.Pow(dGl, 0.8344);
                    Sig0m = (0.171 + 0.6927 * Math.Pow(Ky0t, 0.22)) * SigR0;
                    SigRz = 0.1639 + 0.1815 * Math.Pow(KyZt, 10.206) * Math.Pow(dGl, 0.8344);
                    SigZm = (0.06883 + 0.36032 * Math.Pow(KyZt, 31.0534)) * SigRz;
                    mm = m;
                    tb_ChP.Text = Npr.ToString();
                    tb_KPnP.Text = KyZ.ToString();
                    tb_KPvGL.Text = Ky0.ToString();
 
 
                    //if ((Ky == KyZt) && (Ky0 == Ky0t))
                    //{ goto l6; }
 
                    if (Ky0 < Ky0t)
                    {
                        goto l8;
                    }
                    else
                    {
                        if (KyZ < KyZt)
                        { goto l8; }
                        else
                        {
                            if (sig0 < Sig0m)
                            {
                                //MessageBox.Show("Уплотнение не возможно", "Внимание!", MessageBoxButtons.OK, MessageBoxIcon.Error);
                                goto l7;
                            }
                            else
                            {
                                if (sigz < SigZm)
                                {
                                    //MessageBox.Show("Уплотнение не возможно", "Внимание!", MessageBoxButtons.OK, MessageBoxIcon.Error);
                                    goto l7;
                                }
                                else
                                { goto l10; }
                            }
                        }
                    }
 
                    l8:
 
                    n = n + 1;
 
                }
                l6:
 
                if (Wl == 1)
                { EkWl = 1; }
                else
                {
                    A1 = 25 * (1 + 0.03 * (dGl - 3));
                    A2 = 1.29 - 2.23 * Math.Exp(-0.52 * dGl);
                    EkWl = Math.Pow((A1 * Math.Abs(1 - Wl)), A2);
                }
 
                /*switch (cb_TipVtBand.Text.ToString())
                {
                    case "Гладкий валец":
                        Nv = 2;
                        Nsh = 0;
                        break;
                    case "Блок пневмоколёс":
                        Nv = 1;
                        Nsh = 1;
                        break;
                    case "Ведущие пневмоколёса":
                        Nv = 1;
                        Nsh = 0;
                        break;
                }*/
                n = Math.Ceiling(n * EkWl);
                Npr = n;
                Npr = Math.Ceiling(Npr / (Math.Pow(Nv, 0.8)));
                Npr = Npr + 1 - Nsh;
                if (Npr < 3)
                { Npr = Npr + 1; }
                l7:
                MessageBox.Show("Уплотнение не возможно, т.к мин упл не может быть достигнуто", "Внимание!", MessageBoxButtons.OK, MessageBoxIcon.Error);
                l10:
                int end = 0;
                //end
            }
 
 
        }

@Lexeq · 21.03.2018, 13:55

ДмитрийБатраков, ну вот вы на собственном опыте узнаете, чем плох goto и метод в 300 строчек. Думаю, вы уже сами не очень ориентируетесь в том что за чем там выполняется, а для остальных это вообще полная абракадабра. Единственный нормальный вариант тут - всё переделать, пользуясь благами, которые даёт ООП.

Новые блоги и статьи Все статьи Все блоги /
Как я обхитрил таблицу Word Alexander-7 21.03.2026 Когда мигает курсор у внешнего края таблицы, и нам надо перейти на новую строку, а при нажатии Enter создается новый ряд таблицы с ячейками, то мы вместо нервных нажатий Энтеров мы пишем любые буквы. . .	Krabik - рыболовный бот для WoW 3.3.5a AmbA 21.03.2026 без регистрации и смс. Это не торговля, приложение не содержит рекламы. Выполняет свою непосредственную задачу - автоматизацию рыбалки в WoW - и ничего более. Однако если админы будут против -. . .	Программный отбор значений справочника Maks 21.03.2026 Установка программного отбора значений справочника "Сотрудники" из модуля формы документа. В качестве фильтра для отбора служит предопределенное значение перечислений. Процедура. . .	Переходник USB-CAN-GPIO Eddy_Em 20.03.2026 Достаточно давно на работе возникла необходимость в переходнике CAN-USB с гальваноразвязкой, оный и был разработан. Однако, все меня терзала совесть, что аж 48-ногий МК используется так тупо: просто. . .
Оттенки серого Argus19 18.03.2026 Оттенки серого Нашёл в интернете 3 прекрасных модуля: Модуль класса открытия диалога открытия/ сохранения файла на Win32 API; Модуль класса быстрого перекодирования цветного изображения в оттенки. . .	SDL3 для Desktop (MinGW): Рисуем цветные прямоугольники с помощью рисовальщика SDL3 на Си и C++ 8Observer8 17.03.2026 Содержание блога Финальные проекты на Си и на C++: finish-rectangles-sdl3-c. zip finish-rectangles-sdl3-cpp. zip	Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая ссылка» (hard link),. . .	[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .

@Lexeq 1151 / 743 / 483 Регистрация: 21.01.2014 Сообщений: 1,903
	21.03.2018, 13:55
	ДмитрийБатраков, ну вот вы на собственном опыте узнаете, чем плох goto и метод в 300 строчек. Думаю, вы уже сами не очень ориентируетесь в том что за чем там выполняется, а для остальных это вообще полная абракадабра. Единственный нормальный вариант тут - всё переделать, пользуясь благами, которые даёт ООП. 0