CN1815212B - 金属冲压过程中的诊断方法及其设备 - Google Patents

Abstract

本发明公开了一种对工件在冲压过程中进行诊断的方法，包括确定所述工件的热能和温度分布；提供冲压后所述工件的热能分布图像；以及将所述的热能和温度分布与所述热能分布图像进行比较以判断工件在冲压过程中是否存在缺陷，其中，热能和温度分布与所述热能分布图像的差异表明所述工件在冲压过程中存在缺陷。本发明还提供对工件在冲压过程中进行诊断的设备。利用本发明的方法和设备可对工件冲压过程中的出现的问题进行准确的分析和诊断。

Description

金属冲压过程中的诊断方法及其设备

发明领域

本发明涉及金属加工领域，特别涉及工件冲压过程中的诊断方法和实施该方法的设备。

背景技术

工件的冲压过程是一个复杂的过程，它涉及到工件板料的瞬态弹性变形和塑性变形，以及冲床的静态和动态行为。影响工件冲压过程的因素很多，包括冲床、模具、工件以及冲压条件(如：冲压速度，润滑状态等)。目前，由于对降低生产成本日益增长的需求，金属冲压过程通常处于最大负荷下运行(包括几何形状，冲压速度，冲床功率的充分利用等)。这往往导致各种故障或异常的发生，例如：产品尺寸超限，撕裂，皱褶等，当这些故障发生时，目前上无科学有效的方法找出故障发生的原因。

虽然金属冲压中的应力应变分布可以借助与有限元商业软件包如

和计算得到，但是这些商业软件包仅能够预测金属冲压过程中工件的应力应变的变化，而且这种预测是基于理想冲压条件下计算得到的，不可能包括诸如模具制造误差(包括加工误差，装配误差及表面粗糙度等)，机床的动态特性(如振动)，以及冲压条件(润滑，摩擦等)。因此，这种预测有高达25％的误差。而且，当生产中出现故障时，这种方法无法找出故障发生的根本原因。

本发明正是为解决现有技术中的该问题而提出。

发明内容

本发明第一方面提供一种工件冲压过程中进行诊断的方法，所述方法包括：

提供冲压后所述工件的热能分布图像；以及

根据所述的热能图像判断工件在冲压过程中是否存在缺陷，

其中热能在所述热能图像中的异常分布表明所述工件在冲压过程中存在缺陷。

本发明第二方面提供一种对工件在冲压过程中进行诊断的方法，所述方法包括：

确定所述工件的热能和温度分布；

提供冲压后所述工件的热能分布图像；以及

将所述的热能和温度分布与所述热能分布图像进行比较以判断工件在冲压过程中是否存在缺陷，

其中，热能和温度分布与所述热能分布图像的差异表明所述工件在冲压过程中存在缺陷。

在本发明方法的一实施方案中，所述的确定热能和温度分布包括对所述冲压过程进行建模。

在本发明方法的一另实施方案中，所述的对所述冲压过程进行建模通过有限元分析进行。

在本发明的又一实施方案中，所述的热能分布图像通过红外照相获得，优选地，所述热能分布图像由通过对工件多次红外照相获得，所述多次红外照相为沿所述工件的不同方向进行。

本发明第三方面提供一种对工件在冲压过程中进行诊断的设备，包括提供冲压后确定所述工件的热能分布图像的装置。

在本发明设备的一实施方案中，进一步包括：

提供冲压前所述工件的热能和温度分布的装置；以及

数据处理装置，与所述确定所述工件的热能和温度分布的装置和所述提供冲压后所述工件的热能分布图像的装置电连接，对所述工件的热能和温度分布与所述的热能分布图像进行比较。

在本发明设备的另一实施方案中，所述提供冲压后确定所述工件的热能分布图像的装置包括：

支撑和固定所述工件的支撑装置；

红外照相机，将其摄取的工件红外图像输送到所述数据处理装置。

在本发明设备的另一实施方案中，所述提供冲压后确定所述工件的热能分布图像的装置进一步包括绝热壳体，所述绝热壳体具有：至少一壁以固定红外相机。

在本发明设备的又一实施方案中，所述提供冲压后确定所述工件的热能分布图像的装置进一步包括位于驱动所述支撑装置的驱动装置。

