Branch Switch

From FEAP Wiki
Revision as of 15:22, 13 July 2023 by S g (talk | contribs) (Created page with "Branch switch in FEAP is possible if one use arclength continuation methods. Here is a simple example. Consider a built-in column with an axial load. <pre> feap ** Buckling Path Following with Arclength ** 0 0 0 3 3 8 para e = 1.0 d = 0.5 h = 20 pi = 4.*atan(1.) i = 16.*d*d*d*d/12 pe = pi*pi*e*i/h/h/4 ! Euler load p = 0.1*pe/4/d/d ! Reach euler load a prop value of 1 mate solid elastic moon e 0.4 0.05 finite mixed block cart 4 4 20*2...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Branch switch in FEAP is possible if one use arclength continuation methods. Here is a simple example. Consider a built-in column with an axial load.

feap ** Buckling Path Following with Arclength **
0 0 0 3 3 8

para
 e = 1.0
 d = 0.5
 h = 20
 pi = 4.*atan(1.)
 i  = 16.*d*d*d*d/12
 pe = pi*pi*e*i/h/h/4       ! Euler load
 p  = 0.1*pe/4/d/d          ! Reach euler load a prop value of 1


mate
 solid
  elastic moon e 0.4 0.05
  finite
  mixed

block
 cart 4 4 20*2
 1 -d -d 0
 2  d -d 0
 3  d  d 0
 4 -d  d 0
 5 -d -d h
 6  d -d h
 7  d  d h
 8 -d  d h

eboun         ! Fixed base
 3 0 1 1 1

csur          ! Pressure load (dead)
 surface
 1 -d -d h -p
 2  d -d h -p
 3  d  d h -p
 4 -d  d h -p

end