Ich bin mir nicht sicher, ob ich so die Vorspannkraft simulieren kann oder ob es dafür andere Wege gäbe. Ich versuche derzeit, die Vorspannkraft über eine axiale Normalkraft zu modellieren (siehe Code unten).
Meine Handrechnungen und das Experiment ergeben jedoch größere Abweichungen.
Für ein besseres Verständnis, unser Aufbau war folgender:
- Schwingungsmesser an den Viertelspunkten.
- Querschnitt ist aus dem Code ersichtlich.
- mehrere dünne verpresste Carbonmatten
Das Kapitel auf dem meine Handrechnungen basieren ist aus Dynamik der Baukonstruktionen von Petersen “10.2.2.3.3 Grundstäbe – Stabwerke bei Einwirkung einer konstanten Druckkraft“

wobei omega_j die j-te Eigenform unter der Vorspannungskraft D ist, D die kritische Knicklast der j-ten Eigenform und omega_0,j die j_te Eigenform ohne eine einwirkenge Druckkraft also die j_te Grundform.
Zwischen keiner Vorspannkraft und 40% Vorspannkraft mit umgerechnet 47kN ergab sich bereits eine Differenz von 1 Hz (f_0=78,72 Hz und f_40=77,77Hz), welche ich in Sofistik jedoch noch nicht ansatzweise Nachweisen konnte.
Mein Teddy-Code lautet:
!+!Kapitel Material und Querschnitte
$ ==================================
$ Norm, Material und Querschnitt
+prog aqua urs:1
kopf Norm, Material und Querschnitt
NORM EN 199X-200X
!*!Label Material und Querschnitt
BETO NR 1 ART C FCN 52
!*!Label Querschnitte
QB NR 1 H 15[cm] B 15[cm] HO 10[cm] BO 4[cm] MNR 1
!*!Label globale Variablen
STO#L 1.65[m]
ende
!+!Kapitel SOFIMSHA: System und Geometrie
$ ========================================
+prog sofimsha urs:2
kopf System und Geometrie
!*!Label SYST Koordinatensystem
SYST ART 3D GDIV 1000
$Knoten - Absolutkoordinaten
$Knoten unten
KNOT NR 1 X 0 Y 0 FIX ppmx
KNOT NR 2 X #L Y 0 FIX pypzmx
GRUP NR 1
STAB KA 1 KE 2 QNR 1 TEIL 100 KM 0
ende
!+!Kapitel Belastung
$===================
+prog sofiload urs:3
kopf Belastungen mit Vorspannkraft
LF 1 TYP NONE BEZ ‘PX’
STAB VON GRP BIS 1 TYP PXX PA 47 PE -47
ende
!#!KAPITEL Analysis
+PROG ASE URS:8
KOPF Zustandslastfall für Beuleigenwertermittlung
LF 1000 BEZ Beuleigenwert_PLF
LC 1 $ PX
ENDE
!#!KAPITEL Buckling eigenvalues
+PROG ASE urs:4
KOPF Beuleigenwerte für Spannungen aus LF 100
SYST PLF 1000 FAKV 1.0
EIGE 5
ENDE
!+!Kapitel Schnittgrößen - linear - Vergleichswerte
$ ===============================
+prog ase urs:9
kopf Schnittgrößen linear ohne Berücksichtigung von Imperfektionen
SYST PROB line
LF 100 EGZ 1.35 $Eigengewicht wird nun von Sofistik berechnet, mit 1.35facher berücksichtigung sonst 0
ende
!#!Kapitel Eigenformen - linear
+prog dyna urs:5
kopf Die ersten 5 Eigenformen
SYST ART REST QNR 1 PROB LINE STAT SERV
EIGE NEIG 5
ende
Es freut mich sehr, dass Sie sich dennoch für dieses komplexe Problem interessieren und bin für jede Hilfe dankbar.
Mit freundlichen Grüßen