附图说明

图1是本发明有限元分析构建中某单元的四点四边形示图；

图2是本发明图1在中单元在虚拟空间中的四点四边形示图；

图3为由一步冲压成形的杯子的示图；

图4为图3的杯子根据本发明有限元分析得到的热能分布结果；

图5是对一工件进行连续性冲压得到的视图；

图6是图7工件在连续性冲压过程中的热能图像；

图7是本发明诊断方法的方块图；

图8是本发明诊断系统的装置示意图；

图9为图3冲压杯子的六幅热能图像；

图10为从不同角度重构后的冲压杯子热能分布。

具体实施方式

本发明的基础是能量守恒定律，即在金属冲压过程中，金属吸收能量而发生变形，变形能量最后又转换成热量。因此，工件(或模具，因为热量由工件传递给模具)的温度分布是与工件的变形(应变)和应力(冲压力)有关的。换句话说，工件的温度分布描述了工件在何处以及怎样发生变形的。因此，通过分析成形工件(或模具)的温度分布，我们可以发现模具的发生过度的变形(应力/应变)的位置，从而诊断出故障的根源。

有限元分析已被应用于金属冲压中，商业上也可获得对于金属冲压的有限元分析软件。然而，现有技术中没有对工件冲压中热能分布的有限元分析建模。

本发明人通过对金属冲压热能分布的有限元分析建模，进而运用该结果对金属冲压过程中进行诊断。藉此提供本发明的方法和设备。

图7的方块图表示本发明的方法。简单地说，我们首先模拟工件的热能有限元分析。然后，我们再照数张刚冲压完成的红外图像。对结构简单的工件，可以从得到的图像上判断工件在冲压过程中是否正常。

如果得到的图像很复杂，则简单的二维红外图像不能完全反应工件的每一部分的温度。在此情况下，便利用得到的复杂的二维热红外图像重构其三维的几何与热能分布。这样，我们便能比较三维的热能有限元结果和三维重构的热能分布。则冲压过程中的问题便能分析和诊断。

冲压过程中涉及大量能量的消耗，因而产生明显的温度改变。温度的改变会影响工件的物质性质，并改变最终的能量消耗。由于要补偿这些误差，以下便描述了本发明的建模。

参照图1，首先，我们将工件分解成大量细小的单元。每一单元都是四边形的，并存在四个点1，2，3，4，而每个点均联系了两个坐标，r与z。

接着，我们模拟冲压过程中的其中一个单元为一个四点四边形的单元。Γ_e，1，Γ_e，2，Γ_e，3与Γ_e，4是单元的边界。

物理坐标r与z分别在虚空间对应了两个新的坐标ξ与η。它们描述于图2。

图2中的四个点1’，2’，3’，4’，其对应空间的形函数是：

$N_1=\frac{(1-\mathrm{\ξ})(1-\mathrm{\η})}{4}$

$N_2=\frac{(1+\mathrm{\ξ})(1-\mathrm{\η})}{4}$

$N_3=\frac{(1+\mathrm{\ξ})(1+\mathrm{\η})}{4}$

$N_4=\frac{(1-\mathrm{\ξ})(1+\mathrm{\η})}{4}$

要把(ξ，η)系统转化成(r，z)系统，我们使用以下的公式：

$\frac{{\∂N}_i}{\∂\mathrm{\ξ}}=\frac{{\∂N}_i}{\∂r}\frac{\mathrm{dr}}{\mathrm{d\ξ}}+\frac{{\∂N}_i}{\∂z}\frac{\mathrm{dz}}{\mathrm{d\ξ}}$

$\frac{{\∂N}_i}{\∂\mathrm{\η}}=\frac{{\∂N}_i}{\∂r}\frac{\mathrm{dr}}{\mathrm{d\η}}+\frac{{\∂N}_i}{\∂z}\frac{\mathrm{dz}}{\mathrm{d\η}}$

即

