We often see warnings such as: "such-and-such a formula is true only if theta is 'measured in radians'." For example: "s = r theta" or "d[sin(theta)]/d(theta) = cos(theta)." Unfortunately, the warning is not correctly (or completely) stated. What is meant is: "theta is the (dimensionless) NUMERICAL VALUE of the angle when the angle is expressed in radians," in other words: theta = (angle of interest)/rad.
In order to avoid confusion, let me use theta to represent an angle (dimension A); and theta* to represents the numerical value of theta when theta is expressed in radians: theta* = theta/rad. The above (correct) formulas are: s = r theta* = r theta/rad and d[sin(theta*)]/d(theta*) = cos(theta*) or d[sin(theta/rad)]/d(theta) = cos(theta/rad). Then it doesn't matter what units are used for (the angle) theta.
We need to distinguish between "geometric" trigonometric functions (ratios of lengths of sides in a right-angled triangle) for which the argument is the angle theta (dimension A) and "mathematical" trigonometric functions: solutions of certain differential equations for which the argument is the independent variable, a (dimensionless) number (dimension: 1). Let me use an initial capital for the former. Then:
Sin(theta) = opp/hyp = sin(theta*) = sin(theta/rad)
and so on for the other trigonometric functions. The length ratios and mathematical functions are themselves dimensionless.
In terms of geometric functions of the (dimensional) angle theta, we have, for example:
d[Sin(theta)]/d(theta) = Cos(theta)/rad
easily proved geometrically, using properties of similar triangles (and s = r theta/rad).
Note the dimensional consistency: dim(LHS) = 1/A; dim(RHS) = 1/A. And, since d[Cos(theta)]/d(theta) = –Sin(theta)/rad, if we define y(theta) = Sin(theta) we find that y is the solution of:
d^2y/d(theta)^2 + y/rad^2 = 0
which, with appropriate "initial" conditions, yields: y(theta) = sin(theta/rad).