Multiphysics investigation of Ice Adhesion over PVC surface
Abstract
This work investigates the ice over a PVC (polyvinyl chloride also known as ‘PVC’) surface as a two-layer laminate model. In this study, ice was frozen over a PVC surface and allowed to adhere. The built samples were tested experimentally in a four-point loading setup. The experimental results contain strain data gathered through data acquisition system using LabView® software. The data was collected at the rate of 1 KHz per load step. This model is analysed theoretically using Euler–Bernoulli beam theory and the rule of mixtures. The correlations from Euler–Bernoulli beam theory and the rule of mixtures were coded in MATLAB® script for theoretical analysis. In addition, numerical simulations were performed using ANSYS® Multiphysics. The FEM model of ice and PVC sample was built using solid elements. The mesh was tested for its sensitivity. Finally theoretical results, experimental results and numerical simulation results were compared. A good agreement between the results was observed.
Experimental Setup
Conclusions
• Euler-Bernoulli beam theory can be used to solve the four-point bending problem. The results give the correlation of displacements with load, longitudinal stress and shear stress.
• Problem containing more than two materials requires Rule of Mixtures. This rule helps to derive common variables from two materials, such as Young’s modulus, moment of inertia, and moment of area.
• A good agreement is found between all of three methodologies.
Contact
H. Khawaja
Assoc. Professor,
University of Tromsø, Tromsø, Norway
E-mail:
[email protected]
Theoretical Analysis
H. Xue, H. Khawaja
UiT-The Arctic University of Norway, Tromsø, Norway
MULTIPHYSICS 2015
PVC Sample with Rosette Strain Gauge PVC and Ice sample placed for four-point bending test Delamination of Ice over PVC surface 𝑃𝑃 Ice
PVC
𝐿𝐿1 𝐿𝐿2 𝐿𝐿1
L
𝑥𝑥 𝑦𝑦
𝑃𝑃 2
M V
𝑉𝑉 = − 𝑃𝑃2 𝑀𝑀 = 𝑃𝑃𝐿𝐿21
𝑃𝑃 2
𝐿𝐿1 𝐿𝐿2 𝐿𝐿1
𝑉𝑉 = 𝑃𝑃2
Four-Point Bending Setup Shear Force and Bending Moment Diagram
Multiphysics Numerical Simulation using ANSYS®
1
X Y Z
AUG 29 2015
17:17:10 ELEMENTS
MAT NUM
1
MN
MX
X Y Z
-.843E-03 -.707E-03 -.571E-03
-.436E-03
-.300E-03
-.164E-03
-.283E-04
.107E-03
.243E-03
.379E-03 SEP 6 2015 00:10:19 NODAL SOLUTION
STEP=1 SUB =1 TIME=1
UY (AVG) RSYS=0
DMX =.843E-03 SMN =-.843E-03 SMX =.379E-03
1
MN MX
X Y Z
-.265E+07 -.212E+07 -.160E+07 -.107E+07 -548636 -23573.6 501489 .103E+07 .155E+07 .208E+07 SEP 5 2015 23:29:54 NODAL SOLUTION
STEP=1 SUB =1 TIME=1
SZ (AVG) RSYS=0
DMX =.001865 SMN =-.265E+07 SMX =.208E+07
FEM Model of Four-Point Bending Displacement Contours Stress Contours