Sławomir Czarnecki - Academia.edu (original) (raw)

Uploads

Papers by Sławomir Czarnecki

Research paper thumbnail of Constrained Versions of the Free Material Design Methods and Their Applications in 3D Printing

Advances in Structural and Multidisciplinary Optimization, 2017

The paper deals with the free material design and its constrained versions constructed by imposin... more The paper deals with the free material design and its constrained versions constructed by imposing: (a) cubic symmetry (cubic material design, CMD), (b) isotropy with: (b1) independent bulk and shear moduli (isotropic material design, IMD), and (b2) fixed Poisson’s ratio (Young’s modulus design, YMD). In the latter case the Young modulus is the only design variable. The moduli are viewed as non-negative, thus allowing for the appearance of void domains within the design domain. The paper shows that all these methods (CMD, IMD, YMD) reduce to two mutually dual problems:

Research paper thumbnail of Recovery of the Auxetic Microstructures Appearing in the Least Compliant Continuum Two‐Dimensional Bodies

physica status solidi (b), 2020

The article discusses the 2D problem of manufacturability of the minimum compliance designs of th... more The article discusses the 2D problem of manufacturability of the minimum compliance designs of the structural elements made of inhomogeneous materials of local isotropic or square symmetry properties. The available isotropic material design (IMD) and cubic material design (CMD) methods deliver the optimal distribution of the elastic moduli within the design domain. Within the 2D setting, the cubic symmetry reduces to the symmetry of a square. The varying underlying microstructures corresponding to the optimal designs are recovered by matching the values of the optimal moduli with the values of the effective moduli of the representative volume elements (RVEs) computed by the asymptotic homogenization method for periodic media. The shape of the RVE and its internal symmetries are properly selected providing assumed isotropy or symmetry of a square of the homogenized constitutive tensor. The microstructure topology is described by parametric description of single (or several) fibers in...

Research paper thumbnail of Topology optimization of bone using cubic material design and evolutionary methods based on internal remodeling

Mechanics Research Communications, 2018

Research paper thumbnail of The Isotropic and Cubic Material Designs. Recovery of the Underlying Microstructures Appearing in the Least Compliant Continuum Bodies

Materials, 2017

The paper discusses the problem of manufacturability of the minimum compliance designs of the str... more The paper discusses the problem of manufacturability of the minimum compliance designs of the structural elements made of two kinds of inhomogeneous materials: the isotropic and cubic. In both the cases the unit cost of the design is assumed as equal to the trace of the Hooke tensor. The Isotropic Material Design (IMD) delivers the optimal distribution of the bulk and shear moduli within the design domain. The Cubic Material Design (CMD) leads to the optimal material orientation and optimal distribution of the invariant moduli in the body made of the material of cubic symmetry. The present paper proves that the varying underlying microstructures (i.e., the representative volume elements (RVE) constructed of one or two isotropic materials) corresponding to the optimal designs constructed by IMD and CMD methods can be recovered by matching the values of the optimal moduli with the values of the effective moduli of the RVE computed by the theory of homogenization. The CMD method leads to a larger set of results, i.e., the set of pairs of optimal moduli. Moreover, special attention is focused on proper recovery of the microstructures in the auxetic sub-domains of the optimal designs.

Research paper thumbnail of Pareto Optimal Design of Non-Homogeneous Isotropic Material Properties for the Multiple Loading Conditions

physica status solidi (b), 2017

The stiffest structural elements to resist given systems of loads of prescribed weights are optim... more The stiffest structural elements to resist given systems of loads of prescribed weights are optimally formed from an inhomogeneous isotropic material. The products are cut out from the design domain, thus solving the shape design problem simultaneously. The optimal designs reflect the weights of the loads, see the illustrative design for mutually orthogonal loads of different weights.

Research paper thumbnail of Optimal internal architectures of femoral bone based on relaxation by homogenization and isotropic material design

Mechanics Research Communications, 2016