$\left\{\begin{array}{c}\frac{{\∂N}_i}{\∂\mathrm{\ξ}}\\ \frac{{\∂N}_i}{\∂\mathrm{\η}}\end{array}\right\}=\left[\begin{array}{cc}\frac{\mathrm{dr}}{\mathrm{d\ξ}}& \frac{\mathrm{dz}}{\mathrm{d\ξ}}\\ \frac{\mathrm{dr}}{\mathrm{d\η}}& \frac{\mathrm{dz}}{\mathrm{d\η}}\end{array}\right]\left\{\begin{array}{c}\frac{{\∂N}_i}{\∂r}\\ \frac{{\∂N}_i}{\∂z}\end{array}\right\}$

其(2×2)矩阵

$\left[J_e\right]=\left[\begin{array}{cc}\frac{\mathrm{dr}}{\mathrm{d\ξ}}& \frac{\mathrm{dz}}{\mathrm{d\ξ}}\\ \frac{\mathrm{dr}}{\mathrm{d\η}}& \frac{\mathrm{dz}}{\mathrm{d\η}}\end{array}\right]$

其包含不同的经分解的工件的单元。我们可以从物理空间中的点1、2、3和4表示出单元上的任何一点。

r(ξ，η)＝N₁r₁+N₂r₂+N₃r₃+N₄r₄

z(ξ，η)＝N₁z₁+N₂z₂+N₃z₃+N₄z₄

其r_i与z_i是点i的坐标。

所以[J_e]变为

$\left[J_e\right]=\left[\begin{array}{cccc}\frac{{\mathrm{dN}}_1}{\mathrm{d\ξ}}& \frac{{\mathrm{dN}}_2}{\mathrm{d\ξ}}& \frac{{\mathrm{dN}}_3}{\mathrm{d\ξ}}& \frac{{\mathrm{dN}}_4}{\mathrm{d\ξ}}\\ \frac{{\mathrm{dN}}_1}{\mathrm{d\η}}& \frac{{\mathrm{dN}}_2}{\mathrm{d\η}}& \frac{{\mathrm{dN}}_3}{\mathrm{d\η}}& \frac{{\mathrm{dN}}_4}{\mathrm{d\η}}\end{array}\right]\left[\begin{array}{cc}{r}_1& z_1\\ r_2& z_2\\ r_3& z_3\\ r_4& z_4\end{array}\right]=\left[\begin{array}{cccc}\mathrm{\η}-1& 1-\mathrm{\η}& \mathrm{\η}& -\mathrm{\η}\\ \mathrm{\ξ}-1& -\mathrm{\ξ}& \mathrm{\ξ}& 1-\mathrm{\ξ}\end{array}\right]\left[\begin{array}{cc}{r}_1& z_1\\ r_2& z_2\\ r_3& z_3\\ r_4& z_4\end{array}\right]$

其中

$\left[J_e\right]=\left[\begin{array}{cc}{J}_{11}& J_{12}\\ J_{21}& J_{22}\end{array}\right]$ 和 ${\left[J_e\right]}^{-1}=\left[\begin{array}{cc}{J}_{11}^{-1}& J_{12}^{-1}\\ J_{21}^{-1}& J_{22}^{-1}\end{array}\right]$

[J_e]容许我们建立两者空间的关系并简化其积分，参考下面公式

$\underset{S_{\mathrm{physical}}}{\∫\∫}\mathrm{dS}=\underset{S_{\mathrm{virtual}}}{\∫\∫}\det \left[J_e\right]\mathrm{dS}=\underset{-1}{\overset{1}{\∫}}\underset{-1}{\overset{1}{\∫}}\det \left[J_e\right]\mathrm{d\ξd\η}$ 与 $\det \left[J_e\right]=\left|\begin{array}{cc}\underset{j=1}{\overset{4}{\mathrm{\Σ}}}{r}_j\frac{{\mathrm{dN}}_j}{\mathrm{d\ξ}}& \underset{j=1}{\overset{4}{\mathrm{\Σ}}}{z}_j\frac{{\mathrm{dN}}_j}{\mathrm{d\ξ}}\\ \underset{j=1}{\overset{4}{\mathrm{\Σ}}}{r}_j\frac{{\mathrm{dN}}_j}{\mathrm{d\η}}& \underset{j=1}{\overset{4}{\mathrm{\Σ}}}{z}_j\frac{{\mathrm{dN}}_j}{\mathrm{d\η}}\end{array}\right|.$

根据于热能学第一定律，能量守恒的比率(Rate of Conservation ofEnergy)描述于下面公式

