Mathematical modeling and experimental studies on axial drilling load for rotary ultrasonic drilling of C/SiC composites