The influence of the loading conditions on the trabecular architecture of a femur is investigated... more The influence of the loading conditions on the trabecular architecture of a femur is investigated by using topology optimization methods. The response of the bone to physiological loads results in changes of the internal architecture of bone, reflected by a modification of internal effective density and mechanical properties. The homogenization based optimization model is developed for predicting optimal bone density distribution, wherein bone tissue is assumed to be a composite material consisting of a mixture of material and void. The homogenization scheme treats the geometric parameters of the microstructures and their orientation as design variables and homogenizes the properties in that microstructure, which is generally anisotropic. The penalization of the optimal material density then leads to a classical optimal structure which consists of regions with bone material and regions without bone material. The IMD (Isotropic Material Design) approach is next applied to determine the optimal elasticity tensor in terms of the bulk and shear moduli as well for the present loading applied to the femoral bone sample. IMD is able to provide both the external shape and topology together with the optimal layout of the isotropic moduli. Both topology optimization methods appear to be complementary. Simulations of the internal bone architecture of the human proximal femur results in a density distribution pattern with good consistency with that of the real bone.

Research paper thumbnail of The emergence of auxetic material as a result of optimal isotropic design

physica status solidi (b), 2015

Research paper thumbnail of A stress-based formulation of the free material design problem with the trace constraint and multiple load conditions

Structural and Multidisciplinary Optimization, 2014

The paper deals with minimization of the weighted sum of compliances related to the load cases ap... more The paper deals with minimization of the weighted sum of compliances related to the load cases applied non-simultaneously. The design variables are all components of the Hooke tensor, subject to the isoperimetric condition bounding the integral of the sum of the Kelvin moduli. This free material design problem is reduced to an equilibrium problem-in two formulations-of an effective body with locking. The stress-based formulation reduces to minimization of an integral of a certain norm of stress fields over the stress fields which equilibrate the given loads. The equivalent displacement-based formulation involves a locking locus defined by using a norm being dual to the previous one. The optimal Hooke tensor is determined by using the stress fields solving the auxiliary locking problem. To make the optimal Hooke tensor positive definite one should consider at least 3 load conditions in the 2D case and not less than 6 load conditions in the 3D case.

Research paper thumbnail of On minimum compliance problems of thin elastic plates of varying thickness

Structural and Multidisciplinary Optimization, 2013

The paper deals with two minimum compliance problems of variable thickness plates subject to an i... more The paper deals with two minimum compliance problems of variable thickness plates subject to an in-plane loading or to a transverse loading. The first of this problem (called also the variable thickness sheet problem) is reduced to the locking material problem in its stress-based setting, thus interrelating the stress-based formulation by Allaire (2002) with the kinematic formulation of Golay and Seppecher (Eur J Mech A Solids 20:631-644, 2001). The second problem concerning the Kirchhoff plates of varying thickness is reduced to a non-convex problem in which the integrand of the minimized functional is the square root of the norm of the density energy expressed in terms of the bending moments. This proves that the problem cannot be interpreted as a problem of equilibrium of a locking material. Both formulations discussed need the numerical treatment in which stresses (bending moments) are the main unknowns.

Research paper thumbnail of The stiffest designs of elastic plates: Vector optimization for two loading conditions

Computer Methods in Applied Mechanics and Engineering, 2011

The paper deals with optimal design of linearly elastic plates of the Kelvin moduli being distrib... more The paper deals with optimal design of linearly elastic plates of the Kelvin moduli being distributed according to a given pattern. The case of two loading conditions is discussed. The optimal plate is characterized by the minimum value of the weighted sum of the compliances corresponding to the two kinds of loads. The problem is reduced to the equilibrium problem of a hyperelastic mixture of properties expressed in terms of two stress fields. The stress-based formulation (P) is rearranged to the displacement-based form (P∗). The latter formulation turns out to be well-posed due to convexity of the relevant potential expressed in terms of strains. Due to monotonicity of the stress–strain relations the problem (P∗) is tractable by the finite element method, using special Newton’s solvers. Exemplary numerical results are presented delivering layouts of variation of elastic characteristics for selected values of the weighting factors corresponding to two kinds of loadings.

Research paper thumbnail of Sandwich plates of minimal compliance

Computer Methods in Applied Mechanics and Engineering, 2008

