Radek Kolman
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Radek Kolman has a total of 94 co-authors. Here are the top co-authors by number of shared publications:
Jiří Plešek
7 shared publications- Kwang-Choon Park
Ann and H. J. Smead Aerospace Engineering Sciences University of Colorado Boulder Colorado USA
7 shared publications - José A. González
Escuela Técnica Superior de Ingeniería Universidad de Sevilla Seville Spain
6 shared publications D. Gabriel
5 shared publicationsJ. Plešek
5 shared publications- Ján Kopačka
Institute of Thermomechanics of the Czech Academy of Sciences, Dolejskova 5, Prague, 18200, Czech Republic
5 shared publications
Latest Publications
- Open Access
- 30 October 2026
- Published by Elsevier in Computer Methods in Applied Mechanics and Engineering
- Vol. 464
Localized Lagrange multipliers with weak enforcement of constraints and Mortar method equivalence
- 30 October 2026
- Published by Elsevier in Computer Methods in Applied Mechanics and Engineering
- Vol. 464
Open AccessLocalized Lagrange multipliers with weak enforcement of constraints and Mortar method equivalence
- Open Access
- 5 October 2026
- Published by Elsevier in Journal of Materials Research and Technology
- Vol. 45, 3931-3942
Fundamental investigation of CuCr1Zr – Mar60 horizontal interface processed by PBF-LB
- 5 October 2026
- Published by Elsevier in Journal of Materials Research and Technology
- Vol. 45, 3931-3942
Open AccessFundamental investigation of CuCr1Zr – Mar60 horizontal interface processed by PBF-LB
- Open Access
Experimental and Numerical Study of the Transient Response of a Cantilever Beam with a Piezoelectric Disc Sensor
Online and real-time sensing and monitoring of the health state of complex structures, such as aircraft and critical components of power stations, are essential aspects of research in dynamics. Several types of sensors are used to capture dynamic responses and monitor changes during the operation of critical parts of complex systems. Piezoelectric (PZ) materials belong to a class of electroactive materials that convert mechanical deformation into an electrical response. For example, PZ ceramics or PVDF foils are employed for online sensing of the time history of mechanical deformation. This paper focuses on the dynamical response of a cantilever beam structure equipped with a glued PZ sensor and combines experimental and modelling approaches to achieve accurate and reliable results. The time history of the normal velocity at a point on the beam surface was recorded with a laser vibrometer during transient vibrations of the beam, triggered by the sudden removal of a mass load at the beam's free end. Simultaneously, the output voltage of the PZ sensor was measured with an electronic device. An elastodynamic model of a cantilever beam coupled with a piezoelectric sensor is introduced, along with its discretization using the finite element method. The mathematical model includes additional terms that enforce a floating-potential boundary condition to maintain a constant charge on one of the sensor's electrodes and is presented in an extended form suitable for sensitivity analysis or parameter identification. The model implementation is validated using a numerical example corresponding to the experimental setup. The computed results show good agreement with the experimental data. Furthermore, values of the Rayleigh damping parameters were identified based on the experimental measurements.
Open AccessExperimental and Numerical Study of the Transient Response of a Cantilever Beam with a Piezoelectric Disc Sensor
Online and real-time sensing and monitoring of the health state of complex structures, such as aircraft and critical components of power stations, are essential aspects of research in dynamics. Several types of sensors are used to capture dynamic responses and monitor changes during the operation of critical parts of complex systems. Piezoelectric (PZ) materials belong to a class of electroactive materials that convert mechanical deformation into an electrical response. For example, PZ ceramics or PVDF foils are employed for online sensing of the time history of mechanical deformation. This paper focuses on the dynamical response of a cantilever beam structure equipped with a glued PZ sensor and combines experimental and modelling approaches to achieve accurate and reliable results. The time history of the normal velocity at a point on the beam surface was recorded with a laser vibrometer during transient vibrations of the beam, triggered by the sudden removal of a mass load at the beam's free end. Simultaneously, the output voltage of the PZ sensor was measured with an electronic device. An elastodynamic model of a cantilever beam coupled with a piezoelectric sensor is introduced, along with its discretization using the finite element method. The mathematical model includes additional terms that enforce a floating-potential boundary condition to maintain a constant charge on one of the sensor's electrodes and is presented in an extended form suitable for sensitivity analysis or parameter identification. The model implementation is validated using a numerical example corresponding to the experimental setup. The computed results show good agreement with the experimental data. Furthermore, values of the Rayleigh damping parameters were identified based on the experimental measurements.
- 15 February 2024
- Published by Wiley in International Journal for Numerical Methods in Engineering
- Vol. 125 (11)
Higher‐order inverse mass matrices for the explicit transient analysis of heterogeneous solids
New methods are presented for the direct computation of higher‐order inverse mass matrices (also called reciprocal mass matrices) that are used for explicit transient finite element analysis. The motivation of this work lies in the need of having appropriate sparse inverse mass matrices, which present the same structure as the consistent mass matrix, preserve the total mass, predict suitable frequency spectrum and dictate sufficiently large critical time step sizes. For an efficient evaluation of the reciprocal mass matrix, the projection matrix should be diagonal. This condition can be satisfied by adopting dual shape functions for the momentum field, generated from the same shape functions used for the displacement field. A theoretically consistent derivation of the inverse mass matrix is based on the three‐field Hamilton principle and requires the projection matrix to be evaluated from the integral of these shape functions. Unfortunately, for higher‐order FE shape functions and serendipity FE elements, the projection matrix is not positive definitive and can not be employed. Therefore, we study several lumping procedures for higher order reciprocal mass matrices considering their effect on total‐mass preserving, frequency spectra and accuracy in explicit transient simulations. The article closes with several numerical examples showing suitability of the direct inverse mass matrix in dynamics.
- 15 February 2024
- Published by Wiley in International Journal for Numerical Methods in Engineering
- Vol. 125 (11)
Higher‐order inverse mass matrices for the explicit transient analysis of heterogeneous solids
New methods are presented for the direct computation of higher‐order inverse mass matrices (also called reciprocal mass matrices) that are used for explicit transient finite element analysis. The motivation of this work lies in the need of having appropriate sparse inverse mass matrices, which present the same structure as the consistent mass matrix, preserve the total mass, predict suitable frequency spectrum and dictate sufficiently large critical time step sizes. For an efficient evaluation of the reciprocal mass matrix, the projection matrix should be diagonal. This condition can be satisfied by adopting dual shape functions for the momentum field, generated from the same shape functions used for the displacement field. A theoretically consistent derivation of the inverse mass matrix is based on the three‐field Hamilton principle and requires the projection matrix to be evaluated from the integral of these shape functions. Unfortunately, for higher‐order FE shape functions and serendipity FE elements, the projection matrix is not positive definitive and can not be employed. Therefore, we study several lumping procedures for higher order reciprocal mass matrices considering their effect on total‐mass preserving, frequency spectra and accuracy in explicit transient simulations. The article closes with several numerical examples showing suitability of the direct inverse mass matrix in dynamics.
- Open Access
Corrosion and Electrochemical Properties of Laser-Shock-Peening-Treated Stainless Steel AISI 304L in VVER Primary Water Environment
Laser Shock Peening (LSP) is a surface treatment technique for metallic materials. It induces plastic deformation at the surface of up to around 1 mm in depth. This process introduces residual stresses that lead to strain hardening, and potentially improvements in fatigue, stress corrosion cracking (SCC) and general corrosion behaviour in many, but not all, corrosive media. In this paper, two specimens made of AISI 304L stainless steel, one LSP-treated and one un-treated, were tested at 280 °C and 8 MPa in VVER (or PWR) primary circuit water chemistry using in situ Electrochemical Impedance Spectroscopy (EIS). This experiment serves to qualify the influence of LSP on the changes in corrosion behaviour in high-temperature, high-density water. The residual stress (RS) measurement of the surface showed a compression RS. Before LSP treatment, RS at the surface was 52.2 MPa in the rolling direction 0°RD and 10.42 MPa in the transverse rolling direction 90°RD. After the treatment, surface RS was −175.27 MPa and −183.51 MPa for Scan and TScan directions, respectively. The effect of compressive RS at the surface was studied and showed an increase in corrosion rate. The analysis of oxide layer by SEM revealed differences between LSP-treated and untreated AISI 304L specimens and their connection to corrosion rates.
Open AccessCorrosion and Electrochemical Properties of Laser-Shock-Peening-Treated Stainless Steel AISI 304L in VVER Primary Water Environment
Laser Shock Peening (LSP) is a surface treatment technique for metallic materials. It induces plastic deformation at the surface of up to around 1 mm in depth. This process introduces residual stresses that lead to strain hardening, and potentially improvements in fatigue, stress corrosion cracking (SCC) and general corrosion behaviour in many, but not all, corrosive media. In this paper, two specimens made of AISI 304L stainless steel, one LSP-treated and one un-treated, were tested at 280 °C and 8 MPa in VVER (or PWR) primary circuit water chemistry using in situ Electrochemical Impedance Spectroscopy (EIS). This experiment serves to qualify the influence of LSP on the changes in corrosion behaviour in high-temperature, high-density water. The residual stress (RS) measurement of the surface showed a compression RS. Before LSP treatment, RS at the surface was 52.2 MPa in the rolling direction 0°RD and 10.42 MPa in the transverse rolling direction 90°RD. After the treatment, surface RS was −175.27 MPa and −183.51 MPa for Scan and TScan directions, respectively. The effect of compressive RS at the surface was studied and showed an increase in corrosion rate. The analysis of oxide layer by SEM revealed differences between LSP-treated and untreated AISI 304L specimens and their connection to corrosion rates.
- 2 April 2022
- Published by Elsevier in International Journal of Mechanical Sciences
- Vol. 222, 107222
Meshfree reduced order model for turbomachinery blade flutter analysis
- 2 April 2022
- Published by Elsevier in International Journal of Mechanical Sciences
- Vol. 222, 107222
Meshfree reduced order model for turbomachinery blade flutter analysis
- Open Access
Dynamic Hardening of AISI 304 Steel at a Wide Range of Strain Rates and Its Application to Shot Peening Simulation
The hardening behavior of AISI 304 steel is investigated at various strain rates, from the quasi-static state to ultra-high strain rates, because it is necessary for numerical simulation of high-speed deformation problems. This kind of testing at a wide range of strain rates has not been yet reported in the literature although it is indispensable for accurate numerical analyses where deformation takes place with a wide spectrum of strain rates. AISI 304 steel is a kind of austenitic stainless steel used in various engineering fields, which does not harden by heat treatment, but by cold working such as shot peening. In order to obtain hardening properties at each strain rate, tensile tests were carried out using a universal testing machine of the INSTRON 5583 for the quasi-static state, a high speed material testing machine of a servo-hydraulic type for intermediate strain rates, and a tensile split-Hopkinson bar for high strain rates with a digital image correlation method. The hardening properties at the ultra-high strain-rate region were obtained from Taylor impact test results by calibration with an experimental–numerical hybrid inverse optimization for reliable extrapolation results. Finally, the hardening flow stress curves were obtained at various strain rates from the quasi-static state to ultra-high strain rates by interpolating the data with the extended Lim–Huh model. The result shows that the yield stress of 759 MPa at the quasi-static state increased to 1429 MPa at a strain rate of , which is about 1.9 times of the quasi-static yield stress. As a demonstration example, the dynamic hardening properties obtained were applied to a shot peening simulation that required hardening curves at a wide range of strain rates from the static state to . The simulation result with the dynamic hardening properties is compared to that with the quasi-static properties. The comparison shows a notable difference in the maximum compressive residual stress by 32%, demonstrating that it is important to consider the dynamic hardening properties in such high strain rate simulation as shot peening for accurate simulation results.Open AccessDynamic Hardening of AISI 304 Steel at a Wide Range of Strain Rates and Its Application to Shot Peening Simulation
The hardening behavior of AISI 304 steel is investigated at various strain rates, from the quasi-static state to ultra-high strain rates, because it is necessary for numerical simulation of high-speed deformation problems. This kind of testing at a wide range of strain rates has not been yet reported in the literature although it is indispensable for accurate numerical analyses where deformation takes place with a wide spectrum of strain rates. AISI 304 steel is a kind of austenitic stainless steel used in various engineering fields, which does not harden by heat treatment, but by cold working such as shot peening. In order to obtain hardening properties at each strain rate, tensile tests were carried out using a universal testing machine of the INSTRON 5583 for the quasi-static state, a high speed material testing machine of a servo-hydraulic type for intermediate strain rates, and a tensile split-Hopkinson bar for high strain rates with a digital image correlation method. The hardening properties at the ultra-high strain-rate region were obtained from Taylor impact test results by calibration with an experimental–numerical hybrid inverse optimization for reliable extrapolation results. Finally, the hardening flow stress curves were obtained at various strain rates from the quasi-static state to ultra-high strain rates by interpolating the data with the extended Lim–Huh model. The result shows that the yield stress of 759 MPa at the quasi-static state increased to 1429 MPa at a strain rate of , which is about 1.9 times of the quasi-static yield stress. As a demonstration example, the dynamic hardening properties obtained were applied to a shot peening simulation that required hardening curves at a wide range of strain rates from the static state to . The simulation result with the dynamic hardening properties is compared to that with the quasi-static properties. The comparison shows a notable difference in the maximum compressive residual stress by 32%, demonstrating that it is important to consider the dynamic hardening properties in such high strain rate simulation as shot peening for accurate simulation results. - 1 July 2021
- Published by Elsevier in Applied Acoustics
- Vol. 178
Experimental, analytical, and numerical study of transient elastic waves from a localized source in an aluminium strip
- 1 July 2021
- Published by Elsevier in Applied Acoustics
- Vol. 178
Experimental, analytical, and numerical study of transient elastic waves from a localized source in an aluminium strip
- 8 May 2021
- Published by Wiley in International Journal for Numerical Methods in Engineering
- Vol. 122 (17) , 4609-4636
Partitioned formulation of contact‐impact problems with stabilized contact constraints and reciprocal mass matrices
This work presents an efficient and accuracy‐improved time explicit solution methodology for the simulation of contact‐impact problems with finite elements. The proposed solution process combines four different existent techniques. First, the contact constraints are modeled by a bipenalty contact‐impact formulation that incorporates stiffness and mass penalties preserving the stability limit of contact‐free problems for efficient explicit time integration. Second, a method of localized Lagrange multipliers is employed, which facilitates the partitioned governing equations for each substructure along with the completely localized contact penalty forces pertaining to each free substructure. Third, a method for the direct construction of sparse inverse mass matrices of the free bodies in contact is combined with the localized Lagrange multipliers approach. Finally, an element‐by‐element mass matrix scaling technique that allows the extension of the time integration step is adopted to improve the overall performance of the algorithm. A judicious synthesis of the four numerical techniques has resulted in an increased stable explicit step‐size that boosts the performance of the bipenalty method for contact problems. Classical contact‐impact numerical examples are used to demonstrate the effectiveness of the proposed methodology.
- 8 May 2021
- Published by Wiley in International Journal for Numerical Methods in Engineering
- Vol. 122 (17) , 4609-4636
Partitioned formulation of contact‐impact problems with stabilized contact constraints and reciprocal mass matrices
This work presents an efficient and accuracy‐improved time explicit solution methodology for the simulation of contact‐impact problems with finite elements. The proposed solution process combines four different existent techniques. First, the contact constraints are modeled by a bipenalty contact‐impact formulation that incorporates stiffness and mass penalties preserving the stability limit of contact‐free problems for efficient explicit time integration. Second, a method of localized Lagrange multipliers is employed, which facilitates the partitioned governing equations for each substructure along with the completely localized contact penalty forces pertaining to each free substructure. Third, a method for the direct construction of sparse inverse mass matrices of the free bodies in contact is combined with the localized Lagrange multipliers approach. Finally, an element‐by‐element mass matrix scaling technique that allows the extension of the time integration step is adopted to improve the overall performance of the algorithm. A judicious synthesis of the four numerical techniques has resulted in an increased stable explicit step‐size that boosts the performance of the bipenalty method for contact problems. Classical contact‐impact numerical examples are used to demonstrate the effectiveness of the proposed methodology.