(在V中热能增加的比率)＝

(经S传到V中热能传导的比率)+(在V中热能产生的比率)(1)

其中V是体积，S是包含V的面。

设u代表工件介质特定内部能量，则

其中ρ是工件的密度(kg/m³)。

根据

我们得出

利用傅立叶的热传导理论，我们得出经S传到V中热能传导比率的表达式

$q=-k\left(\mathrm{grad}\left(T\right)\right)n=-k\frac{\∂T}{\∂n}$

其中k是工件的热能传导性(W/m.K)，并在V中假定为常数，

n是面S的垂直方向；

q是n方向中的热流量；

所以我们得出：

假设热能产生在以每体积单位Q的比率发生，这样

其中

是每体积单位热转移的比率，而是因每体积单位的功的热转移，所以其平衡方程可写成

$\underset{V}{\∫}(\mathrm{\ρc}\frac{\∂T}{\∂t}-k\mathrm{div}\left(\mathrm{grad}\left(T\right)\right)-\stackrel{\·}{w})\mathrm{dV}=0$

因体积V在开始时是随意挑选，因此它按照以下公式

$\mathrm{\ρc}\frac{\∂T}{\∂t}=k\mathrm{div}\left(\mathrm{grad}\left(T\right)\right)+\stackrel{\·}{w}={\mathrm{k\Δ}}^{2}T+\stackrel{\·}{w}$

边界条件必须在之前(1)的守恒公式中使用。

在冲压过程中，工件下部的面(Γ₁与其垂直向量n₁)是与模具的接触面。当时间t＝0，T_part＝T_die＝T_air，我们假设工具与摩擦的表面能量消耗是

$-k\frac{\∂T}{{\∂n}_1}=h_{\mathrm{part}/\mathrm{die}}(T-T_{\mathrm{die}})-\frac{b}{b+b_{\mathrm{die}}}{\mathrm{\τv}}_s$

其中h是当模具在温度T_die的热交换参数，b与b_die是渗透率

τ摩擦的应力，而v_s是其模具与工件的相对速度。

那向上的面是与空气接触(Γ₃与其垂直向量n₃)。当t＝0，T_part＝T_air，我们能够设定其辐射与对流近似为

$-k\frac{\∂T}{{\∂n}_3}=h_{\mathrm{conv}}(T-T_{\mathrm{air}})+{\mathrm{\ϵ}}_r{\mathrm{\σ}}_r(T^4-T_0^4)$

其中ε_r是发射率参数，σ_r是Stefan常数(Stefan constant)，T_air是外面温度，而h_conv是空气对流的热交换系数。

工件与周遭空气的热交换能够被传导(当工件与空气的温度差不大于克服空气流动的阻力)建模。因此，引发了自然对流电流。加上，我们忽略辐射的因素而简化了问题，所以，

$-k\frac{\∂T}{{\∂n}_3}=h_{\mathrm{part}/\mathrm{air}}(T-T_{\mathrm{air}})$

假设流经单元的热量为q。这热量流动的速率是不知道的，但我们可以用形函数来估计出来，

$\stackrel{\·}{q}=\underset{j}{\mathrm{\Σ}}{N}_j{\stackrel{\·}{q}}_j=\left\{N_j\right\}\left\{\stackrel{\·}{q}\right\}$

所以，我们可用以下公式来描述这传导

$-k\frac{\∂T}{\∂n}=\stackrel{\·}{q}$

其中n＝{n₂；n₄}

当t＝0，我们假设热流动与比率在工件的任何地方均是零。运用(Galerkin’s method)，我们得出一个另一形态的公式

${\∫N}_i(k\frac{{\∂}^{2}T}{{\∂r}^2}+k\frac{{\∂}^{2}T}{{\∂z}^2}+\frac{k}{r}\frac{\∂T}{\∂r}+\stackrel{\·}{w}-\mathrm{\ρc}\frac{\∂T}{\∂t})\mathrm{rdrd\θdz}=0$

在首两项分部积分

$\text{\∫}{N}_{i}k\frac{{\∂}^{2}T}{{\∂r}^2}\mathrm{rdrd\θdz}=\∫{\mathrm{krN}}_i\frac{\∂T}{\∂r}\mathrm{d\θdz}-\∫k\frac{\∂T}{\∂r}\frac{\∂\left(N_{i}r\right)}{\∂r}\mathrm{drd\θdz}$

