38.8MSMar 20
A Matlab code for analysis and topology optimization with Third Medium ContactAndreas Henrik Frederiksen, Ole Sigmund, Federico Ferrari
We present a Matlab code for modelling and topology optimization of hyperelastic structures, including contact modelled by the Third Medium Contact (TMC) approach. By using the so-called HuHu-regularization we penalize the skew distortion of the bilinear finite elements discretizing void regions, thus promoting convergence of the nonlinear solver. First, we show how this method is implemented in a compact code, allowing to simulate contact and force transfer in hyperelastic structures. We then solve two topology optimization problems for minimum end-compliance of structures exhibiting contact. In the first example, contact happens at the supported boundary, while the second features self-contact. The Matlab scripts that replicate the results are included, and we discuss some possible extensions to more general problems.
MLJun 2, 2024Code
Bayesian Joint Additive Factor Models for Multiview LearningNiccolo Anceschi, Federico Ferrari, David B. Dunson et al.
It is increasingly common to collect data of multiple different types on the same set of samples. Our focus is on studying relationships between such multiview features and responses. A motivating application arises in the context of precision medicine where multi-omics data are collected to correlate with clinical outcomes. It is of interest to infer dependence within and across views while combining multimodal information to improve the prediction of outcomes. The signal-to-noise ratio can vary substantially across views, motivating more nuanced statistical tools beyond standard late and early fusion. This challenge comes with the need to preserve interpretability, select features, and obtain accurate uncertainty quantification. To address these challenges, we introduce two complementary factor regression models. A baseline Joint Factor Regression (\textsc{jfr}) captures combined variation across views via a single factor set, and a more nuanced Joint Additive FActor Regression (\textsc{jafar}) that decomposes variation into shared and view-specific components. For \textsc{jfr}, we use independent cumulative shrinkage process (\textsc{i-cusp}) priors, while for \textsc{jafar} we develop a dependent version (\textsc{d-cusp}) designed to ensure identifiability of the components. We develop Gibbs samplers that exploit the model structure and accommodate flexible feature and outcome distributions. Prediction of time-to-labor onset from immunome, metabolome, and proteome data illustrates performance gains against state-of-the-art competitors. Our open-source software (\texttt{R} package) is available at https://github.com/niccoloanceschi/jafar.
NASep 22, 2018
Revisiting topology optimization with buckling constraintsFederico Ferrari, Ole Sigmund
We review some features of topology optimization with a lower bound on the critical load factor, as computed by linearized buckling analysis. The change of the optimized design, the competition between stiffness and stability requirements and the activation of several buckling modes, depending on the value of such lower bound, are studied. We also discuss some specific issues which are of particular interest for this problem, as the use of non-conforming finite elements for the analysis, the use of inconsistent sensitivities in the optimization and the replacement of the single eigenvalue constraints with an aggregated measure. We discuss the influence of these practices on the optimization result, giving some recommendations.