SOFiSTiK Forum

It is currently Fri Apr 19, 2019 4:40 pm

All times are UTC




Post new topic Reply to topic  [ 4 posts ] 
Author Message
 Post subject: Pendelschlagversuch
PostPosted: Mon Jan 21, 2019 11:47 am 
Offline

Joined: Mon Jan 21, 2019 11:39 am
Posts: 4
Hallo!
Ich brauche dringend Hilfe.
Ich habe das Beispiel aus Teddy -ase - Pendelschlagversuch.
Nun ich erkenne nicht aus dem Beispiel, wo die Fallhöhe reguliert werden kann.

Ich wäre sehr dankbar,wenn mir jemand helfen könnte, da ich leider keine Ansprechperson zu diesem Thema finden kann und bald meine Arbeit abgeben muss... :?

Beispiel:

[color=#008000]Rechenmodell nach DIN 18008-4 Anhang C:
Nachweis der Stoßsicherheit von Glasaufbauten durch Berechnung

C.3 Nachweis über die volldynamisch transiente Simulation des Stoßvorganges

Referenzverglasung mit Abmessungen 876 mm x 1938 mm vierseitig gelagert
8 mm Glasdicke
450 mm Pendelhöhe mit den Referenzwerten

Die Ersatzsteifigkeit des Pendelkörpers wird mit 400 kN/m angegeben.

Bitte im Animator die Massen sichtbar schalten!

</TEXT>

+PROG AQUA URS:1
KOPF
NORM 'DIN' 'en1992-2004'
ECHO MAT VOLL
MATE 1 TYP GLAS
STAH 9 S 355 GAM 0
ROHR 9 D 100 T 1 MNR 9
ENDE

+PROG SOFIMSHA urs:2
KOPF Pendelschlagversuch
UNIT 6 $ mm
SYST RAUM GDIV 1000 POSZ
sto#h 876
sto#b 1938
KNOT NR X Z
101 -#b/2 -#h/2 $ contact springs:
102 -100 -#h/2 $ 203 204 205
103 -50 -#h/2 $ 302 303 304 305 306
104 0 -#h/2 $ 402 403 404 405 406
105 +50 -#h/2 $ 502 503 504 505 506
106 +100 -#h/2 $ 603 604 605
107 +#b/2 -#h/2
TRAN KNOT 101 107 DNR 100 DZ #h/2-100
TRAN KNOT 101 107 DNR 200 DZ #h/2-50
TRAN KNOT 101 107 DNR 300 DZ #h/2
TRAN KNOT 101 107 DNR 400 DZ #h/2+50
TRAN KNOT 101 107 DNR 500 DZ #h/2+100
TRAN KNOT 101 107 DNR 600 DZ #h

TRAN KNOT 202 606 DNR 1000 DY 200[mm]
TRAN KNOT 1404 DNR 1000 DY 10[mm]
TRAN KNOT 1404 DNR 2000 DY 210[mm]
KNOT 1404,2404,3404 FIX ypmm

GRUP 1
let#n1 4 $ center area 5*5 nodes - do not change '4'
let#n2 8 $ above center area - change allowed
let#n3 10 $ beside center area - change allowed
QUAD PROP T 8[mm]
QUAD - 101 201 202 102 M #n2 #n3
QUAD - 102 202 206 106 M #n2 #n1
QUAD - 106 206 207 107 M #n2 #n3

QUAD - 201 601 602 202 M #n2 #n3
QUAD - 202 602 606 206 M #n2 #n1
QUAD - 206 606 607 207 M #n2 #n3

QUAD - 601 701 702 602 M #n2 #n3
QUAD - 602 702 706 606 M #n2 #n1
QUAD - 606 706 707 607 M #n2 #n3

sto#cftot 400 $ Summe der Federsteifigkeiten aller Federn 400 kN\m.
sto#cfi 0.70*#cftot/9 $ 9 springs inside
sto#cfo 0.30*#cftot/12 $ 12 springs outside
GRUP 2
$ 9 inner springs:
let#nodes 9,303,304,305,403,404,405,503,504,505
loop#i #nodes(0)
FEDE #nodes(#i+1) KA #nodes(#i+1) KE #nodes(#i+1)+1000 CP #cfi RISS 0 GAP 5[mm]
if #nodes(#i+1)<>404
KNOT #nodes(#i+1)+1000 FIX KF 1404
endif
endloop
$ 12 outer springs: with a GAP t simulate later contact !
let#nodes 12,203,204,205,302,402,502,306,406,506,603,604,605
loop#i #nodes(0)
FEDE #nodes(#i+1) KA #nodes(#i+1) KE #nodes(#i+1)+1000 CP #cfo RISS 0 GAP 25[mm]
KNOT #nodes(#i+1)+1000 FIX KF 1404
endloop
GRUP 9
FEDE - KA 2404 KE 1404 CP 1E9 $ for plot use only
FACH - KA 3404 KE 2404 QNR 9 $ for animation (without mass - gamma=0)
MASS 3404 MX 50[kg]

$ --------------------- fixings ---------------------------------------------
SBOX A1X -1000 0 0 1000 0 0 A2X 0 -1000 0 0 1000 0 A3X 0 0 -#h/2-1 0 0 -#h/2+1
SBOX A1X -1000 0 0 1000 0 0 A2X 0 -1000 0 0 1000 0 A3X 0 0 +#h/2-1 0 0 +#h/2+1 OPT OR
SBOX A1X -#b/2-1 0 0 -#b/2+1 0 0 A2X 0 -1000 0 0 1000 0 A3X 0 0 -1000 0 0 +1000 OPT OR
SBOX A1X +#b/2-1 0 0 +#b/2+1 0 0 A2X 0 -1000 0 0 1000 0 A3X 0 0 -1000 0 0 +1000 OPT OR
KNOT SBOX FIX PY
FEDE SBOX DY 1 CQ 100
ENDE

+PROG ASE URS:6
KOPF Primärlastfall
LF 1 EGZ 1
ENDE

$ +PROG ASE
KOPF
ECHO VERS,REAK,SCHN,NOST,LAST,GRUP NEIN
ECHO FEDE JA $ aktiviert nur den Ausdruck der Feder- und Seilergebnisse

STEU ITER 0 $ quicker in case of little geometric nonlinear effect
STEP 60 DT 0.001 THE 0.70 LFSM 9901 $ alpha method acc. Hilber-Hughes-Taylor input THE= 1 - |alpha|. example alpha -0.3 --> THE=0.7
SYST PROB TH3 PLF 1
$ grp - ; grp 2 rada 0 radb 0.005 $ stiffness proportional damping in contact springs - not useful
LF 2 EGZ 1
V0 1404,2404,3404 VY -2.97 $ Anfangsgeschwindigkeit 2.97 m/s bei 450 mm Pendelhöhe
ENDE

+PROG ASE urs:3
KOPF Plot after analysis
PLOT LF 1 KNR 3404 RICH AY $ Pendelbeschleunigung
PLOT LF 1 NR 9001 RICH N $ Kontaktkraft gesamt
PLOT LF 1 NR 2404 RICH N $ innerste Feder
PLOT LF 1 NR 2304 RICH N $ innere Feder
PLOT LF 1 NR 2203 RICH N $ äußere Feder (->GAP)
ENDE

<TEXT>
Bitte im Animator die Massen sichtbar schalten
</TEXT>

+PROG RESULTS urs:4
KOPF
DIAG LFTP ZEIT RECH LF $ x-Axis = time
FILT NAME "spri_res.nr" TYP EINZ VAL1 9001 $ select this element no.
$ DBO NAME "spri_res.p" FAC -1 $ turn plot, neg. values upside
LF NR (1 61 1)
STAB TYP N STYP FEDE DARS DIAG
ENDE
$ Übersicht Ergebnisausgaben mit PROG RESULTS siehe ase.dat\deutsch\ase1_uebersicht.dat
[/color]


Top
 Profile  
 
 Post subject: Re: Pendelschlagversuch
PostPosted: Mon Jan 21, 2019 12:05 pm 
Offline

Joined: Tue Jun 20, 2006 9:32 am
Posts: 428
Hallo,

aus der Fallhöhe wird die Anfangsgeschwindigkeit VY = -2.97 m/s berechnet und angesetzt:
v = SQRT(2 * g * h) in m/s
g = 9.81 m/s^2
h = Fallhöhe in Metern!!

_________________
Schöne Grüße

Philipp Loosmann


Top
 Profile  
 
 Post subject: Re: Pendelschlagversuch
PostPosted: Mon Jan 21, 2019 12:12 pm 
Offline

Joined: Mon Jan 21, 2019 11:39 am
Posts: 4
Hallo!
Danke, wollte grade schreiben, dass es sich erledigt hat:)
Danke trotzdem nochmal. :)