$\text{\∫}{N}_{i}k\frac{{\∂}^{2}T}{{\∂r}^2}\mathrm{rdrd\θdz}=\∫{\mathrm{krN}}_i\frac{\∂T}{\∂r}\mathrm{d\θdz}-\∫k\frac{\∂T}{\∂r}(r\frac{{\∂N}_i}{\∂r}+N_i)\mathrm{drd\θdz}$

$\text{\∫}{N}_{i}k\frac{{\∂}^{2}T}{{\∂z}^2}\mathrm{rdrd\θdz}=\∫{\mathrm{krN}}_i\frac{\∂T}{\∂z}\mathrm{d\θdz}-\∫k\frac{\∂T}{\∂r}r\frac{{\∂N}_i}{\∂z}\mathrm{drd\θdz}$

我们得出

$\text{\∫}{\mathrm{krN}}_i\frac{\∂T}{\∂r}\mathrm{d\θdz}+\∫{\mathrm{krN}}_i\frac{\∂T}{\∂z}\mathrm{d\θdr}$

$-\∫(k\frac{{\∂N}_i}{\∂r}\frac{\∂T}{\∂r}+k\frac{{\∂N}_i}{\∂z}\frac{\∂T}{\∂z}+N_i\stackrel{\·}{w}-N_i\mathrm{\ρc}\frac{\∂T}{\∂t})\mathrm{rdrd\θdz}=0$

独立的θ允许我们对公式的简化。

所以，那几何建模的结果是

$r=\underset{j}{\mathrm{\Σ}}{N}_j{r}_j=\left\{N_j\right\}\left\{r\right\}$

$z=\underset{j}{\mathrm{\Σ}}{N}_j{z}_j=\left\{N_j\right\}\left\{z\right\}$

$T=\underset{j}{\mathrm{\Σ}}{N}_j{T}_j=\left\{N_j\right\}\left\{T\right\}$

drdz＝det[J]dξdη

$\frac{{\∂N}_i}{\∂r}=J_{11}^{-1}\frac{{\∂N}_i}{\∂\mathrm{\ξ}}+J_{12}^{-1}\frac{{\∂N}_i}{\∂\mathrm{\η}}$

$\frac{{\∂N}_i}{\∂z}=J_{21}^{-1}\frac{{\∂N}_i}{\∂\mathrm{\ξ}}+J_{22}^{-1}\frac{{\∂N}_i}{\∂\mathrm{\η}}$

$\frac{\∂T}{\∂r}=\left\{\frac{{\∂N}_j}{\∂r}\right\}\left\{T\right\}$

$\frac{\∂T}{\∂z}=\left\{\frac{{\∂N}_j}{\∂z}\right\}\left\{T\right\}$

$\left\{T_e\right\}=\left\{\begin{array}{c}{T}_1\\ T_2\\ T_3\\ T_4\end{array}\right\}$

$\left\{T\right\}=\left\{\begin{array}{c}{T}_1\\ T_2\\ .\\ .\\ T_n\end{array}\right\}$

$\text{\∫}{\mathrm{krN}}_i\frac{\∂T}{\∂r}\mathrm{dz}+\∫{\mathrm{krN}}_i\frac{\∂T}{\∂z}\mathrm{dr}-\∫(k\frac{{\∂N}_i}{\∂r}\frac{\∂T}{\∂r}+k\frac{{\∂N}_i}{\∂z}\frac{\∂T}{\∂z}+N_i\stackrel{\·}{w}-N_i\mathrm{\ρc}\frac{\∂T}{\∂t})r\det \left[J\right]\mathrm{d\ξd\η}=0$

考虑边界条件，我们得出

$-\underset{{\mathrm{\Γ}}_1}{\∫}{\mathrm{rN}}_i{h}_{\mathrm{part}/\mathrm{die}}(T-T_{\mathrm{die}}){\mathrm{d\Γ}}_1-\underset{{\mathrm{\Γ}}_3}{\∫}{\mathrm{rN}}_i{h}_{\mathrm{part}/\mathrm{air}}(T-T_{\mathrm{air}}){\mathrm{d\Γ}}_3-\underset{{\mathrm{\Γ}}_2}{\∫}{\mathrm{rN}}_i{\stackrel{\·}{q}}_2-\underset{{\mathrm{\Γ}}_4}{\∫}{\mathrm{rN}}_i{\stackrel{\·}{q}}_{4}d{\mathrm{\Γ}}_4$