ABSTRACT The subject of the paper is an optimal choice of material parameters characterizing the ... more ABSTRACT The subject of the paper is an optimal choice of material parameters characterizing the core layer of sandwich plates within the framework of the conventional plate theory in which the core layer is treated as soft in the in-plane direction. The mathematical description is similar to the Hencky–Reissner model of plates with transverse shear deformation. Here, however, the bending stiffnesses and the transverse shear stiffnesses can be designed independently. The present paper deals only with optimal design of the core layer to make the plate compliance minimal. Two core materials are at our disposal, which leads to the ill-posed problem. To consider it one should relax this problem by admitting composite domains and characterize their overall properties by the homogenization formulae. The numerical approach is based on this relaxed formulation thus making it mesh-independent. The equilibrium problem is solved by the DSG3 finite element method. The optimization results are found with using the convergent updating schemes of the COC method.

Research paper thumbnail of Topology optimization in structural mechanics

Bulletin of the Polish Academy of Sciences: Technical Sciences, 2013

Optimization of structural topology, called briefly: topology optimization, is a relatively new b... more Optimization of structural topology, called briefly: topology optimization, is a relatively new branch of structural optimization. Its aim is to create optimal structures, instead of correcting the dimensions or changing the shapes of initial designs. For being able to create the structure, one should have a possibility to handle the members of zero stiffness or admit the material of singular constitutive properties, i.e. void. In the present paper, four fundamental problems of topology optimization are discussed: Michell’s structures, two-material layout problem in light of the relaxation by homogenization theory, optimal shape design and the free material design. Their features are disclosed by presenting results for selected problems concerning the same feasible domain, boundary conditions and applied loading. This discussion provides a short introduction into current topics of topology optimization

Research paper thumbnail of An Explicit Construction of the Underlying Laminated Microstructure of the Least Compliant Elastic Bodies

physica status solidi (b), 2018

The paper discusses the problem of recovery of the microstructure of the least compliant bodies o... more The paper discusses the problem of recovery of the microstructure of the least compliant bodies of non‐homogeneous optimal isotropic properties predicted by the Isotropic Material Design method. The three‐dimensional microstructure is assumed as constructed by a subsequent lamination process in which two isotropic materials of ordered properties are used, the process being repeated three or six times, subsequently. The Hooke tensors closest to the optimal ones are found by making use of the analytical Francfort–Murat formulae for effective moduli of third and sixth rank sequential laminates. The ability to model an auxetic behavior within the subdomains where the optimal Poisson's ratio assumes negative values is also shown.

Research paper thumbnail of Constrained Versions of the Free Material Design Methods and Their Applications in 3D Printing

Advances in Structural and Multidisciplinary Optimization, 2017

The paper deals with the free material design and its constrained versions constructed by imposin... more The paper deals with the free material design and its constrained versions constructed by imposing: (a) cubic symmetry (cubic material design, CMD), (b) isotropy with: (b1) independent bulk and shear moduli (isotropic material design, IMD), and (b2) fixed Poisson’s ratio (Young’s modulus design, YMD). In the latter case the Young modulus is the only design variable. The moduli are viewed as non-negative, thus allowing for the appearance of void domains within the design domain. The paper shows that all these methods (CMD, IMD, YMD) reduce to two mutually dual problems:

Research paper thumbnail of Recovery of the Auxetic Microstructures Appearing in the Least Compliant Continuum Two‐Dimensional Bodies

physica status solidi (b), 2020

The article discusses the 2D problem of manufacturability of the minimum compliance designs of th... more The article discusses the 2D problem of manufacturability of the minimum compliance designs of the structural elements made of inhomogeneous materials of local isotropic or square symmetry properties. The available isotropic material design (IMD) and cubic material design (CMD) methods deliver the optimal distribution of the elastic moduli within the design domain. Within the 2D setting, the cubic symmetry reduces to the symmetry of a square. The varying underlying microstructures corresponding to the optimal designs are recovered by matching the values of the optimal moduli with the values of the effective moduli of the representative volume elements (RVEs) computed by the asymptotic homogenization method for periodic media. The shape of the RVE and its internal symmetries are properly selected providing assumed isotropy or symmetry of a square of the homogenized constitutive tensor. The microstructure topology is described by parametric description of single (or several) fibers in...

Research paper thumbnail of Topology optimization of bone using cubic material design and evolutionary methods based on internal remodeling

Mechanics Research Communications, 2018

Research paper thumbnail of The Isotropic and Cubic Material Designs. Recovery of the Underlying Microstructures Appearing in the Least Compliant Continuum Bodies

Materials, 2017

