CNR Author
operator:
and / or
Typology
operator:
and / or
Language
operator:
and / or
Date
operator:
and / or
Rights
operator:
and / or
2012
Contribution to conference
Open Access
Progetto NOSA-ITACA - caso studio: il Voltone, piazza della Repubblica, Livorno
Binante V., Girardi M., Padovani C., Pasquinelli G., Porcelli M., Pellegrini D., Lucchesi M.
NOSA-ITACA project, description of the case studySource: Salone dell'Arte e del Resaturo di Firenze, Firenze, 8-9-10 Novembre 2012
Binante V., Girardi M., Padovani C., Pasquinelli G., Porcelli M., Pellegrini D., Lucchesi M.
NOSA-ITACA project, description of the case studySource: Salone dell'Arte e del Resaturo di Firenze, Firenze, 8-9-10 Novembre 2012
See at:
ISTI Repository 
| CNR ExploRA
2014
Software
Open Access
NOSA-ITACA
Binante V., Girardi M., Padovani C., Pasquinelli G., Pellegrini D., Porcelli M.
NOSA-ITACA is a code for the nonlinear structural analysis of historical masonry constructions. It the result of the integration of the finite element code NOSA into the open-source SALOME platform.
Binante V., Girardi M., Padovani C., Pasquinelli G., Pellegrini D., Porcelli M.
NOSA-ITACA is a code for the nonlinear structural analysis of historical masonry constructions. It the result of the integration of the finite element code NOSA into the open-source SALOME platform.
See at:
ISTI Repository | CNR ExploRA | www.nosaitaca.it
2017
Software
Open Access
NOSA-ITACA 1.1
Binante V., Girardi M., Padovani C., Pasquinelli G., Pellegrini D., Porcelli M., Robol L.
NOSA-ITACA is a code for the nonlinear structural analysis of historical masonry constructions. It the result of the integration of the finite element code NOSA into the open-source SALOME platform.
Binante V., Girardi M., Padovani C., Pasquinelli G., Pellegrini D., Porcelli M., Robol L.
NOSA-ITACA is a code for the nonlinear structural analysis of historical masonry constructions. It the result of the integration of the finite element code NOSA into the open-source SALOME platform.
See at:
ISTI Repository | CNR ExploRA | www.nosaitaca.it
2013
Report
Open Access
Aircraft fuselage sizing with multilevel optimization
Colson B., Porcelli M., Toint P.
In this technical report we describe the activity research carried out in the years 2010-2012 at Department of Mathematics, University of Namur, Namur (Belgium) in collaboration with LMS Samtech, Angleur (Belgium), in the framework of the project "Méthodes de résolution de problèmes d'optimisation de grande taille pour les structures en matériaux composites" (Acronym LARGO "LARge-scale Optimization problems"). LARGO was granted by the Walloon Region and LMS Samtech in the context of the First Program (convention number 916981).Source: ISTI Technical reports, 2013
Colson B., Porcelli M., Toint P.
In this technical report we describe the activity research carried out in the years 2010-2012 at Department of Mathematics, University of Namur, Namur (Belgium) in collaboration with LMS Samtech, Angleur (Belgium), in the framework of the project "Méthodes de résolution de problèmes d'optimisation de grande taille pour les structures en matériaux composites" (Acronym LARGO "LARge-scale Optimization problems"). LARGO was granted by the Walloon Region and LMS Samtech in the context of the First Program (convention number 916981).Source: ISTI Technical reports, 2013
See at:
ISTI Repository | CNR ExploRA
2023
Journal article
Open Access
Numerical solution of a class of quasi-linear matrix equations
Porcelli M., Simoncini V.
Given the matrix equation AX + XB + f (X)C = D in the unknown n * m matrix X, we analyze existence and uniqueness conditions, together with computational solution strategies for [...] being a linear or nonlinear function. We characterize different properties of the matrix equation and of its solution, depending on the considered classes of functions f. Our analysis mainly concerns small dimensional problems, though several considerations also apply to large scale matrix equations.Source: Linear algebra and its applications 664 (2023): 349–368. doi:10.1016/j.laa.2023.01.024
DOI: 10.1016/j.laa.2023.01.024
Metrics:
Porcelli M., Simoncini V.
Given the matrix equation AX + XB + f (X)C = D in the unknown n * m matrix X, we analyze existence and uniqueness conditions, together with computational solution strategies for [...] being a linear or nonlinear function. We characterize different properties of the matrix equation and of its solution, depending on the considered classes of functions f. Our analysis mainly concerns small dimensional problems, though several considerations also apply to large scale matrix equations.Source: Linear algebra and its applications 664 (2023): 349–368. doi:10.1016/j.laa.2023.01.024
DOI: 10.1016/j.laa.2023.01.024
Metrics:
See at:
ISTI Repository | www.sciencedirect.com 
| CNR ExploRA
2012
Contribution to conference
Open Access
Innovative multilevel techniques for structural optimization
Porcelli M., Colson B., Toint P. L.
We address the structural optimization problem of sizing an aircraft fuselage. The problem consists in computing the dimensions of the dierent elements constituting a fuselage mini- mizing the total mass subject to some mechanical constraints. Mathematically, the problem may be formulated as a very large nonlinear optimization problem subjected to several non- linear inequality constraints. We show that this problem possesses a natural hierarchical structure that can be exploited by a multilevel approach. This approach is innovative in the industrial sector and represents a promising alternative to the commonly employed decom- position strategies. Hence, we propose a multilevel procedure which embeds the Recursive Multilevel Trust Region method [1] into an Augmented Lagrangian framework. Some results on both academic and industrial test cases are presented.Source: Congresso Biennale SIMAI MiniSimposio: MSP - Large-Scale Numerical Linear Algebra And Optimization, Torino, 25-28 June 2012
Porcelli M., Colson B., Toint P. L.
We address the structural optimization problem of sizing an aircraft fuselage. The problem consists in computing the dimensions of the dierent elements constituting a fuselage mini- mizing the total mass subject to some mechanical constraints. Mathematically, the problem may be formulated as a very large nonlinear optimization problem subjected to several non- linear inequality constraints. We show that this problem possesses a natural hierarchical structure that can be exploited by a multilevel approach. This approach is innovative in the industrial sector and represents a promising alternative to the commonly employed decom- position strategies. Hence, we propose a multilevel procedure which embeds the Recursive Multilevel Trust Region method [1] into an Augmented Lagrangian framework. Some results on both academic and industrial test cases are presented.Source: Congresso Biennale SIMAI MiniSimposio: MSP - Large-Scale Numerical Linear Algebra And Optimization, Torino, 25-28 June 2012
See at:
ISTI Repository | CNR ExploRA
2012
Contribution to conference
Unknown
Innovative multilevel techniques for structural optimization.
Colson B., Porcelli M., Toint P. L.
We address the structural optimization problem of sizing an aircraft fuselage. The problem consists in computing the dimensions of the dierent elements constituting a fuselage mini- mizing the total mass subject to some mechanical constraints. Mathematically, the problem may be formulated as a very large nonlinear optimization problem subjected to several non- linear inequality constraints. We show that this problem possesses a natural hierarchical structure that can be exploited by a multilevel approach. This approach is innovative in the industrial sector and represents a promising alternative to the commonly employed decom- position strategies. Hence, we propose a multilevel procedure which embeds the Recursive Multilevel Trust Region method [1] into an Augmented Lagrangian framework. Some results on both academic and industrial test cases are presented.Source: Congresso Biennale SIMAI MiniSimposio: MSP - Large-Scale Numerical Linear Algebra And Optimization, Torino, Italy, 25-28 June 2012
Colson B., Porcelli M., Toint P. L.
We address the structural optimization problem of sizing an aircraft fuselage. The problem consists in computing the dimensions of the dierent elements constituting a fuselage mini- mizing the total mass subject to some mechanical constraints. Mathematically, the problem may be formulated as a very large nonlinear optimization problem subjected to several non- linear inequality constraints. We show that this problem possesses a natural hierarchical structure that can be exploited by a multilevel approach. This approach is innovative in the industrial sector and represents a promising alternative to the commonly employed decom- position strategies. Hence, we propose a multilevel procedure which embeds the Recursive Multilevel Trust Region method [1] into an Augmented Lagrangian framework. Some results on both academic and industrial test cases are presented.Source: Congresso Biennale SIMAI MiniSimposio: MSP - Large-Scale Numerical Linear Algebra And Optimization, Torino, Italy, 25-28 June 2012
See at: CNR ExploRA
2013
Report
Unknown
A variable fixing version of the two-block nonlinear constrained Gauss-Seidel algorithm for l1-regularized least-squares
Porcelli M., Rinaldi F.
The problem of finding sparse solutions to underdetermined systems of linear equations is very common in many fields like e.g. signal/image processing and statistics. A standard tool for dealing with sparse recovery is the l1-regularized least-squares approach that has been recently attracting the attention of many researchers. In this paper, we describe a new version of the two-block nonlinear constrained Gauss- Seidel algorithm for solving l1-regularized least-squares that at each step of the iteration process fixes some variables to zero according to a simple rule. We prove the global convergence of the method and we report numerical results on some test problems showing the efficiency of the implemented algorithm.Source: ISTI Technical reports, 2013
Porcelli M., Rinaldi F.
The problem of finding sparse solutions to underdetermined systems of linear equations is very common in many fields like e.g. signal/image processing and statistics. A standard tool for dealing with sparse recovery is the l1-regularized least-squares approach that has been recently attracting the attention of many researchers. In this paper, we describe a new version of the two-block nonlinear constrained Gauss- Seidel algorithm for solving l1-regularized least-squares that at each step of the iteration process fixes some variables to zero according to a simple rule. We prove the global convergence of the method and we report numerical results on some test problems showing the efficiency of the implemented algorithm.Source: ISTI Technical reports, 2013
See at: CNR ExploRA
2013
Contribution to conference
Open Access
New preconditioner updates in Newton-Krylov methods for nonlinear systems
Porcelli M., Bellavia S., Morini B.
We address the problem of preconditioning sequences of large sparse indefinite systems of linear equations arising in the solution of large nonlinear systems via Newton-Krylov methods. We present two new strategies to construct approximate updates of a factorized preconditioner for a reference matrix of the sequence. Both updates are based on the availability of an incomplete factorization for one matrix of the sequence and differ in the approximation of the so-called ideal updates. Furthermore, nearly matrix-free implementations are discussed.Source: 11th EUROPT Workshop on Advances in Continuous Optimization, Firenze, 26-28 Giugno 2013
Porcelli M., Bellavia S., Morini B.
We address the problem of preconditioning sequences of large sparse indefinite systems of linear equations arising in the solution of large nonlinear systems via Newton-Krylov methods. We present two new strategies to construct approximate updates of a factorized preconditioner for a reference matrix of the sequence. Both updates are based on the availability of an incomplete factorization for one matrix of the sequence and differ in the approximation of the so-called ideal updates. Furthermore, nearly matrix-free implementations are discussed.Source: 11th EUROPT Workshop on Advances in Continuous Optimization, Firenze, 26-28 Giugno 2013
See at:
ISTI Repository | CNR ExploRA
2013
Report
Unknown
A solution procedure for constrained eigenvalue problems and its application within the structural finite-element code NOSA-ITACA
Porcelli M., Binante V., Girardi M., Padovani C., Pasquinelli G.
The paper presents an efficient and reliable implementation of numerical methods for constrained generalized eigenvalue problems, specialized for the modal analysis of linear elastic structures in a finite-element setting. The implementation, which takes into account the sparsity of the stiffness and mass matrices and the features of master-slave constraints, is based on open-source packages embedded in the finite-element code NOSA-ITACA. Numerical tests on historical building are performed, with the aims of calculating their vibration frequencies and mode shape vectors, comparing them to the results of a general purpose commercial code and assessing the accuracy of the tool developed.Source: ISTI Technical reports, 2013
Porcelli M., Binante V., Girardi M., Padovani C., Pasquinelli G.
The paper presents an efficient and reliable implementation of numerical methods for constrained generalized eigenvalue problems, specialized for the modal analysis of linear elastic structures in a finite-element setting. The implementation, which takes into account the sparsity of the stiffness and mass matrices and the features of master-slave constraints, is based on open-source packages embedded in the finite-element code NOSA-ITACA. Numerical tests on historical building are performed, with the aims of calculating their vibration frequencies and mode shape vectors, comparing them to the results of a general purpose commercial code and assessing the accuracy of the tool developed.Source: ISTI Technical reports, 2013
See at: CNR ExploRA
2018
Conference article
Open Access
Fea for masonry structures and vibration-based model updating using NOSA-ITACA
Girardi M., Padovani C., Pellegrini D., Porcelli M., Robol L.
NOSA-ITACA is a finite-element code developed by the Mechanics of Materials and Structures Laboratory of ISTI-CNR for the structural analysis of masonry constructions of historical interest via the constitutive equation of masonry-like materials. The latest improvements in the software allow applying model updating techniques to match experimentally measured frequencies in order to fine-tune calculation of the free parameters in the model. The numerical method is briefly presented and applied to two historical buildings in Lucca, the Church of San Francesco and the Clock Tower.Source: 10th International Masonry Conference, pp. 723–735, Milano, Italy, 9-11 July 2018
Girardi M., Padovani C., Pellegrini D., Porcelli M., Robol L.
NOSA-ITACA is a finite-element code developed by the Mechanics of Materials and Structures Laboratory of ISTI-CNR for the structural analysis of masonry constructions of historical interest via the constitutive equation of masonry-like materials. The latest improvements in the software allow applying model updating techniques to match experimentally measured frequencies in order to fine-tune calculation of the free parameters in the model. The numerical method is briefly presented and applied to two historical buildings in Lucca, the Church of San Francesco and the Clock Tower.Source: 10th International Masonry Conference, pp. 723–735, Milano, Italy, 9-11 July 2018
See at:
ISTI Repository | CNR ExploRA
2021
Journal article
Open Access
Solving nonlinear systems of equations via spectral residual methods: stepsize selection and applications
Meli E., Morini B., Porcelli M., Sgattoni C.
Spectral residual methods are derivative-free and low-cost per iteration procedures for solving nonlinear systems of equations. They are generally coupled with a nonmonotone linesearch strategy and compare well with Newton-based methods for large nonlinear systems and sequences of nonlinear systems. The residual vector is used as the search direction and choosing the steplength has a crucial impact on the performance. In this work we address both theoretically and experimentally the steplength selection and provide results on a real application such as a rolling contact problem.Source: Journal of scientific computing 90 (2021). doi:10.1007/s10915-021-01690-x
DOI: 10.1007/s10915-021-01690-x
DOI: 10.48550/arxiv.2005.05851
Metrics:
Meli E., Morini B., Porcelli M., Sgattoni C.
Spectral residual methods are derivative-free and low-cost per iteration procedures for solving nonlinear systems of equations. They are generally coupled with a nonmonotone linesearch strategy and compare well with Newton-based methods for large nonlinear systems and sequences of nonlinear systems. The residual vector is used as the search direction and choosing the steplength has a crucial impact on the performance. In this work we address both theoretically and experimentally the steplength selection and provide results on a real application such as a rolling contact problem.Source: Journal of scientific computing 90 (2021). doi:10.1007/s10915-021-01690-x
DOI: 10.1007/s10915-021-01690-x
DOI: 10.48550/arxiv.2005.05851
Metrics:
See at:
arXiv.org e-Print Archive | Journal of Scientific Computing | ISTI Repository | Journal of Scientific Computing | doi.org | Archivio istituzionale della ricerca - Alma Mater Studiorum Università di Bologna | link.springer.com | CNR ExploRA
2016
Conference article
Open Access
Safety assessment of masonry constructions via numerical tools: the NOSA-ITACA code
Binante V., Girardi M., Lucchesi M., Padovani C., Pellegrini D., Margherita P.
This paper describes the main features of the NOSA-ITACA code, software for the structural analysis of masonry buildings of historical interest resulting from integration of the finite element code NOSA and the open-source platform SALOME. After a short description of the constitutive equation used to model the mechanical behaviour of masonry constructions, some details are given concerning the code's implementation. Then, the result of a static analysis of the "Voltone" in Livorno, performed via the NOSAITACA code, is presented with the aim of highlighting the important role of mathematical models and numerical tools in assessing the safety of historical masonry buildings.Source: 18th ICOMOS General Assembly and Symposium: "Heritage and Landscape as Human Values", pp. 418–424, Florence, Italy, 9-14 November 2014
Binante V., Girardi M., Lucchesi M., Padovani C., Pellegrini D., Margherita P.
This paper describes the main features of the NOSA-ITACA code, software for the structural analysis of masonry buildings of historical interest resulting from integration of the finite element code NOSA and the open-source platform SALOME. After a short description of the constitutive equation used to model the mechanical behaviour of masonry constructions, some details are given concerning the code's implementation. Then, the result of a static analysis of the "Voltone" in Livorno, performed via the NOSAITACA code, is presented with the aim of highlighting the important role of mathematical models and numerical tools in assessing the safety of historical masonry buildings.Source: 18th ICOMOS General Assembly and Symposium: "Heritage and Landscape as Human Values", pp. 418–424, Florence, Italy, 9-14 November 2014
See at:
ISTI Repository | CNR ExploRA
2022
Journal article
Open Access
Exploiting problem structure in derivative free optimization
Porcelli M., Toint P. L.
A structured version of derivative-free random pattern search optimization algorithms is introduced, which is able to exploit coordinate partially separable structure (typically associated with sparsity) often present in unconstrained and bound-constrained optimization problems. This technique improves performance by orders of magnitude and makes it possible to solve large problems that otherwise are totally intractable by other derivative-free methods. A library of interpolation-based modelling tools is also described, which can be associated with the structured or unstructured versions of the initial pattern search algorithm. The use of the library further enhances performance, especially when associated with structure. The significant gains in performance associated with these two techniques are illustrated using a new freely-available release of the Brute Force Optimizer (BFO) package firstly introduced in [Porcelli and Toint 2017], which incorporates them. An interesting conclusion of the numerical results presented is that providing global structural information on a problem can result in significantly less evaluations of the objective function than attempting to building local Taylor-like models.Source: ACM transactions on mathematical software (Online) 48 (2022). doi:10.1145/3474054
DOI: 10.1145/3474054
DOI: 10.48550/arxiv.2001.04801
Metrics:
Porcelli M., Toint P. L.
A structured version of derivative-free random pattern search optimization algorithms is introduced, which is able to exploit coordinate partially separable structure (typically associated with sparsity) often present in unconstrained and bound-constrained optimization problems. This technique improves performance by orders of magnitude and makes it possible to solve large problems that otherwise are totally intractable by other derivative-free methods. A library of interpolation-based modelling tools is also described, which can be associated with the structured or unstructured versions of the initial pattern search algorithm. The use of the library further enhances performance, especially when associated with structure. The significant gains in performance associated with these two techniques are illustrated using a new freely-available release of the Brute Force Optimizer (BFO) package firstly introduced in [Porcelli and Toint 2017], which incorporates them. An interesting conclusion of the numerical results presented is that providing global structural information on a problem can result in significantly less evaluations of the objective function than attempting to building local Taylor-like models.Source: ACM transactions on mathematical software (Online) 48 (2022). doi:10.1145/3474054
DOI: 10.1145/3474054
DOI: 10.48550/arxiv.2001.04801
Metrics:
See at:
arXiv.org e-Print Archive | ACM Transactions on Mathematical Software | ISTI Repository | dl.acm.org | ACM Transactions on Mathematical Software | doi.org | Archivio istituzionale della ricerca - Alma Mater Studiorum Università di Bologna | CNR ExploRA
2014
Journal article
Restricted
New updates of incomplete LU factorizations and applications to large nonlinear systems
Bellavia S., Morini B., Porcelli M.
In this paper, we address the problem of preconditioning sequences of large sparse indefinite systems of linear equations and present two new strategies to construct approximate updates of factorized preconditioners. Both updates are based on the availability of an incomplete factorization for one matrix of the sequence and differ in the approximation of the so-called ideal update. For a general treatment, an incomplete LU (ILU) factorization is considered, but the proposed approaches apply to incomplete factorizations of symmetric matrices as well. The first strategy is an approximate diagonal update of the ILU factorization; the second strategy relies on banded approximations of the factors in the ideal update. The efficiency and reliability of the proposed preconditioners are shown in the solution of nonlinear systems of equations by preconditioned Newton-Krylov methods. Nearly matrix-free implementations of the updating strategy are provided, and numerical experiments are carried out on application problems.Source: Optimization methods & software (Print) 29 (2014): 321–340. doi:10.1080/10556788.2012.762517
DOI: 10.1080/10556788.2012.762517
Metrics:
Bellavia S., Morini B., Porcelli M.
In this paper, we address the problem of preconditioning sequences of large sparse indefinite systems of linear equations and present two new strategies to construct approximate updates of factorized preconditioners. Both updates are based on the availability of an incomplete factorization for one matrix of the sequence and differ in the approximation of the so-called ideal update. For a general treatment, an incomplete LU (ILU) factorization is considered, but the proposed approaches apply to incomplete factorizations of symmetric matrices as well. The first strategy is an approximate diagonal update of the ILU factorization; the second strategy relies on banded approximations of the factors in the ideal update. The efficiency and reliability of the proposed preconditioners are shown in the solution of nonlinear systems of equations by preconditioned Newton-Krylov methods. Nearly matrix-free implementations of the updating strategy are provided, and numerical experiments are carried out on application problems.Source: Optimization methods & software (Print) 29 (2014): 321–340. doi:10.1080/10556788.2012.762517
DOI: 10.1080/10556788.2012.762517
Metrics:
See at:
Optimization Methods and Software | Archivio istituzionale della ricerca - Alma Mater Studiorum Università di Bologna | www.tandfonline.com | CNR ExploRA
2019
Journal article
Open Access
Finite element model updating for structural applications
Girardi M., Padovani C., Pellegrini D., Porcelli M., Robol L.
A novel method for performing model updating on finite element models is presented. The approach is particularly tailored to modal analyses of buildings, by which the lowest frequencies, obtained by using sensors and system identification approaches, need to be matched to the numerical ones predicted by the model. This is done by optimizing some unknown material parameters (such as mass density and Young's modulus) of the materials and/or the boundary conditions, which are often known only approximately. In particular, this is the case when considering historical buildings. The straightforward application of a general-purpose optimizer can be impractical, given the large size of the model involved. In the paper, we show that, by slightly modifying the projection scheme used to compute the eigenvalues at the lowest end of the spectrum one can obtain local parametric reduced order models that, embedded in a trust-region scheme, form the basis for a reliable and efficient specialized algorithm. We describe an optimization strategy based on this approach, and we provide numerical experiments that confirm its effectiveness and accuracy.Source: Journal of computational and applied mathematics 370 (2019). doi:10.1016/j.cam.2019.112675
DOI: 10.1016/j.cam.2019.112675
DOI: 10.48550/arxiv.1801.09122
Metrics:
Girardi M., Padovani C., Pellegrini D., Porcelli M., Robol L.
A novel method for performing model updating on finite element models is presented. The approach is particularly tailored to modal analyses of buildings, by which the lowest frequencies, obtained by using sensors and system identification approaches, need to be matched to the numerical ones predicted by the model. This is done by optimizing some unknown material parameters (such as mass density and Young's modulus) of the materials and/or the boundary conditions, which are often known only approximately. In particular, this is the case when considering historical buildings. The straightforward application of a general-purpose optimizer can be impractical, given the large size of the model involved. In the paper, we show that, by slightly modifying the projection scheme used to compute the eigenvalues at the lowest end of the spectrum one can obtain local parametric reduced order models that, embedded in a trust-region scheme, form the basis for a reliable and efficient specialized algorithm. We describe an optimization strategy based on this approach, and we provide numerical experiments that confirm its effectiveness and accuracy.Source: Journal of computational and applied mathematics 370 (2019). doi:10.1016/j.cam.2019.112675
DOI: 10.1016/j.cam.2019.112675
DOI: 10.48550/arxiv.1801.09122
Metrics:
See at:
arXiv.org e-Print Archive | Journal of Computational and Applied Mathematics | ISTI Repository | Journal of Computational and Applied Mathematics | doi.org | Archivio istituzionale della ricerca - Alma Mater Studiorum Università di Bologna | www.sciencedirect.com | CNR ExploRA
2022
Journal article
Open Access
A semidefinite programming approach for the projection onto the cone of negative semidefinite symmetric tensors with applications to solid mechanics
Padovani C., Porcelli M.
We propose an algorithm for computing the projection of a symmetric second-order tensor onto the cone of negative semidefinite symmetric tensors with respect to the inner product defined by an assigned positive definite symmetric fourth-order tensor C. The projection problem is written as a semidefinite programming problem and an algorithm based on a primal-dual path-following interior point method coupled with a Mehrotra's predictor-corrector approach is proposed. Implementations based on well-known symmetrization schemes and on direct methods are theoretically and numerically investigated taking into account tensors C arising in the modelling of masonry-like materials. For these special cases, indications on the preferable symmetrization scheme that take into account the conditioning of the arising linear systems are given.Source: Calcolo (Online) 59 (2022). doi:10.1007/s10092-022-00478-1
DOI: 10.1007/s10092-022-00478-1
Metrics:
Padovani C., Porcelli M.
We propose an algorithm for computing the projection of a symmetric second-order tensor onto the cone of negative semidefinite symmetric tensors with respect to the inner product defined by an assigned positive definite symmetric fourth-order tensor C. The projection problem is written as a semidefinite programming problem and an algorithm based on a primal-dual path-following interior point method coupled with a Mehrotra's predictor-corrector approach is proposed. Implementations based on well-known symmetrization schemes and on direct methods are theoretically and numerically investigated taking into account tensors C arising in the modelling of masonry-like materials. For these special cases, indications on the preferable symmetrization scheme that take into account the conditioning of the arising linear systems are given.Source: Calcolo (Online) 59 (2022). doi:10.1007/s10092-022-00478-1
DOI: 10.1007/s10092-022-00478-1
Metrics:
See at:
link.springer.com | ISTI Repository | CNR ExploRA
2021
Journal article
Open Access
Improved penalty algorithm for mixed integer PDE constrained optimization problems
Garmatter D., Porcelli M., Rinaldi F., Stoll M.
Optimal control problems including partial differential equation (PDE) as well as integer constraints merge the combinatorial difficulties of integer programming and the challenges related to large-scale systems resulting from discretized PDEs. So far, the branch-and-bound framework has been the most common solution strategy for such problems. In order to provide an alternative solution approach, especially in a large-scale context, this article investigates penalization techniques. Taking inspiration from a well-known family of existing exact penalty algorithms, a novel improved penalty algorithm is derived, whose key ingredients are a basin hopping strategy and an interior point method, both of which are specialized for the problem class. A thorough numerical investigation is carried out for a standard stationary test problem. Extensions to a convection-diffusion as well as a nonlinear test problem finally demonstrate the versatility of the approach.Source: Computers & mathematics with applications (1987) 116 (2021): 2–14. doi:10.1016/j.camwa.2021.11.004
DOI: 10.1016/j.camwa.2021.11.004
DOI: 10.48550/arxiv.1907.06462
Metrics:
Garmatter D., Porcelli M., Rinaldi F., Stoll M.
Optimal control problems including partial differential equation (PDE) as well as integer constraints merge the combinatorial difficulties of integer programming and the challenges related to large-scale systems resulting from discretized PDEs. So far, the branch-and-bound framework has been the most common solution strategy for such problems. In order to provide an alternative solution approach, especially in a large-scale context, this article investigates penalization techniques. Taking inspiration from a well-known family of existing exact penalty algorithms, a novel improved penalty algorithm is derived, whose key ingredients are a basin hopping strategy and an interior point method, both of which are specialized for the problem class. A thorough numerical investigation is carried out for a standard stationary test problem. Extensions to a convection-diffusion as well as a nonlinear test problem finally demonstrate the versatility of the approach.Source: Computers & mathematics with applications (1987) 116 (2021): 2–14. doi:10.1016/j.camwa.2021.11.004
DOI: 10.1016/j.camwa.2021.11.004
DOI: 10.48550/arxiv.1907.06462
Metrics:
See at:
arXiv.org e-Print Archive | Computers & Mathematics with Applications | ISTI Repository | Computers & Mathematics with Applications | doi.org | www.sciencedirect.com | CNR ExploRA
2022
Journal article
Open Access
An improved penalty algorithm using model order reduction for MIPDECO problems with partial observations
Garmatter D., Porcelli M., Rinaldi F., Stoll M.
This work addresses optimal control problems governed by a linear time-dependent partial differential equation (PDE) as well as integer constraints on the control. Moreover, partial observations are assumed in the objective function. The resulting problem poses several numerical challenges due to the mixture of combinatorial aspects, induced by integer variables, and large scale linear algebra issues, arising from the PDE discretization. Since classical solution approaches such as the branch-and-bound framework are typically overwhelmed by such large-scale problems, this work extends an improved penalty algorithm proposed by the authors, to the time-dependent setting. The main contribution is a novel combination of an interior point method, preconditioning, and model order reduction yielding a tailored local optimization solver at the heart of the overall solution procedure. A thorough numerical investigation is carried out both for the heat equation as well as a convection-diffusion problem demonstrating the versatility of the approach.Source: Computational optimization and applications (Dordr., Online) (2022). doi:10.1007/s10589-022-00386-8
DOI: 10.1007/s10589-022-00386-8
Metrics:
Garmatter D., Porcelli M., Rinaldi F., Stoll M.
This work addresses optimal control problems governed by a linear time-dependent partial differential equation (PDE) as well as integer constraints on the control. Moreover, partial observations are assumed in the objective function. The resulting problem poses several numerical challenges due to the mixture of combinatorial aspects, induced by integer variables, and large scale linear algebra issues, arising from the PDE discretization. Since classical solution approaches such as the branch-and-bound framework are typically overwhelmed by such large-scale problems, this work extends an improved penalty algorithm proposed by the authors, to the time-dependent setting. The main contribution is a novel combination of an interior point method, preconditioning, and model order reduction yielding a tailored local optimization solver at the heart of the overall solution procedure. A thorough numerical investigation is carried out both for the heat equation as well as a convection-diffusion problem demonstrating the versatility of the approach.Source: Computational optimization and applications (Dordr., Online) (2022). doi:10.1007/s10589-022-00386-8
DOI: 10.1007/s10589-022-00386-8
Metrics:
See at:
ISTI Repository | rdcu.be | CNR ExploRA
2023
Contribution to book
Open Access
Numerical modelling of historical masonry structures with the finite element code NOSA-ITACA
Girardi M., Padovani C., Pellegrini D., Porcelli M., Robol L.
This chapter presents the finite element code NOSA-ITACA for static and modal analyses of masonry structures of architectural interest. NOSA-ITACA adopts the constitutive equation of masonrylike materials, which considers masonry a non-linear elastic material with zero tensile strength. The capability of modelling restoration and consolidation operations makes the code a helpful tool for maintaining historical buildings. In recent years, long-term vibration monitoring turned out to be an effective non-destructive technique to investigate the dynamic behaviour and check the health status of historical buildings. Changes in their dynamic properties, such as natural frequencies, can represent effective damage indicators. The latest NOSA-ITACA developments are oriented towards structural health monitoring. The availability of the experimental modal properties of a structure makes it possible to calibrate its finite element model via model updating procedures. In particular, the unknown structure's characteristics, such as materials' properties and boundary conditions, can be determined by solving a minimum problem whose objective function is expressed as the discrepancy between experimental frequencies and mode shapes and their numerical counterparts. Several case studies are presented to show the main features of NOSA-ITACA and its effectiveness in the conservation of architectural heritage.Source: Mathematical Modeling in Cultural Heritage. MACH 2021, edited by Bretti G., Cavaterra C., Solci M., Spagnuolo M., pp. 133–152, 2023
DOI: 10.1007/978-981-99-3679-3_9
Metrics:
Girardi M., Padovani C., Pellegrini D., Porcelli M., Robol L.
This chapter presents the finite element code NOSA-ITACA for static and modal analyses of masonry structures of architectural interest. NOSA-ITACA adopts the constitutive equation of masonrylike materials, which considers masonry a non-linear elastic material with zero tensile strength. The capability of modelling restoration and consolidation operations makes the code a helpful tool for maintaining historical buildings. In recent years, long-term vibration monitoring turned out to be an effective non-destructive technique to investigate the dynamic behaviour and check the health status of historical buildings. Changes in their dynamic properties, such as natural frequencies, can represent effective damage indicators. The latest NOSA-ITACA developments are oriented towards structural health monitoring. The availability of the experimental modal properties of a structure makes it possible to calibrate its finite element model via model updating procedures. In particular, the unknown structure's characteristics, such as materials' properties and boundary conditions, can be determined by solving a minimum problem whose objective function is expressed as the discrepancy between experimental frequencies and mode shapes and their numerical counterparts. Several case studies are presented to show the main features of NOSA-ITACA and its effectiveness in the conservation of architectural heritage.Source: Mathematical Modeling in Cultural Heritage. MACH 2021, edited by Bretti G., Cavaterra C., Solci M., Spagnuolo M., pp. 133–152, 2023
DOI: 10.1007/978-981-99-3679-3_9
Metrics:
See at:
ISTI Repository | link.springer.com | CNR ExploRA