ΑΡΧΙΤΕΚΤΟΝΙΚΗ ΚΑΙ ΑΝΤΙΣΕΙΣΜΙΚΗ ΘΩΡΑΚΙΣΗ ΕΡΕΥΝΑ

[seismic1]

Γεια σου ρε Γιαννάκη παίχτουρα...επιστημονικέ παίχτουρα...

Από ότι φαίνεται εκ των πραγμάτων εσείς είσαστε παίχτουρες γιατί....το έχετε ξεχειλώσει το πουλί σας.
Αν και το ξέρετε δεν παραδέχεστε ότι η πάκτωση και η προ ένταση αυξάνουν την σταθερότητα και την δυναμική των κατασκευών.
Είναι σαν να μην παραδέχεστε ότι η γη γυρίζει, και όταν η γη κουνιέται και εμείς κουνιόμαστε.

πακτώθηκε η μαλακία στο μυαλό σου και την κάνεις προένταση στους άλλους...

άστοχος...

Θα δούμε την δυναμική προσομοίωση.
# Ενημερωμένος κώδικας για την προσομοίωση με πάκτωση και προένταση

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp

# Σταθερές του προβλήματος
m = 155000 # kg (μάζα φρεατίου)
k = 3e6 # N/m (δυσκαμψία)
c = 5000 # Ns/m (απόσβεση)
g = 9.81 # m/s² (επιτάχυνση βαρύτητας)

# Γεωμετρικά χαρακτηριστικά
h = 10 # ύψος φρεατίου (m)
b = 3 # πλάτος φρεατίου (m)
L1 = 3 # μήκος φρεατίου (m)
t = 0.40 # πάχος τοιχίων (m)

b2 = 4 # πλάτος βάσης (m)
L = 4 # μήκος βάσης (m)
h2 = 1 # ύψος βάσης (m)

# Ιδιότητες τενόντων
num_tendons = 6
prestress_force = 200 * 1e3 # 200 τόνοι σε Newton (N)
anchorage_force = 1000 * 1e3 # 1000 τόνοι σε Newton (N)

# Προσαρμογή δυσκαμψίας λόγω προέντασης και πάκτωσης
k_effective = k + (num_tendons * prestress_force / h)

# Σεισμική επιτάχυνση (συνθετικός σεισμός)
def seismic_acceleration(t):
return 2.41 * g * np.sin(2 * np.pi * 1 * t) # 1Hz ταλάντωση

# Εξισώσεις κίνησης (Newmark-β)
def equation_of_motion(t, y):
u, v = y # u: μετατόπιση, v: ταχύτητα
a = (seismic_acceleration(t) - c * v - k_effective * u) / m # επιτάχυνση
return [v, a]

# Χρονικά βήματα προσομοίωσης
T = 5 # συνολικός χρόνος σε sec
dt = 0.01 # χρονικό βήμα
time_span = (0, T)
initial_conditions = [0, 0] # αρχικές συνθήκες (μηδενική μετατόπιση και ταχύτητα)
time_eval = np.arange(0, T, dt)

# Αριθμητική επίλυση
solution = solve_ivp(equation_of_motion, time_span, initial_conditions, t_eval=time_eval, method='RK45')

# Αποτελέσματα
u_max = np.max(np.abs(solution.y[0])) # Μέγιστη απόκριση μετατόπισης

# Γραφική παράσταση
plt.figure(figsize=(8, 5))
plt.plot(solution.t, solution.y[0], label="Μετατόπιση (m)", color='b')
plt.xlabel("Χρόνος (s)")
plt.ylabel("Μετατόπιση (m)")
plt.title("Δυναμική Απόκριση Κατασκευής με Προένταση και Πάκτωση")
plt.legend()
plt.grid()
plt.show()

# Εκτύπωση αποτελέσματος
u_max

[Bill Hicks]

import nump(t)y as np

fixed

[taxalata xalasa]

Από ότι φαίνεται εκ των πραγμάτων εσείς είσαστε παίχτουρες γιατί....το έχετε ξεχειλώσει το πουλί σας.
Αν και το ξέρετε δεν παραδέχεστε ότι η πάκτωση και η προ ένταση αυξάνουν την σταθερότητα και την δυναμική των κατασκευών.
Είναι σαν να μην παραδέχεστε ότι η γη γυρίζει, και όταν η γη κουνιέται και εμείς κουνιόμαστε.

