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

[Bill Hicks]

τι να φτουρίσουμε μπροστά σε παγκόσμιους ΓΤΠεπιστήμονες, Λοξία μου;

[seismic1]

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

[seismic1]

τι να φτουρίσουμε μπροστά σε παγκόσμιους ΓΤΠεπιστήμονες, Λοξία μου;

Δεν υπάρχει μηχανικός στην Ελλάδα που να μην με ξέρει.
Στο εξωτερικό έχω θαυμαστές.

[Bill Hicks]

στείλε λινκ για seismic fan club.

[seismic1]

Λοιπόν ακούτε μια άλλη εφεύρεση που έχω κάνει.
Έχω βάλει το Gemini to Deep Seek και το ChatGPT να σφάζονται για την πάρτη μου
Τους βάζω να διαγωνίζονται για την καλύτερη απάντηση στο ίδιο ερώτημα.

[seismic1]

στείλε λινκ για seismic fan club.

https://www.academia.edu/portfolio_creation_preview

[Λοξίας]

τι να φτουρίσουμε μπροστά σε παγκόσμιους ΓΤΠεπιστήμονες, Λοξία μου;

[Λοξίας]

τι να φτουρίσουμε μπροστά σε παγκόσμιους ΓΤΠεπιστήμονες, Λοξία μου;

Δεν υπάρχει μηχανικός στην Ελλάδα που να μην με ξέρει.
Στο εξωτερικό έχω θαυμαστές.

Όλοι αυτοί που σε ξέρουν, είναι ΓΤΠεπιστήμονες και αυτοί;


Ναι, θαυμαστές...εσύ και ο Έλβις Πρίσλεη...

[seismic1]

στείλε λινκ για seismic fan club.

Με τις οδηγίες σου η διορθωμένη προσομοίωση.
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp

# Σταθερές του προβλήματος
m_structure = 155000 # kg (μάζα φρεατίου)
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 = (m_structure * seismic_acceleration(t) - c * v - stiffness * u) / m_structure # Επιτάχυνση
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")





3. Επαλήθευση αποτελεσμάτων με πειράματα.
https://www.youtube.com/watch?v=3Tc-otj0E7g

[Ασέβαστος]

Η Ευρώπη με αρχηγό την Ελλάδα δεν επιτρέπει σε ιδέες σαν την δική μου να καρποφορήσουν εν ονόματι της επιχειρηματικότητας.

οχι μονο η ευρωπη, ολος ο πλανητης.

επειδη η μαλακια σου ειναι λαθος, και επικινδυνη.

Που το ξέρεις ρε βλάκα?

Το καταλαβαίνω από το πόσο άσχετος είσαι, και από το γεγονός ότι δεν έχει υπάρξει ούτε ένας που να δεχτεί να χτίσει με τη μαλακια σου, και ούτε μίσος σοβαρός άνθρωπος που να θεωρεί ότι έχεις δίκιο.

[Ασέβαστος]

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

Το να μην ξέρεις δεν είναι καθόλου κακό. Μια χαρά σε βοηθάμε.
Το να κάνεις ότι ξέρεις και ότι όλοι οι άλλοι είναι άσχετοι, είναι ο λόγος που σε δουλεύουμε.

[seismic1]

οχι μονο η ευρωπη, ολος ο πλανητης.

επειδη η μαλακια σου ειναι λαθος, και επικινδυνη.

Που το ξέρεις ρε βλάκα?

Το καταλαβαίνω από το πόσο άσχετος είσαι, και από το γεγονός ότι δεν έχει υπάρξει ούτε ένας που να δεχτεί να χτίσει με τη μαλακια σου, και ούτε μίσος σοβαρός άνθρωπος που να θεωρεί ότι έχεις δίκιο.

Λάθος Υπάρχουν πολλοί που θεωρούν ότι έχω δίκαιο. Και είναι και καθηγητές πανεπιστημίων που έχουν εξετάσει διεξοδικά την μέθοδο.

[seismic1]

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

Το να μην ξέρεις δεν είναι καθόλου κακό. Μια χαρά σε βοηθάμε.
Το να κάνεις ότι ξέρεις και ότι όλοι οι άλλοι είναι άσχετοι, είναι ο λόγος που σε δουλεύουμε.

