Citation for published version (APA):
Menken, C. M. (1982). Interactieve buigingsknik van een eenvoudig discreet model. (DCT rapporten; Vol. 1982.004). Technische Hogeschool Eindhoven.
Document status and date: Gepubliceerd: 01/01/1982
Document Version:
Uitgevers PDF, ook bekend als Version of Record
Please check the document version of this publication:
• A submitted manuscript is the version of the article upon submission and before peer-review. There can be important differences between the submitted version and the official published version of record. People interested in the research are advised to contact the author for the final version of the publication, or visit the DOI to the publisher's website.
• The final author version and the galley proof are versions of the publication after peer review.
• The final published version features the final layout of the paper including the volume, issue and page numbers.
Link to publication
General rights
Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain
• You may freely distribute the URL identifying the publication in the public portal.
If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license above, please follow below link for the End User Agreement:
www.tue.nl/taverne Take down policy
If you believe that this document breaches copyright please contact us at: openaccess@tue.nl
providing details and we will investigate your claim.
INTERACTIEVE BUIGINGSKNIK VAN EEN EENVOUDIG DISCRZET MODEL C. M. Menken 1. Gebruikte symbolen : a, b C e 1 U x9 z C E F P a 7r 7r2 hulpgrootheden breedte
indrukking veer, vervorming lengte axiale indrukking coördinaten hulpgrootheid (veer-)stijfheid drukkracht op veer drukkracht op constructie hoekverdraaiing
totale potentiële energie
kwadratische toename potentiële energie
Indices :
O initieel
k knik, critisch
L lokaal
2. Inleiding
Bij het onderzoek naar fenomenen die zich zouden kunnen voordoen bij door dwarskracht belaste geëxtrudeerde Al-profielen wordt o.a. gedacht aan een interactie van lokale knik en kip (inter- actieve kip).
Aan het fenomeen interactieve kpik van axiaal gedrukte kolommen, verstijfde plaatpanelen en cirkelcilinders werd door verschillende onderzoekers veel aandacht besteed
C l l .
Het bleek dat de oudste analysemethode, waarbij de lokale knik ontkoppeld geanalyseerd werd, waarna het daarbijhorende na-knik gedrag als een geredueeer- de stijfheid bij de constructie in rekening werd gebracht, door meer verfijnde methoden in grote lijnen bevestigd werd[21.
Van der Neut C31 was in 1969 de eerste die een continu model op deze wijze analyseerde en welde interactieve knik van een kolom. Hij besteedde daarbij o.a. aandacht aan de invloed van de model- lengte.
Alvorens de interactieve kip te analyseren leek het nuttig de bij knik gehanteerde werkwijzen te exploreren aan een discreet model,
íx we1 het m o d e l van Shanley L41. Thm~ps~:: C . S . [ S I grhriiikte d i t model om kwalitatief inzicht te krijgen in de knik van verstijfde platen. Bij deze studie was de lokale knikbelasting PL variabel, terwijl de modellengte 1 constant gehouden werd. Om kwalitatieve aansluiting bij het continue model van van der Neut te krijgen zal daarom hier ook bij het discrete model van Thompson aandacht aan de invloed van de lengte besteed worden. Bovendien zal aandacht geschonken worden aan het gemak van de werkwijze, zodat daarvan ook bij een komende analyse van de interactieve kip gebruik gemaakt kan worden. Alleen het perfecte model zal geanalyseerd worden.
3. Modelbeschrijving
Het mechanisch model is in fig.
1
gegeven. Bij dit model is zowel axiale indrukking als buiging mogelijk. Verondersteld is dat het model knikt bij een belasting Pk'Fig. 2 toont de karakteristieken van de ondersteunende veren. Deze representeren een lokaal constructiedeel dat kan knikken onder een belasting 'L/2 en dat daarna toch nog een (gemiddelde) gereduceerde stijfheid ET heeft.
/-
Fig.
1
Het bi-lineaire gedrag van de veren laat nu drie knikmogelijkheden toe.
1 .
Pk < PL ; beide veren hebben tijdens het knikken de stijfheid E. 2. PK > PL ; beide veren hebben tijdens het knikken de stijfheid ET. 3. Het grensgeval Pk = P L Y d.w.z. de totale knikbelasting is preciesgelijk aan de lokale. Bij het uitknikken zal de ene veer ontlas- ten en de stijfheid E behouden, terwijl de andere veer iets ver- der belast wordt en de gereduceerde stijfheid ET krijgt.
Het laatste, meest algemene, geval zal allereerst beschouwd worden. De beide andere kunnen daaruit verkregen worden. Het model zal op twee manieren geanalyseerd worden :
a) We gaan uit van een fictief model met lineaire veren. In plaats van de bi-lineaire veer nemen we een lineaire met initiële indrukking eo en stijfheid ET. Bovendien beschou- wen we de verplaatsingen van het "zwaartepunt van de dwars- doorsnede", punt Z in fig.
1 .
b) We gaan uit van het werkelijke model en beschouwen alleen de toename van de potentiële energie bij uirknikken. De verplaatsingen worden gerelateerd aan het "drukpunt" N met de eigenschap dat een vertikale belasting door N géén
rotatie veroorzaakt.
4 . Analyse van het fictieve model.
We veronderstellen dat de constructie bij het knikken rechtsom draait ( a > o). In het fictieve model heeft de linker veer dan de stijfheid E en de rechter de stijfheid ET, samen met een initiële indrukking eo. De grootheden die in onbelaste toestand een rol spelen zijn in fig. 3 gedefinieerd en de bij de belaste toestand horende in fig. 4 .
€!,*=O
-
Wij kiezen u (z = o) en a als vrijheidsgraden. Het model houdt op fictief te zijn voor a 2 o. Het verband tussen de fictieve voorvervorming e20 en de lokale kniklast PL volgt uit fig. 2 :
*
pL
e2 = - 2E PL BC = - ETDe bij het fictieve model horende potentiële energie is :
1 2
K = 1Eel2
+
yET(e2-
e20)-
P(u - u - 1 c o s a+
1
cos ao)O
De vervormingen el en e2 drukken we uit in de verplaatsingen :
C
2
e l = u - - a C 2 e 2 = u + - aBovendien nemen we aan dat u en a zo klein zijn dat we ons, na
v n - 1 Lccn,,LLLw,kkeling, "nntTT; t ~ t kvadratische termen kimmen beperken.
De potentiële energie wordt dan :
( 4 . 2 )
De evenwichtspaden volgen uit :
C C
- o = E(u - -a)
+
ET(u + - a - e20) - Pau 2 2
ar
_ -
( 4 . 3 )
ar C C C C
Eliminatie van u geeft :
E -E -
P - Pla = o
Deze niet lineaire vergelijking schrijven we als :
aa - b - CP
-
Pla = omet :
2EET
b : =
(2)
CE+ET
e20 < o (zie4.1)
(4.4)
(4.5)
We gaan ons nu een beeld vormen van de verbanden tussen P, 1 en a in de buurt van a = o en wel voor vaste P
L'
b
C L
p = - - = p (4.9)
Dit was ook ons uitgangspunt.
Indien hieraan ook voldaan is voor a t o, dan levert ( 4 . 5 ) ons een eigenwaardeprobleem :
(a - P 1)a = o, L
met als critische lengte :
a
1
= -pL
(4.
IO)Bij deze lengte hebben we dus met zogenaamd neutraal evenwicht te maken.
Met (4.5) kunnen we ook het verband tussen P en a in de buurt van
a = o en bij constante 1 onderzoeken. Uit (4.5) volgt :
a - P 1
dP L
C
<da>
=ol=0
We kunnen nu twee gevallen onderscheiden :
1. 1 < lky dus (a
-
PL1) > o,
dP
dan geldt
(z)
> o,
fig. 5a=0
2. i > lk, dus (a - PLl) < o
,
dP dol
dan geldt
(-1
< o , fig. 6 a=0o
-ct
I
(4.11)
-a
Fig. 5 Verband tussen P en a Fig. 6 Verband tussen P en a
voor 1 > lk k
We houden ons aan het geldigheidsgebied a
2
o . Zou de constructie linksom knikken, dan moeten we de figuren spiegelen om de P-as.De
---
stabiliteit bij P = PL beoordelen we aan de hand van het verbandtussen P en de daarbijhorende verplaatsing u. Na eliminatie van a
uit het oorspronkelijke stelsel vergelijkingen ( 4 . 3 ) vinden we :
dP a= O (ZÜ) a
-
P,1 L-
-
-P 1 L-
E+ETZoals hiervoor onderscheiden we weer twee gevallen :
1.
1
< lk9 dus (a-
PLi) > o : < o,
fig.7
O
-U
Fig.
7
Verband tussen P en u voor1
< lkFig.
8
Verband tussen P en u voor1
> lkVolgens fig. 5 geldt alleen het bij toenemende P horende deel van de vertakking. Het gedrag is dus stabiel.
2. 1 > lk, dus (a
-
PLl) < o :> o fig. 8 dP
+
a=0
Volgens fig. 4 . 6 geldt alleen het bij dalende P horende deel van de vertakking. Het gedrag is dus labiel.
Hiermee is het grensgeval
P
(zie pag. 3 ) zijn hieruit te verkrijgen. De bij het grensgeval horende critische belasting en/of lengte volgt immers uit (zie ( 4 . 6 ) en ( 4 . 1 0 ) ) :
= PL behandeld. De beide andere gevallen k
C 2 4EET
(Pl), = a =
(y)
- E+ET ( 4 . 1 2 )De critische belastingen horende bij P < P zijn hieruit te verkrijgen door
E
te vervangen doorE
:k L T 2 k < L' = iC E voor P ( 4 . 1 3 ) vervangen we in ( 4 . 1 2 ) E door E
-
k'
L' T 'Voor het geval P
( 4 . 1 4 ) 2
L (Pi), = iC ET voor Pk > P
De drie gevallen, samen met de mate van stabiliteit zijn verzameld in fig. 9. Kwalitatief komt deze figuur geheel overeen met door van der Neut C31 voor een continu model verkregen resultaten!
Imperfectie gevoeligheid :
dP Voor 1 > lk zien we uit ( 4 . 1 1 ) dat
tiever wordt. Uit fig. 4 . 6 is dan te concluderen dat de gevoeligheid voor imperfecties daarmee toe zal nemen. Het is te verwachten dat de grootste gevoeligheid zal optreden als de helling dP/da het grootste
P
I
I
2 C2E C EET-
2 P L PL(E+ET) 2% -Fig. 9 Het stabiliteitsgedrag van de constructie.
is. Dit treedt op als de kniklast van de constructie zonder lokale knik gelijk is aan de lokale kniklast (punt A in fig. 7 ) .
5. Vereenvoudigde werkwijzen.
Uit ( 4 . 1 2 ) , ( 4 . 1 3 ) en ( 4 . 1 4 ) blijkt dat ter bepaling van het gedrag van de perfecte constructie (fig. 9 ) uitdrukking ( 4 . 1 2 ) belangrijk is. Omdat het om een vertakkingspunt gaat kan volstaan worden met de kwa- dratische uitdrukking van de toename van de potentiële energie bij verstoring of uitknikken. Daarnaast zal het handiger zijn indien niet met de verplaatsingen van het hart van de constructie gewerkt wordt,
doch met die van het drukpunt van het fictieve model, d.w.z. dat punt waar de werklijn van de belasting doorheen moet gaan om te zorgen dat er beneden de critieke belasting géén rotatie optreedt.
Fig. 10
We nemen aan dat dit punt op een afstand d links van de oorsprong ligt. De indrukking van de veren is dan :
C = u
-
(7-
d)a e 1 C 2 2 e = u+
( - + d)aDe toename van de potentiële energie bij uitknikken is :
2
2
C 2 C
ïi2 = ~E{u
-
(7
-
d)a)+
4ET{u+
(-+
d)a) - 1P laUit de eis dat hierin géén producttermen ua voorkomen volgt de ligging van het drukpunt :
C E-ET 2 E+ET d = - -
Hiermee wordt n2 :
Een stationaire waarde treedt op als :
an2 - = o = (E+E )u au T an2 C2 EET - - aa - o = {
-
P11
a=o E+ET LUit de eerste voorwaarde volgt, zoals te verwachten [ S I , u=o. Uit de tweede voorwaarde volgt opnieuw (4.12).
LITERATUUR
1. Tvergaard, V., Buckling Behaviour of Plate and Shell Structures, Proc. 14th IUTAM Congress, Delft, 232
-
247.2. Budiansky B,, Hutchinson J.W., Buckling : progress and challenge, Proc. of the Symposium dedicated to the 65th Birthday o f W.T. Koiter. Sijthoff & Noordhoff International Publishers, 1979.
3 . Van der Neut A., The interaction of local buckling and column failure of thin-walled compression members, Proc. of the Twelfth Intern. Congr. on Appl. Mech., Stanford Univ. 1968.
4.
Shanley, F.R., J. Aero. Sci.-
14,
261, 1947.5. Thompson J.M.T., Tul1 J.D., Walker A.C., An experimental study of imperfection-sensitivity in the interactive buckling of stiffened plates. Buckling of Structures, IUTAM Symposium, Cambridge USA, 1974.
