The following matlab project contains the source code and matlab examples used for jacobian free newton krylov jfnk method. Also download the constants file here, and the ground motion file. Reliabilitybased optimal design software for earthquake. Opensees opensees modelling of tall buildings in fire. Opensees is an objectoriented framework for finite element analysis. Can anyone help me to solve this problem in open sees. In this study, genetic algorithm and particle swarm optimization algorithm are used for truss optimization.
For strongly nonlinear problems, on the one hand, the elements must consider the material and geometric nonlinearity simultaneously. Tidriri the international journal of high performance computing applications 2016 14. Stko scientific toolkit for opensees is a powerful pre and postprocessor and cuttingedge data visualization tool for opensees. Modified newton algorithm uses the tangent at the first iteration to iterate to convergence. This analysis also used a krylovnewton solution algorithm. The website provides an excellent wiki as well as several example scripts outlining the steps needed to create a model and perform the different types of analysis.
Find file copy path fetching contributors cannot retrieve contributors at this time. How to tackle the convergence problems in opensees. Opensees the open system for earthquake engineering simulation is a software framework for simulating the seismic response of structural and geotechnical systems and it is primarily developed and maintained by peer pacific earthquake engineering research center at. Solving nonlinear equations with newtonkrylov method based on automatic di. This function solves a system of nonlinear equations using the jacobianfree newtonkrylov jfnk method.
Novel structural nvstructural is the advanced nonlinear structural analysis and earthquake engineering software. The software is not commercialized yet, and it is still an opensource program to date, it has mostly been used in academia by several research teams around the world. The work extends the opensees software, which is rapidly gaining users in the earthquake engineering community. Solving nonlinear equations with newtonkrylov method based. Opensees is an objectoriented framework for building models of structural and geotechnical systems, performing nonlinear analysis with the model, and processing the response results. As for material models, the bilinear model, which is named steel01 in opensees, is adopted for steel rebars. Graphical representation of the rst two krylov newton iterations software implementation within the framework of opensees mckenna et al. Two n concentrated loads were distributed over three nodes each at a distance of about 0. Regarding iterative integration and solution, this paper employs krylov newton algorithm to facilitate the convergence of program computation. Dynamic analysis of multihazardresistant bridge piers having. Thus, the linear algorithm will be used to solve the system of equations.
This command is used to construct a solutionalgorithm object, which determines the sequence of steps taken to solve the nonlinear equation. It is still a big challenge to calculate the structural reliability index. At high levels of the framework, new classes handle meshing and interface detection of. Interior point solver with krylov newton system solver. Printed on 1 july, 2007 open system for earthquake engineering simulation opensees opensees command language manual silvia mazzoni, frank mckenna, michael h. Evaluation of deformationbased seismic performance of recc. For this purpose, we implement and apply a stateoftheart algorithm for reliabilitybased design optimization. Solving nonlinear equations with newtonkrylov method. Opensees, as an opensource finite element software, is now widely used owing to its high transparency and freedom mckenna, scott and fenves 2009.
Focus is on seismic design, in which nonlinear structural response must be considered. To ensure the convergence, multiple iterative algorithms were preset, such as newtonraphson method, modified newton. Dynamic behavior prediction of a reinforced concrete building. The code for the opensees software is open source, or freely available to. Recommendations for numerical modelling of concentrically. In this study, the open source fea software opensees is used to apply incremental dynamic analysis. This command is used to construct a krylovnewton algorithm object which uses a krylov subspace accelerator to accelerate the convergence. Some numerical applications has been performed, and the results obtained by using opensees code have been compared with the ones provided by the software c. I tried krylovnewton algorithm, and checked that all the node are been connected. Nonlinear finiteelement analysis software architecture using object composition. Introduccion opensees civil engineering structural. The interpreter is an extension of the tcl page 14 scripting language. In opensees, the residual is referred to as the unbalance.
As the algorithm approaches an optimal solution, this system of equations becomes very poorly conditioned, making iterative methods a poor choice. The most popular one is the substructuring method in which the response of the foundation soil and its interaction with the pile foundation and the abutment system are represented by a set of onedimensional springs and dashpots. Opensees intended users are in the research community. Krylovnewton algorithm, newton with line search algorithm, broyden algorithm. Due to the size of the numerical models, a parallel version of opensees openseesmp version 3. Reliability assessment of reinforce concrete structures. To solve the nonlinear equilibrium equation, the krylov newton algorithm provided by opensees was used. The opensees interpreter adds commands to tcl for finite element analysis. It allows users to create both serial and parallel finite element computer. Asdea software will take part in the conference which aims to entice new users and bring together the opensees european community. It has advanced capabilities for modeling and analyzing the nonlinear response of systems using a wide range of material models, elements, and solution algorithms. Comparison of genetic algorithm ga and particle swarm. Open sees command language manual june 2006 free ebook download as pdf file. Make sure the numpy, matplotlib and eqsig packages are installed in your python distribution.
Homebrew is a free and opensource software package management system that simplifies the. The schematic diagram of analytical procedure is shown in figure 3. Globalized newtonkrylovschwarz algorithms and software. This interpreter is an extension of the tcltk language for use with opensees.
The source code is shown below, which can be downloaded here. The performance of resistance progressive collapse analysis. The performance of resistance progressive collapse. Opensees abstract the opensees nite element software framework is extended for simulating uid structure interaction fsi by the particle nite element method pfem. The program will use displacementcontrol as the algorithm. Originally, opensees was conceived as a framework, this is apparent from the main page of the wiki. A single rc frame, an ecc frame and an eccrc composite frame were extracted from the reference building respectively and modeled with 2d elements for simplicity.
The type of solution algorithm created and the additional arguments required depends on the algorithmtype. A finite element model fem of frameshear structure was constructed using opensees program based on the nonlinear flexibility theory and multiverticalline theory that considered bendingshear coupling, and its progressive collapse resistance under abnormal conditions was analyzed. How can imply cyclic displacement in opensees software. In this study, the finite element model of polonggou cablestayed bridge is established using opensees software considering multiple nonlinear effects. Contribute to nassermarafiopenseesapi development by creating an account on github. Schellenberg department of civil and environmental engineering university of california, berkeley stephen a. Analytical modeling of the contribution of transverse. To analyse and determine the stresses and displacements, opensees software is used and linked with the codes of genetic algorithm and particle swarm optimization algorithm provided in the matlab software environment. An accelerated newton algorithm based on krylov subspaces is applied to solving nonlinear. Thus, a software framework for experimental testing should ideally. Line search algorithm modified newton algorithm krylovnewton. Dynamic behavior prediction of a reinforced concrete. Open sees command language manual june 2006 plasticity. In short, you can employ them freely assuming you cite the original source and the relevant publication but if you want to build upon, extend or redistribute them, then the derivative software products will also have to be covered under the gpl i.
Evaluation of deformationbased seismic performance of. The combination of all of these features in one software makes stko a revolutionary visualization toolkit for opensees. Although first order reliability method form is effective in calculating. First european conference on opensees opensees day europe 2017 the first european conference on opensees days will take place on 1920 june 2017 at the faculty of engineering of the university of porto, portugal. More options for krylov newton nassermarafiopenseesapi. This class is called acceleratednewton, and is composed of an instance of the accelerator class, as shown in figure 2. The newtonraphson method is the most widely used and most robust method for solving nonlinear algebraic equations.
After the first analysis was complete, the shear load was applied. The resulting natural period of the first mode was 5. A key feature of opensees is the interchangeability of components and the ability to integrate existing libraries and new components into the framework not just new element classes without the need to change the existing code. Threedimensional nonlinear analysis of dynamic soilpile. Flexibilitybased method for modeling shear wall finite element and multiverticalline element sfimvlem. Is has been shown to successfully reach convergence on problems where newton and modified newton. Advanced implementation of hybrid simulation by stephen mahin. Reliabilitybased structural design optimization for. Application to dynamic progressive collapse simulation of frames. Coleman2 1 department of mathematics, tongji university shanghai, china, 200092 2 department of combinatorics and optimization university of waterloo.
Opensees open system for earthquake engineering simulation is an opensource software framework developed at university of california berkeley to simulate the performance of structural and geotechnical systems subjected to earthquakes opensees wiki. This software use the beam model for the masonry walls, named sam method, that is actually considered the most accurate modeling technique of masonry buildings in the literature. Open system for earthquake engineering simulation opensees is a finite element program developed by mckenna and fenves 2001 at university of california, berkeley. This command is used to construct a krylovnewton algorithm object which uses a modified newton method with krylov subspace acceleration to advance to the next time step.
Patran is an excellent preprocessing commercial package for finite element analysis and offers the. In addition, the krylov newton method is selected as the nonlinear solution algorithm since it is proved to be effective. I tried to imply cyclic displacement in top node of concrete frame in opensees software, can anyone help me. Inexact newton methods, newtonkrylov methods, and extensions. The software is designed for parallel computing to allow scalable simulations on highend computers or for parametric studies. Performance of inelastic structures on mitigated and. If it is ease of development, your choices are julia and python, but python is not performant although you could boost its pe. The command algorithm linear construct a linear algorithm object which takes one iteration to solve the system of equations. Parallel scalability and meshindependence for the linear solve.
The goal for opensees is to support a wide range of simulation applications in earthquake engineering. This document is intended to outline the basic commands currently available with the opensees interpreter. If youre not sure which to choose, learn more about installing packages. I tried krylov newton algorithm, and checked that all the node are been connected. Lacking of powerful preprocessor, the preparing of input scripts for opensees is time consuming and errorprone for large and complex models.
Coleman2 1 department of mathematics, tongji university shanghai, china, 200092 2 department of combinatorics and optimization university of waterloo waterloo, on. This command is used to construct a krylovnewton algorithm object which uses a modified newton method with krylov subspace. A software framework is a reusable design for a system or subsystem and defines the overall architecture of such system, meaning its fundamental components as well as the relationships among them. The convergence criterion is the norm of the residual less than. Opensees has more than 160 element types, 220 material types, 15 solution algorithms and 30 solver types. Is your solver using significantly more time or memory than the physics. If your requirement is performance, the best choices are fortran, c and julia. The vibration properties of columns c5 and c6 were similar to the one of column c4.
Opensees performancebased structural engineering for. Opensees has advanced capabilities for modeling and analyzing the nonlinear response of systems using a wide range of material models, elements, and solution algorithms. Bridge designers have adopted simple approximate methods to take into account soilstructureinteraction ssi in dynamic analysis of bridge systems. The software framework opensees open system for earthquake engineering simulation provides a platform for structural and geotechnical engineers.
Linear algorithm uses the solution at the first iteration and continues. However, the shell elements and explicit algorithm in opensees still require further improvements. Overview introduction seismic design according to eurocode 8 description of the software optimization algorithms case studies final comments 3. This command is used to construct a krylovnewton algorithm object which uses a krylov subspace accelerator to accelerate the convergence of the modified newton method. Therefore, in this work, a triangular shell element, nldkgt, and an explicit algorithm are proposed and implemented in opensees. Graphical representation of the first two krylov newton iterations. Opensees, the open system for earthquake engineering simulation, is an objectoriented, open source software framework. Globalized newtonkrylovschwarz algorithms and software for parallel implicit cfd william gropp, david keyes, lois curfman mcinnes, and m. Software framework for parameter updating and finite element response sensitivity analysis. The mpcorecorder is a new recorder class for opensees that allows opensees to write results data in the format required by stko. The software is designed for parallel computing to allow scalable simulations on highend computers or for parameter studies. With its userfriendly interface, you can model complex structures in just minutes and get the analysis results and structural response you want to inspect, quickly and. Linear algorithm uses the solution at the first iteration and continues newton algorithm page 329 uses the tangent at the current iteration to iterate to convergence modified newton algorithm uses the tangent at the first iteration to iterate to convergence newton with line search algorithm krylov newton algorithm bfgs algorithm. Newton algorithm algorithm newton, secantfalse, initialfalse, initialthencurrentfalse.
Static and cyclic triaxial tests were performed to determine the parameters required for the dynamic response analysis of the existing pond ash embankment. Implementation of opensource structural engineering. A krylov subspace accelerated newton algorithm citeseerx. Sep 15, 2014 for linear integrated with the krylovnewton algorithm 37 to solve the non systems, in addition to using inputoutput data for model updating, linear equation of motion in each case with a time step of 0. Newton method in site response analysis not converging in. This command is used to construct a solutionalgorithm object, which. The krylov newton algorithm decomposes uk 1 into two components as follows. Opensees provides beamcolumn and continuum elements for structureal and geotechnical models, and contains a wide range of uniaxial materials and section models. By using stko, a user can create a tool command language tcl input file for opensees and an output file readable by its graphic interface. Krylov subspace accelerated newton algorithm oregon state. Since the frame is elastic the system of equation can be solved in one iteration. The concreted material is adopted to simulate the behavior of concrete, as it is. Contribute to opensees opensees development by creating an account on github.
Can anyone help with a problem with opensee program. The open system for earthquake engineering simulation is a software framework for simulating the static and seismic response of structural and geotechnical systems. This paper presents the dynamic response analysis of talcher pond ash embankment in india, considering both full saturation and existing water table condition subjected to earthquake excitation. Find file copy path opensees src analysis algorithm solutionalgorithm. Hybrid formsampling simulation method for finding design. Stko scientific toolkit for opensees asdea software. Opensees is an objectoriented opensource finite element software platform suitable for advanced simulation of structural and geotechnical systems.
Matlab code for krylov subspace acceleration algorithm. Using the newtonraphson algorithm, convergence is quadratic for both the residual and the displacement increment. Change the line 7 below to set the right path where the openseespy library located. The software has advanced capabilities for modeling and analyzing the nonlinear response of systems using a wide range of material models, elements, and solution algorithms. A leading group of earthquake engineering researchers has produced an opensource program, opensees, to simulate structural behavior under dynamic excitation, incorporating advanced material models and nonlinear dynamic analysis. This solution algorithm tests convergence on the norm of the displacement increment vector with a tolerance set to 1e5 and a maximum number of iterations set to 10. Syllabus lecture date 1 march introduction to opensees and tcl topic introduction to opensees.
Is has been shown to successfully reach convergence on problems where newton and modified. Transformation method is used for boundary condition processing, while sparsesym is utilized for solving dynamic analysis equations 21. Algorithm and particle swarm optimization algorithm are used for truss optimization. Inexact newton methods, newton krylov methods, and extensions for largescale underdetermined systems homer walker doe o ce of advanced scienti c computing. Development of preprocessor for opensees based on patran. This analysis also used a krylov newton solution algorithm. The opensource finite element software, opensees, is widely used in the earthquake engineering community. Newton algorithm page 329 uses the tangent at the current iteration to iterate to.
320 691 206 1120 1538 747 22 896 1062 411 480 298 252 123 327 443 44 782 782 795 1013 798 1021 883 979 901 1531 953 394 1484 1244 1121 602 356 357 1096 730 1476