강 파이프의 고주파 전기저항 용접공정에 대한 수치 해석

2.1 전자기장 해석 모델

HF-ERW 공정에서는 대강 소재에 고주파 교류 전류가 인가되며 이로 인해 와전류가 생성되는 전자기장 현상이 발생하게 되고, 전류의 전기저항에 의해 발생하는 열을 용접 열원으로 이용하게 된다. 이때의 전자기장 현상을 묘사하는 모델은 맥스웰 방정식으로부터 유도하여 다음과 같은 방정식으로 나타낼 수 있다[11].

여기서, μ는 투자율, $$ \overrightarrow{A}$$는 자기 벡터 포텐셜, σ는 전기전도도, $$ {\overrightarrow{J}}_s$$는 전원 전류밀도, $$ \overrightarrow{J_e}$$는 와전류밀도, $$ \overrightarrow{E}$$는 전기장의 세기이다.

식(1)에서 $$ \overrightarrow{A}$$를 구하고, 이를 이용하여 식(2)에서 와전류를 구할 수 있다. 따라서, 소재에 인가된 전원 전류와 이로 인해 생성된 와전류에서 소재에 발생하는 줄열(Joule Heat), q는 식(3)과 같이 나타낼 수 있다.

방정식(1)에 대한 해는 유한요소해석 등 수치해석적 방법으로 구할 수 있고, 식(3)의 줄열을 구하면 이는 열전달 해석에서의 입력 열원인 발열량이 된다.

한편, HF-ERW에서 활용하는 300 kHz 이상의 고주파 전원과 강 소재의 경우, 침투깊이는 대략 0.0065 mm 이하 크기로 예측된다[12]. 따라서 표면을 제외한 소재 내부 영역에서는 전자기장이 존재하지 않는다고 가정할 수 있고, 줄열은 대부분이 소재 표면에서 발생한다고 판단할 수 있다. 이러한 문제에서는 표면 임피던스 경계조건(Surface Impedance Boundary Condition, SIBC)을 적용하면 시간과 용량을 크게 줄여 해석할 수 있다[12,13].

Fig. 2는 해석에 사용된 강 파이프 모델로서 대강 소재와 두 전극으로 구성된 기하학적 형상을 보이고 있다. 대강의 두 모서리 표면이 서로 만나는 V-수렴점과 용접 위치는 동일하다고 가정하고 전자기 현상은 용접부 주위에 국한하여 나타나므로, 전체 길이가 400 mm인 파이프를 대상으로 파이프의 중심에서 용접선을 기준으로 ±30^o 이내를 해석 영역으로 설정하였다.

Fig. 2

Geometric model for electromagnetic analysis

파이프 끝에서 전극 중심까지는 80 mm, 전극에서 용접 위치까지는 130 mm로, 전극 단면은 가로 22 mm, 세로 30 mm의 사각형으로, 전극 중심 간은 58 mm로, 파이프와 전극 사이 간극은 0.1 mm로 전극의 크기와 상대 위치를 각각 설정하였다. 한편, 파이프 모델로부터 일정 거리로 이격된 공간 위치까지는 전자기장의 영향이 발생한다고 가정하고 그 내부는 공기층으로 설정하였다.

용접 공정을 모사하기 위한 전기 회로는 전압이 공급되는 좌측 전극, 좌측 전극과 파이프 표면 사이의 좌측 공극, 파이프, 우측 공극, 우측 전극 등의 요소로 구성되고, 전극의 상부 표면에 주파수 320 kHz의 전압이 인가된다. 고주파수의 특성을 이용하여, 전극과 파이프는 표면 임피던스 경계를 가진 전도체로 가정하였다. 강재의 비투자율은 300, 전기 저항은 0.5 × 10^-6 Ω로, 전극의 비투자율은 1, 전기 저항은 1.72 × 10^-8 Ω로 각각 일정하다고 가정하였다. 전극이 파이프와 접촉되는 공극에서의 접촉 저항은 10^-4 Ω로 가정하였다.

2.2 온도 분포 해석 모델