The paper discusses the problem of manufacturability of the minimum compliance designs of the str... more The paper discusses the problem of manufacturability of the minimum compliance designs of the structural elements made of two kinds of inhomogeneous materials: the isotropic and cubic. In both the cases the unit cost of the design is assumed as equal to the trace of the Hooke tensor. The Isotropic Material Design (IMD) delivers the optimal distribution of the bulk and shear moduli within the design domain. The Cubic Material Design (CMD) leads to the optimal material orientation and optimal distribution of the invariant moduli in the body made of the material of cubic symmetry. The present paper proves that the varying underlying microstructures (i.e., the representative volume elements (RVE) constructed of one or two isotropic materials) corresponding to the optimal designs constructed by IMD and CMD methods can be recovered by matching the values of the optimal moduli with the values of the effective moduli of the RVE computed by the theory of homogenization. The CMD method leads to a larger set of results, i.e., the set of pairs of optimal moduli. Moreover, special attention is focused on proper recovery of the microstructures in the auxetic sub-domains of the optimal designs.

Research paper thumbnail of Pareto Optimal Design of Non-Homogeneous Isotropic Material Properties for the Multiple Loading Conditions

physica status solidi (b), 2017

The stiffest structural elements to resist given systems of loads of prescribed weights are optim... more The stiffest structural elements to resist given systems of loads of prescribed weights are optimally formed from an inhomogeneous isotropic material. The products are cut out from the design domain, thus solving the shape design problem simultaneously. The optimal designs reflect the weights of the loads, see the illustrative design for mutually orthogonal loads of different weights.

Research paper thumbnail of Optimal internal architectures of femoral bone based on relaxation by homogenization and isotropic material design

Mechanics Research Communications, 2016

The influence of the loading conditions on the trabecular architecture of a femur is investigated... more The influence of the loading conditions on the trabecular architecture of a femur is investigated by using topology optimization methods. The response of the bone to physiological loads results in changes of the internal architecture of bone, reflected by a modification of internal effective density and mechanical properties. The homogenization based optimization model is developed for predicting optimal bone density distribution, wherein bone tissue is assumed to be a composite material consisting of a mixture of material and void. The homogenization scheme treats the geometric parameters of the microstructures and their orientation as design variables and homogenizes the properties in that microstructure, which is generally anisotropic. The penalization of the optimal material density then leads to a classical optimal structure which consists of regions with bone material and regions without bone material. The IMD (Isotropic Material Design) approach is next applied to determine the optimal elasticity tensor in terms of the bulk and shear moduli as well for the present loading applied to the femoral bone sample. IMD is able to provide both the external shape and topology together with the optimal layout of the isotropic moduli. Both topology optimization methods appear to be complementary. Simulations of the internal bone architecture of the human proximal femur results in a density distribution pattern with good consistency with that of the real bone.

Research paper thumbnail of The emergence of auxetic material as a result of optimal isotropic design

physica status solidi (b), 2015

Research paper thumbnail of A stress-based formulation of the free material design problem with the trace constraint and multiple load conditions

Structural and Multidisciplinary Optimization, 2014

The paper deals with minimization of the weighted sum of compliances related to the load cases ap... more The paper deals with minimization of the weighted sum of compliances related to the load cases applied non-simultaneously. The design variables are all components of the Hooke tensor, subject to the isoperimetric condition bounding the integral of the sum of the Kelvin moduli. This free material design problem is reduced to an equilibrium problem-in two formulations-of an effective body with locking. The stress-based formulation reduces to minimization of an integral of a certain norm of stress fields over the stress fields which equilibrate the given loads. The equivalent displacement-based formulation involves a locking locus defined by using a norm being dual to the previous one. The optimal Hooke tensor is determined by using the stress fields solving the auxiliary locking problem. To make the optimal Hooke tensor positive definite one should consider at least 3 load conditions in the 2D case and not less than 6 load conditions in the 3D case.

Research paper thumbnail of On minimum compliance problems of thin elastic plates of varying thickness

Structural and Multidisciplinary Optimization, 2013

