CN105500712B - 一种三自由度三维打印机的控制方法 - Google Patents
一种三自由度三维打印机的控制方法 Download PDFInfo
- Publication number
- CN105500712B CN105500712B CN201610010142.5A CN201610010142A CN105500712B CN 105500712 B CN105500712 B CN 105500712B CN 201610010142 A CN201610010142 A CN 201610010142A CN 105500712 B CN105500712 B CN 105500712B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mtr
- mtd
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B22—CASTING; POWDER METALLURGY
- B22F—WORKING METALLIC POWDER; MANUFACTURE OF ARTICLES FROM METALLIC POWDER; MAKING METALLIC POWDER; APPARATUS OR DEVICES SPECIALLY ADAPTED FOR METALLIC POWDER
- B22F3/00—Manufacture of workpieces or articles from metallic powder characterised by the manner of compacting or sintering; Apparatus specially adapted therefor ; Presses and furnaces
- B22F3/115—Manufacture of workpieces or articles from metallic powder characterised by the manner of compacting or sintering; Apparatus specially adapted therefor ; Presses and furnaces by spraying molten metal, i.e. spray sintering, spray casting
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Mechanical Engineering (AREA)
Abstract
本发明公开了一种三自由度三维打印机的控制方法，属于机械工具技术领域。具体实施例中的三自由度三维打印机，包括送料装置，所述送料装置包括绕线盘、送丝机构、固定轴，绕线盘活动连接在固定轴上，耗材丝一端缠绕在绕线盘上，另一端与送丝机构相连接，送丝机构包括压线盒、送料轮、压线机构、电机，压线盒设置有凹槽，用以放置送料轮和压线机构，压线盒一端与电机相连接，电机轴穿过压线盒与送料轮固定连接。本发明结构较为简单，可靠性高，利用NURBS曲线进行轨迹规划具有较高的位置精度。
Description
技术领域
本发明涉及一种三自由度三维打印机的控制方法，属于机械工具技术领域。
背景技术
三维打印技术是快速成型制造技术的一种，它是一种以数字模型文件为基础，运用粉末状金属或者塑料等可粘合材料，将一定厚度的材料逐层反复打印在平台上，循环往复，直到生成整个成型件。其工艺工程是先利用三维造型软件创建三维实体造型，再将设计出的实体造型通过快速成型设备的处理软件进行“分区”成逐层的截面切片，然后将处理过的数据输入设备进行制造，最后还需要进行一定的后处理以得到最终的成品。国内开发的三维立体打印快速成型设备一般设置可升降的工作台，喷嘴活动连接在工作台上，根据需要控制喷嘴的运动，喷嘴的位置调整比较复杂，一般通过步进电机牵引同步带，同步带拉动丝杠转动来带动喷嘴移动，这种开环控制加上机构间隙等因素，并不能使喷嘴位置具有较高精度，采用的送丝机构结构也比较复杂，无法有效送丝，不能满足三维打印的要求。针对目前现有技术中存在的上述缺陷，实有必要进行研究，解决现有技术中存在的缺陷。
发明内容
为解决上述问题，本发明的目的在于提供一种结构简单、打印精度高三自由度三维打印机，解决现有技术的三自由度三维打印机送料装置机构复杂，不能准确有效送丝问题，无法满足高精度的三维打印要求。
为实现上述目的，本发明的技术方案为：
一种三自由度三维打印机，包括送料装置，所述送料装置包括绕线盘、送丝机构、固定轴，绕线盘活动连接在固定轴上，耗材丝一端缠绕在绕线盘上，另一端与送丝机构相连接，送丝机构包括压线盒、送料轮、压线机构、电机，压线盒设置有凹槽，用以放置送料轮和压线机构，压线盒一端与电机相连接，电机轴穿过压线盒与送料轮固定连接。压线盒设置有不规则凹槽，根据压线机构的形状设置凹槽的形状，这样凹槽能够有效固定压线机构和送料轮。
进一步地，所述压线机构为月牙形中间位置安装有轴承，压线机构安装在压线盒凹槽内与送料轮保持一定的间隔空隙，使得耗材丝穿过并且对耗材丝形成一定的挤压摩擦力，当送料轮旋转时，送料轮在压线机构的配合下，带动耗材丝往一个方向运动。压线机构通过螺栓安装在压线盒凹槽内，压线机构并不完全固定，能够相对移动，用于调整压线机构与送料轮的间隙。
进一步地，所述压线盒靠近所述压线机构的一端面开有螺栓孔，并安装有弹簧，弹簧一端在压线盒内与压线机构相接触，另一端在压线盒侧面的螺栓孔内并处于压缩状态，压线盒装配有调节螺母、圆垫片，调节螺母从压线盒的螺栓孔内旋进推动圆垫片，圆垫片压缩弹簧，进而弹簧被压缩产生弹性力推动压线机构移动，进行实现调节压线机构的目的。
进一步地，所述送料轮设置有紧定螺钉，通过拧紧紧定螺钉能够使得电机轴与送料轮同步旋转。
进一步地，所述送料轮设置有齿结构，且齿轮的中间有浅凹槽，使得耗材丝容易定位。由于送料轮是圆柱形，耗材丝一般也是圆柱形，二者之间的摩擦力比较小，因此送料轮设置带齿结构并且中间设置有浅凹槽，使得耗材丝能够在浅凹槽中穿过，防止耗材丝滑动。
进一步地，在压线盒靠近送料轮一端设置有辅助豁口，方便调节紧定螺钉。
进一步地，还包括机械臂装置，机械臂装置设置有控制机构，控制机构包括控制芯片和显示屏，控制芯片根据需要打印的位置信息，控制三个电机旋转，进而实现调整末端执行器位姿的目的。利用机械臂控制三自由度三维打印机的喷嘴能够实现三维空间快速准确定位，并且其结构紧凑能够放置在狭窄的空间内。优选地，机械臂装置为三自由度机械臂。
进一步地，还包括打印执行装置，打印执行装置包括连接架、加热机构、散热机构、金属散热管，连接架通过螺钉固定在末端执行器上，同时与加热机构固定连接，机械臂装置带动连接架进行运动。所述末端执行器中间位置设置有通孔，金属散热管一端连接有保护软管并穿过末端执行器的通孔，另一端通过连接螺柱固定在加热机构上。所述连接螺柱为中空外螺纹结构，一端通过螺纹配合连接金属散热管另一端通过螺纹配合连接喷嘴。
一种三自由度三维打印机的控制方法，包括以下步骤：
第10步，采集轮廓离散点信息包括X、Y、X坐标信息，输入到控制芯片；
第20步，控制芯片对打印轮廓利用NURBS曲线进行描述；
第30步，由于三维打印过程，喷嘴需要匀速运动，利用NURBS曲线理论知识，求出对应的下一插补点参数u值，计算公式如下所示：
Ts是采样周期，v喷嘴的运动速度,ui为上一插补点的参数值，X'(u)为NURBS曲线X坐标的一阶导数，Y'(u)为NURBS曲线Y坐标的一阶导数，Z'(u)为NURBS曲线Z坐标的一阶导数，X”(u)为NURBS曲线X坐标的二阶导数，Y”(u)为NURBS曲线Y坐标的二阶导数，Z”(u)为NURBS曲线Z坐标的二阶导数；
第40步，根据计算出的参数u，计算下一插补点坐标，即控制芯片控制打印机喷嘴需要运动到坐标位置，插补点c坐标计算公式如下：
其中：pk+i为控制点、wk+i为权重因子、Ni-p,p为基函数；
第50步，根据NURBS曲线插补算法，得到插补点坐标，控制芯片根据每个采样周期内的插补点坐标计算出需要对每个电机的发送脉冲数，控制电机旋转，实现三维打印机喷嘴的移动。
进一步地，第20步，包括以下步骤：
第201步，控制芯片确定节点向量，已知m+p+1个控制点p为所构造的基函数次数，根据控制点之间的距离远近，计算出合适的节点向量值u，
先计算控制点之间距离dj，公式为：
其中xj表示控制点Pj的横轴坐标，yj表示控制点Pj的纵轴坐标；
哈德利—贾德方法递推公式进行变形，变形后公式为：
其中ui表示第i个节点，进而可得节点向量U递推公式为：
为了使得控制点首末端点和通过点的首末端点相重合，一般节点向量首末节点都取p+1重，求出的节点向量形式为：
第202步，控制芯片根据节点向量计算样条基函数，利用递推公式求基函数：给定参数U轴一个剖分：用下列递推方式确定的基函数称为相应于U的P次，即p+1阶B样条基函数：
ui称为节点，Ni,k(u)的双下标中第一下标i表示序号，第二个下标p表示次数，u为样条参数变量；
第203步，进而可得第k段NURBS曲线Pk(u),表达式为：
u∈[uk+p,uk+p+1],k＝1,2...m+1
其中：pk+i为控制点、wk+i为权重因子、Ni-p,p为基函数；
第204步，反算控制点，由于控制点比轮廓离散点多2个，假设控制点与轮廓离散点的始末端点重合，这样未知量就减少2个，就可以直接利用关系矩阵反算出控制点；
第205步，计算NURBS曲线的导矢:一阶导数，二阶导数。
与现有技术相比，本发明具有以下有益效果：
本发明利用送料轮与压线机构的摩擦力实现耗材丝的输送，结构简单，可靠性高；本发明的压线盒设置有不规则凹槽，根据压线机构的形状设置凹槽的形状，这样凹槽能够有效固定压线机构和送料轮。
本发明利用机械臂控制三自由度三维打印机的喷嘴能够实现三维空间快速准确定位，并且其结构紧凑能够放置在狭窄的空间内。
本发明结构较为简单，可靠性高，在已有的控制算法下，具有较高的位置精度。
附图说明
图1为本发明具体实施例的三自由度三维打印机装配图；
图2为本发明具体实施例的三自由度三维打印机爆炸视图；
图3为本发明具体实施例的三自由度三维打印机局部爆炸视图；
图4为本发明具体实施例的送丝机构装配图；
图5为本发明具体实施例的送丝机构爆炸视图；
图6为本发明具体实施例的控制流程图。
其中：100-送料装置，200-机械臂装置，300-打印执行装置，1-绕线盘，2-固定轴，3-支撑套筒，4-第一轴承，5-保护软管，6-送丝机构，7-支撑竖架，9-末端执行器，10-耗材丝，11-金属散热管，12-第一螺钉，13-连接螺柱，14-连接架，15-第二螺钉，16-加热机构，17-喷嘴，18-散热机构，19-机械臂底座，20-固定底板，601-电机，602-调节螺母，603-弹簧，604-圆垫片，605-辅助进丝件，606-第二轴承，607-压线机构，610-固定螺栓，611-紧定螺钉，612-送料轮，613-压线盒，614-辅助豁口。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
相反，本发明涵盖任何由权利要求定义的在本发明的精髓和范围上做的替代、修改、等效方法以及方案。进一步，为了使公众对本发明有更好的了解，在下文对本发明的细节描述中，详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。
如图1、2所示的具体实施例中，一种三维打印机，包括送料装置100、机械臂装置200、打印执行装置300，送料装置100和机械臂装置200都固定在固定底板20上，打印执行装置300与机械臂装置200的末端执行器9相连接，送料装置100通过保护软管5与打印执行装置300相连接。送料装置100包括绕线盘1、送丝机构6，耗材丝10一端缠绕在绕线盘1上，另一端通过送丝机构6送到打印执行装置300，送丝机构6与打印执行装置300之间连接有保护软管5一方面用以保护耗材丝10，另一方面限制耗材丝10的运动路径。打印执行装置300设置有加热机构16，使得耗材丝10熔化通过喷嘴17进行三维打印，同时在机械臂装置200带动下，打印执行装置300能够灵活的选择打印位置，进行立体成型。
具体实施例中送料装置100包括绕线盘1、送丝机构6、支撑套筒3、第一轴承4、固定轴2，固定底板20一端面固定有支撑竖架7，支撑竖架7上端设置有固定轴2，固定轴2与支撑套筒3通过第一轴承4相连接，使得支撑套筒3一方面能够在固定轴2上旋转另一方面同绕线盘1固定连接，使得绕线盘1能够在固定轴2上旋转。
具体实施例中机械臂装置200包括机械臂本体、控制机构，控制机构包括控制芯片和显示屏，控制芯片根据需要打印的位置信息，控制三个直流伺服电机601旋转，进而实现调整机械臂本体的末端执行器9位姿的目的。
如图3所示具体实施例中打印执行装置300包括连接架14、加热机构16、散热机构18、金属散热管11。连接架14通过第一螺钉12固定在末端执行器9上，同时与加热机构16固定连接，机械臂装置200带动连接架14进行运动。末端执行器9中间位置设置有通孔，金属散热管11一端连接有保护软管5并穿过末端执行器9的通孔，另一端通过连接螺柱13固定在加热机构16上。连接螺柱13为中空外螺纹结构，一端通过螺纹配合连接金属散热管11另一端通过螺纹配合连接喷嘴17。耗材丝10通过保护软管5进入金属散热管11，加热机构16对其加热，使得耗材丝10熔化，通过喷嘴17流出到规划位置，逐步层叠形成三维实体。连接架14为L型结构，外侧有三个螺纹孔，散热机构18通过第二螺钉15固定连接架14外端面上，用于整个机构散热。
如图4、5所示具体实施例中，送丝机构6包括压线盒613、送料轮612、压线机构607、辅助进丝件605、电机601。压线盒613设置有凹槽，用以放置送料轮612和压线机构607，压线盒613一端与电机601相连接，电机601轴穿过压线盒613与送料轮612固定连接，为保证电机601轴与送料轮612的有效连接，送料轮612设置有紧定螺钉611，通过拧紧紧定螺钉611能够使得电机601轴与送料轮612同步旋转，电机601带动送料轮612旋转。送料轮612设置有齿结构，且齿轮的中间有浅凹槽，使得耗材丝10容易定位。
压线机构607为月牙形中间位置安装有轴承，压线机构607固定在压线盒613凹槽内与送料轮612保持一定的间隔空隙，使得耗材丝10穿过并且对耗材丝10形成一定的挤压摩擦力，当送料轮612旋转时，送料轮612在压线机构607的配合下，带动耗材丝10往一个方向运动。压线机构607中间位置安装有第二轴承606，并通过螺栓固定在压线盒613凹槽内，第二轴承606用以同送料轮612相配合。为调节送料轮612与压线机构607的间隙，在压线盒613靠近压线机构607的一端面开有螺栓孔，并安装有弹簧603，弹簧603一端在压线盒613内与压线机构607相接触，另一端在压线盒613侧面的螺栓孔内并处于压缩状态，压线盒613装配有调节螺母602、圆垫片604，调节螺母602从压线盒613的螺栓孔内旋进推动圆垫片604，圆垫片604压缩弹簧603，进而弹簧603被压缩产生弹性力推动压线机构607移动，进行实现调节压线机构607的目的。通过调节螺母602改变送料轮612与压线机构607之间的间隙，来适用不同线径的耗材丝10。同时在压线盒613用以耗材丝10通过的端口处安装有辅助进丝件605，方便耗材丝10通过。由于紧定螺钉611安装在凹槽内，为方便调节紧定螺钉611，在压线盒613靠近送料轮612一端设置有辅助豁口614。
如图1、2、3、4、5所示具体实施例中，电机601带动送料轮612旋转，送料轮612靠摩擦力带动耗材丝10从辅助进丝孔进入，穿过送料轮612与压线机构607的间隙，从另一个辅助进丝孔穿出，进入保护软管5，保护软管5一端连接在送丝机构6的辅助进丝孔另一端与金属散热管11相连接，耗材丝10在送料轮612的推动下，从保护软管5穿过，进入金属散热管11，加热机构16对其加热，耗材丝10在高温条件熔化从喷嘴17喷出。同时机械臂装置200在控制芯片的控制下按照规划方案控制末端执行器9运动，末端执行器9同连接架14连接，喷嘴17固定在连接架14上，进而喷嘴17随着末端执行器9运动而运动，喷嘴17在运动的同时喷出耗材丝10液体，最终实现立体打印、快速成型的目的。
如图6所示的具体实施例中，一种三自由度三维打印机的控制方法，包括以下步骤：
第10步，采集轮廓离散点信息包括X、Y、X坐标信息，输入到控制芯片；
第20步，控制芯片对打印轮廓利用NURBS曲线进行描述；
第30步，由于三维打印过程，喷嘴需要匀速运动，利用NURBS曲线理论知识，求出对应的下一插补点参数u值，计算公式如下所示：
Ts是采样周期，v喷嘴的运动速度,ui为上一插补点的参数值，X'(u)为NURBS曲线X坐标的一阶导数，Y'(u)为NURBS曲线Y坐标的一阶导数，Z'(u)为NURBS曲线Z坐标的一阶导数，X”(u)为NURBS曲线X坐标的二阶导数，Y”(u)为NURBS曲线Y坐标的二阶导数，Z”(u)为NURBS曲线Z坐标的二阶导数；
第40步，根据计算出的参数u，计算下一插补点坐标，即控制芯片控制打印机喷嘴需要运动到坐标位置，插补点c坐标计算公式如下：
其中：pk+i为控制点、wk+i为权重因子、Ni-p,p为基函数；
第50步，根据NURBS曲线插补算法，得到插补点坐标，控制芯片根据每个采样周期内的插补点坐标计算出需要对每个电机的发送脉冲数，控制电机旋转，实现三维打印机喷嘴的移动。
具体实施例中，第20步，进一步包括以下步骤：
第201步，控制芯片确定节点向量，已知m+p+1个控制点p为所构造的基函数次数，根据控制点之间的距离远近，计算出合适的节点向量值u，
先计算控制点之间距离dj，公式为：
其中xj表示控制点Pj的横轴坐标，yj表示控制点Pj的纵轴坐标；
哈德利—贾德方法递推公式进行变形，变形后公式为：
其中ui表示第i个节点，进而可得节点向量U递推公式为：
为了使得控制点首末端点和通过点的首末端点相重合，一般节点向量首末节点都取p+1重，求出的节点向量形式为：
第202步，控制芯片根据节点向量计算样条基函数，利用递推公式求基函数：给定参数U轴一个剖分：用下列递推方式确定的基函数称为相应于U的P次，即p+1阶B样条基函数：
ui称为节点，Ni,k(u)的双下标中第一下标i表示序号，第二个下标p表示次数，u为样条参数变量；
第203步，进而可得第k段NURBS曲线Pk(u),表达式为：
u∈[uk+p,uk+p+1],k＝1,2...m+1
其中：pk+i为控制点、wk+i为权重因子、Ni-p,p为基函数；
第204步，反算控制点，由于控制点比轮廓离散点多2个，假设控制点与轮廓离散点的始末端点重合，这样未知量就减少2个，就可以直接利用关系矩阵反算出控制点；
第205步，计算NURBS曲线的导矢:一阶导数，二阶导数。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (2)
1.一种三自由度三维打印机的控制方法，其特征在于，所述三自由度三维打印机包括送料装置(100)，所述送料装置(100)包括绕线盘(1)、送丝机构(6)、固定轴(2)，绕线盘(1)活动连接在固定轴(2)上，耗材丝(10)一端缠绕在绕线盘(1)上，另一端与送丝机构(6)相连接，送丝机构(6)包括压线盒(613)、送料轮(612)、压线机构(607)、电机(601)，压线盒(613)设置有凹槽，用以放置送料轮(612)和压线机构(607)，压线盒(613)一端与电机(601)相连接，电机(601)轴穿过压线盒(613)与送料轮(612)固定连接，所述三自由度三维打印机的控制方法包括：
第10步，采集轮廓离散点信息包括X、Y、X坐标信息，输入到控制芯片；
第20步，控制芯片对打印轮廓利用NURBS曲线进行描述；
第30步，由于三维打印过程，喷嘴需要匀速运动，利用NURBS曲线理论知识，求出对应的下一插补点参数u值，计算公式如下所示：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mi>T</mi>
<mi>s</mi>
</mrow>
<msqrt>
<mrow>
<msup>
<mi>X</mi>
<mrow>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mrow>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mo>+</mo>
<msup>
<mi>Z</mi>
<mrow>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mi>T</mi>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>X</mi>
<mo>,</mo>
</msup>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>X</mi>
<mrow>
<mo>,</mo>
<mo>,</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mo>,</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>Y</mi>
<mrow>
<mo>,</mo>
<mo>,</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mo>+</mo>
<msup>
<mi>Z</mi>
<mo>,</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>Z</mi>
<mrow>
<mo>,</mo>
<mo>,</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>X</mi>
<mrow>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mrow>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
<mo>+</mo>
<mo>+</mo>
<msup>
<mi>Z</mi>
<mrow>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
Ts是采样周期，v喷嘴的运动速度,ui为上一插补点的参数值，X'(u)为NURBS曲线X坐标的一阶导数，Y'(u)为NURBS曲线Y坐标的一阶导数，Z'(u)为NURBS曲线Z坐标的一阶导数，X”(u)为NURBS曲线X坐标的二阶导数，Y”(u)为NURBS曲线Y坐标的二阶导数，Z”(u)为NURBS曲线Z坐标的二阶导数；
第40步，根据计算出的参数u，计算下一插补点坐标，即控制芯片控制打印机喷嘴需要运动到坐标位置，插补点c坐标计算公式如下：
<mrow>
<mi>c</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>p</mi>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>N</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mi>p</mi>
<mo>,</mo>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>p</mi>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>N</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mi>p</mi>
<mo>,</mo>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中：pk+i为控制点、wk+i为权重因子、Ni-p,p为基函数；
第50步，根据NURBS曲线插补算法，得到插补点坐标，控制芯片根据每个采样周期内的插补点坐标，计算出需要对每个电机的发送脉冲数，控制电机旋转，实现三维打印机喷嘴的移动。
2.如权利要求1所述的三自由度三维打印机的控制方法，其特征在于，第20步，进一步包括以下步骤：
第201步，控制芯片确定节点向量，已知m+p+1个控制点p为所构造的基函数次数，根据控制点之间的距离远近，计算出合适的节点向量值u，
先计算控制点之间距离dj，公式为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3...</mn>
<mi>m</mi>
<mo>+</mo>
<mi>p</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
其中xj表示控制点Pj的横轴坐标，yj表示控制点Pj的纵轴坐标；
哈德利—贾德方法递推公式进行变形，变形后公式为：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mi>i</mi>
<mo>-</mo>
<mi>p</mi>
</mrow>
<mrow>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>d</mi>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mi>i</mi>
<mo>-</mo>
<mi>p</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>d</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>p</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mi>m</mi>
<mo>+</mo>
<mi>p</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中ui表示第i个节点，进而可得节点向量U递推公式为：
<mrow>
<mtable>
<mtr>
<mtd>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>i</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>p</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mi>m</mi>
<mo>+</mo>
<mi>p</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
为了使得控制点首末端点和通过点的首末端点相重合，一般节点向量首末节点都取p+1重，求出的节点向量形式为：
第202步，控制芯片根据节点向量计算样条基函数，利用递推公式求基函数：给定参数U轴一个剖分：用下列递推方式确定的基函数称为相应于U的P次，即p+1阶B样条基函数：
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>u</mi>
<mo>-</mo>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mi>p</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>N</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>p</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mi>u</mi>
</mrow>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<msub>
<mi>N</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>p</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>p</mi>
<mo>&GreaterEqual;</mo>
<mn>2</mn>
</mrow>
ui称为节点，Ni,k(u)的双下标中第一下标i表示序号，第二个下标p表示次数，u为样条参数变量；
第203步，进而可得第k段NURBS曲线Pk(u),表达式为：
<mrow>
<msub>
<mi>P</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>p</mi>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>N</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mi>p</mi>
<mo>,</mo>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>p</mi>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>N</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mi>p</mi>
<mo>,</mo>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
u∈[uk+p,uk+p+1],k＝1,2...m+1
其中：pk+i为控制点、wk+i为权重因子、Ni-p,p为基函数；
第204步，反算控制点，由于控制点比轮廓离散点多2个，假设控制点与轮廓离散点的始末端点重合，这样未知量就减少2个，就可以直接利用关系矩阵反算出控制点；
第205步，计算NURBS曲线的导矢:一阶导数，二阶导数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610010142.5A CN105500712B (zh) | 2016-01-08 | 2016-01-08 | 一种三自由度三维打印机的控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610010142.5A CN105500712B (zh) | 2016-01-08 | 2016-01-08 | 一种三自由度三维打印机的控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105500712A CN105500712A (zh) | 2016-04-20 |
| CN105500712B true CN105500712B (zh) | 2018-01-30 |
Family
ID=55709190
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610010142.5A Active CN105500712B (zh) | 2016-01-08 | 2016-01-08 | 一种三自由度三维打印机的控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105500712B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023133534A1 (en) * | 2022-01-06 | 2023-07-13 | Augmenta Inc. | Techniques for generating composite structures that combine metal and polymer compositions |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107234808B (zh) * | 2017-07-21 | 2023-02-24 | 苏州大学 | 一种多材料喷嘴沉积系统及沉积方法 |
| CN114269670A (zh) | 2019-08-16 | 2022-04-01 | 科思创(荷兰)有限公司 | 用于支撑长丝的线轴和长丝线轴 |
| CN114889128A (zh) * | 2022-04-16 | 2022-08-12 | 大连辰邦三维科技有限公司 | 一种高温打印机喷头 |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN203485449U (zh) * | 2013-09-27 | 2014-03-19 | 东莞市粤船机械技术有限公司 | 一种用于fdm成型技术的3d打印机送丝装置 |
| CN203726795U (zh) * | 2013-09-23 | 2014-07-23 | 南京融会咨策投资管理有限公司 | 一种基于熔融层积成型技术3d打印机的螺杆挤出装置 |
| CN104097326A (zh) * | 2014-07-09 | 2014-10-15 | 西安交通大学 | 一种纤维增强复合材料多自由度3d打印机及其打印方法 |
| CN203876233U (zh) * | 2014-04-18 | 2014-10-15 | 江西环彩三维科技有限公司 | 一种fdm的3d扫描打印装置送料机构 |
| CN104441645A (zh) * | 2013-09-15 | 2015-03-25 | 南京大五教育科技有限公司 | 单喷头熔融挤压式三维打印机 |
| CN104626581A (zh) * | 2015-01-19 | 2015-05-20 | 西安交通大学 | 一种面向太空复杂环境的多自由度3d打印机及打印方法 |
| CN204451227U (zh) * | 2015-02-27 | 2015-07-08 | 三维天工（北京）科技有限公司 | 一种改良的fdm熔融式3d打印机供料组件 |
| CN204658947U (zh) * | 2014-12-24 | 2015-09-23 | 珠海天威飞马打印耗材有限公司 | 检测装置及三维打印机 |
| CN204712470U (zh) * | 2015-06-09 | 2015-10-21 | 珠海天威飞马打印耗材有限公司 | 打印头及三维打印机 |
| CN204936224U (zh) * | 2015-09-08 | 2016-01-06 | 北京理工大学珠海学院 | 一种机械手式彩色3d打印机 |
| CN205326297U (zh) * | 2016-01-08 | 2016-06-22 | 浙江理工大学 | 一种机械臂式三维打印机 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20130094109A (ko) * | 2012-02-15 | 2013-08-23 | 목포대학교산학협력단 | 넙스(ｎｕｒｂｓ)를 이용한 프로펠러 블레이드 형상 설계방법 |
| CN102722141B (zh) * | 2012-06-28 | 2014-05-07 | 华南理工大学 | 一种基于pc-fpga的nurbs曲线实时插补方法 |
| CN103048953B (zh) * | 2013-01-05 | 2014-11-12 | 福建工程学院 | 一种基于arm9嵌入式系统和fpga的nurbs曲线分段插补方法 |
| CN104238457B (zh) * | 2013-06-08 | 2016-12-28 | 沈阳高精数控智能技术股份有限公司 | 一种计算复杂度自适应的nurbs曲线插补方法 |
| CN103345201A (zh) * | 2013-07-01 | 2013-10-09 | 浙江工业大学 | Nurbs曲线降阶插补的运动控制系统 |
| CN104281099B (zh) * | 2014-10-28 | 2017-04-05 | 湘潭大学 | 一种考虑加工特性的nurbs直接插补方法及装置 |
-
2016
- 2016-01-08 CN CN201610010142.5A patent/CN105500712B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104441645A (zh) * | 2013-09-15 | 2015-03-25 | 南京大五教育科技有限公司 | 单喷头熔融挤压式三维打印机 |
| CN203726795U (zh) * | 2013-09-23 | 2014-07-23 | 南京融会咨策投资管理有限公司 | 一种基于熔融层积成型技术3d打印机的螺杆挤出装置 |
| CN203485449U (zh) * | 2013-09-27 | 2014-03-19 | 东莞市粤船机械技术有限公司 | 一种用于fdm成型技术的3d打印机送丝装置 |
| CN203876233U (zh) * | 2014-04-18 | 2014-10-15 | 江西环彩三维科技有限公司 | 一种fdm的3d扫描打印装置送料机构 |
| CN104097326A (zh) * | 2014-07-09 | 2014-10-15 | 西安交通大学 | 一种纤维增强复合材料多自由度3d打印机及其打印方法 |
| CN204658947U (zh) * | 2014-12-24 | 2015-09-23 | 珠海天威飞马打印耗材有限公司 | 检测装置及三维打印机 |
| CN104626581A (zh) * | 2015-01-19 | 2015-05-20 | 西安交通大学 | 一种面向太空复杂环境的多自由度3d打印机及打印方法 |
| CN204451227U (zh) * | 2015-02-27 | 2015-07-08 | 三维天工（北京）科技有限公司 | 一种改良的fdm熔融式3d打印机供料组件 |
| CN204712470U (zh) * | 2015-06-09 | 2015-10-21 | 珠海天威飞马打印耗材有限公司 | 打印头及三维打印机 |
| CN204936224U (zh) * | 2015-09-08 | 2016-01-06 | 北京理工大学珠海学院 | 一种机械手式彩色3d打印机 |
| CN205326297U (zh) * | 2016-01-08 | 2016-06-22 | 浙江理工大学 | 一种机械臂式三维打印机 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023133534A1 (en) * | 2022-01-06 | 2023-07-13 | Augmenta Inc. | Techniques for generating composite structures that combine metal and polymer compositions |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105500712A (zh) | 2016-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105500712B (zh) | 一种三自由度三维打印机的控制方法 | |
| US20210016495A1 (en) | Print assembly for additive manufacturing system, and methods of use thereof | |
| KR100384135B1 (ko) | 선형 열절단 시스템을 이용한 단속적 재료 공급식가변적층 쾌속조형 공정 및 장치 | |
| CN105711095B (zh) | 三维打印用的支撑装置以及三维打印方法 | |
| CN205129792U (zh) | 一种3d打印机工作平台调平装置 | |
| CN103639249B (zh) | 一种板料双点渐进成形装置及成形方法 | |
| CN106003728A (zh) | 用于圆筒壁成型的3d打印机 | |
| CN105414981A (zh) | 一种电弧增材和铣削加工装置 | |
| CN108788446A (zh) | 一种机械臂差速搅拌摩擦焊 | |
| CN205326297U (zh) | 一种机械臂式三维打印机 | |
| CN111002298A (zh) | 一种3d打印机喷头用夹持机械臂 | |
| CN105666867B (zh) | 一种三维打印机以及控制方法 | |
| CN205341922U (zh) | 一种金属熔融快速成型极坐标3d打印设备 | |
| CN205467378U (zh) | 一种三维打印机 | |
| CN110126267B (zh) | 一种基于电涡流场补偿加热的三维打印装置及三维打印头 | |
| CN106608047A (zh) | 一种桌面式3d打印机 | |
| CN204914595U (zh) | 滚珠丝杠驱动式并联臂三维打印机的机架安装结构 | |
| WO2020071933A1 (en) | 3d printer and a method of operating the same | |
| CN204914600U (zh) | 半导体制冷式三维打印机 | |
| CN203679274U (zh) | 带中心孔的金属坯3d打印装置 | |
| CN103752829B (zh) | 带中心孔的金属坯3d打印装置 | |
| CN109351966B (zh) | 一种阶梯孔和斜面的金属坯3d打印装置 | |
| CN104589647A (zh) | 一种3d打印平移台装置 | |
| CN103737003B (zh) | 3d打印金属坯中心孔调节杆 | |
| CN203726816U (zh) | 热缩管扩张真空及外径自动调节装置、热缩管扩张设备 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210526 Address after: 710000 Shaanxi Xi'an economic and Technological Development Zone, Fengcheng four road, Ming Guang road, southeast corner of the new business building 1808 Patentee after: SHAANXI ZHUANYI INTELLECTUAL PROPERTY OPERATION Co.,Ltd. Address before: 310018 No. 2, Xiasha Higher Education Zone, Hangzhou, Zhejiang, Jianggan District Patentee before: Zhejiang University of Technology |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211109 Address after: 257000 Room 308, building 3, Dongying Software Park, No. 228, Nanyi Road, development zone, Dongying City, Shandong Province Patentee after: Dongkai Shuke (Shandong) Industrial Park Co.,Ltd. Address before: 710000 Shaanxi Xi'an economic and Technological Development Zone, Fengcheng four road, Ming Guang road, southeast corner of the new business building 1808 Patentee before: SHAANXI ZHUANYI INTELLECTUAL PROPERTY OPERATION CO.,LTD. |
|
| TR01 | Transfer of patent right |