Full text
Simulations of superelastic lattice materials manufactured by additive manufacturing using a hypoelastic material model M.M. Schasching1, O. ˇ Cervinek2, D. Koutn´y2, H.E. Pettermann1, and M. Todt1 1Institute of Lightweight Design and Structural Biomechanics, TU Wien, Austria 2Institute of Machine and Industrial Design, Brno University of Technology, Technick´a 2896/2, Brno, Czechia The funding of the project ”Building Actions in Smart Aviation with Environmental Gains” by the European Union Programme Horizon Europe under grant agreement no. 101079091 is gratefully acknowledged. Session: S04 Magdeburg, GAMM2024 20.03.2024
1 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Outline Motivation FEM implementation of constitutive law Methodology Applications Summary
1 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Outline Motivation FEM implementation of constitutive law Methodology Applications Summary
2 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Motivation Figure: Taken and modified from [1,2]. •Superelastic lattice material - large deformations while ensuring reversibility •Predictions of large-scale lattice structures 1Deformation mechanisms of internal architecture •Beam-based modeling 2Material response of parent material •Constitutive law to easily take custom-based material data into account [1] Jamshidi, P., et al., Development, characterisation, and modelling of processability of nitinol stents using laser powder bed fusion, J. Alloys Compd., vol. 909, 164681, 2022 [2] Biffi, C.A., et al., Microstructure and Martensitic Transformation of Selective Laser Melted NiTi Shape Memory Alloy Parts, II Int. Conf. on Sim. for Add. Man. - Sim-AM 2019
2 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Motivation Figure: Taken and modified from [1,2]. •Superelastic lattice material - large deformations while ensuring reversibility •Predictions of large-scale lattice structures 1Deformation mechanisms of internal architecture •Beam-based modeling 2Material response of parent material •Constitutive law to easily take custom-based material data into account [1] Jamshidi, P., et al., Development, characterisation, and modelling of processability of nitinol stents using laser powder bed fusion, J. Alloys Compd., vol. 909, 164681, 2022 [2] Biffi, C.A., et al., Microstructure and Martensitic Transformation of Selective Laser Melted NiTi Shape Memory Alloy Parts, II Int. Conf. on Sim. for Add. Man. - Sim-AM 2019
2 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Motivation Figure: Taken and modified from [1,2]. •Superelastic lattice material - large deformations while ensuring reversibility •Predictions of large-scale lattice structures 1Deformation mechanisms of internal architecture •Beam-based modeling 2Material response of parent material •Constitutive law to easily take custom-based material data into account [1] Jamshidi, P., et al., Development, characterisation, and modelling of processability of nitinol stents using laser powder bed fusion, J. Alloys Compd., vol. 909, 164681, 2022 [2] Biffi, C.A., et al., Microstructure and Martensitic Transformation of Selective Laser Melted NiTi Shape Memory Alloy Parts, II Int. Conf. on Sim. for Add. Man. - Sim-AM 2019
2 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Motivation Figure: Taken and modified from [1,2]. •Superelastic lattice material - large deformations while ensuring reversibility •Predictions of large-scale lattice structures 1Deformation mechanisms of internal architecture •Beam-based modeling 2Material response of parent material •Constitutive law to easily take custom-based material data into account [1] Jamshidi, P., et al., Development, characterisation, and modelling of processability of nitinol stents using laser powder bed fusion, J. Alloys Compd., vol. 909, 164681, 2022 [2] Biffi, C.A., et al., Microstructure and Martensitic Transformation of Selective Laser Melted NiTi Shape Memory Alloy Parts, II Int. Conf. on Sim. for Add. Man. - Sim-AM 2019
2 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Motivation Figure: Taken and modified from [1,2]. •Superelastic lattice material - large deformations while ensuring reversibility •Predictions of large-scale lattice structures 1Deformation mechanisms of internal architecture •Beam-based modeling 2Material response of parent material •Constitutive law to easily take custom-based material data into account [1] Jamshidi, P., et al., Development, characterisation, and modelling of processability of nitinol stents using laser powder bed fusion, J. Alloys Compd., vol. 909, 164681, 2022 [2] Biffi, C.A., et al., Microstructure and Martensitic Transformation of Selective Laser Melted NiTi Shape Memory Alloy Parts, II Int. Conf. on Sim. for Add. Man. - Sim-AM 2019
2 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Motivation Figure: Taken and modified from [1,2]. •Superelastic lattice material - large deformations while ensuring reversibility •Predictions of large-scale lattice structures 1Deformation mechanisms of internal architecture •Beam-based modeling 2Material response of parent material •Constitutive law to easily take custom-based material data into account [1] Jamshidi, P., et al., Development, characterisation, and modelling of processability of nitinol stents using laser powder bed fusion, J. Alloys Compd., vol. 909, 164681, 2022 [2] Biffi, C.A., et al., Microstructure and Martensitic Transformation of Selective Laser Melted NiTi Shape Memory Alloy Parts, II Int. Conf. on Sim. for Add. Man. - Sim-AM 2019
3 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB FEM implementation of constitutive law •Hypoelastic user-interface UHYPEL of ABAQUS •Tangential elastic modulus, ET(ε) •Poisson ratio, ν(ε)→ν •Stress-strain curve of experimental data of uniaxial compression tests of additive manufactured NiTi •Assumptions of stress-strain relations •Transformation (loading) and re-transformation processes (unloading) ⇒σ(ε) = aε3+bε2+cε •Fully transformed material ⇒σ(ε) = cε •Coefficients a, b, c obtained by curve fitting •Tangential elastic modulus (conservative systems) ET(ε) = ∂σ ∂ε = 3aε2+ 2bε +c •Stress-strain curve of a general loading sequence
3 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Outline Motivation FEM implementation of constitutive law Methodology Applications Summary
4 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Models and material parameter calibration •Investigation of superelastic lattice structures •Three different models 1bUHYP, beam-based model using UHYPEL 2bASM, beam-based model using ABAQUS’ superelasticity model (ASM) [3,4] 3cASM, high-fidelity model using continuum elements with ASM •Material parameter calibration, coefficients a,b, and cfor bUHYP a b c EtL;T 2.5791426E+7 −1.790820E+6 4.5940E+4 EM,Et;T 6.7159E+4 EtU;T 2.3729021E+7 −1.082987E+6 1.7937E+4 •Input parameter for cASM and bASM EAνAEMνMεLσS tL σE tL σS tU σE tU σS cL 1E+5 3E−1 7E+4 3E−1 3.4E−2 2.51E+2 5.74E+2 2.14E+2 3.6E+1 σS tL [3] F. Auricchio, and R. L. Taylor, Comput. Methods Appl. Mech. Eng., 143(1–2), 175–194 (1997) [4] F. Auricchio, et al., Comput. Methods Appl. Mech. Eng., 146(3–4), 281–312 (1997)
4 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Models and material parameter calibration •Investigation of superelastic lattice structures •Three different models 1bUHYP, beam-based model using UHYPEL 2bASM, beam-based model using ABAQUS’ superelasticity model (ASM) [3,4] 3cASM, high-fidelity model using continuum elements with ASM •Material parameter calibration, coefficients a,b, and cfor bUHYP a b c EtL;T 2.5791426E+7 −1.790820E+6 4.5940E+4 EM,Et;T 6.7159E+4 EtU;T 2.3729021E+7 −1.082987E+6 1.7937E+4 •Input parameter for cASM and bASM EAνAEMνMεLσS tL σE tL σS tU σE tU σS cL 1E+5 3E−1 7E+4 3E−1 3.4E−2 2.51E+2 5.74E+2 2.14E+2 3.6E+1 σS tL [3] F. Auricchio, and R. L. Taylor, Comput. Methods Appl. Mech. Eng., 143(1–2), 175–194 (1997) [4] F. Auricchio, et al., Comput. Methods Appl. Mech. Eng., 146(3–4), 281–312 (1997)
4 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Models and material parameter calibration •Investigation of superelastic lattice structures •Three different models 1bUHYP, beam-based model using UHYPEL 2bASM, beam-based model using ABAQUS’ superelasticity model (ASM) [3,4] 3cASM, high-fidelity model using continuum elements with ASM •Material parameter calibration, coefficients a,b, and cfor bUHYP a b c EtL;T 2.5791426E+7 −1.790820E+6 4.5940E+4 EM,Et;T 6.7159E+4 EtU;T 2.3729021E+7 −1.082987E+6 1.7937E+4 •Input parameter for cASM and bASM EAνAEMνMεLσS tL σE tL σS tU σE tU σS cL 1E+5 3E−1 7E+4 3E−1 3.4E−2 2.51E+2 5.74E+2 2.14E+2 3.6E+1 σS tL [3] F. Auricchio, and R. L. Taylor, Comput. Methods Appl. Mech. Eng., 143(1–2), 175–194 (1997) [4] F. Auricchio, et al., Comput. Methods Appl. Mech. Eng., 146(3–4), 281–312 (1997)
4 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Models and material parameter calibration •Investigation of superelastic lattice structures •Three different models 1bUHYP, beam-based model using UHYPEL 2bASM, beam-based model using ABAQUS’ superelasticity model (ASM) [3,4] 3cASM, high-fidelity model using continuum elements with ASM 0.00 0.01 0.02 0.03 0.04 0.05 εin - 0 200 400 600 800 1000 σin MPa exp. data curve fit loading curve fit unloading martensite input ASM •Material parameter calibration, coefficients a,b, and cfor bUHYP a b c EtL;T 2.5791426E+7 −1.790820E+6 4.5940E+4 EM,Et;T 6.7159E+4 EtU;T 2.3729021E+7 −1.082987E+6 1.7937E+4 •Input parameter for cASM and bASM EAνAEMνMεLσS tL σE tL σS tU σE tU σS cL 1E+5 3E−1 7E+4 3E−1 3.4E−2 2.51E+2 5.74E+2 2.14E+2 3.6E+1 σS tL [3] F. Auricchio, and R. L. Taylor, Comput. Methods Appl. Mech. Eng., 143(1–2), 175–194 (1997) [4] F. Auricchio, et al., Comput. Methods Appl. Mech. Eng., 146(3–4), 281–312 (1997)
5 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Comparison studies •Verification of 1beam-based modeling approach 2bUHYP via a cantilever subjected to pure bending and transverse shear bending •Lattice structures,stiff and compliant lattices •Infinite lattices •Finite lattices •Plane stress assumption for all models •Geometrically nonlinear simulations to account for large deformations using Newton-Raphson solution procedure
5 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Comparison studies •Verification of 1beam-based modeling approach 2bUHYP via a cantilever subjected to pure bending and transverse shear bending •Lattice structures,stiff and compliant lattices •Infinite lattices •Finite lattices •Plane stress assumption for all models •Geometrically nonlinear simulations to account for large deformations using Newton-Raphson solution procedure
5 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Comparison studies •Verification of 1beam-based modeling approach 2bUHYP via a cantilever subjected to pure bending and transverse shear bending •Lattice structures,stiff and compliant lattices •Infinite lattices •Finite lattices •Plane stress assumption for all models •Geometrically nonlinear simulations to account for large deformations using Newton-Raphson solution procedure
5 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Comparison studies •Verification of 1beam-based modeling approach 2bUHYP via a cantilever subjected to pure bending and transverse shear bending •Lattice structures,stiff and compliant lattices •Infinite lattices •Finite lattices •Plane stress assumption for all models •Geometrically nonlinear simulations to account for large deformations using Newton-Raphson solution procedure
6 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Verification - Description •Cantilever subjected to pure bending (left) and transverse shear bending (right) •Dimensions of the cantilever •L1= 10,t= 1, and h= 1 for the out-of-plane thickness. •Pure bending •left end is clamped without constraining the transverse deformations •right end a rotation of UR3=−0.6is prescribed •Transverse shear bending •left end is clamped without constraining the transverse deformations •right end a displacement of U2=−3is applied
7 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Verification - Results •Nodal moment-rotation curve of pure bending (left) and force-displacement curve of transverse shear bending (right) •Equivalent transformation strains of cASM (top) and bASM (bottom) pure bending transverse shear bending ⇒Beam-based modeling approach is suitable
7 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Verification - Results •Nodal moment-rotation curve of pure bending (left) and force-displacement curve of transverse shear bending (right) •Equivalent transformation strains of cASM (top) and bASM (bottom) pure bending transverse shear bending ⇒Beam-based modeling approach is suitable
7 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Verification - Results •Nodal moment-rotation curve of pure bending (left) and force-displacement curve of transverse shear bending (right) •Equivalent transformation strains of cASM (top) and bASM (bottom) pure bending transverse shear bending ⇒Beam-based modeling approach is suitable
8 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Infinite lattices - Description •Infinite lattices via periodic boundary conditions (PBC) NxN Li=L t A =Lh - in mm in mm in mm2 1x1 10 0.5 10 •Compression (i) and pure shear (ii) LCs •Compression •U2=−L/10 is prescribed at node NW •Pure shear •U1=L/10 and U2=L/10 are prescribed at nodes NW and SE, respectively •Bifurcation problem for the stiff lattice is expected
8 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Infinite lattices - Description •Infinite lattices via periodic boundary conditions (PBC) NxN Li=L t A =Lh - in mm in mm in mm2 1x1 10 0.5 10 •Compression (i) and pure shear (ii) LCs •Compression •U2=−L/10 is prescribed at node NW •Pure shear •U1=L/10 and U2=L/10 are prescribed at nodes NW and SE, respectively •Bifurcation problem for the stiff lattice is expected
8 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Infinite lattices - Description •Infinite lattices via periodic boundary conditions (PBC) NxN Li=L t A =Lh - in mm in mm in mm2 1x1 10 0.5 10 •Compression (i) and pure shear (ii) LCs •Compression •U2=−L/10 is prescribed at node NW •Pure shear •U1=L/10 and U2=L/10 are prescribed at nodes NW and SE, respectively •Bifurcation problem for the stiff lattice is expected
8 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Infinite lattices - Description •Infinite lattices via periodic boundary conditions (PBC) NxN Li=L t A =Lh - in mm in mm in mm2 1x1 10 0.5 10 •Compression (i) and pure shear (ii) LCs •Compression •U2=−L/10 is prescribed at node NW •Pure shear •U1=L/10 and U2=L/10 are prescribed at nodes NW and SE, respectively •Bifurcation problem for the stiff lattice is expected
9 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Infinite lattices - Results •Normal stress-strain curves of compression LC (left) and shear stress-strain curves of pure shear LC (right) stiff compliant stiff compliant 0.0 0.2 0.4 −22 in - 0 1 2 3 4 −σ22 in MPa 0.0 0.2 0.4 −22 in - 0 1 2 3 4 0.00 0.05 0.10 γ12 in - 0 20 40 60 σ12 in MPa 0.00 0.05 0.10 γ12 in - 0 5 10 15 0.00 0.01 0.02 0.03 0.04 0.05 12 in - 0 20 40 σ12 in N/mm2 0.00 0.01 0.02 0.03 0.04 0.05 12 in - 0 5 10 15 bASM stiff bUHYP stiff cASM stiff bASM compliant bUHYP compliant cASM compliant •Equivalent transformation strains of cASM for the compression LC (Avg: 75%) TEEQ +0.000e+00 +5.667e−03 +1.133e−02 +1.700e−02 +2.267e−02 +2.833e−02 +3.400e−02 −9.088e−04 +3.675e−02 ⇒Very good agreement of all models
9 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Infinite lattices - Results •Normal stress-strain curves of compression LC (left) and shear stress-strain curves of pure shear LC (right) stiff compliant stiff compliant 0.0 0.2 0.4 −22 in - 0 1 2 3 4 −σ22 in MPa 0.0 0.2 0.4 −22 in - 0 1 2 3 4 0.00 0.05 0.10 γ12 in - 0 20 40 60 σ12 in MPa 0.00 0.05 0.10 γ12 in - 0 5 10 15 0.00 0.01 0.02 0.03 0.04 0.05 12 in - 0 20 40 σ12 in N/mm2 0.00 0.01 0.02 0.03 0.04 0.05 12 in - 0 5 10 15 bASM stiff bUHYP stiff cASM stiff bASM compliant bUHYP compliant cASM compliant •Equivalent transformation strains of cASM for the compression LC (Avg: 75%) TEEQ +0.000e+00 +5.667e−03 +1.133e−02 +1.700e−02 +2.267e−02 +2.833e−02 +3.400e−02 −9.088e−04 +3.675e−02 ⇒Very good agreement of all models
11 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Finite lattices - Results •Force-displacement curves at the load application point •Equivalent transformation strains at maximum loading based on bASM 0 10 20 30 40 U2in mm 0 10 20 30 40 50 RF2in N bASM stiff bUHYP stiff bASM compliant bUHYP compliant (Avg: 75%) TEEQ Bottom +0.000e+00 +5.667e−03 +1.133e−02 +1.700e−02 +2.267e−02 +2.833e−02 +3.400e−02 ⇒More realistic response in the lower loading regime with bUHYP than of bASM
11 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Outline Motivation FEM implementation of constitutive law Methodology Applications Summary
12 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Summary •Uniaxial hypoelastic material model for beam elements via UHYPEL interface of ABAQUS •Verification of beam-based modeling approach and bUHYP •Predictions of the mechanical response of infinite lattices •Beam-based modeling allows us to simulate large-scale (finite) lattices •More realistic response in the lower loading regime with our model than with ASM In combination with much higher numerical efficiency, we provide a tool suitable for large-scale lattice simulations including non-uniform loading scenarios Outlook: Validation of model with experiments of superelastic lattice materials
12 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Summary •Uniaxial hypoelastic material model for beam elements via UHYPEL interface of ABAQUS •Verification of beam-based modeling approach and bUHYP •Predictions of the mechanical response of infinite lattices •Beam-based modeling allows us to simulate large-scale (finite) lattices •More realistic response in the lower loading regime with our model than with ASM In combination with much higher numerical efficiency, we provide a tool suitable for large-scale lattice simulations including non-uniform loading scenarios Outlook: Validation of model with experiments of superelastic lattice materials
12 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Summary •Uniaxial hypoelastic material model for beam elements via UHYPEL interface of ABAQUS •Verification of beam-based modeling approach and bUHYP •Predictions of the mechanical response of infinite lattices •Beam-based modeling allows us to simulate large-scale (finite) lattices •More realistic response in the lower loading regime with our model than with ASM In combination with much higher numerical efficiency, we provide a tool suitable for large-scale lattice simulations including non-uniform loading scenarios Outlook: Validation of model with experiments of superelastic lattice materials
12 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Summary •Uniaxial hypoelastic material model for beam elements via UHYPEL interface of ABAQUS •Verification of beam-based modeling approach and bUHYP •Predictions of the mechanical response of infinite lattices •Beam-based modeling allows us to simulate large-scale (finite) lattices •More realistic response in the lower loading regime with our model than with ASM In combination with much higher numerical efficiency, we provide a tool suitable for large-scale lattice simulations including non-uniform loading scenarios Outlook: Validation of model with experiments of superelastic lattice materials
12 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Summary •Uniaxial hypoelastic material model for beam elements via UHYPEL interface of ABAQUS •Verification of beam-based modeling approach and bUHYP •Predictions of the mechanical response of infinite lattices •Beam-based modeling allows us to simulate large-scale (finite) lattices •More realistic response in the lower loading regime with our model than with ASM In combination with much higher numerical efficiency, we provide a tool suitable for large-scale lattice simulations including non-uniform loading scenarios Outlook: Validation of model with experiments of superelastic lattice materials
12 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Summary •Uniaxial hypoelastic material model for beam elements via UHYPEL interface of ABAQUS •Verification of beam-based modeling approach and bUHYP •Predictions of the mechanical response of infinite lattices •Beam-based modeling allows us to simulate large-scale (finite) lattices •More realistic response in the lower loading regime with our model than with ASM In combination with much higher numerical efficiency, we provide a tool suitable for large-scale lattice simulations including non-uniform loading scenarios Outlook: Validation of model with experiments of superelastic lattice materials
12 / 12 Motivation FEM implementation of constitutive law Methodology Applications Summary ILSB Summary •Uniaxial hypoelastic material model for beam elements via UHYPEL interface of ABAQUS •Verification of beam-based modeling approach and bUHYP •Predictions of the mechanical response of infinite lattices •Beam-based modeling allows us to simulate large-scale (finite) lattices •More realistic response in the lower loading regime with our model than with ASM In combination with much higher numerical efficiency, we provide a tool suitable for large-scale lattice simulations including non-uniform loading scenarios Outlook: Validation of model with experiments of superelastic lattice materials
Simulations of superelastic lattice materials manufactured by additive manufacturing using a hypoelastic material model M.M. Schasching1, O. ˇ Cervinek2, D. Koutn´y2, H.E. Pettermann1, and M. Todt1 1Institute of Lightweight Design and Structural Biomechanics, TU Wien, Austria 2Institute of Machine and Industrial Design, Brno University of Technology, Technick´a 2896/2, Brno, Czechia The funding of the project ”Building Actions in Smart Aviation with Environmental Gains” by the European Union Programme Horizon Europe under grant agreement no. 101079091 is gratefully acknowledged. Session: S04 Magdeburg, GAMM2024 20.03.2024