CN106513969B - 一种焊接电流实时测量计算方法 - Google Patents
一种焊接电流实时测量计算方法 Download PDFInfo
- Publication number
- CN106513969B CN106513969B CN201611163840.5A CN201611163840A CN106513969B CN 106513969 B CN106513969 B CN 106513969B CN 201611163840 A CN201611163840 A CN 201611163840A CN 106513969 B CN106513969 B CN 106513969B
- Authority
- CN
- China
- Prior art keywords
- current
- value
- angle
- welding
- welding current
- 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
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23K—SOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
- B23K11/00—Resistance welding; Severing by resistance heating
- B23K11/24—Electric supply or control circuits therefor
- B23K11/25—Monitoring devices
- B23K11/251—Monitoring devices using analog means
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Resistance Welding (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及检测技术与信号处理领域,具体涉及一种焊接电流实时测量计算方法,计算步骤如下:以晶闸管触发角为输入参数,通过模拟‑数字采集装置采集焊接电流,搜索电流序列最大值所在的采集时间点;计算一个周期内电流的导通角、功率因数角;分别计算依赖模型部分有效值和不依赖模型部分有效值,最后联合两者求得整体焊接电流有效值。本发明可以大大提高电阻点焊设备操作中的电流有效值计算效率,满足了点焊过程实时控制需求。可为电阻点焊控制系统在线参数控制、在线质量评估、在线安全检测等做出贡献。
Description
技术领域
本发明涉及检测技术与信号处理领域,具体涉及一种焊接电流实时测量计算方法,可应用于单相交流电阻点焊机电流测量计算及类似电路的环境。
背景技术
在单相交流电阻点焊机设备中,流经金属焊片的电流值是一个重要的控制参数和过程变量。目前大部分的电阻焊机都是以电流值作为控制变量输入的,也即根据电流值在调整焊接电源的输入。同时,电流值还可以作为重要的计算参数,用以在线质量估计模型的建立。另外,电流值还可以用以动态电阻值的计算,以用于焊接过程分析以及控制策略设计,焊接安全性操作的在线检测等。这些应用大部分需要用到实时数据计算,而当前的电阻点焊电流值测量,通常使用均方根值来对每一个采集值进行计算,需要采集的实时数据量很多,计算涉及到的数学算法步骤复杂,计算量大,对于以微处理器系统为主的中央处理系统造成很大的压力。
在单相交流电阻焊机中,唯一的输入参数为每一个控制周期的晶闸管触发角α,主要决定了有多少能量进入焊接系统,当能量进入后即在每一个控制周期内产生了持续一定时间的焊接电流,在传统的方法中,是根据平方和开方-也即离散点逐次积分的方法求得,具体是对在数据采集的每一个点,均需要做平方运算,然后相加,最后除以总数N(一个周期内采集到的电流值的个数)以及开平方,工作量很大。通常在以微处理器为主要工具的点焊操作中,计算空间非常有限,因此这样的运算给其他运算以及整体任务的完成造成很大压力,必须用更加有效快速的方法予以改进。
方云等人通过理论分析,建立了不同功率因数角的电流有效值系数神经网络计算模型,进而推算出在1/2周期内估计出电流有效值的计算方法,但是这种方法与具体建立模型的焊接过程有关,使用范围非常有限,而且计算过程也较为复杂。贡亮等人建立了电流过零微分值与电流有效值之间的数学关系,进而同样利用神经网络建立了数学模型,但是神经网络过程过于复杂,且分析过程同样与单一实验环境有很大的关系,导致误差较大,也只能利用在比较简单的应用中。
由于电阻点焊的应用越来越广泛,实时计算越来越复杂,对微处理器的计算速度要求也越来越高,亟需发明一种计算更为简单、计算方面更为通用的电流测量方法。
发明内容
本发明的目的在于克服现有技术的缺点和不足,根据电阻点焊机的电力系统特性,从电流有效值计算的基本原理和方法出发,提供通用性强,计算工作量小的一种焊接电流实时测量计算方法。
本发明的目的通过下述技术方案实现:
(1)以晶闸管触发角α为单相交流电阻点焊机唯一输入参数,控制输入焊接系统能量;
(2)通过模拟-数字采集装置采集焊接电流id;
(3)通过采集到的焊接电流id,搜索电流序列最大值所在的采集时间点γ;
(4)计算一个周期内电流的导通角θ;
(5)根据触发角α,导通角θ与功率因数角之间的数学关系求解功率因数角;
(6)计算依赖模型部分有效值IRMS1;
(7)利用神经网络方法计算不依赖模型部分有效值IRMS2;
(8)计算最终的焊接电流有效值IRMS。
进一步,步骤(4)中的导通角θ按以下公式计算:
式中N为一个控制周期内电流采集值的总个数,fs为电流数据采样频率,Tc为控制周期,在工频为50Hz的情况下,一般为10毫秒。
进一步,步骤(5)中的功率因数角按以下公式计算:
式中Ts为电流数据采样周期,M为到达每一个控制周期内电流峰值时的电流采集值的个数,γ为电流峰值,f为工频频率,即50Hz,ω为角频率。
进一步,步骤(6)中的依赖模型部分有效值IRMS1按以下公式计算:
式中,U为电压和Z为系统等效阻抗。
进一步,步骤(7)中的不依赖模型部分有效值IRMS2计算采用神经网络方法计算可以大大节省运算时间,神经网络架构使用2×5×1形式可以取得较好效果。以P=[α,θ]T为输入值,离线训练获得神经网络输入层与隐层之间的权值WI,隐层与输出层之间的权值WII;输入层与隐层之间的偏离值BI和隐层与输出层之间偏离值BII,控制系统实时计算过程如下:
T1=WI·P+BI (4)
IRMS2=WII×T2+BII (6)
式中T1,T2均为5行1列的矩阵,T1,i,T2,i表示矩阵中每一个具体的元素。
进一步,步骤(8)中的最终的焊接电流有效值IRMS按以下公式计算:
1RMS=IRMS1×IRMS2 (7)
本发明与现有技术相比,具有如下优点和有益效果:
本发明主要根据单相交流电阻点焊机电路系统的基本原理,结合数学积分运算,推理出每一个焊接控制周期有效值计算的公式,将有效值计算分解为依赖模型部分和不依赖与模型部分,对于依赖模型部分采用峰值直接除法的方式,而对于不依赖模型部分采用神经网络建模。更为重要的是,不依赖模型部分对于所有的单相交流电阻点焊操作而言都是一样的,完全可以事先确定好,在具体需要计算的时候再代入计算,增强了本方法的普遍适用性。这样可以实现从原理上直接求解电流有效值,而且计算过程简单,计算量小,可以大大改进当前的依赖于均方根计算逐次积分有效值的方法。
附图说明
图1.本发明一种焊接电流实时测量计算方法流程图;
图2.单相交流电阻点焊机电路系统原理图;
图3.本发明实施例触发角为104.4°时,两种方法电流有效值曲线对比图;
图4.本发明实施例触发角为104.4°时,电流计算误差曲线图;
图5.本发明实施例触发角为109.8°时,两种方法电流有效值曲线对比图;
图6.本发明实施例触发角为109.8°时,电流计算误差曲线图;
图7.本发明实施例触发角为113.4°时,两种方法电流有效值曲线对比图;
图8.本发明实施例触发角为113.4°时,电流计算误差曲线图;
图9.本发明实施例触发角为117°时,两种方法电流有效值曲线对比图;
图10.本发明实施例触发角为117°时,电流计算误差曲线图;
图11.本发明实施例触发角为120.6°时,两种方法电流有效值曲线对比图;
图12.本发明实施例触发角为120.6°时,电流计算误差曲线图。
附图说明:1-焊接电源;2-晶闸管1;3-晶闸管2;4-焊接变压器;5-次级线圈的等效电感;6-次级线圈的等效电阻;7-负载;8-主线圈的等效电感;9-主线圈的等效电阻。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅说明书图1:
图1为本发明一种焊接电流实时测量计算方法的流程图,具体包括如下步骤:
(1)以晶闸管触发角α为单相交流电阻点焊机唯一输入参数,控制输入焊接系统能量;
(2)通过模拟-数字采集装置采集焊接电流id;
(3)通过采集到的焊接电流id,搜索电流序列最大值所在的采集时间点γ;
(4)计算一个周期内电流的导通角θ;
(5)根据触发角α,导通角θ与功率因数角之间的数学关系求解功率因数角;
(6)计算依赖模型部分有效值IRMS1;
(7)利用神经网络方法计算不依赖模型部分有效值IRMS2;
(8)计算最终的焊接电流有效值IRMS。
请参阅说明书图2:
在图2单相交流电阻点焊机电路系统原理图中,焊接电源1通过晶闸管2、3的触发角α为系统输入参数,控制输入能量;焊接变压器4实现降压,之后提供能量给负载7,焊接电流id为流经负载7的电流数字采集值。5为次级线圈电路的等效电感;6为次级线圈电路的等效电阻;8为主线圈电路的等效电感;9为主线圈电路的等效电阻。
本发明以普通单相交流电阻点焊操作为例,实验平台包括电阻点焊机,电流测量仪器罗氏线圈,工控机,控制及数据采集电路板等,其中控制及数据采集电路板为自行研制而成,微处理器采用DSP进行设计,以常规的连续逐次积分方法为基准进行计算。电压为380V的标准电压,焊接有效控制时间均为20个控制周期(每个控制周期为10毫秒)。
选择晶闸管触发角度分别为104.4°,109.8°,113.4°,117°和120.6°这5种情况,分别采用本发明一种焊接电流实时测量计算方法与常规的连续积分方法进行计算,并对计算结果进行对比。每个触发角度时两种方法电流有效值对比曲线图和电流计算误差曲线图见说明书图3至图12。可以看出两种方法计算误差很小,基本在-200A至400A之间振动,平均误差不超过1.6%,最大误差不超过3.5%。
表1为不同触发角情况下两种方法误差与CPU处理时间的对比。
表1.不同触发角情况下两种方法误差与CPU处理时间对比
结果可以看出,误差很小,但是本发明CPU工作完成时间远远低于传统方法CPU工作时间,且传统处理方法在前三组触发角情况下,根本无法完成在线计算。
本实施例所述一种焊接电流实时测量计算方法应用过程如下:晶闸管触发角是单相交流电阻点焊机的唯一一个输入,每一个控制周期需要一个晶闸管触发角来决定有多少能量进入焊接系统中,进而熔化固态金属。当能量进入焊接系统后,产生电流,通过罗氏线圈,根据合适的采集频率(本实施例中使用的采集频率为6.4kHz)采集电流数值,即为id,当采集完成后,可以相应的得到电流峰值γ,计算功率因数角以及电流作用域范围导通角θ。根据这几个数值可以计算依赖模型部分有效值IRMS1,构建计算非模型依赖部分的神经网络系统,之后计算不依赖模型部分有效值IRMS2。由于神经网络模型是基于不依赖模型部分,因此可以适用于任何相似电路中,也即构建一次之后,以后可以无限次的使用。分别计算出电流有效值的两个部分IRMS1和IRMS2后,最终计算总的电流有效值IRMS。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (4)
1.一种焊接电流实时测量计算方法,具体步骤如下:
(1)以晶闸管触发角α为单相交流电阻点焊机唯一输入参数,控制输入焊接系统能量;
(2)通过模拟-数字采集装置采集焊接电流id;
(3)通过采集到的焊接电流id,搜索电流序列最大值所在的采集时间点γ;
(4)计算一个周期内电流的导通角θ;
(5)根据触发角α,导通角θ与功率因数角之间的数学关系求解功率因数角
(6)将焊接电流有效值计算公式中与电流采集量有关的部分设置为电流有效值依赖模型部分IRMS1,该部分数值可通过电流序列最大值及其采集时间点γ,结合触发角α与功率因数角计算而得到,依赖模型部分有效值IRMS1按以下公式计算:
式中,U为电压和Z为系统等效阻抗;
(7)将焊接电流有效值计算公式中仅与触发角α,导通角θ与功率因数角有关的部分设置为电流有效值不依赖模型部分IRMS2,该部分数值可通过这三个角度值,利用神经网络计算而得到;
(8)计算最终的焊接电流有效值IRMS,该值由电流有效值依赖模型部分IRMS1与电流有效值不依赖模型部分IRMS2的乘积计算而得到,公式如下:
IRMS=IRMS1×IRMS2。
2.根据权利要求1所述的一种焊接电流实时测量计算方法,其特征在于:步骤(4)中的导通角θ按以下公式计算:
式中N为一个控制周期内电流采集值的总个数,fs为电流数据采样频率,Tc为控制周期,在工频为50赫兹的情况下,一般为10毫秒。
3.根据权利要求1所述的一种焊接电流实时测量计算方法,其特征在于:步骤(5)中的功率因数角按以下公式计算:
式中Ts为电流数据采样周期,f为工频频率,即50Hz,M为到达每一个控制周期内电流峰值时的电流采集值的个数,γ为每一个控制周期内的电流峰值,ω为角频率。
4.根据权利要求1所述的一种焊接电流实时测量计算方法,其特征在于:步骤(7)中的不依赖模型部分有效值IRMS2计算采用神经网络方法计算,神经网络架构使用2×5×1形式可以取得较好效果,以P=[α,θ]T为输入值,离线训练获得神经网络输入层与隐层之间权值WI,隐层与输出层之间权值WII;输入层与隐层之间偏离值BI和隐层与输出层之间偏离值BII,控制系统实时计算过程如下:
T1=WI·P+BI
IRMS2=WⅡ×T2+BⅡ
式中T1,T2均为5行1列的矩阵,T1,i,T2,i表示矩阵中每一个具体的元素。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611163840.5A CN106513969B (zh) | 2016-12-09 | 2016-12-09 | 一种焊接电流实时测量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611163840.5A CN106513969B (zh) | 2016-12-09 | 2016-12-09 | 一种焊接电流实时测量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106513969A CN106513969A (zh) | 2017-03-22 |
CN106513969B true CN106513969B (zh) | 2019-05-14 |
Family
ID=58339894
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611163840.5A Active CN106513969B (zh) | 2016-12-09 | 2016-12-09 | 一种焊接电流实时测量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106513969B (zh) |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2742545B2 (ja) * | 1994-02-25 | 1998-04-22 | ミヤチテクノス株式会社 | 抵抗溶接制御方法 |
RU2253552C1 (ru) * | 2003-12-01 | 2005-06-10 | Открытое акционерное общество "АВТОВАЗ" (ОАО "АВТОВАЗ") | Способ управления сварочным током при контактной точечной сварке |
CN100575967C (zh) * | 2007-09-20 | 2009-12-30 | 上海交通大学 | 交流点焊动态电阻实时测量装置与方法 |
CN101733540B (zh) * | 2009-11-23 | 2011-11-16 | 上海交通大学 | 交流电阻点焊中实时获得触发角和熄灭角数值的方法 |
CN102451951A (zh) * | 2010-10-27 | 2012-05-16 | 西安扩力机电科技有限公司 | 一种电阻焊焊接设备用电流检测装置 |
CN105522271A (zh) * | 2014-10-27 | 2016-04-27 | 西安越度机电科技有限公司 | 一种电阻焊焊接的电流检测方法 |
-
2016
- 2016-12-09 CN CN201611163840.5A patent/CN106513969B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106513969A (zh) | 2017-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nascimento et al. | Harmonic distortion monitoring for nonlinear loads using neural-network-method | |
CN109217371B (zh) | 考虑锁相环影响的电压源型变流器并网系统稳定性分析方法、装置及系统 | |
CN102253266B (zh) | 一种基于线性滤波器的三相四线制谐波电流检测方法 | |
CN104578134B (zh) | 一种最大功率点跟踪方法和系统 | |
CN102565523A (zh) | 一种电流谐波检测系统及工作方法 | |
CN102005762B (zh) | 有源滤波器中直流电压在线辨识方法 | |
CN104833852A (zh) | 一种基于人工神经网络的电力系统谐波信号估计测量方法 | |
CN107863774A (zh) | 谐波指令电流的获取方法、装置、系统及可读存储介质 | |
CN106513969B (zh) | 一种焊接电流实时测量计算方法 | |
CN202564962U (zh) | 三相三线不平衡负荷补偿装置 | |
CN104143829B (zh) | Pwm整流器的直接电流控制方法及装置 | |
CN101963787B (zh) | 电力有源滤波器的双反馈环控制方法 | |
Gupta et al. | Neural network-based shunt active filter with direct current control for power quality conditioning | |
CN101079582B (zh) | 一种逆变系统的控制方法 | |
CN109142857A (zh) | 列车用带相位自动校准功能的能量计算系统和方法 | |
Gong et al. | Model-based real-time dynamic power factor measurement in AC resistance spot welding with an embedded ANN | |
CN103094909A (zh) | 基于静止无功发生器的高精度单相数字锁相环的控制方法 | |
CN109307800A (zh) | 一种电网总谐波检测方法 | |
CN105811419B (zh) | 一种精确抑制电网谐波的控制方法 | |
CN103488891A (zh) | 基于风险的大电网安全稳定预防控制性能代价指标计算方法 | |
CN107611987A (zh) | 基于gmn算法的静止同步补偿器的控制方法 | |
Sebasthirani et al. | Design of shunt active power filter with fuzzy logic control for mitigating harmonics | |
Tan et al. | A simulation research on single phase bridge full control resistive load rectifying circuit based on MATLAB | |
CN108268856B (zh) | 基于l2范数和真正跟踪误差的变步长自适应谐波检测方法 | |
CN109407543B (zh) | 一种风电机组电气模型电压响应特性的验证方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20190326 Address after: 510635 No. 576 Longkou West Road, Tianhe District, Guangzhou City, Guangdong Province Applicant after: Guangdong Normal University of Technology Applicant after: BEIJING INSTITUTE OF TECHNOLOGY Address before: 510665 West 293 Zhongshan Avenue, Tianhe District, Guangzhou City, Guangdong Province Applicant before: Guangdong Technical Normal College |
|
GR01 | Patent grant | ||
GR01 | Patent grant |