Corrigendum to “Macroscale Chirality in Planar Elasticity: Predictions of Granular Micromechanics based Extended Micropolar Model Confirmed via Classical Micro Finite Element Simulations.” [Journal of the Mechanics and Physics of Solids, 214 (2026) 106664.]
Other Scholarly Work
Yilmaz, N, Placidi, L, Misra, A. (2026). Corrigendum to “Macroscale Chirality in Planar Elasticity: Predictions of Granular Micromechanics based Extended Micropolar Model Confirmed via Classical Micro Finite Element Simulations.” [Journal of the Mechanics and Physics of Solids, 214 (2026) 106664.]
. JOURNAL OF THE MECHANICS AND PHYSICS OF SOLIDS, 10.1016/j.jmps.2026.106713
Yilmaz, N, Placidi, L, Misra, A. (2026). Corrigendum to “Macroscale Chirality in Planar Elasticity: Predictions of Granular Micromechanics based Extended Micropolar Model Confirmed via Classical Micro Finite Element Simulations.” [Journal of the Mechanics and Physics of Solids, 214 (2026) 106664.]
. JOURNAL OF THE MECHANICS AND PHYSICS OF SOLIDS, 10.1016/j.jmps.2026.106713
The authors regret that in the final published version of the above paper, the numerical values of micropolar continuum model parameters should be twice as large as those given in the original tables 1 and 2. This error are due to certain peculiarities in the implementation of the weak form expression in COMSOL Multiphysics. Such a situation arises in those cases in which the unit of measure of the kinematic variables are not same. In the case considered in this paper, the kinematic variables are the vector displacement field uᵢ(X,t) (in meter) and the pseudo-scalar micro-rotation ψ(X,t) (dimensionless). Therefore, in the strain energy density w (given in Eqs. 18 and 19), the unit of measure of the strains, and consequently, the associated stiffness parameters are different. In COMSOL implementation, one can either create two weak form PDEs (one with two displacement variables and the second with a single micro-rotation variable), each addressing the variables with different units of measure, or a single weak form PDE applying appropriate normalization. When two weak form PDEs implementation is considered in COMSOL, each weak form PDE independently contributes to terms along the boundary. Consequently, the reaction force contribution from the strain energy density associated with each PDEs is considered, and therefore the Lagrange multiplier associated with the boundary conditions measures twice the actual reaction. To avoid this peculiarity for dependent variables having different units of measure in COMSOL, we created one weak form PDE such that all the dependent variables have same unit of measure (in meter). This is achieved by multiplying the microrotation with an arbitrary length (say 1m), such that micro-rotation gradients are dimensionless. Further, stiffness is treated accordingly to achieve the same unit of measure for all stiffness components. With this implementation, the double counting of boundary reactions is avoided. As a result, the reactions given in original Fig. 5 are reduced by half. To obtain the correct reactions, the stiffness parameters must be doubled. The corrected stiffness parameters are given in new Table 1 and 2. The authors would like to apologise for any inconvenience caused.