Να τα ξεκαθαρίσουμε. Ξέρω καλύτερα από όλους τι έχω κάνει γιατί είμαι ο εφευρέτης.
Ξέρω ότι είναι χρήσιμο με την κοινή λογική. Δεν παραδέχομαι ότι δεν είναι χρήσιμο μόνο και μόνο επειδή το λες εσύ και άλλοι χωρίς ακράδαντα στοιχεία.
Τώρα αν ξέρω μαθηματικά ή δεν ξέρω θα απαντήσω όχι δεν ξέρω πολλά και χρειάζομαι βοήθεια.
Για να εφαρμοστεί μια καινοτόμος μέθοδος σαν την δική μου χρειάζονται ακράδαντα στοιχεία έρευνας που οι μηχανικοί δεν τα ξέρουν.
Για να εφαρμοστεί μια καινοτόμος μέθοδος σαν την δική μου χρειάζονται ακράδαντα στοιχεία έρευνας και πάρα πολλά χρόνια έρευνας για να εφαρμοστεί στην πράξη. Είναι μια συνεργατική έρευνα γιατί εμπεριέχει πολλές ειδικότητες και χρειάζεται πολλή χρήμα το οποίο σαν ιδιότης δεν μπορώ να διαθέτω. Είναι δουλειά των κρατικών φορέων και χρειάζονται επιδοτήσεις τις οποίες δίνουν μόνο σε κρατικούς φορείς έρευνας.
Ότι κάνω το κάνω μόνος μου με τα μέσα και χρήματα που έχω έως τώρα.
Το δεν ξέρεις και δεν έχει εφαρμοστεί είναι γενικόλογα.

[seismic1]

Τώρα τροποποιώ τον κώδικα προσομοίωσης βάζοντας μέσα το επιταχυνσιογράφημα του μεγαλύτερου σεισμού στον κόσμο.
Το πρόγραμμα τώρα χρησιμοποιεί πραγματικά δεδομένα σεισμικής επιτάχυνσης από τον σεισμό του Tōhoku, αντικαθιστώντας την ημιτονοειδή διέγερση.
Οι σεισμικές επιταχύνσεις που καταγράφηκαν στον σταθμό Tsukidate στην Ιαπωνία κατά τον σεισμό του Μαρτίου 2011 (σεισμός Tōhoku). Συγκεκριμένα, παρουσιάζονται οι μέγιστες επιταχύνσεις που καταγράφηκαν:

Μέγιστη οριζόντια επιτάχυνση σε μία κατεύθυνση: 2,700 cm/sec² (ή 2.75g).
Μέγιστη οριζόντια επιτάχυνση (διάνυσμα αθροίσματος των δύο κατευθύνσεων): 2,983 cm/sec² (ή 3.04g).
Μέγιστη συνολική επιτάχυνση (3-D): 3,526 cm/sec² (ή 3.59g).

Επίσης, καταγράφηκε:

Μέγιστη κατακόρυφη επιτάχυνση: 1,880 cm/sec² (ή 1.92g).

# Φόρτωση δεδομένων σεισμικής επιτάχυνσης από το αρχείο
seismic_data = np.loadtxt(file_path, skiprows=1) # Παράλειψη της επικεφαλίδας

# Απόσπαση των στηλών δεδομένων
time_seismic = seismic_data[:, 0] # Χρόνος (s)
acc_seismic = seismic_data[:, 1] * 0.01 # Οριζόντια επιτάχυνση σε m/s²

# Συνάρτηση σεισμικής επιτάχυνσης μέσω παρεμβολής
from scipy.interpolate import interp1d
seismic_acceleration = interp1d(time_seismic, acc_seismic, fill_value="extrapolate")

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

# Χρονικά βήματα προσομοίωσης
T = time_seismic[-1] # Συνολικός χρόνος από το σεισμικό αρχείο
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="Conventional Structure", color='r') # σε mm
plt.plot(solution_innovative.t, solution_innovative.y[0] * 1000, label="Innovative Structure", color='b') # σε mm
plt.xlabel("Time (s)")
plt.ylabel("Displacement (mm)")
plt.title("Dynamic Response Comparison - Tōhoku Earthquake (Tsukidate)")
plt.legend()
plt.grid()
plt.show()

# Εκτύπωση αποτελεσμάτων
u_max_conventional, u_max_innovative

[Ασέβαστος]

Που το ξέρεις ρε βλάκα?

Το καταλαβαίνω από το πόσο άσχετος είσαι, και από το γεγονός ότι δεν έχει υπάρξει ούτε ένας που να δεχτεί να χτίσει με τη μαλακια σου, και ούτε μίσος σοβαρός άνθρωπος που να θεωρεί ότι έχεις δίκιο.

Λάθος Υπάρχουν πολλοί που θεωρούν ότι έχω δίκαιο. Και είναι και καθηγητές πανεπιστημίων που έχουν εξετάσει διεξοδικά την μέθοδο.

Μηδέν. Κανένας. Σε δουλεύουνε.