.. image:: https://badges.gitter.im/mousai_hb/Lobby.svg :alt: Join the chat at https://gitter.im/mousai_hb/Lobby :target: https://gitter.im/mousai_hb/Lobby?utm_source=badge&utm_medium=badge&utm_campaign=pr-badge&utm_content=badge
.. image:: https://badge.fury.io/py/mousai.png/ :target: http://badge.fury.io/py/mousai
.. image:: https://travis-ci.org/josephcslater/mousai.svg?branch=master :target: https://travis-ci.org/josephcslater/mousai
.. image:: https://img.shields.io/badge/Say%20Thanks-!-1EAEDB.svg :target: https://saythanks.io/to/josephcslater
.. image:: https://img.shields.io/badge/patreon-donate-yellow.svg :target: https://www.patreon.com/josephcslater
.. image:: http://pepy.tech/badge/mousai :target: http://pepy.tech/project/mousai :alt: PyPi Download stats
.. image:: https://mybinder.org/badge.svg :target: https://mybinder.org/v2/gh/josephcslater/mousai/binder
Python-based Harmonic Balance solvers and relevant tools.
Please see the manual <https://josephcslater.github.io/mousai/>__ for usage and installation, and a nice overview from my Dayton Engineering Sciences Symposium 2017 presentation <https://josephcslater.github.io/mousai/_static/Theory_and_Examples.slides.html#/>_. You can play around with it by selecting the launch binder <https://mybinder.org/v2/gh/josephcslater/mousai/binder> button above.
Mousai can solve sets of first-order and second-order ordinary differential equations written in state-space form (solved for acceleration for second-order form) subject to a harmonic excitation. All you need to provide is the name of a Python function which may simply be a wrapper to an external code. The function must accept:
The function must return the state derivatives (acceleration for second order form). It is expected that in most cases this will simply be a python wrapper function to call an external finite element, electro-magnetic, molecular dynamic, computational fluid dynamic, or other simulation code.
No traditional numerical integration is performed in harmonic balance, so issues of stability, reliability, etc. are separate concerns. An inability to determine a solution is reported if failure occurs. A good initial guess is always helpful, sometimes essential, and may simply have the form of a nearby solution. Sweeping through frequencies or increasing amplitudes of excitation is the best way to ensure swift convergence.
Quasi-linear models can (not yet implemented) obtain good low amplitude solutions that can be used as starting guesses for solutions (appropriate scaling may also help).
The manual is still under development, but the Tutorial provides two examples of solutions to the Duffing Oscillator, one in first order form, the other in second order form.
The Reference section of the manual illustrates supporting code for building out more time-refined solutions, finding velocities from displacements, and specifics on function calls.
Please see the manual <https://josephcslater.github.io/mousai/>__ for usage and installation.
Please see release notes <https://github.com/josephcslater/mousai/tree/master/release_notes>_ for updates, fixes, and plans.
If you use this for any work, please cite it as:
Slater, Joseph C., "Mousai: An Open Source Harmonic Balance Solver for Nonlinear Systems," 13th ASME Dayton Engineering Sciences Symposium, Oct. 30, 2017, Dayton, OH. <https://josephcslater.github.io/DESS2017/#/>_
Like this module, buy me a coffee! <https://www.buymeacoffee.com/s6BCSuEiU>_