πακτώθηκε η μαλακία στο μυαλό σου και την κάνεις προένταση στους άλλους...

άστοχος...

Θα δούμε την δυναμική προσομοίωση.

Spoiler

# Ενημερωμένος κώδικας για την προσομοίωση με πάκτωση και προένταση

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp

# Σταθερές του προβλήματος
m = 155000 # kg (μάζα φρεατίου)
k = 3e6 # N/m (δυσκαμψία)
c = 5000 # Ns/m (απόσβεση)
g = 9.81 # m/s² (επιτάχυνση βαρύτητας)

# Γεωμετρικά χαρακτηριστικά
h = 10 # ύψος φρεατίου (m)
b = 3 # πλάτος φρεατίου (m)
L1 = 3 # μήκος φρεατίου (m)
t = 0.40 # πάχος τοιχίων (m)

b2 = 4 # πλάτος βάσης (m)
L = 4 # μήκος βάσης (m)
h2 = 1 # ύψος βάσης (m)

# Ιδιότητες τενόντων
num_tendons = 6
prestress_force = 200 * 1e3 # 200 τόνοι σε Newton (N)
anchorage_force = 1000 * 1e3 # 1000 τόνοι σε Newton (N)

# Προσαρμογή δυσκαμψίας λόγω προέντασης και πάκτωσης
k_effective = k + (num_tendons * prestress_force / h)

# Σεισμική επιτάχυνση (συνθετικός σεισμός)
def seismic_acceleration(t):
return 2.41 * g * np.sin(2 * np.pi * 1 * t) # 1Hz ταλάντωση

# Εξισώσεις κίνησης (Newmark-β)
def equation_of_motion(t, y):
u, v = y # u: μετατόπιση, v: ταχύτητα
a = (seismic_acceleration(t) - c * v - k_effective * u) / m # επιτάχυνση
return [v, a]

# Χρονικά βήματα προσομοίωσης
T = 5 # συνολικός χρόνος σε sec
dt = 0.01 # χρονικό βήμα
time_span = (0, T)
initial_conditions = [0, 0] # αρχικές συνθήκες (μηδενική μετατόπιση και ταχύτητα)
time_eval = np.arange(0, T, dt)

# Αριθμητική επίλυση
solution = solve_ivp(equation_of_motion, time_span, initial_conditions, t_eval=time_eval, method='RK45')

# Αποτελέσματα
u_max = np.max(np.abs(solution.y[0])) # Μέγιστη απόκριση μετατόπισης

# Γραφική παράσταση
plt.figure(figsize=(8, 5))
plt.plot(solution.t, solution.y[0], label="Μετατόπιση (m)", color='b')
plt.xlabel("Χρόνος (s)")
plt.ylabel("Μετατόπιση (m)")
plt.title("Δυναμική Απόκριση Κατασκευής με Προένταση και Πάκτωση")
plt.legend()
plt.grid()
plt.show()

# Εκτύπωση αποτελέσματος
u_max

τι θα δουμε ρε τσιμεντοκονιάρη;
μονος σου τα αντιγραφεις μονος σου τα βλέψ...

δεν ασχολειται κάποιος άλλος με την εγκεφαλική σου μαλακία...

σε έχουν στείλει αδιάβαστο κι ακόμα δεν εχεις χαμπεριάσει πως εισαι zombie...

[Ασέβαστος]

Αν και το ξέρετε δεν παραδέχεστε ότι η πάκτωση

σταματα να χρησιμοποιεις λεξεις που δεν ξερεις τι σημαινουν.

[seismic1]

πακτώθηκε η μαλακία στο μυαλό σου και την κάνεις προένταση στους άλλους...

άστοχος...

Θα δούμε την δυναμική προσομοίωση.

Spoiler

# Ενημερωμένος κώδικας για την προσομοίωση με πάκτωση και προένταση

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp

# Σταθερές του προβλήματος
m = 155000 # kg (μάζα φρεατίου)
k = 3e6 # N/m (δυσκαμψία)
c = 5000 # Ns/m (απόσβεση)
g = 9.81 # m/s² (επιτάχυνση βαρύτητας)

# Γεωμετρικά χαρακτηριστικά
h = 10 # ύψος φρεατίου (m)
b = 3 # πλάτος φρεατίου (m)
L1 = 3 # μήκος φρεατίου (m)
t = 0.40 # πάχος τοιχίων (m)