전자기장 해석에서 산출된 발열량 분포를 열원으로 하여, HF-ERW 공정 중 파이프의 온도 분포를 해석하였다. 이때, 파이프의 이송 속력이 열전도보다 매우 빠르다면 파이프의 이송 방향(z축 방향)으로 열전도가 없다고 가정할 수 있다. 이에 따라, Fig. 3과 같이 시간의 변화에 따라 변화되는 해석 영역이 나타나 있듯이, x−y축 해석 단면을 파이프의 진행 속도에 의한 시간 변화(t₁, t₂, ..., t_n)에 따라 변형하면서 z축 방향 위치의 해당 발열량을 열원으로 하는 2차원 과도상태로 공정을 모델화하고 온도분포 해석을 수행할 수 있다. 2차원 해석영역에서 시간 변화에 따른 과도상태의 열전달 지배방정식은 다음과 같다.

Fig. 3

Solution domain in temperature-distribution analysis model

여기서, T는 온도, k는 열전도 계수, Q_i는 발열량, ρ는 밀도, c_i는 비열, t는 시간이고, 초기조건은 상온으로 설정할 수 있다.

경계조건은 다음과 같은 표면에서의 대류 조건을 사용하였다.

여기서, n은 경계면에 수직한 방향 벡터, h는 대류 계수, T_a는 대기 온도, T_s는 경계면 온도이다.

파이프의 온도 분포 해석에서는 전자기장 해석과 동일한 기하학적 모델을 적용하였다. 용접 속도는 파이프 이송 속도인 220 mm/s로 설정하였고, 밀도, 비열, 열전도도 등 파이프 소재의 열적 물성치는 온도변화에 따라 변화되는 특성을 고려하였다[14].

2.3 유한요소해석

두께 15.4 mm의 X65 급 대강 소재로써 외경 609.6 mm, 내경 578.8 mm의 파이프를 제작하는 공정에 대한 전자기장 및 온도 분포에 관한 모델을 유한요소해석 모델로 구현화하였다.

전자기장 유한요소해석 모델의 해를 구하기 위해서 상용소프트웨어인 FLUX3D를 사용하여 3차원 정상상태 AC 해석을 수행하였다[15]. 해석영역은 용접부 주위는 크기를 줄여 분할하여, 특히 V-수렴점 주위에서 0.01 mm 크기로 사면체 요소 약 450,000개의 요소로 분할하였다.

공급 전압, 전극 거리(전극에서 V-수렴점 위치까지의 거리), 대강 V 각도 등의 변화에 따른 공정 현상을 해석하고자 하였고, Table 1은 입력 변수의 변화를 보여 주고 있다. 공급 전압의 변화에 대한 해석에서는 전극 거리와 대강 각도를 각각 130 mm와 6^o로, 거리 변화의 해석에서는 전압과 각도를 각각 240 V와 6^o로, 각도 변화의 해석에서는 전압과 거리를 각각 240 V와 130 mm로 일정하게 두고 해석하였다.

Table 1

Changes of welding variables

전자기장 해석에서는 전류밀도, 자속밀도, 발열량 분포를 계산하여 용접부에서의 전류 밀집 효과를 분석하고, 소재 표면에 작용하는 자속밀도의 크기를 파악하며 온도 분포 해석을 위한 열원을 도출한다.

HF-ERW에서 파이프의 온도 분포에 대한 유한요소해석 모델의 해는 상용소프트웨어인 ABAQUS를 활용하여 구하였다[16]. 해석 영역의 분할은 사각형 요소를 사용하였으며 접합부 부근에는 조밀하게 하여 전체 27,480개 요소로 분할하였다.

강관의 3차원 정상상태 전자기장 해석으로 먼저 발열량을 구하고, 각 위치에서의 발열량을 이용하여 파이프 길이 방향 각 2차원 단면에서의 천이 과정 온도 분포 해석을 수행하였다. 이를 위해, 해석영역 내의 파이프 단면 형상을 시간에 따라 변형하여 파이프 이송에 따라 모서리 간격이 줄어들어, 모서리가 접합되는 V-수렴점 위치의 형상에 도달하게 될 때까지의 열전달 과정을 모사하여 해석하였다. 온도 해석 결과 분석을 위해, 용접부에서 용융 온도 이상에 도달하는 크기를 추출하였으며 두께 중앙부의 폭을 용융 폭으로 나타내었다.

3.1 전자기장 해석

