umu.sePublications
Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
On using a zero lower bound on the physical density in material distribution topology optimization
Umeå University, Faculty of Science and Technology, Department of Computing Science.ORCID iD: 0000-0001-9019-2795
Umeå University, Faculty of Science and Technology, Department of Computing Science.
2020 (English)In: Computer Methods in Applied Mechanics and Engineering, ISSN 0045-7825, E-ISSN 1879-2138, Vol. 359, article id 112669Article in journal (Refereed) Published
Abstract [en]

The current paper studies the possibility of allowing a zero lower bound on the physical density in material distribution based topology optimization. We limit our attention to the standard test problem of minimizing the compliance of a linearly elastic structure subject to a constant forcing. First order tensor product Finite Elements discretize the problem. An elementwise constant material indicator function defines the discretized, elementwise constant, physical density by using filtering and penalization. To alleviate the ill-conditioning of the stiffness matrix, due to the variation of the elementwise constant physical density, we precondition the system. We provide a specific spectral analysis for large matrix sizes for the one-dimensional problem with Dirichlet-Neumann conditions in detail, even if most of the mathematical tools apply also in a d-dimensional setting, d >= 2. It is easy to find an elementwise constant material indicator function so that the resulting preconditioned system matrix is singular when allowing the vanishing physical densities. However, for a large class of material indicator functions, the corresponding preconditioned system matrix has a condition number of the same order as the system matrix for the case when the physical density is one in all elements. Finally, we critically report and illustrate results from numerical experiments: as a conclusion, it is indeed possible to solve large-scale topology optimization problems, allowing a vanishing physical density, without encountering ill-conditioned system matrices during the optimization.

Place, publisher, year, edition, pages
Elsevier, 2020. Vol. 359, article id 112669
Keywords [en]
Topology optimization, Conditioning, Preconditioning, Large-scale problems
National Category
Computational Mathematics
Identifiers
URN: urn:nbn:se:umu:diva-167344DOI: 10.1016/j.cma.2019.112669ISI: 000505216600028OAI: oai:DiVA.org:umu-167344DiVA, id: diva2:1391459
Available from: 2020-02-04 Created: 2020-02-04 Last updated: 2020-02-04Bibliographically approved

Open Access in DiVA

No full text in DiVA

Other links

Publisher's full text

Authority records BETA

Wadbro, Eddie

Search in DiVA

By author/editor
Nguyen, Quoc KhanhWadbro, Eddie
By organisation
Department of Computing Science
In the same journal
Computer Methods in Applied Mechanics and Engineering
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 80 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf