Λοξίας έγραψε: 19 Φεβ 2025, 21:38
seismic1 έγραψε: 19 Φεβ 2025, 21:31
Λοξίας έγραψε: 19 Φεβ 2025, 21:28
Έτσι μπράβο! Και να θυμάσαι: ΜΗΝ ξαναπείς ψέματα πως έχεις πτυχίο και δικαίωμα μελέτης, διότι θα στο ξαναφέρω αυτό το ποστ...
Και γενικότερα, κόψε τα ψέματα. Δεν θα γίνεις σπουδαίος με ψέματα. Ο περίγελως του ίντερνετ θα γίνεις.
..Έχω πτυχίο και δικαίωμα μελέτης..Έχω πτυχίο και δικαίωμα μελέτης..Έχω πτυχίο και δικαίωμα μελέτης..
Και είσαι μαλάκας και εγώ επιστήμονας.
Πως σε τσούζει όταν ξεμπροστιάζω τα ψέματά σου, ε;
Πάνω που έλεγες ότι θα το χάψουν πως έχεις πτυχίο...φααααάπ! να σου ο Λοξίας να σου ρίξει τη φάπα και να σου χαλάσει την παραμύθα...
Βέλαξε λιγάκι ακόμα...
Ρε όλους σας έχω για τα αρχίδια μου.
Τρέξε εσύ και οι άλλοι αυτόν τον κώδικα που έκανα Μπορείτε? Όχι. Μόνο εγώ μπορώ. Τότε τι στο διάολο επιστήμονες είσαστε?
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp
from scipy.interpolate import interp1d
# Φόρτωση σεισμικών δεδομένων (Tohoku 2011)
file_path = "simulated_tohoku_acceleration.txt" # Αντικατέστησε με το σωστό αρχείο
seismic_data = np.loadtxt(file_path, skiprows=1)
# Απόσπαση δεδομένων
time_seismic = seismic_data[:, 0] # Χρόνος (s)
acc_seismic = seismic_data[:, 1] * 0.01 # Επιτάχυνση (m/s²)
# Παρεμβολή σεισμικής επιτάχυνσης
seismic_acceleration = interp1d(time_seismic, acc_seismic, fill_value="extrapolate")
# Σταθερές του προβλήματος
g = 9.81 # Επιτάχυνση βαρύτητας (m/s²)
m_elastic = 260000 # kg (μάζα ελαστικού φορέα)
m_shear_wall = 155000 # kg (μάζα φρεατίου)
H_structure = 18 # Ύψος κατασκευής (m)
k_elastic = 3e6 # N/m (Δυσκαμψία ελαστικού φορέα)
k_fixed = 10 * k_elastic # Ακαμψία πάκτωσης
k_post_tension = 5 * k_elastic # Ακαμψία προέντασης
k_shear_wall = k_fixed + k_post_tension # Συνολική ακαμψία φρεατίου
c = 5000 # Ns/m (Απόσβεση)
# Σεισμικοί αρμοί που αυξάνονται σταδιακά καθ' ύψος
seismic_joints = np.linspace(0.05, 0.3, num=int(H_structure/3)) # Από 5 cm στη βάση έως 30 cm στην κορυφή
# Εξισώσεις κίνησης
def equation_of_motion(t, y, mass, stiffness):
u, v = y # u: μετατόπιση, v: ταχύτητα
a = (mass * seismic_acceleration(t) - c * v - stiffness * u) / mass
return [v, a]
# Χρονικά βήματα προσομοίωσης
T = time_seismic[-1]
dt = 0.01
time_span = (0, T)
time_eval = np.arange(0, T, dt)
initial_conditions = [0, 0]
# Αριθμητική επίλυση για ελαστικό φορέα
solution_elastic = solve_ivp(equation_of_motion, time_span, initial_conditions, t_eval=time_eval,
args=(m_elastic, k_elastic), method='RK45')
# Αριθμητική επίλυση για φρεάτιο
solution_shear_wall = solve_ivp(equation_of_motion, time_span, initial_conditions, t_eval=time_eval,
args=(m_shear_wall, k_shear_wall), method='RK45')
# Υπολογισμός σύγκρουσης μεταξύ των δύο φορέων μέσω σεισμικών αρμών
collision_force = np.where(np.abs(solution_elastic.y[0] - solution_shear_wall.y[0]) > seismic_joints[-1],
k_elastic * (np.abs(solution_elastic.y[0] - solution_shear_wall.y[0]) - seismic_joints[-1]), 0)
# Υπολογισμός μέγιστης μετατόπισης
u_max_elastic = np.max(np.abs(solution_elastic.y[0])) * 1000 # mm
u_max_shear_wall = np.max(np.abs(solution_shear_wall.y[0])) * 1000 # mm
# Δεδομένα προέντασης και πακτώσεων
Tendon_force = 750 * 3 * 1e3 # Συνολική εφελκυστική δύναμη από 3 τένοντες (N)
Anchorage_force = 1000 * 3 * 1e3 # Συνολική αντοχή πακτώσεων σε έλξη (N)
lever_arm_tendons = 3.5 # m (μοχλοβραχίονας τένοντα)
lever_arm_anchorage = 3.6 / 2 # m (μοχλοβραχίονας πάκτωσης)
# Υπολογισμός πρόσθετης αντίρροπης ροπής λόγω προέντασης και πακτώσεων
M_resisting_tendons = Tendon_force * lever_arm_tendons
M_resisting_anchorage = Anchorage_force * lever_arm_anchorage
M_total_resisting = M_resisting_weight + M_resisting_tendons + M_resisting_anchorage
# Συντελεστής ασφάλειας
FS_shear_wall = M_total_resisting / M_overturning_shear_wall
# Γραφική παράσταση
plt.figure(figsize=(8, 5))
plt.plot(solution_elastic.t, solution_elastic.y[0] * 1000, label="Elastic Frame", color='g')
plt.plot(solution_shear_wall.t, solution_shear_wall.y[0] * 1000, label="Shear Wall", color='b')
plt.xlabel("Time (s)")
plt.ylabel("Displacement (mm)")
plt.title("Dynamic Response with Seismic Joints")
plt.legend()
plt.grid()
plt.show()
# Εκτύπωση αποτελεσμάτων
print(f"Μέγιστη μετατόπιση (ελαστικός φορέας): {u_max_elastic:.2f} mm")
print(f"Μέγιστη μετατόπιση (φρεάτιο): {u_max_shear_wall:.2f} mm")
print(f"Ροπή ανατροπής φρεατίου: {M_overturning_shear_wall:.2f} MN·m")
print(f"Συνολική αντίρροπη ροπή: {M_total_resisting:.2f} MN·m")
print(f"Συντελεστής ασφάλειας φρεατίου: {FS_shear_wall:.2f}")