Fig. 4는 파이프의 전류밀도(Current Density)에 대한 해석 결과를 보여 주고 있다. Fig. 4(a)는 두 전극과 대강 소재 표면에 나타난 전류밀도 분포를 보이고 있다. 소재 모서리 표면에서 전류밀도가 소재 두께에 걸쳐 대칭적 분포를 보이며, 용접 위치에 근접할수록 크기가 증가하는 것을 알 수 있다.

Fig. 4

Current density distribution on the workpiece

Fig. 4(b)는 대강 소재의 한쪽 끝에서 다른 쪽 끝까지 상부 모서리를 따르는 경로 그리고 하부 모서리를 따르는 경로의 각각에 대하여 전류밀도를 나타내고 있다. 주위에 좌측 전극이 존재하고 있는 80 mm 근방에서 전류밀도가 크게 증가하고, 용접 위치라고 할 수 있는 V-수렴점 근방에서는 전류밀도가 급격히 증가하고 있다. 우측 전극이 존재하는 340 mm 근방에서부터 전류밀도는 감소하게 된다. 용접부 주위에서의 급격한 증가는 고주파 전류가 반대 방향으로 흐르는 두 표면이 매우 근접해, 표피효과에 더하여 근접효과가 뚜렷하게 나타나는 결과로 판단된다. 모서리 간격이 어느 정도 벌어져 있는 V-영역에 비하여 용접 위치에서 전류가 표피로 2배 이상 더 밀집되고 있다. 전극은 소재의 외부 표면에서 위치하고 있으나, 전류밀도 분포는 상부와 하부모서리에서 유사하게 나타나고 있다. 전류밀도의 분포로 인하여 자속밀도 분포와 발열량 분포도 유사한 형태로 나타나게 된다.

Fig. 5는 전극에 인가하는 입력 전압의 크기를 220 V에서 300 V까지 변화하는 경우, V-수렴점에서, 소재 두께 방향에 분포하는 전류밀도, 자속밀도(Magnetic Flux Density), 발열량(Generated Heat Density)의 변화를 나타내고 있다. 여기서 두께 방향은 파이프 소재의 외면에서 내면까지의 반경 방향을 나타낸다. 모든 결과에서 두께 중앙에서는 대체로 일정한 분포를 보이나, 소재의 표면과 이면에 해당하는 가장자리에서 크기가 급격히 증가하고 있다. 전압이 증가할수록 전류밀도, 자속밀도, 발열량은 비례적으로 증가하고 있으며 가장자리에서의 비례 폭이 중앙부보다 더 크게 나타나고 있다.

Fig. 5

Results of the analysis along the thickness of workpiece at V-convergence point with change of input source voltage

발열량은 소재에 흐르는 전류의 제곱에 비례하는 크기를 가지므로 발열량의 분포는 전류밀도의 분포와 유사하게 된다. 발열량 분포로서는 용접 위치에서 소재의 두께 방향 중심부보다 외면과 내면에서 가열이 크게 발생하게 되는 것을 예측할 수 있다.

공급 전압에 의해 양 전극 사이에 흐르는 전류는 전압의 변화에 따라 Table 2와 같이 계산되었다. 이는 전압에 비례하며 소재 표면에서의 임피던스와 와전류에 의해 결정된 크기로 판단된다.

Table 2

Current flowing through the pipe with input voltage

Fig. 6은 전극과 용접 위치 사이의 거리 변화에 따라 V-수렴점에서 소재 두께 방향에 분포하는 전류밀도, 자속밀도, 발열량의 변화를 보여 주고 있다. 두께의 중앙부에서는 일정한 크기를 나타내고, 가장자리에서는 크기가 급격하게 증가하고 있다. 거리가 근접할수록 전류밀도가 크게 증가하고 있고, 자속밀도와 발열량도 이에 비례하는 형태로 나타나고 있다. 이는 거리가 증가할수록 전기저항이 증가하고, 이로 인해 전류가 감소하고, 자속의 감소와 이로 인한 와전류의 감소 등에 기인한 결과로 판단된다. 특히 용접 위치 부근에서의 감소가 크게 나타나고 있다. 한편, 거리가 130 mm 이상에서는 거리에 의한 효과가 크게 줄어들고 있다.

Fig. 6

Results of the analysis with change of distance between contactor and V-convergence point

