Arc-Length

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Post Reply
imancivil
Posts: 55
Joined: Mon Feb 11, 2008 11:15 pm
Location: Kerman

Arc-Length

Post by imancivil »

Hi all

In nonlinear finite element analysis references,such as Crisfield, Arc-Length is one of solution algorithm for nonlinear equations, but in Opensees Arc-Length is one of integrator methods. Is there any difference between them? I want to use Arc-Length method for algorithm command, and Newmark method for integrator command. Is it possible?

Thanks.
Last edited by imancivil on Mon Dec 01, 2008 8:40 am, edited 1 time in total.
Iman Mansouri
PhD Candidate
University of Kerman
fmk
Site Admin
Posts: 5884
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley
Contact:

Post by fmk »

it is not possible. Arc-Length is a continuation method attempting to solve (9.15 in Crisfield)

g(P,lambda) = Qi(P) - lambda *Qref, Qi the internal forces a function of fixed external forces P and displacements, and Qref the reference loads. arclength it is attempting to determine the lambda. the loadfactor to be applied to external reference loads to solve this
condition and the other conditions of the continuation method, in arc-length case this condition is based on the arc-length.

In OpenSees we call these static integrators, and we use solution algorithms such as Newmark to iterate to the solution.
imancivil
Posts: 55
Joined: Mon Feb 11, 2008 11:15 pm
Location: Kerman

Post by imancivil »

Dear Frank

Thank you.
Last edited by imancivil on Mon Dec 01, 2008 8:42 am, edited 1 time in total.
Iman Mansouri
PhD Candidate
University of Kerman
TWI
Posts: 12
Joined: Thu Aug 14, 2008 3:11 pm

Post by TWI »

fmk wrote:it is not possible. Arc-Length is a continuation method attempting to solve (9.15 in Crisfield)

g(P,lambda) = Qi(P) - lambda *Qref, Qi the internal forces a function of fixed external forces P and displacements, and Qref the reference loads. arclength it is attempting to determine the lambda. the loadfactor to be applied to external reference loads to solve this
condition and the other conditions of the continuation method, in arc-length case this condition is based on the arc-length.

In OpenSees we call these static integrators, and we use solution algorithms such as Newmark to iterate to the solution.
A follow-up question on this: Newmark method is a direct integration method which finally gives a nonlinear equations with effective stiffness and load vector, so in order to solve it, an advancing scheme still need to be specified, say 'arc-length',' loadcontrol' or 'displacement control'. But opensees did not ask for it when performing the transient analysis, so is there a default setting regarding this issue?
fmk
Site Admin
Posts: 5884
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley
Contact:

Post by fmk »

ArcLength, DisplacementConrol and others are Static integrators. Newmark, HHT are dynamic integrators. The integrators just define what the equilibrium equation is at each point in time.

In OpenSees you specify in addition o others the Integrator and the Algorithm. It is the solution algorithm that handles solvuing the nonlinear equation. Example algorithms are Linear, NewtonRaphson, ModifiedNewton, etc.
TWI
Posts: 12
Joined: Thu Aug 14, 2008 3:11 pm

Post by TWI »

fmk wrote:ArcLength, DisplacementConrol and others are Static integrators. Newmark, HHT are dynamic integrators. The integrators just define what the equilibrium equation is at each point in time.

In OpenSees you specify in addition o others the Integrator and the Algorithm. It is the solution algorithm that handles solvuing the nonlinear equation. Example algorithms are Linear, NewtonRaphson, ModifiedNewton, etc.
Thanks for the reply, Frank. My understanding was Newton-family algorithms provide a corrector phase when solving the nonlinear equation, however, certain advacing scheme (let the solution moves forward) has to be provided too. Maybe Opensees used a different way of handling this issue.
susan
Posts: 110
Joined: Sun Mar 02, 2008 12:10 am

Post by susan »

Hello,

Dr. Frank McKenna suggested alpha=0.0
would you please tell me what the reason is?

generally, how should we assign $arclength $alpha?

Thank you in advance.
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

do you really have a snap-back/snap-through problem?
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Post Reply