Source code for streng.codes.greek.kanepe.raw.ch7.rotation

# import math


[docs]def θycalc(φy, Ls, av, z, h, db, fy, fc, units='m-MPa'): """ Στροφή διαρροής .. highlight:: python .. code-block:: python :linenos: theta = rotation.θy(φy=0.00638, Ls=2.525, av=1, z=0.506, h=0.506, db=0.018, fy=500., fc=16, units='m-MPa') Args: φy (float): Καμπυλότητα διαρροής [m-1] Ls (float): Μήκος διάτμησης [m] av (float): συντελεστής που ισούται με 1 εάν η τέμνουσα VRc που προκαλεί λοξή ρηγμάτωση του στοιχείου, υπολείπεται της τιμής της τέμνουσας κατά την καμπτική διαρροή VMu=My/Ls, και με 0 αν είναι μεγαλύτερη z (float): Ο μοχλοβραχίονας των εσωτερικών δυνάμεων h (float): το ύψος της διατομής [m] db (float): η διάμετρος των διαμήκων ράβδων [m] fy (float): το όριο διαρροής του χάλυβα [MPa] fc (float): η θλιπτική αντοχή του σκυροδέματος [MPa] Returns: float: Given using the expression: .. math:: θ_y=φ_y\dfrac{L_s+a_vz}{3}+0.0014 \Big(1+1.5\dfrac{h}{L_s} \Big)+\dfrac{φ_yd_bf_y}{8\sqrt{f_c}} """ return φy * (Ls + av * z) / 3 + 0.0014 * (1 + 1.5 * h / Ls) + φy * db * fy / (8 * fc**0.5)
[docs]def θum(ν, ωtot, ω2, αs, α, ρs, ρd, fc, fyw, units='MPa'): """ Στροφή αστοχίας θum .. highlight:: python .. code-block:: python :linenos: theta_u = rotation.θum(ν=0.093, ωtot=0.351, ω2=0.132, αs=3.78, α=0.4737, ρs=0.00279, ρd=0., fc=16., fyw=500., units='MPa') Args: ν (float): ανηγμένο αξονικό ωtot (float): Συνολικό μηχανικό ποσοστό οπλισμού ω2 (float): Μηχανικό ποσοστό θλιβόμενου οπλισμού αs (float): Λόγος διάτμησης ρs (float): Γεωμετρικό ποσοστό εγκάρσιου οπλισμού ρd (float): Γεωμετρικό ποσοστό δισδιαγώνιου οπλισμού fc (float): η θλιπτική αντοχή του σκυροδέματος [MPa] fyw (float): το όριο διαρροής του χάλυβα ου εγκάρσιου οπλισμού [MPa] Returns: float: Given using the expression: .. math:: θ_{um}=0.016\cdot(0.3^ν)\cdot\Bigg[\dfrac{max(0.01,ω')}{max(0.01, (ω_{tot}-ω'))} \cdot f_c\Bigg]^{0.225}\cdot (α_s)^{0.35} \cdot 25^{(α\cdot ρ_s \cdot \dfrac{f_{yw}}{f_c})}\cdot (1.25^{100ρ_d}) """ part1 = (fc * max(0.01, ω2) / max(0.01, (ωtot-ω2)))**0.225 part2 = αs**0.35 part3 = 25.0**(α * ρs * fyw / fc) part4 = 1.25**(100 * ρd) return 0.016 * 0.3**ν * part1 * part2 * part3 * part4
[docs]def αcalc(sh, bc, hc, Σbi2): α = (1-sh/(2*bc))*(1-sh/(2*hc))*(1-Σbi2/(6*bc*hc)) return max(0.0, α)