A computational framework is proposed to perform parameter continuation of periodic solutions of nonlinear, distributed-parameter systems represented by partial differential equations with time-dependent coefficients and excitations. The path-following procedure, encoded in the general-purpose Matlab-based computational continuation core (referred to below as coco), employs only the evaluation of the vector field of an appropriate spatial discretization; for example as formulated through an explicit finite-element discretization or through reliance on a black-box discretization. An original contribution of this paper is a systematic treatment of the coupling of coco with Comsolmultiphysics, demonstrating the great flexibility afforded by this computational framework. Comsolmultiphysics provides embedded discretization algorithms capable of accommodating a great variety of mechanical/physical assumptions and multiphysics interactions. Within this framework, it is shown that a concurrent bifurcation analysis may be carried out together with parameter continuation of the corresponding monodromy matrices. As a case study, we consider a nonlinear beam, subject to a harmonic, transverse direct excitation for two different sets of boundary conditions and demonstrate how the proposed approach may be able to generate results for a variety of structural models with great ease. The numerical results include primary-resonance, frequency-response curves together with their stability and two-parameter analysis of multistability regions bounded by the loci of fold bifurcations that occur along the resonance curves. In addition, the results of comsol are validated for the Mettler model of slender beams against an in-house constructed finite-element discretization scheme, the convergence of which is assessed for increasing number of finite elements.

References

References
1.
Krauskopf
,
B.
, and
Osinga
,
H. M.
,
2007
,
Numerical Continuation Methods for Dynamical Systems
,
J.
Galan-Vioque
, ed.,
Springer and Canopus Publishing Limited
,
Bristol, UK
.
2.
Govaerts
,
W. J. F.
,
2000
,
Numerical Methods for Bifurcations of Dynamical Equilibria
,
Society for Industrial and Applied Mathematics
,
Philadelphia, PA
.
3.
Salinger
,
A. G.
,
Burroughs
,
E. A.
,
Pawlowski
,
R. P.
,
Phipps
,
E. T.
, and
Romero
,
L. A.
,
2005
, “
Bifurcation Tracking Algorithms and Software for Large Scale Applications
,”
J. Bifur. Chaos Appl. Sci. Eng.
,
15
(
3
), pp.
1015
1032
.10.1142/S0218127405012508
4.
Bank
,
R. E.
,
PLTMG
,
1998
,
A Software Package for Solving Elliptic Partial Differential Equations, Users’ Guide 8.0
,
Society for Industrial and Applied Mathematics
,
Philadelphia, PA
.
5.
Arena
,
A.
,
Formica
,
G.
,
Lacarbonara
,
W.
, and
Dankowicz
,
H.
,
2011
, “
Nonlinear Finite Element-Based Path Following of Periodic Orbits
,”
Proceedings of the ASME 2011 IDETC MSNDC 2011
,
Washington, DC
, Aug. 29-31, Paper No. DETC2011-48673.
6.
Zimmerman
,
W. B. J.
,
2006
, “
Comsol Multiphysics and the Basics of Numerical Analysis
,”
Multiphysics Modeling With Finite Element Methods, Series on Stability Vibration and Control of Systems
, vol.
18
,
World Scientific Publishing Company
,
Singapore
.
7.
Formica
,
G.
,
Arena
,
A.
,
Lacarbonara
,
W.
, and
Dankowicz
,
H.
,
2011
, “
General-Purpose Finite Element-Based Path Following of Nonlinear Dynamical Systems
,”
Proceedings of the XX AIMETA Congress
,
Bologna
, Sept. 12–15.
8.
Huang
,
Z.
,
Conway
,
P. P.
,
Thomson
,
R. C.
,
Dinsdale
,
A. T.
, and
Robinson
,
J. A. J.
,
2008
, “
A Computational Interface for Thermodynamic Calculations Software MTDATA
CALPHAD: Comput. Coupling Phase Diagrams Thermochem.
,
32
, pp.
129
134
.10.1016/j.calphad.2007.07.003
9.
Dankowicz
,
H.
, and
Schilder
,
F.
,
2011
, “
An Extended Continuation Problem for Bifurcation Analysis in the Presence of Constraints
,”
J. Comput. Nonlinear Dyn.
,
6
(3), art. no.
031003
.10.1115/1.4002684
10.
Ambrose
,
D. M.
, and
Wilkening
,
J.
,
2010
, “
Computation of Time-Periodic Solutions of the Benjamin-Ono Equation
,”
J. Nonlinear Sci.
,
20
(
3
), pp.
277
308
.10.1007/s00332-009-9058-x
11.
Lacarbonara
,
W.
, and
Yabuno
,
H.
,
2006
, “
Refined Models of Elastic Beams Undergoing Large In-Plane Motions: Theory and Experiment
,”
Int. J. Solids Struct.
,
43
, pp.
5066
5084
.10.1016/j.ijsolstr.2005.07.018
12.
Lacarbonara
,
W.
,
2013
,
Nonlinear Structural Mechanics. Theory, Dynamical Phenomena and Modeling
,
Springer
,
New York
.
13.
Hughes
,
T. J. R.
,
2000
,
The Finite Element Method: Linear Static and Dynamic Finite Element Analysis
,
Dover Publications
,
New York
.
14.
Doedel
,
E. J.
,
Paffenroth
,
R. C.
,
Champneys
,
A. R.
,
Fairgrieve
,
T. F.
,
Kuznetsov
,
Yu. A.
,
Sandstede
,
B.
, and
Wang
,
X.
,
2001
, “
AUTO 2000: Continuation and Bifurcation Software for Ordinary Differential Equations (With HomCont)
,”
Technical Report
, Caltech.
15.
Lacarbonara
,
W.
,
Arafat
,
H.
, and
Nayfeh
,
A. H.
,
2005
, “
Nonlinear Interactions in Imperfect Beams at Veering
,”
Int. J. Nonlinear Mech.
,
40
, pp.
987
1003
.10.1016/j.ijnonlinmec.2004.10.006
16.
COMSOL AB
,
2008
,
COMSOL Multiphysics User Guide and Model Library
, version 3.5a,
Sweden
.
You do not currently have access to this content.