b2 = 4 # πλάτος βάσης (m)
L = 4 # μήκος βάσης (m)
h2 = 1 # ύψος βάσης (m)

# Ιδιότητες τενόντων
num_tendons = 6
prestress_force = 200 * 1e3 # 200 τόνοι σε Newton (N)
anchorage_force = 1000 * 1e3 # 1000 τόνοι σε Newton (N)

# Προσαρμογή δυσκαμψίας λόγω προέντασης και πάκτωσης
k_effective = k + (num_tendons * prestress_force / h)

# Σεισμική επιτάχυνση (συνθετικός σεισμός)
def seismic_acceleration(t):
return 2.41 * g * np.sin(2 * np.pi * 1 * t) # 1Hz ταλάντωση

# Εξισώσεις κίνησης (Newmark-β)
def equation_of_motion(t, y):
u, v = y # u: μετατόπιση, v: ταχύτητα
a = (seismic_acceleration(t) - c * v - k_effective * u) / m # επιτάχυνση
return [v, a]

# Χρονικά βήματα προσομοίωσης
T = 5 # συνολικός χρόνος σε sec
dt = 0.01 # χρονικό βήμα
time_span = (0, T)
initial_conditions = [0, 0] # αρχικές συνθήκες (μηδενική μετατόπιση και ταχύτητα)
time_eval = np.arange(0, T, dt)

# Αριθμητική επίλυση
solution = solve_ivp(equation_of_motion, time_span, initial_conditions, t_eval=time_eval, method='RK45')

# Αποτελέσματα
u_max = np.max(np.abs(solution.y[0])) # Μέγιστη απόκριση μετατόπισης

# Γραφική παράσταση
plt.figure(figsize=(8, 5))
plt.plot(solution.t, solution.y[0], label="Μετατόπιση (m)", color='b')
plt.xlabel("Χρόνος (s)")
plt.ylabel("Μετατόπιση (m)")
plt.title("Δυναμική Απόκριση Κατασκευής με Προένταση και Πάκτωση")
plt.legend()
plt.grid()
plt.show()

# Εκτύπωση αποτελέσματος
u_max

τι θα δουμε ρε τσιμεντοκονιάρη;
μονος σου τα αντιγραφεις μονος σου τα βλέψ...

δεν ασχολειται κάποιος άλλος με την εγκεφαλική σου μαλακία...

σε έχουν στείλει αδιάβαστο κι ακόμα δεν εχεις χαμπεριάσει πως εισαι zombie...

Δεν περιμένω μπράβο από εσάς.
Να σας πηδήξω θέλω.

[seismic1]

Αν και το ξέρετε δεν παραδέχεστε ότι η πάκτωση

σταματα να χρησιμοποιεις λεξεις που δεν ξερεις τι σημαινουν.

Όταν σου τον χώσω θα καταλάβεις τι είναι πάκτωση.

[Στύξ]

Αν και το ξέρετε δεν παραδέχεστε ότι η πάκτωση

σταματα να χρησιμοποιεις λεξεις που δεν ξερεις τι σημαινουν.

Όταν σου τον χώσω θα καταλάβεις τι είναι πάκτωση.

ΗΜΑΡΤΟΝ

[taxalata xalasa]

Θα δούμε την δυναμική προσομοίωση.

Spoiler

# Ενημερωμένος κώδικας για την προσομοίωση με πάκτωση και προένταση

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp

# Σταθερές του προβλήματος
m = 155000 # kg (μάζα φρεατίου)
k = 3e6 # N/m (δυσκαμψία)
c = 5000 # Ns/m (απόσβεση)
g = 9.81 # m/s² (επιτάχυνση βαρύτητας)

# Γεωμετρικά χαρακτηριστικά
h = 10 # ύψος φρεατίου (m)
b = 3 # πλάτος φρεατίου (m)
L1 = 3 # μήκος φρεατίου (m)
t = 0.40 # πάχος τοιχίων (m)

b2 = 4 # πλάτος βάσης (m)
L = 4 # μήκος βάσης (m)
h2 = 1 # ύψος βάσης (m)

# Ιδιότητες τενόντων
num_tendons = 6
prestress_force = 200 * 1e3 # 200 τόνοι σε Newton (N)
anchorage_force = 1000 * 1e3 # 1000 τόνοι σε Newton (N)