The paper deals with two minimum compliance problems of variable thickness plates subject to an i... more The paper deals with two minimum compliance problems of variable thickness plates subject to an in-plane loading or to a transverse loading. The first of this problem (called also the variable thickness sheet problem) is reduced to the locking material problem in its stress-based setting, thus interrelating the stress-based formulation by Allaire (2002) with the kinematic formulation of Golay and Seppecher (Eur J Mech A Solids 20:631-644, 2001). The second problem concerning the Kirchhoff plates of varying thickness is reduced to a non-convex problem in which the integrand of the minimized functional is the square root of the norm of the density energy expressed in terms of the bending moments. This proves that the problem cannot be interpreted as a problem of equilibrium of a locking material. Both formulations discussed need the numerical treatment in which stresses (bending moments) are the main unknowns.

Research paper thumbnail of The stiffest designs of elastic plates: Vector optimization for two loading conditions

Computer Methods in Applied Mechanics and Engineering, 2011

The paper deals with optimal design of linearly elastic plates of the Kelvin moduli being distrib... more The paper deals with optimal design of linearly elastic plates of the Kelvin moduli being distributed according to a given pattern. The case of two loading conditions is discussed. The optimal plate is characterized by the minimum value of the weighted sum of the compliances corresponding to the two kinds of loads. The problem is reduced to the equilibrium problem of a hyperelastic mixture of properties expressed in terms of two stress fields. The stress-based formulation (P) is rearranged to the displacement-based form (P∗). The latter formulation turns out to be well-posed due to convexity of the relevant potential expressed in terms of strains. Due to monotonicity of the stress–strain relations the problem (P∗) is tractable by the finite element method, using special Newton’s solvers. Exemplary numerical results are presented delivering layouts of variation of elastic characteristics for selected values of the weighting factors corresponding to two kinds of loadings.

Research paper thumbnail of Sandwich plates of minimal compliance

Computer Methods in Applied Mechanics and Engineering, 2008

ABSTRACT The subject of the paper is an optimal choice of material parameters characterizing the ... more ABSTRACT The subject of the paper is an optimal choice of material parameters characterizing the core layer of sandwich plates within the framework of the conventional plate theory in which the core layer is treated as soft in the in-plane direction. The mathematical description is similar to the Hencky–Reissner model of plates with transverse shear deformation. Here, however, the bending stiffnesses and the transverse shear stiffnesses can be designed independently. The present paper deals only with optimal design of the core layer to make the plate compliance minimal. Two core materials are at our disposal, which leads to the ill-posed problem. To consider it one should relax this problem by admitting composite domains and characterize their overall properties by the homogenization formulae. The numerical approach is based on this relaxed formulation thus making it mesh-independent. The equilibrium problem is solved by the DSG3 finite element method. The optimization results are found with using the convergent updating schemes of the COC method.

Research paper thumbnail of Topology optimization in structural mechanics

Bulletin of the Polish Academy of Sciences: Technical Sciences, 2013

Optimization of structural topology, called briefly: topology optimization, is a relatively new b... more Optimization of structural topology, called briefly: topology optimization, is a relatively new branch of structural optimization. Its aim is to create optimal structures, instead of correcting the dimensions or changing the shapes of initial designs. For being able to create the structure, one should have a possibility to handle the members of zero stiffness or admit the material of singular constitutive properties, i.e. void. In the present paper, four fundamental problems of topology optimization are discussed: Michell’s structures, two-material layout problem in light of the relaxation by homogenization theory, optimal shape design and the free material design. Their features are disclosed by presenting results for selected problems concerning the same feasible domain, boundary conditions and applied loading. This discussion provides a short introduction into current topics of topology optimization

Research paper thumbnail of An Explicit Construction of the Underlying Laminated Microstructure of the Least Compliant Elastic Bodies

physica status solidi (b), 2018

The paper discusses the problem of recovery of the microstructure of the least compliant bodies o... more The paper discusses the problem of recovery of the microstructure of the least compliant bodies of non‐homogeneous optimal isotropic properties predicted by the Isotropic Material Design method. The three‐dimensional microstructure is assumed as constructed by a subsequent lamination process in which two isotropic materials of ordered properties are used, the process being repeated three or six times, subsequently. The Hooke tensors closest to the optimal ones are found by making use of the analytical Francfort–Murat formulae for effective moduli of third and sixth rank sequential laminates. The ability to model an auxetic behavior within the subdomains where the optimal Poisson's ratio assumes negative values is also shown.