编写NPT 1/2内螺纹尺寸的数控程序，可以参考以下步骤和示例代码：
设置坐标系和螺纹参数
```
G90 ； 绝对坐标
G54 ； 坐标系选择
G40 ； 取消半径补偿
G80 ； 取消模态钻孔
```
指定螺纹参数
```
M06 T1 ； 刀具1换刀
G92 S3000 ； 主轴转速设定为每分钟3000转
G94 ； 进给率单位设定为每分钟
```
定义起始点和螺距
```
G00 X50 Y50 ； 将刀具移动到起始点
```
定义螺纹切削参数
NPT 1/2" 内螺纹的参数如下：
大径：21.223 mm
小径：18.321 mm
螺纹深度：18 mm
螺纹高度：1.451 mm
螺纹起点（X坐标）：21.223 - 18/16 = 20.098 mm
螺纹刀定位在Z5，计算R值：
R = （18 + 5） * 1/16 / 2 = 0.719 mm
毛胚孔口直径 = 基面上的小径 = 18.321 mm
毛胚孔底直径 = 基面上的小径 - 孔深 * 1/16 = 18.321 - 1.125 = 17.196 mm
编写螺纹加工程序
```
G0 X15 Z5 ； 螺纹刀定位在Z5
G76 P020260 Q50 R0.05 ； 定义螺纹参数
G76 X20.098 Z-18 R0.719 ； 定义螺纹切削终点
P1451 Q 400 I14 ； 定义螺纹的导程和牙型
```
循环切削直到完成
```
； 循环切削直到完成，重复上述步骤直到所有螺纹部分切削完毕
```
结束螺纹切削
```
G00 Z100 ； 撤离螺纹部分
G80 ； 取消模态钻孔
```
请注意，以上代码示例仅供参考，实际编程时可能需要根据具体的机床型号和编程软件进行调整。在实际操作中，建议先进行充分的测试和验证，以确保生成的螺纹符合设计要求。