강성행렬 요소

왜 배우는가: 강성행렬 요소는 구조물 전체 강성행렬을 조립하는 출발점이므로 유한요소해석과 강성행렬법의 기본 원리를 이해하는 데 필수적이다. 트러스 요소 강성행렬의 유도 과정은 기술사 시험에서 반복적으로 출제되며, 강성행렬 각 원소의 물리적 의미를 묻는 용어형 문제와 연결된다.

출제 이력: 제122회 3교시 6번, 제110회 1교시 9번

1. 한 줄 정의와 핵심 키워드

요소강성행렬은 요소 절점변위와 절점력의 관계를 나타내는 행렬로, 부재의 탄성계수·단면적·길이와 방향에 따라 결정된다.

핵심 키워드

• 요소강성행렬(Element Stiffness Matrix)
• direct stiffness method(직접강성법)
• 방향코사인(direction cosine, λ=cosθ, μ=sinθ)
• 절점변위와 절점력
• 좌표변환행렬(transformation matrix)
• 축강성(axial stiffness, EA/L)
• 전체좌표계(global coordinate system)와 국부좌표계(local coordinate system)

2. 개념 풀이

쉽게 말하면 요소강성행렬은 부재 하나가 끝점에서 얼마만큼 힘을 버티는지를 숫자 표로 적은 것이다. 끝점을 조금 밀면 힘이 얼마나 생기는지를 한 칸씩 기록한다.

정식으로는 요소의 절점력 벡터 {F}와 절점변위 벡터 {u} 사이의 관계를 {F}=[K]{u}로 나타내는 행렬 [K]를 요소강성행렬이라 한다. 트러스 요소는 축방향 힘만 전달하므로 국부좌표계에서 축강성 행렬은 EA/L 을 이용해 쉽게 쓸 수 있다.

요소강성행렬은 탄성 거동을 가정한 선형 관계식이다. 행렬의 각 항은 한 자유도의 단위 변위가 다른 자유도에 주는 힘을 뜻한다.

좌표변환행렬은 국부좌표계의 변위와 전체좌표계의 변위를 연결한다. 2차원 트러스에서는 방향코사인 λ와 μ로 구성된다.

방향코사인은 부재축이 전체 X축, Y축과 이루는 각의 코사인과 사인 값이다. 부재의 기울기를 나타내는 수치이다.

직접강성법은 모든 요소의 강성행렬을 절점 번호에 맞춰 더해 전체 강성행렬을 만드는 방법이다. 요소강성행렬은 이 과정의 최소 단위가 된다.

3. 등장 배경과 필요성

초기 골조 해석은 부재 단위의 처짐각법이나 모멘트분배법으로 수작업 처리되었다. 구조물의 자유도가 많아지면서 부재 강성을 행렬로 체계화하고 컴퓨터로 연산하는 방식이 필요해졌다.

유한요소법과 강성행렬법은 이러한 요구에 따라 발전하였다. 요소를 독립적으로 정식화하고 전체 구조에 조립하는 절차가 표준화되었다.

트러스 해석은 부재가 축력만 받는다는 단순 가정으로 시작하기 좋다. 따라서 요소강성행렬 유도의 기본 예제로 널리 쓰인다.

좌표변환을 도입하면 임의 방향의 부재를 같은 식으로 처리할 수 있다. 이 점이 실무 프로그램에서 요소강성행렬을 공통 형식으로 쓰는 이유이다.

4. 구성요소·원리·절차

강성행렬 요소는 한 부재의 단부 힘과 단부 변위를 연결하는 행렬이다. 트러스 요소는 축방향 힘만 전달하므로 부재 축 방향 강성 하나로 요소 강성을 표현한다. 이 값을 전체 좌표계로 변환한 것이 요소 강성행렬이며, 이후 전체 강성행렬 조립의 기본 단위가 된다.