$\text{-\∫}(k\frac{{\∂N}_i}{\∂r}\frac{\∂T}{\∂r}+k\frac{{\∂N}_i}{\∂z}\frac{\∂T}{\∂z}+N_i\stackrel{\·}{w}-N_i\mathrm{\ρc}\frac{\∂T}{\∂t})r\det \left[J\right]\mathrm{d\ξd\η}=0$

由于Γ_i是代表线性单元(线性单元内插值是线性的)的两个节点，我们可以计算出dΓ_i是

${\mathrm{d\Γ}}_i=\sqrt{{\left(\frac{\mathrm{dr}}{\mathrm{d\ξ}}\right)}^2+{\left(\frac{\mathrm{dz}}{\mathrm{d\ξ}}\right)}^2}\mathrm{d\ξ},$ 对i＝1，η＝-1，而对i＝3，η＝1；

${\mathrm{d\Γ}}_i=\sqrt{{\left(\frac{\mathrm{dr}}{\mathrm{d\η}}\right)}^2+{\left(\frac{\mathrm{dz}}{\mathrm{d\η}}\right)}^2}\mathrm{d\η},$ 对i＝2，ξ＝1，而对i＝4，ξ＝-1。

即是

${\mathrm{d\Γ}}_1=\frac{1}{2}\sqrt{{(r_1-r_2)}^2+{(z_1-z_2)}^2}\mathrm{d\ξ}$

${\mathrm{d\Γ}}_2=\frac{1}{2}\sqrt{{(r_2-r_3)}^2+{(z_2-z_3)}^2}\mathrm{d\η}$

${\mathrm{d\Γ}}_3=\frac{1}{2}\sqrt{{(r_3-r_4)}^2+{(z_3-z_4)}^2}\mathrm{d\ξ}$

${\mathrm{d\Γ}}_4=\frac{1}{2}\sqrt{{(r_4-r_1)}^2+{(z_4-z_1)}^2}\mathrm{d\η}$

这些公式可以写成以下形式

$M_e\frac{d\left\{T_e\right\}}{\mathrm{dt}}+K_e\left\{T_e\right\}=f_e+{\widehat{\stackrel{\·}{q}}}_e$

其中

M_e，ij＝∫ρcN_iN_j{N_k}{r}det[J]dξdη

$K_{e,\mathrm{ij}}=-\∫h_{\mathrm{part}/\mathrm{die}}\left\{N_k\right\}\left\{r\right\}{N}_i{N}_j{\mathrm{d\Γ}}_1-\∫h_{\mathrm{part}/\mathrm{air}}\left\{N_k\right\}\left\{r\right\}{N}_i{N}_j{\mathrm{d\Γ}}_3$

$-\∫k\left\{N_k\right\}\left\{r\right\}(\frac{{\∂N}_i}{\∂r}\frac{{\∂N}_j}{\∂r}+\frac{{\∂N}_i}{\∂r}\frac{{\∂N}_j}{\∂z})\det \left[J\right]\mathrm{d\ξd\η}$

${\widehat{\stackrel{\·}{q}}}_{e,i}=\underset{{\mathrm{\Γ}}_2}{\text{\∫}}r{N}_i{\stackrel{\·}{q}}_2{\mathrm{d\Γ}}_2+\underset{{\mathrm{\Γ}}_4}{\∫}{\mathrm{rN}}_i{\stackrel{\·}{q}}_4{\mathrm{d\Γ}}_4$

用数值法来描述热流动率，我们得出

${\widehat{\stackrel{\·}{q}}}_{e,i}=\underset{k=1}{\overset{{\mathrm{Nb}}_{\mathrm{nodes}}}{\mathrm{\Σ}}}(\left(\underset{{\mathrm{\Γ}}_2}{\∫}{\mathrm{rN}}_i{N}_k{\mathrm{d\Γ}}_2\right){\stackrel{\·}{q}}_{2,k}+\left(\underset{{\mathrm{\Γ}}_4}{\∫}{\mathrm{rN}}_i{N}_k{\mathrm{d\Γ}}_4\right){\stackrel{\·}{q}}_{4,k})$