# Προσαρμογή δυσκαμψίας λόγω προέντασης και πάκτωσης
k_effective = k + (num_tendons * prestress_force / h)

# Σεισμική επιτάχυνση (συνθετικός σεισμός)
def seismic_acceleration(t):
return 2.41 * g * np.sin(2 * np.pi * 1 * t) # 1Hz ταλάντωση

# Εξισώσεις κίνησης (Newmark-β)
def equation_of_motion(t, y):
u, v = y # u: μετατόπιση, v: ταχύτητα
a = (seismic_acceleration(t) - c * v - k_effective * u) / m # επιτάχυνση
return [v, a]

# Χρονικά βήματα προσομοίωσης
T = 5 # συνολικός χρόνος σε sec
dt = 0.01 # χρονικό βήμα
time_span = (0, T)
initial_conditions = [0, 0] # αρχικές συνθήκες (μηδενική μετατόπιση και ταχύτητα)
time_eval = np.arange(0, T, dt)

# Αριθμητική επίλυση
solution = solve_ivp(equation_of_motion, time_span, initial_conditions, t_eval=time_eval, method='RK45')

# Αποτελέσματα
u_max = np.max(np.abs(solution.y[0])) # Μέγιστη απόκριση μετατόπισης

# Γραφική παράσταση
plt.figure(figsize=(8, 5))
plt.plot(solution.t, solution.y[0], label="Μετατόπιση (m)", color='b')
plt.xlabel("Χρόνος (s)")
plt.ylabel("Μετατόπιση (m)")
plt.title("Δυναμική Απόκριση Κατασκευής με Προένταση και Πάκτωση")
plt.legend()
plt.grid()
plt.show()

# Εκτύπωση αποτελέσματος
u_max

τι θα δουμε ρε τσιμεντοκονιάρη;
μονος σου τα αντιγραφεις μονος σου τα βλέψ...

δεν ασχολειται κάποιος άλλος με την εγκεφαλική σου μαλακία...

σε έχουν στείλει αδιάβαστο κι ακόμα δεν εχεις χαμπεριάσει πως εισαι zombie...

Δεν περιμένω μπράβο από εσάς.
Να σας πηδήξω θέλω.

θέλε...

[Λοξίας]

Γεια σου ρε Γιαννάκη παίχτουρα...επιστημονικέ παίχτουρα...

Από ότι φαίνεται εκ των πραγμάτων εσείς είσαστε παίχτουρες γιατί....το έχετε ξεχειλώσει το πουλί σας.

Γιαννάκη κακομοίρη...το ποιος έχει εχει ξεχειλώσει το πουλί του, φαίνεται εκ των πραγμάτων όντως. Και από...

Εδώ μπαίνω για να τον παίζω.

...την ίδια του την δήλωση.

Γιαννάκη μαλακοπαίχτουρα...

[Daje]

Να σας σωσουμε, δεν θελετε.
Να σας γαμησουμε, παλι δεν θελετε.
Τι στον πεουλα θελετε???

[seismic1]

Να σας σωσουμε, δεν θελετε.
Να σας γαμησουμε, παλι δεν θελετε.
Τι στον πεουλα θελετε???

Χέσιμο γερό.

[seismic1]

Άκυρο.

[seismic1]

Μια χώρα χωρίς οικονομία δεν είναι χώρα είναι χώρος.
Η οικονομία χωρίς την καινοτομία δεν είναι ισχυρή.
Η Ευρώπη με αρχηγό την Ελλάδα δεν επιτρέπει σε ιδέες σαν την δική μου να καρποφορήσουν εν ονόματι της επιχειρηματικότητας.
Ένα τσούρμο αργόσχολοι υπογράφουν διατάξεις για να σταματήσουν την επιχειρηματικότητα.
Αποτέλεσμα. Ο Τραμπ και ο Πούτιν να πετάξουν την Ευρώπη έξω από την μοιρασιά δημιουργώντας μόνοι τους
μια νέα Γιάλτα.
Υπογεννητικότητα και τα λιγοστά ποια μορφωμένα παιδιά μας μεταναστεύουν και αντικαθίστανται από αλλοδαπούς.
Αποτέλεσμα όλων αυτών οι φόροι φόροι φόροι...της παγκοσμιοποίησης.
Η υπερφορολόγηση χωρίς παραγωγή είναι σαν να προσπαθείς να αντλήσεις νερό από ένα άδειο πηγάδι. Αν δεν υπάρχουν επιχειρήσεις, καινοτομία και πραγματική οικονομική δραστηριότητα, το κράτος δεν έχει από πού να εισπράξει.

