Difference between revisions of "Branch Switch"

From FEAP Wiki
Jump to navigation Jump to search
(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...")
(No difference)

Revision as of 15:22, 13 July 2023

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