其中是从点k中边界Γj未知的热流动率，已知在Γ₂，ξ＝1与Γ₄，ξ＝-1。由于流动率是从边向外的垂直方向，得出，

$f_{e,i}=\∫N_i\stackrel{\·}{w}\left\{N_k\right\}\left\{r\right\}\det \left[J\right]\mathrm{d\ξd\η}-\∫h_{\mathrm{part}/\mathrm{die}}\left\{N_k\right\}\left\{r\right\}{N}_i{T}_{\mathrm{die}}{\mathrm{d\Γ}}_1$

$-\∫h_{\mathrm{part}/\mathrm{air}}\left\{N_k\right\}\left\{r\right\}{N}_i{T}_{\mathrm{air}}{\mathrm{d\Γ}}_3$

如果是以Ω空间的拉格朗日插值(Lagrange interpolation)来表达，那近似值

便能用来简化向量{f}的评估。

之后，我们得出

与

的关系

${\stackrel{\·}{w}}_j=\mathrm{\Σ}{\mathrm{\σ}}_{i,j}{\stackrel{\·}{\mathrm{\ϵ}}}_{i,j}.$

加上，

$f_{e,i}=\underset{j=1}{\overset{4}{\mathrm{\Σ}}}\left({\∫N}_j{N}_i\right\{N_k\left\}\right\{r\left\}\det \right[J\left]\mathrm{d\ξd\η}\right){\stackrel{\·}{w}}_j$

$-\∫h_{\mathrm{part}/\mathrm{die}}\left\{N_k\right\}\left\{r\right\}{N}_i{T}_{\mathrm{die}}d{\mathrm{\Γ}}_1-{\∫h}_{\mathrm{part}/\mathrm{air}}\left\{N_k\right\}\left\{r\right\}{N}_i{T}_{\mathrm{air}}d{\mathrm{\Γ}}_3$

以及σ_j与