Top
 Profile  
 
 Post subject: Re: Pendelschlagversuch
PostPosted: Wed Mar 20, 2019 9:30 am 
Offline

Joined: Mon Jan 21, 2019 11:39 am
Posts: 4
Hallo!
Ich bräuchte nochmal gaanz dringend Hilfe :?
Und zwar..wenn man die Werte für die Kanten, also n2 und n3 ändert, so dass n2=10 und n3=8 wird..ändern sich meine Ergebnisse bzw. die Beschleunigung im Knoten 404 (also am Glas) um ein vielfaches und entspricht auch den Messungen aus Versuchen.
So wie ich das verstanden hatte, wird mit dem QUAD "M" die Feinheit bestimmt..Die Abweichungen sollten dementsprechend nicht so groß sein, jedoch sind die Werte um merh als das doppelte höher...
Ich wäre sehr dankbar, wenn mir das jemand erklären könnte..

Liebe Grüße
Zulai


GRUP 1
let#n1 4 $ center area 5*5 nodes - do not change '4'
let#n2 8 $ above center area - change allowed
let#n3 10 $ beside center area - change allowed
QUAD PROP T 8[mm]
QUAD - 101 201 202 102 M #n2 #n3
QUAD - 102 202 206 106 M #n2 #n1
QUAD - 106 206 207 107 M #n2 #n3


Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 4 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 7 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group