요소 강성 의 핵심은 축방향 강성 k = \frac{EA}{L} 이다. 여기서 E 는 탄성계수(MPa), A 는 단면적(㎟), L 은 부재길이(mm)이다. 국부좌표계에서 축방향 2자유도 요소의 강성행렬은 \begin{bmatrix}k & -k\\ -k & k\end{bmatrix} 이다.

방향여현 은 부재축이 전체 X축과 이루는 각 \theta 에 대해 \lambda = \cos\theta, \mu = \sin\theta 로 정의한다.

원리는 다음과 같다. 국부좌표계의 축력 F_{axial} = k(u_{2}-u_{1}) 를 전체좌표계 성분으로 분해하면 힘의 방향이 \lambda, \mu 로 나뉜다. 따라서 요소 강성행렬은 [k]=\frac{EA}{L}\,[T]^T[k_{local}][T] 형태로 얻는다. 변환행렬 [T] 는 방향여현으로 구성된다.

요소 강성행렬은 국부강성을 방향여현으로 좌표변환하여 얻는다

전개 결과는 다음 4×4 행렬이다.

$[K]=\frac{EA}{L}\begin{bmatrix}\lambda^2 & \lambda\mu & -\lambda^2 & -\lambda\mu\\ \lambda\mu & \mu^2 & -\lambda\mu & -\mu^2\\ -\lambda^2 & -\lambda\mu & \lambda^2 & \lambda\mu\\ -\lambda\mu & -\mu^2 & \lambda\mu & \mu^2\end{bmatrix}$

이 행렬은 대칭이며, 행과 열의 합이 0이므로 강체이동에 대해 힘이 발생하지 않는다. 대각 항은 직접강성, 비대각 항은 연성강성이다.

절차 는 다음과 같다.

• 각 부재의 절점좌표로 길이 L 과 방향여현 \lambda, \mu 를 구한다.
• 재료·단면 값으로 \frac{EA}{L} 를 계산한다.
• 위 식으로 부재별 4×4 요소 강성행렬 [K] 를 만든다.
• 절점의 전체 자유도 번호에 맞춰 전체 강성행렬에 직접강성법 으로 조립한다.
• 지점조건을 적용하여 [K]\{u\}=\{F\} 를 풀고 부재력을 구한다.

5. 관련 법령·기준·표준

강성행렬 요소는 법령 조문이 아니라 구조역학과 유한요소법의 이론 체계에 속한다. 설계 실무에서는 해석 결과를 설계기준의 허용응력·처짐 한계와 대조한다. 국내 구조설계 기준은 「건축구조기준(KDS)」 체계를 따르며, 해석 방법 자체를 규정하는 수치는 아래 표의 값을 기준으로 삼는다.

강성행렬법은 1950년대 이후 항공 구조해석에서 발전했다. 국내에서는 구조해석 프로그램의 보편화와 함께 실무 표준 방법이 되었다. 최근에는 KDS 구조설계기준의 해석 조항에서 선형 탄성 해석의 기본 방법으로 다뤄진다.

6. 유형과 비교

강성행렬 요소는 요소의 종류와 좌표계에 따라 형태가 달라진다. 트러스 요소는 축방향 힘만 전달하므로 한 방향 강성만 갖고, 보 요소는 축력과 휨을 함께 받아 절점당 변위 자유도가 더 많다. 두 요소를 비교하면 요소 강성의 구조가 어떻게 달라지는지 드러난다.

위 표에서 트러스 요소는 축강성 \frac{EA}{L} 한 값이 강성의 뼈대가 된다. 보 요소는 같은 골격에 휨 항이 더해져 행렬 크기가 커진다.

트러스 요소는 전역 좌표계로 옮기기 위해 방향여현 \lambda=\cos\theta, \mu=\sin\theta를 쓴다. 이 점은 2차원과 3차원 트러스 모두 같은 패턴을 따른다. 3차원 트러스에서는 \lambda, \mu에 \nu=\cos\gamma가 더해져 6×6 행렬이 된다.