应力与应变都能从点j得出。那主要成分(

其中i＝{1，2，3}的总和便是应力与应变的主要方向。

根据本发明的方法，可以诊断简单和复杂工件或模具在冲压过程中的情况。

如果工件是对称，但得出的热能分布不是，这样便是有问题。如工件的形状很复杂，便需要进一步的分析。它涉及最初的热能分布的有限元分析，然后比较和分析获得的热能图像。这样，我们便能准确地诊断出问题的原因。如果工件再复杂，我们便用刚冲压后获得的数张红外热能图像来三维重构，以展示它的几何和热能分布

图3展示了一由一步冲压成形的杯子。根据本发明的有限元分析，其热能分布结果示于图4中。从图4可知其温度分布情况。通过比较图3和4，便可知道杯子的理想热能分布和冲压过程中出现的问题如何。

图6展示一个由连续性冲压过程制成的传统小部件。在冲压制造中，并没有发生严重的问题，除了偶尔部件超出容许的尺寸范围。

虽然视觉检测能描述容许的尺寸范围，但却不能反映出什么导致超出其要求。我们的在线监控系统能检测出问题，但却不能指出问题的原因。然而，根据本发明的方法，通过分析连续性冲压过程该小部件的热能分布的有限元分析结果，我们可以看出金属片中不同地方的不同温度分布。

本发明中，热能图像的获得可以通过现有技术中任何已知的方法进行。例如，红外照相等手段。

但是热能图像摄取与分析是一项复杂的工作。首先，电磁波幅射覆盖了非常阔的光谱并互相干扰。换句话说，热能图像会被不同的噪音影响，例如，环境光线、反射、操作员的身体热能和机器。

因此，本发明提供的对工件在冲压过程中进行诊断的设备，进一步包括一绝热壳体10(请参考图8)。壳体10通常由绝热材料制成，将处理工件封闭在其中，并且将其构型为尽量避免热损失的形状，使得摄取到的工件热能图像最大程度反应其实际情况。在一实施方案中，将其制成内壁光滑的球形结构。显然，本发明的方法优选对刚冲压完成的工件进行，以减少对流的热能流失。

如上所述，本发明的方法，对于简单的工件，如描述于图5中的工件，根据得到热能图像的结果已足够确认出潜在的问题。图6示出了其红外线照相机(制造商：Guide，型号：IR913)拍摄到的工件的热能分布。在图6中，600处有一高温点(白色)，它显示了工件的非典型变形。这很可能是导致问题的原因：模具在那地方制造得不准确。那定位销的弯曲导致额外的摩擦力，因而产生不同的温度。结果，工件便被推成超出了尺寸。

图8展示了本发明的诊断设备的一实施方案。绝热壳体10具有至少一壁11，红外相机400固定在壁11上。在壳体内具有下部部分12，和底部部分13。支撑工件100的转盘200装配在壳体的下部部分12。驱动转盘200转动的电机300固定在壳体的底部部分13。

红外相机和一计算机500电连接。根据本发明方法获得的工件100热能和温度分析的数据贮存在计算机500内。

将工件100置于转盘200上，转动转盘通过红外相机的连续照相可获得工件不同角度的热能图像，其图像数据传到安装了模拟和分析软件的计算机500进行处理。

根据本发明的方法，可实现一些冲压工件的三维热能重构。为提高准确度，先我们可将诊断系统校验，得出转盘和红外线照相机的外在关系以及红外线照相机的内在参数。接着，我们便从所摄得的图像中抽取工件的轮廓，并用空间雕刻的方法重构其三维热能分布。

我们用图3中的冲压杯子，拍摄了十八幅热能图像来取得其几何形状和热能分布。图9展示了十八幅冲压杯子的其中六幅热能图像。图10展示了从不同角度重构后的冲压杯子热能分布。它是用来与图4的热能有限元分析模型比较。从比较中，我们便可以诊断冲压过程中潜在的问题。

本发明与现有技术相比，根据工件在冲压后的真实应变和热能分布的相似性，通过考察热能分布来诊断冲压过程。它不同于有限元分析，有限元分析只能模拟一个假设的状况，它并不能包含模具表面处理、润滑等因素。因而，本发明提高了准确性。

本发明可用在高速冲压和热能辅助冲压中，其中温度控制是一重要的因素。

虽然本发明已用上述说明和具体实施例进行了详细的描述，但本领域技术人员显然可以通过本发明的公开的教导，对本发明进行进一步的变化和修饰，这些不脱离本发明精髓的变化和修饰应属于本发明的范围。

Claims (7)

1.一种对工件在冲压过程中进行诊断的方法，所述方法包括：

通过有限元分析为所述冲压过程建立模型；

根据所建立的模型确定所述工件冲压前的热能和温度分布；

获取所述工件刚冲压后的三维热能分布图像；以及

将所述的热能和温度分布与所述三维热能分布图像进行比较以判断工件是否存在缺陷，

其中，所述热能和温度分布与所述热能分布图像的差异表明所述工件存在缺陷。

2.根据权利要求1所述的方法，其中，所述热能分布图像通过红外照相获取。

3.根据权利要求2所述的方法，其中，所述的三维热能分布图像由通过对工件数次红外照相获得，所述的数次红外照相为沿所述工件的不同方向进行。

4.一种对工件在冲压过程中进行诊断的设备，包括：

根据通过有限元分析为所述冲压过程建立的模型，确定所述工件冲压前的热能和温度分布的装置；

确定所述工件刚冲压后的热能分布图像的装置；以及

数据处理装置，与确定所述工件的热能和温度分布的装置和确定所述工件冲压后的热能分布图像的装置电连接，对所述热能和温度分布与所述热能分布图像进行比较。

5.根据权利要求4所述的设备，其中，确定所述工件冲压后的热能分布图像的装置包括：

支撑和固定所述工件的支撑装置；以及

红外照相机，将其摄取的工件红外图像输送到所述数据处理装置。

6.根据权利要求5所述的设备，其中，确定所述工件冲压后的热能分布图像的装置进一步包括绝热壳体，所述绝热壳体具有至少一个壁以固定所述红外照相机。

7.根据权利要求6所述的设备，其中，确定所述工件冲压后的热能分布图像的装置进一步包括用于驱动所述支撑装置的驱动装置。

Images