Γι’ αυτό και τα πιο επιτυχημένα κράτη δεν βασίζονται μόνο στη φορολογία, αλλά στη δημιουργία πλούτου. Δίνουν κίνητρα για παραγωγή, προσελκύουν επενδύσεις και ενισχύουν την επιχειρηματικότητα.

Η Ευρώπη όμως μοιάζει να πιστεύει ότι μπορεί να διατηρήσει το βιοτικό της επίπεδο μόνο με φόρους, χωρίς να επενδύει σοβαρά στην παραγωγή.
Θα το συνειδητοποιήσει έγκαιρα ή θα βρεθεί αντιμέτωπη με μια βαθύτερη κρίση?

[seismic1]

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp

# Σταθερές του προβλήματος
m_structure = 155000 # kg (μάζα φρεατίου)
m_soil = 200000 # kg (μάζα εδάφους που συμμετέχει ενεργά)
m_total = m_structure + m_soil # Συνολική μάζα συστήματος
k_conventional = 3e6 # N/m (δυσκαμψία συμβατικής κατασκευής)
c = 5000 # Ns/m (απόσβεση)
g = 9.81 # m/s² (επιτάχυνση βαρύτητας)

# Ενίσχυση δυσκαμψίας λόγω πάκτωσης και μετά-τάνυσης
k_fixed = 10 * k_conventional # Μεγάλη δυσκαμψία λόγω πάκτωσης
k_post_tension = 5 * k_conventional # Πρόσθετη δυσκαμψία λόγω μετά-τάνυσης
k_innovative = k_fixed + k_post_tension # Συνολική δυσκαμψία

# Σεισμική επιτάχυνση (ημιτονοειδής διέγερση)
def seismic_acceleration(t):
return 0.36 * g * np.sin(2 * np.pi * 1 * t) # 1Hz ταλάντωση, 0.36g επιτάχυνση

# Εξισώσεις κίνησης
def equation_of_motion(t, y, stiffness):
u, v = y # u: μετατόπιση, v: ταχύτητα
a = (seismic_acceleration(t) - c * v - stiffness * u) / m_total # Επιτάχυνση
return [v, a]

# Χρονικά βήματα προσομοίωσης
T = 5 # Συνολικός χρόνος σε sec
dt = 0.01 # Χρονικό βήμα
time_span = (0, T)
initial_conditions = [0, 0] # Αρχικές συνθήκες (μηδενική μετατόπιση και ταχύτητα)
time_eval = np.arange(0, T, dt)

# Αριθμητική επίλυση για συμβατική κατασκευή
solution_conventional = solve_ivp(equation_of_motion, time_span, initial_conditions, t_eval=time_eval,
args=(k_conventional,), method='RK45')

# Αριθμητική επίλυση για καινοτόμο κατασκευή (πάκτωση + μετά-τάνυση)
solution_innovative = solve_ivp(equation_of_motion, time_span, initial_conditions, t_eval=time_eval,
args=(k_innovative,), method='RK45')

# Υπολογισμός μέγιστης μετατόπισης
u_max_conventional = np.max(np.abs(solution_conventional.y[0])) * 1000 # σε mm
u_max_innovative = np.max(np.abs(solution_innovative.y[0])) * 1000 # σε mm

# Γραφική παράσταση
plt.figure(figsize=(8, 5))
plt.plot(solution_conventional.t, solution_conventional.y[0] * 1000, label="Συμβατική Κατασκευή", color='r') # σε mm
plt.plot(solution_innovative.t, solution_innovative.y[0] * 1000, label="Καινοτόμος Κατασκευή", color='b') # σε mm
plt.xlabel("Χρόνος (s)")
plt.ylabel("Μετατόπιση (mm)")
plt.title("Σύγκριση Δυναμικής Απόκρισης")
plt.legend()
plt.grid()
plt.show()

# Εκτύπωση αποτελεσμάτων
print(f"Μέγιστη μετατόπιση (συμβατική): {u_max_conventional:.2f} mm")
print(f"Μέγιστη μετατόπιση (καινοτόμος): {u_max_innovative:.2f} mm")

[taxalata xalasa]

είπαν του τρελού να χέσει κι αυτός εξεκωλιάστηκε...