00001 #ifndef NORMMAT_H 00002 #define NORMMAT_H 00003 00004 #include "vector.h" 00005 #include "matrix.h" 00006 #include "iotools.h" 00007 00008 /** 00009 class contains material based on Eurocode 00010 00011 components of array other 00012 00013 00014 */ 00015 class normmat 00016 { 00017 public: 00018 00019 normmat (void); 00020 ~normmat (void); 00021 00022 void read (XFILE *in); 00023 00024 00025 //modul pruznosti betonu v tlaku prislusny 0,6 fc v kPa 00026 void compute_ecnom (); 00027 00028 //pretvoreni oceli pro stanoveni modulu pruznosti e0cn 00029 void compute_epsuk (); 00030 00031 //modul pruznosti odpovidajici zplastizovani oceli v kPa 00032 void compute_e0cn (); 00033 00034 //soucinitel pro definovani pracovniho diagramu betonu 00035 void compute_k (); 00036 00037 //plocha prurezu 00038 void compute_a (); 00039 00040 //ucinna vyska prurezu 00041 void compute_d (); 00042 00043 //pevnost betonu v tlaku uvazovana pro vypocet 00044 void compute_fc (); 00045 00046 //svisla koncova sila ve stycniku i 00047 void compute_z1 (long ipp,double lp); 00048 00049 //momentova koncova sila ve stycniku i 00050 void compute_z2 (long ipp,double lp); 00051 00052 //vodorovna pomocna koncova sila ve stycniku i 00053 void compute_z3p (long ipp,double lp); 00054 00055 //vodorovna koncova sila ve stycniku i 00056 void compute_z3 (); 00057 00058 //vyztuz pro tahovou oblast od ucinku ohyboveho momentu as1 00059 double compute_as1 (double lp); 00060 00061 //vyztuz pro tlakovou oblast od ucinku ohyboveho momentu as2 00062 double compute_as2 (double lp); 00063 00064 //ohybovy moment na dannem elementu v polovine intervalu 00065 void compute_m (double lp); 00066 00067 //absolutni hodnota ohyboveho momentu na dannem elementu v polovine intervalu 00068 void compute_mpp (double lp); 00069 00070 //normalova sila na dannem elementu v polovine intervalu 00071 double compute_n0 (double lp); 00072 00073 //pocet iteraci 00074 long compute_iter (); 00075 00076 //funkce ktera rusi vliv betonove casti prurezu 00077 void compute_tcv (); 00078 00079 //VLIV NORMALOVE SILY 00080 00081 //pretvoreni od ucinku normalove sily 00082 double compute_epsn (); 00083 00084 //maximalni tlakova normalova sila, kterou lze zatizit prurez 00085 double compute_nmax (); 00086 00087 //prvni clen pro linearizaci 00088 double compute_an1 (); 00089 00090 //druhy clen pro linearizaci 00091 double compute_an2 (); 00092 00093 //prace na betonove casti prurezu pro prubeh normalove sily (-z3 + fn x); pres element delky lp 00094 double compute_in1 (double lp); 00095 00096 //prace na ocelove casti prurezu pro prubeh normalove sily (-z3 + fn x); pres element delky lp 00097 double compute_in2 (double lp); 00098 00099 //celkova prace na prurezu pro prubeh normalove sily (-z3 + fn x); pres element delky lp 00100 double compute_kinh (); 00101 00102 //pomocna prace na prurezu pro prubeh normalove sily (-z3 + fn x); pres element delky lp 00103 double compute_kind (double lp); 00104 00105 //pomer prace vnitrnich sil a pomocne prace vnitrnich sil na dannem elementu od ucinku normalove sily 00106 double compute_kin (); 00107 00108 //TLACENA VYZTUZ + NORMALOVA SILA + x = h 00109 00110 //pretvoreni betonu, kdy tlacena plocha x=h 00111 double compute_epschr1 (); 00112 00113 //hodnota ohyboveho momentu, kdy tlacena plocha x=h a nedochazi k teceni oceli As1 a As2 00114 double compute_mhr1 (); 00115 00116 //TLACENA VYZTUZ + NORMALOVA SILA + TECENI V As1 + x < h 00117 00118 //pretvoreni betonu, kdy dochazi k teceni vyztuze As1, ale netece As2 a ani nedochazi k drceni betonu 00119 double compute_epschr2 (); 00120 00121 //hodnota ohyboveho momentu, kdy dochazi k teceni vyztuze As1, ale netece As2 a ani nedochazi k drceni betonu 00122 double compute_mhr2 (); 00123 00124 //TLACENA VYZTUZ + NORMALOVA SILA + TECENI V As2 + x < h 00125 00126 //pretvoreni betonu, kdy dochazi k teceni vyztuze As2, ale netece As1 a ani nedochazi k drceni betonu, x < h 00127 double compute_epschr3 (); 00128 00129 //hodnota ohyboveho momentu, kdy dochazi k teceni vyztuze As2, ale netece As1 a ani nedochazi k drceni betonu, x < h 00130 double compute_mhr3 (); 00131 00132 //TLACENA VYZTUZ + NORMALOVA SILA + TECENI V As2 + x > h 00133 00134 //pretvoreni betonu, kdy dochazi k teceni vyztuze As2, ale netece As1 a ani nedochazi k drceni betonu, x > h 00135 double compute_epschr4 (); 00136 00137 //maximalni pretvoreni betonu od ucinku ohyboveho momentu 00138 double compute_epscmax1 (); 00139 00140 //pretvoreni betonu, kdy uz tece vyztuz As2, netece As1 a ani nedochazi k drceni betonu, ale dochazi ke zmene tlacene plochy x = h 00141 double compute_epschr4b (); 00142 00143 //pomocne pretvoreni betonu 00144 double compute_epschrk1a (); 00145 00146 //hodnota ohyboveho momentu, kdy dochazi k teceni vyztuze As2, ale netece As1 a ani nedochazi k drceni betonu, x > h 00147 double compute_mhr4 (); 00148 00149 //hodnota ohyboveho momentu, kdy uz tece vyztuz As2, netece As1 a ani nedochazi k drceni betonu, ale dochazi ke zmene tlacene plochy x = h 00150 double compute_mhr4b (); 00151 00152 //pretvoreni betonu, ktere odpovida max. dosazenemu ohybovemu momentu, kdy uz tece vyztuz As2, netece As1 a dochazi k drceni betonu, x > h 00153 double compute_epscmax1p (); 00154 00155 //max. hodnota ohyboveho momentu, kdy uz tece vyztuz As2, netece As1 a dochazi k drceni betonu, x > h 00156 double compute_mmax1a (); 00157 00158 //hranicni ohybovy moment 00159 double compute_mh1 (); 00160 00161 //hranicni pretvoreni betonu 00162 double compute_epsch1 (); 00163 00164 //hodnota ohyboveho momentu pro pretvoreni na polovine danneho intervalu, x > h 00165 double compute_mh1hrk1a (); 00166 00167 //pomocna hodnota ohyboveho momentu 00168 double compute_mhrk1a (); 00169 00170 //hranicni ohybovy moment 00171 double compute_mh3 (); 00172 00173 //hranicni pretvoreni betonu 00174 double compute_epsch3 (); 00175 00176 //hranicni ohybovy moment 00177 double compute_mh2 (); 00178 00179 //hranicni pretvoreni betonu 00180 double compute_epsch2 (); 00181 00182 //pomocne pretvoreni betonu 00183 double compute_epscder2 (); 00184 00185 //hodnota maximalniho ohyboveho momentu, kdy uz tece vyztuz As2, netece As1 a dochazi k drceni betonu, x > h 00186 double compute_mmax1 (); 00187 00188 //pomocna hodnota ohyboveho momentu 00189 double compute_mhrk1 (); 00190 00191 //pomocne pretvoreni betonu 00192 double compute_epschrk1 (); 00193 00194 //hodnota ohyboveho momentu pro pretvoreni na polovine danneho intervalu, x > h 00195 double compute_mh1hrk1 (); 00196 00197 //TLACENA VYZTUZ + NORMALOVA SILA + TECENI V As1 + TECENI V As2 + x < h 00198 00199 //PRVNI TEKLO As1 00200 00201 //pretvoreni betonu, kdy dochazi k teceni vyztuze As2, uz teklo As1 a nedochazi k drceni betonu, x < h 00202 double compute_epschr5 (); 00203 00204 //pretvoreni betonu po dosazeni mezniho pretvoreni oceli As1, kdy jeste nedoslo k teceni vyztuze As2 a nedochazi k drceni betonu, x < h 00205 double compute_epscs11 (); 00206 00207 //hodnota ohyboveho momentu, kdy dochazi k teceni vyztuze As2, uz teklo As1 a nedochazi k drceni betonu, x < h 00208 double compute_mhr5 (); 00209 00210 //hodnota ohyboveho momentu po dosazeni mezniho pretvoreni betonu v tlaku, tece As1 a jeste nedoslo k teceni vyztuze As2 00211 double compute_mmax2 (); 00212 00213 //hodnota ohyboveho momentu po dosazeni mezniho pretvoreni oceli As1 pri jejim tazeni, kdy jeste nedoslo k teceni vyztuze As2 a nedochazi k drceni betonu, x < h 00214 double compute_ms11 (); 00215 00216 //pomocna hodnota ohyboveho momentu 00217 double compute_mh5a (); 00218 00219 //pomocne pretvoreni betonu 00220 double compute_epsch5a (); 00221 00222 //pomocny hranicni ohybovy moment 00223 double compute_mh5p (); 00224 00225 //pomocne hranicni pretvoreni betonu 00226 double compute_epsch5p (); 00227 00228 //hodnota ohyboveho momentu pro pretvoreni na polovine danneho intervalu, x < h 00229 double compute_mh25a (); 00230 00231 //mezni pretvorení betonu v tlaku po dosazeni mezni únosnosti pruezu, která se blizí meznímu pretvoreni oceli As1 pri jejim tazeni, kdy jeste nedoslo k teceni vyztuze As2 a nedochazi k drceni betonu, x < h 00232 double compute_epscder4 (); 00233 00234 //hodnota ohyboveho momentu po dosazeni mezni únosnosti pruezu, která se blizí meznímu pretvoreni oceli As1 pri jejim tazeni, kdy jeste nedoslo k teceni vyztuze As2 a nedochazi k drceni betonu, x < h 00235 double compute_mmax5 (); 00236 00237 //hranicni ohybovy moment 00238 double compute_mh5 (); 00239 00240 //hranicni pretvoreni betonu 00241 double compute_epsch5 (); 00242 00243 //PRVNI TEKLO As2 00244 00245 //pomocne pretvoreni betonu 00246 double compute_epscp (); 00247 00248 //pretvoreni betonu v tlaku, kdy dochazi k teceni vyztuze As1, uz teklo As2 a nedochazi k drceni betonu, x < h 00249 double compute_epschr6 (); 00250 00251 //hodnota ohyboveho momentu, kdy dochazi k teceni vyztuze As1, uz teklo As2 a nedochazi k drceni betonu, x < h 00252 double compute_mhr6 (); 00253 00254 //hodnota ohyboveho momentu po dosazeni mezniho pretvoreni betonu v tlaku, tece As2 a jeste nedoslo k teceni vyztuze As1, x < h 00255 double compute_mmax3 (); 00256 00257 //pomocna hodnota ohyboveho momentu 00258 double compute_mh6a (); 00259 00260 //pomocna hodnota pretvoreni betonu v tlaku 00261 double compute_epsch6a (); 00262 00263 //hranicni ohybovy moment 00264 double compute_mh7 (); 00265 00266 //hranicni pretvoreni betonu v tlaku 00267 double compute_epsch7 (); 00268 00269 //pretvoreni betonu v tlaku na polovine danneho intervalu 00270 double compute_epsch12 (); 00271 00272 //pretvorení vyztuze As1, po dosazeni pretvorení betonu v tlaku epsch1 00273 double compute_epss1h1 (); 00274 00275 //pretvorení vyztuze As1, po dosazeni pretvorení betonu v tlaku epsch2 00276 double compute_epss1h2 (); 00277 00278 //hodnota ohyboveho momentu pro pretvoreni na polovine danneho intervalu 00279 double compute_mh01 (); 00280 00281 //hodnota ohyboveho momentu pro pretvoreni na polovine danneho intervalu 00282 double compute_mh12 (); 00283 00284 //hodnota ohyboveho momentu pro pretvoreni na polovine danneho intervalu 00285 double compute_mh25 (); 00286 00287 //hodnota ohyboveho momentu pro pretvoreni na polovine danneho intervalu 00288 double compute_mh36a (); 00289 00290 //pomocne pretvoreni betonu 00291 double compute_epscder1 (); 00292 00293 //hranicni ohybovy moment 00294 double compute_mh6 (); 00295 00296 //hranicni pretvoreni betonu v tlaku 00297 double compute_epsch6 (); 00298 00299 //hodnota ohyboveho momentu pro pretvoreni na polovine danneho intervalu 00300 double compute_mh36 (); 00301 00302 //maximalni ohybovy moment, tece As2 a jeste nedoslo k teceni vyztuze As1 a k dosazeni mezniho pretvoreni betonu, x < h 00303 double compute_mmax4 (); 00304 00305 //hodnota ohyboveho momentu pro pretvoreni na polovine danneho intervalu 00306 double compute_mh7maxa (); 00307 00308 //maximalní mozny ohybovy moment dosazeny na dannem elementu pri danne normalove sile, bez zohlednení epscder3 00309 double compute_mhmaxa (); 00310 00311 //pomocna hodnota pretvoreni betonu v tlaku 00312 double compute_epscder3 (); 00313 00314 //maximalní mozny ohybovy moment dosazeny na dannem elementu pri danne normalove sile, pokud epscder3 lezí v danném definicním oboru 00315 double compute_mhmax (); 00316 00317 //pomocna hodnota ohyboveho momentu 00318 double compute_mh8p (); 00319 00320 //pomocna hodnota pretvoreni betonu v tlaku 00321 double compute_epsch8p (); 00322 00323 //pomocna hodnota ohyboveho momentu 00324 double compute_mh8 (); 00325 00326 //pomocna hodnota pretvoreni betonu v tlaku 00327 double compute_epsch8 (); 00328 00329 //hodnota ohyboveho momentu pro pretvoreni na polovine danneho intervalu 00330 double compute_mh78 (); 00331 00332 //ohybovy moment pri kterem uz nespolupusobi beton a dochazi k teceni vyztuze As2, netece As1 00333 double compute_mhoc1 (); 00334 00335 //ohybovy moment pri kterem uz nespolupusobi beton a dochazi k meznimu teceni vyztuze As2, netece As1 00336 double compute_mmaxoc1 (); 00337 00338 //pretvoreni oceli As1 pri kterem uz nespolupusobi beton a dochazi k meznimu teceni vyztuze As2, netece As1 00339 double compute_epss1oc1 (); 00340 00341 //ohybovy moment pri kterem uz nespolupusobi beton a dochazi k teceni vyztuze As1 a uz teklo As2 00342 double compute_mhoc2 (); 00343 00344 //pretvoreni oceli As1 pri kterem uz nespolupusobi beton a dochazi k teceni vyztuze As1 a uz teklo As2 00345 double compute_epss1oc2 (); 00346 00347 //ohybovy moment pri kterem uz nespolupusobi beton a dochazi k meznimu teceni vyztuze As1 a prekroceni mezniho pretvoreni As2 00348 double compute_mmaxoc2 (); 00349 00350 //hodnota ohyboveho momentu pro prvni bod zlomu na pracovnim diagramu prurezu 00351 double compute_mht1 (); 00352 00353 //hodnota ohyboveho momentu pro druhy bod zlomu na pracovnim diagramu prurezu 00354 double compute_mht2 (); 00355 00356 //hodnota ohyboveho momentu pro treti bod zlomu na pracovnim diagramu prurezu 00357 double compute_mht3 (); 00358 00359 //hodnota ohyboveho momentu pro ctvrty bod zlomu na pracovnim diagramu prurezu 00360 double compute_mht4 (); 00361 00362 //maximalni hodnota ohyboveho momentu na pracovnim diagramu prurezu 00363 double compute_mmax (); 00364 00365 //hodnota pretvoreni betonu v tlaku pro prvni bod zlomu na pracovnim diagramu prurezu 00366 double compute_epscht1 (); 00367 00368 //hodnota pretvoreni betonu v tlaku pro druhy bod zlomu na pracovnim diagramu prurezu 00369 double compute_epscht2 (); 00370 00371 //hodnota pretvoreni betonu v tlaku pro treti bod zlomu na pracovnim diagramu prurezu 00372 double compute_epscht3 (); 00373 00374 //hodnota pretvoreni betonu v tlaku pro ctvrty bod zlomu na pracovnim diagramu prurezu 00375 double compute_epscht4 (); 00376 00377 //maximalni hodnota pretvoreni betonu v tlaku na pracovnim diagramu prurezu 00378 double compute_epscmax (); 00379 00380 //absolutni uvazovana hodnota ohyboveho momentu na dannem elementu v polovine intervalu 00381 void compute_mp (); 00382 00383 //funkce ktera rusi vliv betonove casti prurezu 00384 void compute_tcn (); 00385 00386 //hodnota pretvoreni betonu v tlaku odpovidajici i-temu vnitrnimu uzlu 00387 void compute_epscsp (); 00388 00389 //hodnota pretvoreni oceli As1 v tahu, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00390 void compute_epss1sp (); 00391 00392 //hodnota prvni konstanty pro vypocet ohyboveho momentu v zavislosti na pretvoreni betonu v tlaku, odpovidajici i-temu vnitrnimu uzlu 00393 void compute_a1sp (); 00394 00395 //hodnota druhe pomocne konstanty pro vypocet ohyboveho momentu v zavislosti na pretvoreni betonu v tlaku, odpovidajici i-temu vnitrnimu uzlu 00396 void compute_a2spp (); 00397 00398 //hodnota druhe konstanty pro vypocet ohyboveho momentu v zavislosti na pretvoreni betonu v tlaku, odpovidajici i-temu vnitrnimu uzlu 00399 void compute_a2sp (); 00400 00401 //hodnota prvni konstanty pro vypocet pretvoreni oceli As1 v tahu, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00402 void compute_b1sp (); 00403 00404 //hodnota druhe pomocne konstanty pro vypocet pretvoreni oceli As1 v tahu, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00405 void compute_b2spp (); 00406 00407 //hodnota druhe konstanty pro vypocet ohyboveho momntu v zavislosti na pretvoreni betonu v tlaku, odpovidajici i-temu vnitrnimu uzlu 00408 void compute_b2sp (); 00409 00410 //hodnota ohyboveho momentu pro vypocet prvni hodnoty funkce prace vnitrnich sil,odpovidajici i-temu vnitrnimu uzlu 00411 void compute_mr (); 00412 00413 //hodnota pretvoreni betonu v tlaku pro vypocet prvni hodnoty funkce prace vnitrnich sil, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00414 void compute_epscr (); 00415 00416 //hodnota pretvoreni oceli As1 v tahu pro vypocet prvni hodnoty funkce prace vnitrnich sil, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00417 void compute_epss1r (); 00418 00419 //hodnota pretvoreni oceli As1 v tahu pro vypocet prvni hodnoty funkce pomocne prace vnitrnich sil, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00420 void compute_epss1spp (); 00421 00422 //velikost pomocne tlacene plochy prurezu 00423 void compute_xr (); 00424 00425 //velikost tlacene plochy prurezu 00426 void compute_xsp (); 00427 00428 //prvni integracni hodnota udavajici pocatek funkce prace vnitrnich sil pripadajici betonu, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00429 void compute_ir1 (); 00430 00431 //druha integracni hodnota funkce prace vnitrnich sil pripadajici betonu, od ucinku ohyboveho momentu,odpovidajici i-temu vnitrnimu uzlu 00432 void compute_ims1 (); 00433 00434 //prvni integracni hodnota udavajici pocatek funkce prace vnitrnich sil pripadajici tazene oceli As1, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00435 void compute_ir2 (); 00436 00437 //druha integracni hodnota funkce prace vnitrnich sil pripadajici tazene oceli As1, od ucinku ohyboveho momentu,odpovidajici i-temu vnitrnimu uzlu 00438 void compute_ims2 (); 00439 00440 //prvni integracni hodnota udavajici pocatek funkce prace vnitrnich sil pripadajici tlacene oceli As2, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00441 void compute_ir3 (); 00442 00443 //druha integracni hodnota prace vnitrnich sil pripadajici tlacene oceli As2, od ucinku ohyboveho momentu,odpovidajici i-temu vnitrnimu uzlu 00444 void compute_ims3 (); 00445 00446 //prvni integracni hodnota udavajici pocatek funkce prace vnitrnich sil, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00447 void compute_ic1 (); 00448 00449 //druha integracni hodnota prace vnitrnich sil, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00450 void compute_ic2 (); 00451 00452 //pomocna integracni hodnota prace vnitrnich sil, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00453 void compute_icp (); 00454 00455 //prvni clen integracni funkce prace vnitrnich sil, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00456 void compute_c1 (); 00457 00458 //druhy clen integracni funkce prace vnitrnich sil, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00459 void compute_c2 (); 00460 00461 //prvni clen pomocne integracni funkce prace vnitrnich sil, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00462 void compute_c1p (); 00463 00464 //druhy clen pomocne integracni funkce prace vnitrnich sil, od ucinku ohyboveho momentu, odpovidajici i-temu vnitrnimu uzlu 00465 void compute_c2p (); 00466 00467 //integral prace vnitrnich sil pres dany element 00468 void compute_kih (double lp); 00469 00470 //integral pomocne prace vnitrnich sil pres dany element 00471 void compute_kid (double lp); 00472 00473 //pomer prace vnitrnich sil a pomocne prace vnitrnich sil na dannem elementu 00474 void compute_ki (); 00475 00476 //clen matice tuhosti 00477 double compute_k11 (double lp); 00478 00479 //clen matice tuhosti 00480 double compute_k12 (double lp); 00481 00482 //clen matice tuhosti 00483 double compute_k13 (double lp); 00484 00485 //clen matice tuhosti 00486 double compute_k14 (double lp); 00487 00488 //zatizeni od ucinku rovnomerneho spojiteho zatizeni (prava strana); 00489 double compute_f1 (double lp); 00490 00491 //clen matice tuhosti 00492 double compute_k21 (double lp); 00493 00494 //clen matice tuhosti 00495 double compute_k22 (double lp); 00496 00497 //clen matice tuhosti 00498 double compute_k23 (double lp); 00499 00500 //clen matice tuhosti 00501 double compute_k24 (double lp); 00502 00503 //zatizeni od ucinku rovnomerneho spojiteho zatizeni (prava strana); 00504 double compute_f2 (double lp); 00505 00506 //clen matice tuhosti 00507 double compute_k31 (double lp); 00508 00509 //clen matice tuhosti 00510 double compute_k32 (double lp); 00511 00512 //clen matice tuhosti 00513 double compute_k33 (double lp); 00514 00515 //clen matice tuhosti 00516 double compute_k34 (double lp); 00517 00518 //zatizeni od ucinku rovnomerneho spojiteho zatizeni (prava strana); 00519 double compute_f3 (double lp); 00520 00521 //clen matice tuhosti 00522 double compute_k41 (double lp); 00523 00524 //clen matice tuhosti 00525 double compute_k42 (double lp); 00526 00527 //clen matice tuhosti 00528 double compute_k43 (double lp); 00529 00530 //clen matice tuhosti 00531 double compute_k44 (double lp); 00532 00533 //zatizeni od ucinku rovnomerneho spojiteho zatizeni (prava strana); 00534 double compute_f4 (double lp); 00535 00536 //matice tuhosti pro stanoveni posouvajici sily a ohyboveho momentu na elementu delky lp 00537 void compute_kp (double lp); 00538 00539 //prava strana (+zatizeni od ucinku spojiteho zatizeni); pro stanoveni posouvajici sily a ohyboveho momentu na elementu delky lp 00540 void compute_fkp (double lp); 00541 00542 //matice tuhosti pro stanoveni normalove sily na elementu delky lp 00543 void compute_kn (double lp); 00544 00545 //prava strana (+zatizeni od ucinku spojiteho zatizeni); pro stanoveni normalove sily na elementu delky lp 00546 void compute_fkn (double lp); 00547 00548 00549 void compute_all (long ipp,double lp); 00550 00551 00552 void stiffness_matrix (matrix &sm,long ipp,double lp); 00553 00554 void internal_forces (vector &ifor,long ipp,double lp); 00555 00556 00557 00558 00559 00560 // charakteristicka pevnost betonu v tlaku v [kPa] např. fck = -25000 kPa 00561 double fck; 00562 00563 // pevnost oceli v tahu v [kPa] např. fy = 490000 kPa 00564 double fy; 00565 00566 // pevnost oceli v tahu v [kPa] odpovidajici pretvorevi 0.05 např. ft= 550000 kPa 00567 double ft; 00568 00569 // pretvoreni betonu v tlaku jemuz odpovida max. hodnota napeti fck napr. epsc1 = -0.0022 00570 double epsc1; 00571 00572 // mezni pretvoreni betonu v tlaku napr. epscm = -0.0035 00573 double epscm; 00574 00575 // mezni pretvoreni oceli v tahu napr. epss1lim = 0.01 00576 double epss1lim; 00577 00578 // plocha vyztuze umistene u hornich vlaken prurezu v [m2] napr. ash = 0.0032168 m2 (konvence dle deformacni metody) 00579 double ash; 00580 00581 // plocha vyztuze umistene u dolnich vlaken prurezu v [m2] napr. asd = 0.0032168 m2 (konvence dle deformacni metody) 00582 double asd; 00583 00584 // kryti vyztuze umistene u hornich vlaken prurezu + polomer vyztuze v [m] napr. ch = 0.041 m (konvence dle deformacni metody) 00585 double ch; 00586 00587 // kryti vyztuze umistene u dolnich vlaken prurezu + polomer vyztuze v [m] napr. cd = 0.041 m (konvence dle deformacni metody) 00588 double cd; 00589 00590 // modul pruznosti oceli pred plastickym stavem v [kPa] napr. e0c = 200000000 kPa 00591 double e0c; 00592 00593 // vyska prurezu v [m] napr. h = 0.6 m 00594 double h; 00595 00596 // sirka prurezu v [m] napr. b = 0.3 m 00597 double b; 00598 00599 // hodnota rovnomerneho spojiteho zatizeni kolmo na strednici prutu v [kN/m] (u pruvlaku kladne hodnoty ve smeru gravitace, u sloupu z leva do prava) 00600 double fp; 00601 00602 // hodnota rovnomerneho spojiteho zatizeni ve smeru strednice prutu v [kN/m] (u sloupu kladne hodnoty ve smeru gravitace u pruvlaku z prava do leva) 00603 double fn; 00604 00605 00606 00607 double fc,ecnom,epsuk,e0cn,k,a,d,as1,as2,z1,z2,z3p,z3,n0,epsn,nmax,an1,an2,in1,in2,kinh,kind,kin; 00608 00609 double epschr1,mhr1,epschr2,mhr2,epschr3,mhr3,epschr4,epscmax1,epschr4b,epschrk1a,mhr4,mhr4b,epscmax1p,mmax1a,mh1,epsch1; 00610 00611 double mh1hrk1a,mhrk1a,mh3,epsch3,mh2,epsch2,epscder2,mmax1,mhrk1,epschrk1,mh1hrk1,epschr5,epscs11,mhr5,mmax2,ms11; 00612 00613 double mh5a,epsch5a,mh5p,epsch5p,mh25a,epscder4,mmax5,mh5,epsch5,epscp,epschr6,mhr6,mmax3,mh6a,epsch6a,mh7,epsch7,epsch12; 00614 00615 double epss1h1,epss1h2,mh01,mh12,mh25,mh36a,epscder1,mh6,epsch6,mh36,mmax4,mh7maxa,mhmaxa,epscder3,mhmax,mh8p,epsch8p,mh8; 00616 00617 double epsch8,mh78,mhoc1,mmaxoc1,epss1oc1,mhoc2,epss1oc2,mmaxoc2; 00618 00619 double mht1,mht2,mht3,mht4,mmax,epscht1,epscht2,epscht3,epscht4,epscmax; 00620 00621 double *m,*mpp,*tcv; 00622 00623 double *mp,*tcn,*epscsp,*epss1,*epss1sp,*a1sp,*a2spp,*a2sp,*b1sp,*b2spp,*b2sp,*mr,*epscr,*epss1r,*epss1spp,*xr,*xsp; 00624 00625 double *ir1,*ims1,*ims2,*ir2,*ir3,*ims3,*ic1,*ic2,*icp,*c1,*c2,*c1p,*c2p,*kih,*kid,*ki; 00626 00627 double k11,k12,k13,k14,f1,k21,k22,k23,k24,f2,k31,k32,k33,k34,f3,k41,k42,k43,k44,f4; 00628 00629 00630 vector fkn,fkp; 00631 00632 matrix kp,kn; 00633 00634 // pocet iteraci 00635 long iter; 00636 00637 double zero; 00638 }; 00639 00640 #endif