Fig. 7은 대강 소재가 용접 위치로 진입하는 각도의 변화에 따라, 용접 위치에서의 소재 두께 방향으로 나타낸 전류밀도, 자속밀도, 발열량 분포를 보여주고 있다. 전류밀도는 소재의 두께에 걸쳐 대체로 일정하게 나타나나, 외면과 내면 주위에서는 매우 크게 나타나고 있고 두께 중심을 기준으로 외부와 내부가 대칭적 분포를 보이고 있다. 대강의 진입각도가 작아질수록 전류밀도가 크게 나타나고 있다. 자속밀도와 발열량의 분포는 두께 방향으로의 분포, 그리고 대강의 진입각도 등의 변화에 대하여 전류밀도의 분포와 유사한 형태를 나타내고 있다.

Fig. 7

Results of the analysis with change of V-angle of skelp

3.2 온도 분포 해석

파이프의 단면을 파이프의 이송에 따라 변형하면서 2차원 과도상태 온도 분포 해석을 수행한 결과, 파이프 소재 두께 방향에 나타나는 온도 분포를 시간변화에 따라 얻을 수 있었다. Fig. 8은 파이프의 이송 속력을 고려한 시간의 경과에 따라 단면의 형상이 변화되고 단면 내의 온도 분포가 변화되는 현상을 보여주고 있다. 대강의 모서리는 용접 위치로 근접될수록 간격이 좁아지고 용접 위치에서는 두 모서리가 접합된 상태를 보여 주고 있다. 대강이 용접 위치로 근접할수록 단면 내의 온도가 높아지는 분포 상태를 나타내고 있다.

Fig. 8

Temperature distribution around the edge of the skelp during welding with elapse of time

Fig. 9는 용접된 소재의 두께 방향에서 실제 용접부의 한 쪽과 이에 해당하는 해석의 결과를 동시에 보여 주고 있다. 공급전압이 240 V, 전극과 용접 위치 사이 거리가 130 mm, 대강의 진입 각도를 6^o, 파이프의 이송 속도를 220 mm/s로 수행한 용접결과를 비교하였다.

Fig. 9

Comparison between actual welded surface and predicted temperature distribution of the model surface

결과에서 위는 파이프 소재의 외부 표면, 아래는 내부 표면에 해당한다. 두께의 중앙부에서 실제 용융 반폭은 0.3 mm이고, 해석 결과의 용융 반폭은 0.2 mm로 유사하게 나타나고 있다. 그러나 두께의 가장자리 부분인 파이프의 외면과 내면에서는 실제와 예측의 결과에서 차이가 발생하고 있다. 실제 용접부는 압축롤에 의한 업셋을 기계 가공하여 제거한 후의 결과이고, 해석에서는 전자기력에 의한 소재 밀착 현상과 업셋에 의한 용융부의 유동을 고려하지 않은 결과이다. 이로써 가장자리 부근에서의 결과에서는 다소 차이가 발생한 것으로 판단된다.

Fig. 10은 유한요소해석으로 예측한 용접 조건에 따라 V-수렴점에서 두께 중앙부 용융 폭의 변화를 나타내고 있다. Fig. 10(a)는 전극에 인가하는 입력 전압을 220 V 조건에서 300 V 조건까지 변화할 때 용융 폭의 변화로, 0-0.8 mm까지 비례적으로 증가하고 있음을 보여 주고 있다.

Fig. 10

Fused width at mid-thickness position according to the welding variables

Fig. 10(b)는 전극의 위치가 용접 위치로부터 70 mm 이격된 조건에서 140 mm 이격된 조건까지 변화한 경우, 용융 폭의 변화가 1.08 mm에서 0 mm로 나타나고 있다. 특히 70-130 mm까지는 용융 폭은 거리 변화에 비례적으로 감소하고, 거리가 150 mm인 경우는 용융부가 나타나지 않았다. 거리에 따른 용융 폭의 변화는 전자기장 해석에서 나타난 거리에 따른 전류밀도와 발열량의 분포로써 경향성을 유추할 수 있다.

Fig. 10(c)는 용접 위치로 진입하는 대강의 진입각도가 3^o인 조건에서 7^o인 조건까지 변화할 때 용용지 폭은 0.46 mm에서 0.12 mm로 비례적으로 축소되고 있다. 이는 전자기장해석의 결과에서 예측된 바와 같이 대강 모서리의 진입 각도가 커질수록 전류밀도와 이와 인한 발열량이 작아지는 결과에 기인한다고 판단할 수 있다.