KEC plant pressure vessel training manual 4 of 5 권 PDF

KOREA ENGINEERING CONSULTING, INC.
PRESSURE VESSEL BASIC CONCEPT

PROGRAM : NISA(NISA SOFTWARE, INC. USA) FEA THEORETICAL COURSE
유한요소법의 이론 및 개념 4-1
Pressure Vessel Basic Concept

Training Manual(4권)
(FEM Theoretical/Calculation)
Korea Engineering Consulting, Inc.

Tel) 02-3486-9324/9377, Fax) 02-3486-9325, hkecc@hkecc.co.kr, http://www.hkecc.co.kr
KOREA ENGINEERING CONSULTING, INC. PRESSURE VESSEL BASIC CONCEPT
==== Contents ===
[1권Book] : Theoretical Strength Calculation & FEA

1.Stress-Strain Curve(Experimental) ............................................................................................ 6
2. Simple Stress(σ) & Strain(ε) Calculation ................................................................................. 8
3. Thin Cylinders & Shells............................................................................................................ 9
4. Thick Cylinders & Shells ........................................................................................................ 16
5. Theories of Elastic Failure ...................................................................................................... 23
6. Torsion Theory........................................................................................................................ 24
7. Basic Lug Calculation ............................................................................................................. 26
8. Basic Skirt Calculation.................................................................................................... 28 ~ 31
8-1. Skirt Stress Calculation.................................................................................................... 31
8-2. Skirt Anchor Bolt Stress Calculation ............................................................................. 32
9. Buckling Equation & Analysis Approach ............................................................................... 32
9-1. Buckling Equation & Analysis Approach........................................................................ 35
10. Circular Flat Plate and Head with Uniform Loading ............................................................ 36
10-1 Roark’s Formulas for Stress & Strain ............................................................................. 38
10-2 ASME Code Sec.VIII, Div.2, Appendix4....................................................................... 40
10-3 Perforated Plate Bending Stress Calculation .................................................................. 44
10-4 Doom Head(Dished/Ellipsoidal/Spherical Heads) under External Calculation.............. 46
10-5 Perforated Plate Bending Stress Calculation(H/E Tubesheet) ........................................ 52
11. Blind Flange with Uniform Loading ..................................................................................... 56
12. Internal Components Support Ring under Packing & Pressure Drop Loading............. 58 ~ 62
12-1 Support Ring under Packing & Pressure Drop Loading ................................................. 59
12-2 Intermediate Head Weldment Calculation...................................................................... 63
12-3 Partition Plate Thickness Calculation ............................................................................. 69
12-4 Nozzle Load Application ................................................................................................ 71
12-5 Stress Indices on Nozzle per ASME, Div.2 Para AD-560.7........................................... 74
12-6 Flange Torque Calculation for Sealing Flanges.............................................................. 91
12-7 Bolt Thread Shear/Bending Stress & Nut Length Calculation ....................................... 94
12-8 Girth Flange Calculation................................................................................................. 96
12-9 Buckling(Wrinkling) of Ellipsoidal Head(Internal Pressure) ....................................... 104
12-9-1 Compressive Pressure Calculation Internal Pressure) ........................................... 104
12-9-2 Head Wrinkling Thickness Calculation(ASME Codecase 2261).......................... 110
12-9-3 Stiffening Ring Calculation under External Pressure(ASME UG29) ................... 114
12-9-4 Internal Support Beam Calculation(Pressure Drop + Static Weight) .................... 117
12-10 Linearization of Stress Results and Classification ...................................................... 121
13. ASME Allowable Stress Criteria Concept ...................................................................... 126
14. Vessel Equipment Analysis Approach ................................................................................ 127
14-1. Stress Analysis Approach(based on ASME) ............................................................... 129
14-2. Stress Evaluation Sample 1(based on ASME) ........................................................... 133
14-3. FE Stress/Fatigue Model Sample 2(2D Axisymmetric Solid) ................................... 139
14-4. Stress Analysis Approach(based on ASME) ............................................................... 143
14-5. FE Stress/Fatigue Model Sample 3(3D Shell) ............................................................. 144
14-6. Fatigue Analysis Approach(based on ASME) ............................................................. 147
14-7. Crack/Fracture Analysis Approach(ASME/PD/API) .................................................. 147
14-8. Nozzle Local Stress Analysis....................................................................................... 153
14-9. Agitator Rotating Equipment Analysis ........................................................................ 154
[[2권Book] : Theoretical Strength Calculation & FEA

Appendix ″A″................................................................................................................................ 6
″Additional Vessel/HE Stress Calculation & Information analyzed by NISA″
App 1) 실린더에 외압을 받는 경우 최적두께 계산............................................................ 7
App 2) 고온에서 허용응력이 ASME, Sec.II, Part D에 없는 경우 ..................................... 11
App 3) Internal Support Beam의 Bending Stress & Displacement ........................................... 16
App 4) Spread Beam Calculation due to Rigging Condition ...................................................... 18
App 5) Vibration Analysis & Design of Tall Towers ................................................................. 28
1) Sample 1 : Zorrila’s Function..................................................................................... 33
2) Sample 2 : Zorrila’s Function..................................................................................... 35
3) Sample 3 : ASME STS-2006 Code ............................................................................ 38
App 6) Seismic Spectrum Analysis(UBC/ASCE/IBC Code) ..................................................... 50
App 7) Tubesheet Stress Analysis............................................................................................... 57
App 8) Conical Section Calculation under Nozzle Load ............................................................ 66
App 9) Theoretical procedures of Expansion Bellows Calculation ............................................ 68
App 10) Expansion Bellows Stress & Fatigue Calculation......................................................... 83
App 11) Seismic Responses Spectrum Analysis....................................................................... 100
App 12) Reaction Load on Multi-Saddle Supports ................................................................... 106
App 13) Stress Linerization Analysis........................................................................................ 110
App 14) Imperfection Shape on Equipments ............................................................................ 160
App 15) Spring Constants Calculation(Perforated Plate).......................................................... 172
App 16) Thickness & Deflection Calculation of Impingement................................................. 178
[3권Book] : Theoretical Strength Calculation & FEA

Appendix ″B″ ................................................................................................................................ 6
″Additional Vessel/HE Stress Calculation & Information analyzed by NISA″
App 1) Stability Check for Girth Flange/Bolt under Erection Conditions.................................... 7
App 2) Stability Check for Girth Flange/Bolt under Erection Conditions.................................. 17
App 3) Stability Check for Tube/Header Pipe under Thermal Expansion .................................. 22
App 4) Wind Guide Design for Bending Moment due to Wind Load ........................................ 25
App 5) Knee Brace Support Buckling Stability Analysis ........................................................... 32
App 6) Saddle Bending Checking under Sea Motion Load ........................................................ 42
App 7) Tube to Tubesheet Joint Welding Checking & Fatigue .................................................. 54
App 8) High Temp Creep Damage Procedure & Calculation..................................................... 59
App 9) Bending/Buckling Checking of Spread Beam under Erection........................................ 73
App 10) Bending Checking of Circular Plate subjected to Fluid Load....................................... 81
App 11) Tubesheet Calculation(Miller’s Method) .................................................................... 85
App 12) Thickness Calculation under External Pressure ............................................................ 93
App 12-1) Cylinder Thickness Calculation under External Pressure.......................................... 96
App 12-2) Spherical Thickness Calculation under External Pressure....................................... 103
App 12-3) Cylinder Thickness with Stiffener under External Pressure .................................... 107
App 13) Heating Up Force(Air/Saturated Water) in Isolated Space......................................... 120
App 14) App 14) Main Shell Thickness Calculation(ASME Code 2006 & 2007) ................. 138
[4권Book] : FEA(Finite Element Analysis) Theoretical Course

Appendix ″C″ ................................................................................................................................ 4
″Finite Element Theoretical & Analysis by NISA″
z Theoretical Course ″A″
1. 유한요소법(Finite Element Method)의 개념 ..................................................... 1-1
2. 유한요소법의 단계별 적용 .......................................................................... 2-1
3. 요소의 형상 함수..................................................................................... 3-1
4. 요소 형상 및 특징 ................................................................................... 4-1
5. 요소의 생성기법 및 결과 분석 방법 ............................................................. 5-1
z Theoretical Course ″B″

1-1. 선형용수철 요소(Linear Spring Element) ..................................................... 1-1
1-2. 축 방향 보 요소(Axial Bar Element) .......................................................... 2-1
1-3. 굽힘 보 요소(Bending Beam Element) ....................................................... 3-1
1-4. 평면 구조 요소(Plane Frame Element) ....................................................... 4-1
1-5. 평면 격자 요소(Plane Lattice Element) ...................................................... 5-1
1-6. 평면 요소(Plane Element)....................................................................... 6-1
z Theoretical Course ″C″

[예제1] 스프링 계로 이루어진 구조의 유한요소해석 ........................................... C-1
[예제2] 트러스(Truss) 구조물의 유한요소해석 ................................................... C-2
[예제3] 빕(Beam) 구조물의 유한요소해석 ........................................................ C-3
[5권Book] : Crack Propagation & Facture Mechanics

Appendix ″D″................................................................................................................................ 7
1. Crack Propagation ..................................................................................................................... 7

1.1 Theoretical Descriptions ..................................................................................................... 7
1) Paris Law .......................................................................................................................... 7
2) Forman Law ...................................................................................................................... 8
3) Elber Law.......................................................................................................................... 9
4) Walker Law..................................................................................................................... 10
2. Fracture Mechanics ................................................................................................................. 37
2.1 Theoretical Descriptions ................................................................................................... 37
2.2 Stress Intensity Calculation ..................................................................................... 80 ~ 228
Appendix ″C″
″Finite Element Theoretical & Analysis by NISA″
Theoretical Course ″A″

1. 유한요소법(Finite Element Method)의 개념
2. 유한요소법의 단계별 적용
3. 요소의 형상 함수
4. 요소형상 및 특징
5. 요소의 생성 및 결과 분석
This page internally left blank

1. 유한요소법(Finite Element Method)의 개념
1.1 서론
유한요소법(Finite Element Method, 有限要素法)은 1930년대에 개발된 행렬(Matrix)의

이론을 사용하여 1960년대에 영국과 미국에서 개발되었다. 이론의 개발은 일찍이 되었으나
실제 엔지니어가 사용한 시기는 전산기의 발전이 이루어진 1970년대 이후이다.
유한요소법은 초기에는 공학 및 물리학에서 많이 사용되었으나 현재에는 공학의 여러 분

야에서 전산기의 급속한 발전과 더불어 널리 사용되고 있다.
유한요소법은 초기에는 고전적인 이론식으로 해석할 수 없는 복잡한 형상의 구조물의 응

력해석을 위해 개발되었으나 현재에는 이론의 발전으로 연속체역학(Continuum Mechanics)
등의 모든 분야에서 널리 사용되고 있다.
기술의 발전과 더불어 더욱 복잡한 형상이나 하중 상태 및 물성치 들을 포함하는 문제를

해결할 필요성이 점점 커지면서 더욱 중요한 하나의 구조물의 해석 절차를 반영한 설계에
적용하는 Tool로 인식되게 되었다.
좀더 유한요소법을 설명하기 전에 처음 입문자를 위하여 전산기를 이용해 현재 공학에 적

용하고 있는 CAD, CAM 및 CAE의 차이점 및 적용 분야에 대하여 설명하면 다음과 같다.
단) 현재 CAD, CAM 및 CAE가 복합적으로 통합화 된 프로그램도 있고, CAD 및 CAE만

접목이 된 프로그램도 있기 때문에 본 내용에서는 그 프로그램이 추구하고 있는 큰 방향에
서 접근하여 예를 들었다.
[CAD] : ″Computer Aided Design″

일반적으로 보편화 되어있는 구조물의 도면 작성 단계를 기존의 수 작업에 의해 도면화
단계를 전산기(컴퓨터) 이용한 설계를 CAD(Computer Aided Design)라 한다. 이는 단순히
구조물의 차원을 표시하는 방법으로 최근에는 구조물의 형상을 3차원적인 차원(Solid
Model)에서 접근하는 기법도 도입되고 있다.
-. MICROSTATION, AUTOCAD, CADKEY, CATIA, I-DEAS, PRO/ENGINEER,
SOLID/WORK, SOLID/EDGE,… 등
[CAM] : ″Computer Aided Manufacturing″

컴퓨터를 이용한 구조물 형상을 3차원 Surface 데이터(x, y, z)로 처리하여 NC Machine
에 적용하여 형상을 가공하는 단계를 CAM(Computer Aided Manufacturing)이라 한다. 현
재에 상용화된 CAM 소프트웨어는 CAD기능을 대부분 동시에 포함하고 있다.
-. MDT, I-CAM, CATIA, PRO/ENGINEER, CADAM, Surface CAM,… 등
[CAE] : ″Computer Aided Engineering″

유한요소법 및 여러가지의 알고리즘이 적용되어 구조물의 응력, 변형, 유동 및 기타 역학
해석을 통하여 구조물의 안정성 및 최적화 구조물에 적용할 수 있는 단계의 적용을
CAE(Computer Aided Engineering)이라 한다.
적용되는 분야는 모든 공학 분야에 적용된다고 보아도 무리는 아니다. 특히 공과 대학에
서 숙지하고 있는 역학의 분야, 정역학(Linear/Nonlinear Static), 진동(Dynamic), 열전달
(Heat Transfer), 유체역학(Fluid Flow Dynamic), 유압공학(Pheumerical , 동력학, 의공학,
전자장 등의 구조 및 비구조 분야에 적용하여 구조의 최적화를 구현하는데 정보를 제공하
는 Tool이다.
유한요소법은 단어 자체의 의미에서 알 수 있듯이 복잡한 형상의 구조물을 이론에서 제공
하는 요소(선, 면, 입체의 형상)로 분할하고 요소사이에 공유하는 점, 선 및 면에서 연속이
론을 이용하여 수식을 만들고 해를 구하는 것이다.
더 확장이 되어 구조 및 비 구조에 적용하는 알고리즘은 유한요소법 외에 유한차분법
(FDM, Finite Differential Method), 유한체적법(FVM, Finite Volume Method), P-
Method(High Order Adaptive Method), 경계요소법(Boundary Element Method)등이 있다.
이러한 알고리즘의 적용범위는 결과의 정도 와 결과를 도출하기 까지의 시간에 영향을

받는데 유한요소법에 대하여 이론적인 접근 전에 이들의 대략적인 개념은 다음과 같다.
[FEM] : 유한요소법(Finite Element Method)

일반적인 구조 분야(건축, 토목, 자동차, 항공, 조선, 기계구조, 의공학 등), 열(온도)과
비 구조 분야[열전달, 유체의 유동(내부 및 외부 유동, 압축성, 비압축성), 전자장] 등의 모
든 분야에 적용되는 기본이 되며 오랜 기간 동안 검증 및 일반화된 Tool이라 할 수 있다.
단 유체의 유동 문제는 많은 요소 및 내부 알고리즘의 문제로 많은 시간이 소요되기 때문

에 엔지니어의 견해에서는 선호도가 떨어진다. 반면 어떤 유동의 형상도 모델이 가능하기
때문에 강점도 있다. 예를 들면 해석의 범위에 따라 상이하지만 구조에서 요소가 50,000개
정도면 많은 편이나 Windows NT Version으로 결과를 얻는 시간이 대략 20분 정도 소요되
나 유체의 유동 문제에 있어서는 유동의 특성이 따라 다르지만 대략 1주일이 소요되기 때
문에 선호도가 적은 이유가 여기에 있다. 결과의 정확한 정도는 모델에 따라 상이하여 비교
는 누구도 장담할 수 있는 부분은 아니다.
앞서 언급하였듯이 유한요소법은 거의 모든 분야에 적용이 가능하기 때문에 적용분야 별

로 한계를 짓기는 어렵지만 국내에서 많이 사용하는 분야별로 구분하면 다음과 같다.
그리고 아래의 소프트웨어는 대부분 외국에서 개발된 프로그램이며 국내에서 개발된 프로

그램은 유일하게 MIDAS 뿐이다.
-. 기계 구조물, 자동차, 항공, 조선 등의 입체적인 구조물(Solid) 등
예) NISA, NASTRAN(PATRAN, MARC), ANSYS, ABAQUS, FERMAS, FEMAP,
ADAMS(Kinematic), DADS(Kinematic),… 등
-. 건축, 토목 구조물, 터널 등
예), NISA, MAIDAS(국내개발), ADINA, DIANA, SAP2000, RUSAS, COSMOS, SBETA,
STRAND 7, STADIII,… 등
[FVM] : 유한체적법(Finite Volume Method)

일반적으로 유체의 유동 문제나 화학반응 등의 유체 관련 분야에 많이 적용되고 있으며
국내에서는 유한요소법에 비해 90% 정도 사용빈도가 높다. 또한 사출, 주물 등의 유동 현
상에 적용되기도 하지만 유한차분법 보다는 적용 빈도가 높은 편은 아니다.
반면 구조물의 역학적인 분석에서는 불가능하다.
-. 내부 및 외부 유동, 사출, 주물, 화학반응 등의 압축 및 비압축성 유체

예) CFD++, FLUENT, RAMPANT, STAR-CD, PHOENIX, FIDAP(FEM), PAM-FLOW,
FLOW3D, CDF2000, FLOTHERM, FLOTRAN, N3S, ICEM,… 등
[FDM] : 유한차분법(Finite Differential Method)

일반 유체의 유체유동 분야에도 적용은 하고 있지만 특히 주물, 사출, 쇳물 등의 비 뉴토
니안 유체(Non-Newtonian Fluid)에 많이 적용하고 있다. 뉴토니안 유체는 속도 구배와는
무관한 점성계수를 갖는 공기, 물 등이며, 혈액, 프라스틱, 오일, 수은 비 뉴토니안 유체
(Non-Newtonian Fluid)이며 속도 구배에 따라 점성계수가 상이한 유체를 통칭한다.
-. Moldflow, Flow-3D, C-Mold, MAGMA(주물),… 등

1.2 유한유소법의 장, 단점
일반적으로 유한요소법의 장점과 단점은 다음과 같다. 다음에 간략하게 설명하는 장단점
은 유한요소법을 이용한 소프트웨어를 사용하려고 하는 사용자는 항상 머리 속에 인식하고
있어야 한다.
[장점] :
z 기하학적 형상, 하중 및 구속조건등에 대한 제한이 없다.
유한요소법은 구조물을 작은 요소를 사용하여 표현하므로 형상등에 대한 제한조건이 없다.
다만 복잡한 형상인 경우 자세하게 표현하려면 요소의 수가 기하급수적으로 늘어나므로 계
산시간이 많이 소요된다. 초보자와 경험자의 차이는 실 구조물을 어느 정도 단순화하느냐에
달려있다. 본 내용에 수록된 많은 이론적인 접근 및 실 문제를 분석해 보면 같은 구조물 이
라도 엔지니어마다 그리고 같은 엔지니어라도 모델링 방법(요소/절점)에 따라 다소 상이한
결과가 도출이 될 수 있기 때문에 많은 경험 과 결과분석 능력이 중요하리라 본다.
z 여러 가지 복합재료로 이루어진 연속체에도 사용이 가능하다.

구조물을 요소를 이용하여 분할하므로 구조물이 여러 가지의 재료로 이루어진 경우에도
아무런 제한조건 없이 사용할 수 있다. 자동차에서 사용하는 플라스틱과 강(鋼)으로 구성된
부품에 적용하는 데에도 아무런 어려움이 없다.
z 재료의 물성과 구조물 거동의 비선형에도 적용이 가능하다.

유한요소법에서는 재료의 물성치는 시간이나 온도 및 하중에 따라 변하지 않고 일정하다
고 가정하여 해석을 한다. 그러나 구조물을 이루는 재료들은 작용하는 응력과 하중의 경로
에 따라 탄성계수 및 탄성한도가 변한다(응력경화, 재료의 항복). 혹은 온도에 따라 탄성계
수, 열팽창계수 및 열전달계수(Thermal Conductivity)가 변하게 된다. 응력해석의 경우, 하
중이 작용하는 동안 구조물에 작용하는 응력에 따라 사용하는 탄성계수가 변하게 되며 정확
한 결과를 구하기 위해서는 응력에 따라 산출된 탄성계수를 사용하여야 한다.
이러한 분야를 재료의 비선형(Material Nonlinear)이라 한다. 다른 한가지는 일반적으로
고전적인 이론이나 유한요소법은 구조물의 크기에 비해 변형의 양이 상당히 작다고 가정을
하고 이론을 개발하였다. 예를 들어 보의 경우 변형의 양이 길이의 1/10보다 작은 경우일
때 이론적인 계산이 잘 맞게 된다. 그러나 변형이 기준보다 크게 되면 이론적인 계산과 실
제값은 일치하지 않게 된다.
그러므로 변형이 큰 경우에는 기존의 이론식과 다른 식과 방법을 사용하여야 한다. 이러
한 경우를 기하학적 비선형(Geometry Nonlinear)이라 한다. 이러한 비선형의 경우 간단한
구조물의 경우라도 수계산으로는 구하는 데에 많은 시간이 소요된다. 유한요소법은 이러한
비선형의 경우 여러 가지 방법을 사용하여 정확한 해를 구할 수 있다.
z 응력, 좌굴, 진동, 열전달, 유동/열유동, 전자기장, 운동학 및 피로 및 해석이 가능하다.

유한요소법은 연속체역학에 속하는 모든 공학분야에 적용할 수 있다.
[단점] :
z 초보자도 사용하기가 쉽지만 결과의 신뢰는 공학지식에 비례한다.
위의 장점들은 반대로 단점이 될 수 있다. 선무당이 사람을 잡듯이 장점만을 기억해 사용
하는 경우, 일반적으로는 문제가 없지만 특이한 경우에는 전혀 다른 결과를 얻게 된다.
일반적으로 유한요소법은 소프트웨어를 사용해야 한다. 소프트웨어는 잘 알다시피 주어진
데이터에 대해 정확한 계산을 해 줄 수 있다. 또한 개발자가 경험을 이용하여 사용자가 실
수하기 쉬운 부분은 검토하는 기능을 포함하고 있다. 소프트웨어가 모든 가능성에 대한 검
토기능을 포함할 수는 없다.
만일 사용자가 탄성계수를 207000 MPa을 20700 MPa으로 입력을 잘못하는 경우 잘못된
결과를 출력하게 된다. 이 경우 결과의 신뢰는 사용자의 지식과 경험에 의해 좌우가 된다.
예를 들어 자동차를 움직이는 것은 남자인 경우 몇 가지만 배우면 가능하지만 운전을 잘 하
는 것은 아니다. 유한요소법도 동일한 개념으로 생각하면 된다.
z 국부응력해석에 적용할 시에는 많은 경험이 필요하다.

전체구조물을 해석하는 경우 모델링에 많은 시간이 소요되지만 그렇게 많은 경험과 지식
이 필요하지는 않는다. 전체구조물에서 일부분만을 해석하는 경우에는 구조물이 아무리 단
순하더라도 쉽게 생각하지 않아야 한다.
z 실구조물을 완벽하게 유한요소로 치환하는 것은 어려우므로 구조물을 단순화 한다.

이때 많은 경험과 지식이 필요하다. 이 부분은 다음 장에서 자세히 설명하겠다.
z 대형구조물인 경우 높은 성능의 전산기가 필요하다.

앞에서 설명하였듯이 요소의 수가 많아지면 형상을 정확하게 표현할 수 있지만 반대로
공유하는 점, 선 및 면은 많아지게 된다. 이런 경우 선형대수식의 크기가 커져 해를 구하는
데에 많은 시간이 소요되어 전산기의 성능이 중요한 변수가 된다.
위에서 설명한 내용은 반드시 읽어보아야 하고 기억하여야 한다. 아무리 많은 경험자일지

라도 실수를 하게 마련이다. 저자의 경우에도 시간이 촉박한 프로젝트를 수행하면서 급한
마음에 탄성계수를 잘못 입력하고, 결과검토를 하지 않고 보고서를 제출하여 추후에 많은
고생을 한적이 있다.

2. 유한요소법의 단계별 적용
유한요소법(FEM)은 응력, 진동, 피로 및 파괴 해석 등의 구조 분야 와 열, 유체 및 전자

기장 등의 비 구조 문제를 푸는 방법으로 널리 이용되고 있다.
유한요소법은 구조 해석을 위해서 처음 고안되었고, 오늘날 구조 해석 분야에 보다 효과
적으로 이용되고 있는 관계로 구조 문제에 대해서 유한요소법(FEM)으로 수식화하고 그 해
를 구하는 일련의 과정은 다음과 같다.
1단계, 기하학적 형상 및 연속체를 정립한다(전처리 단계, Pre/Processor)
응력, 변형 등을 분석하고자 하는 구조형상에서 기하학적 형상을 기반으로 해석 영역 부

분을 도출하고, 이를 간단한 기하학적 연속체(선, 삼각형, 사각형, 사면체 등)의 유한요소로
나누는 과정이다. 해석 전체의 과정 중 첫째 단계를 전처리 기능(Pre-Processor)이라 한다.
외부의 CAD 데이터를 이용하

거나 DISPLAY3/4의 기능을 이 점, 선, 면, 입체
용한 2, 3차원 형상 모델링
1단계
Point Mass, Line Element
전처리 과정
삼각형, 사각형, 입체요소 절점, 요소생성
(Pre Processor)
(1차,2차, 2차 요소 생성)
형상특성(Shell, Beam 등)
기계적 물성치(탄성계수 등) 특성치 입력
구조 및 비 구조 경계조건
정력학(선형,비선형), 2단계
동적해석(고유치, 진동 등)
열전달(전도, 대류, 열복사) 연산 과정(Solver)
기구동력학, 최적화해석 연산과정(Solver)
피로,크랙 및 파괴해석
유체해석, 전자장해석 등
3단계
후처리 과정
그래픽 및 Text로 결과 분석 그래픽 결과분석
(Post Processor)
그림[2.1] 유한요소해석의 단계별 기능

1) 기하학적 형상 생성(Geometric Modeling)
모든 형상을 생성하는 명령어는 소, 대문자 구분 없이 2차 및 3차원의 공간을 사용한다.
(1) 점(GRID)
기하학적형상을 구현하기 위한 3차원의 좌표 값을 보유하는 Point이다.
-. 예, 두 점의 좌표 값이 각각 GRID 1(X=0, Y=0) 과 GRID 2(X=0, Y=10) 일때;
• GRD, ADD, 1, 0/0/0 ; ″1″(점의 번호로 생략가능, GRD,ADD,,0/0/0)
• GRD, ADD, 2, 0/10/0
• •
G1(0,0) G2(10,0) 그림[2.2] Grid 생성하기
(2) 선(LINE)
점을 연결하는 연속선으로 2차원의 폐곡선이 존재할 경우 요소의 생성이 가능한 곡선의
형상이다.
-. 예, 두 점(Grid)을 연결하는 선을 생성 하고자 할 때;
• LIN, 2GD, 1, 2 ; ″1″ 과 ″2″[각각 점(GRID)의 번호]
• • • •
G1(0,0) G2(10,0) L1
그림[2.3] Line 생성하기
(3) 면(Patch)
점 및 선의 연결성에 의해 구현이 되는 평면 및 곡면의 형상이다.
-. 예, 두 선을 연결하는 면을 생성 하고자 하 때;
• PAT, 2LN, 1, 2 ; ″1″ 과 ″2″[각각 선(LINE)의 번호]
L2 P1
P1 P1
L1
그림[2.4] Patch 생성하기
(4) 입체(HYPERPATCH)
점, 선, 면 등의 기하학적 기본형상을 이용한 하나의 3차원 입체를 구현하는 형상이다.
-. 예, 두 면 및 세 면 이상을 연결하는 공간상의 입체를 생성하고자 할 때,
• HYP,2PAT, 1, 2 or, HYP,3PAT, 1, 2, 3 ; ″1″, ″2″, ″3″(각각 면의 번호)
P2
P3
H1 P2
H1 H2
P1
P1
그림[2.5] Hyperpatch 생성하기

2) 요소생성(Mesh Generation)
앞의 기하학적 형상(Grid, Line, Patch, Hyperpatch)을 이용하여 연속체로 분할하는데 절

점 및 요소가 생성된다.
(1) 절점(Node)을 이용한 요소생성 : Mass(Point) Element

구조물에 일정한 Mass 및 관성 모멘트가 있는 경우의 Point Mass에 적용된다.
-. 예, 한 절점에 Point Mass를 생성 할 때
• FEG, -----------
• •
N1 Mass 그림[2.6] Point Mass 생성하기
(2) 선(Line)을 이용한 요소생성 : Beam, Bar, Truss, Spring, Cable, Gap,…등
철 구조물과 같이 Line의 형상으로 이상화가 가능한 구조에 적용된다.
-. 예, 선에서 3차원 Beam의 요소를 생성 할 때(Line 1개에 무한대 요소 생성 가능);
• FEG, BAR, 1T4,,, 2, 12 ; 1T4(1 to 4 선의 정의), 2(요소개수),12(3차원beam NKTP)
L3 E6 E5
E7
L4 E4
L2
E8
E3
L1 E1 E2
그림[2.7] Line 요소 생성하기
(3) 면(Patch)을 이용한 요소생성 :

2차원 형상 : Plane Stress, Plane Strain, 2D Axisymmetric,…등
3차원 형상 : Thin Shell, Thick Shell, Membrane, Composite,…등
축 대칭 형상이나 2차원적인 해석이 가능한 구조물 및 3차원 Shell 구조에 적용된다.
-. 예, 축대칭 원통형의 경우 및 3차원 쉘의 구조에서 요소를 생성할 때;
• FEG, QUA, 1,,,1/6,3 ; ″1″(Patch 1),1/6(① 방향으로 1개/② 방향으로 6개 요소생성),
″3″(2차원 축 대칭 Solid NKTP)
• FEG, QUA, 1,,,2/2, 20 ; ″1″(Patch 1), 2/2(① 방향 2개/② 방향 2개 요소생성),
″20″(3차원 General Shell NKTP)
E6
E5
E4 E3
E4 P1 ②
②
E3 E1 E2
E2
①
E1
① 그림[2.8] 2/3차 Shell 요소 생성하기
(4) 입체(Hyperpatch)를 이용한 요소생성 : Solid, Composite,

3차원 입체의 체적을 유지하는 형상으로 2차원 축 대칭 및 3차원 Shell의 형상을 유지할
수 없는 경우의 구조에 적용된다.
-. 예, 3차원 입체에 직육면체요소를 생성할 때;
• FEG, HEX, 1,,,3/1/3,4 ; ″1″(Hyperpatch 1), ″3/1/3″(① 방향 3개/② 방향 1개
/③ 방향 3개), ″4″(입체요소 NKTP)
H1
②
그림[2.9] Hyperpatch 생성하기
②
• FEG, HEX, 1,,,1/3/1,4 ; ″1″(Hyperpatch 1), ″1/3/1″(① 방향 1개/② 방향 3개
/③ 방향 1개), ″4″(입체요소 NKTP)
H1 E3
P3 E2
③ E1
P1
① 그림[2.10] Hyperpatch 생성하기

3) 특성치 입력(Specific Properties)
앞의 기하학적 형상(Geometric Model)을 기반으로 다양한 2차원 및 3차원 요소를 생성한

다음 구조물에 대한 기계적 물성치(Mechnical Properties)나 외력적이 하중(Loading
Condition)을 포함한 경계조건(Boundary Conditions)등이 적용되는 전처리 단계의 마지막
단계이다.
단, 기계적 물성치는 기하학적 모델링 및 요소 생성 전에 입력이 가능하다. 이 의미는 요
소생성과 는 무관한 단계의 형식에서 입력이 되기 때문이다.
또한 하중 및 경계조건은 기하학적 형상 및 요소 생성이 되어있는 구조 형상에 입력을 해
야 하기 때문에 전처리(Pre Processor)입력 조건에서 최종 단계가 된다고 볼 수 있다.
(1) 기계적 물성치(Material Properties)

구조 및 비 구조해석에 적용되는 기계적 물성치를 간력하게 설명하면 다음 표[1.1]와
같다.
(2) 단면특성(Section Properties)

입체요소에는 적용이 되지않고 빔(Beam), 스프링(Spring), 케이블(Cable) 등 및 쉘
(shell) 구조에 입력이 되는 특성이다. 쉘의 구조에는 두께(Thickness)가 입력이 되고,
Beam의 구조에는 단면형상, 단면적, 관성모멘트 등이 입력이 된다.
(3) 경계조건(Boundary Conditions)

기하학적 모델링 및 요소의 생성이 끝나면 각 구조물에 존재하는 구속적인 조건의 입
력을 해야 하는데 이는 각 요소가 보유하고 있는 자유도(Degree of Freedom)로 구속
이 이루어 진다. 예를 들면 3차원 입체요소의 경우는 3자유도(UX, UY, UZ)이고, 3차
원 빔(Beam)의 경우는 6개의 자유도(UX, UY, UZ), RX, RY, RZ)를 보유하고 있다.
각 요소에 대한 자유도 및 해석 범위는 표[1.2] 나타나 있다.
(4) 하중조건(Loading Conditions)

경계조건의 입력 시점과 마찬가지로 모델링 및 요소가 생성된 다음 구조 부분에서는
균일분포하중(Uniform Pressure Loading), 집중하중(Concentrated Nodal Force), 모멘
트(Moment), 빔 하중(Beam Lading), 수두압(Hydrostatic Loading), 절점열하중(Nodal
Specific Temperature) 등 과 비구조해석에서는 온도(Temperature), 열량(Nodal,
Face Flux), 열발생(Heat Generation), 대류(Convection), 복사(Radiation)등을
GUI(Graphic User Interface)를 이용하여 적용한다.
표[1.1] 해석에서 요구되어지는 기계적 물성치 특성

요소 재질특성
해석범위 비 고
등방성(Isotropic Material)
탄성계수(Elastic Modulus, Ex) 각 1개씩 입력해도 내부적으
선형구조해석
포이송비(Poisson’s Ratio, νxy) 로는 모든 방향에 대해 같은
(Linear Static)
Density(ρ) : 자중이 외력으로 작용 시 값으로 해석 처리 됨
비선형구조해석 재료 비선형 기하학적 비선형 9 Stress-Strain Curve
(Nonlinear) (Material Nonlinear) (Geometric Nonlinear) -. Elastic-Perfectly Plastic
탄성계수(Elastic Modulus, Ex) -. Elastic-Linear Harding
포이송비(Poisson’s Ratio, νxy) -. Ramberg-Osgood
Density(ρ) -. Piecewise Linear Harding
소성응력(Yield Stress, σy) 9 Yield Criteria
응력-변형률 곡선(Stress-Strain Curve) -. Von-Mises(Equivalent)
소성이론(Yield Criteria) -. Tresca(Maximum Shear)
Work Harding Parameter(H) -. Mohr-Coulomb
Parameter for Harding rules(K) -. Drucker-Prager
초탄성물질(Hyperelastic) 9 Neo-Hooken
9 Classical Mooney-Rivlin
9 Generalized Mooney-R.
9 Swanson
9 Blatz-Ko
Potential or Strain-Energy Function 9 IHT
9 Biderman
9 Hart-Smith(Exponential)
9 Alexander(Exponential)
Creep 9 Stress Function Type
-. Norton
-. McVetty
-. Soderberg
-. Dorn
-. Johnson
Stress Function Type[F1(σ1)]
-. Garofaro
Time Function[F2(t)]
-. ORNL Law
Temperature Function[F3(t)]
9 Time Function
-. Andrade
-. Bailey
-. Mc Vetty
-. Graham & Walles
2D & 3D Contacts
Node Contact(Node to Node)
선형 및 비선형에 적용
Line Contact(Line to Line)
Surface Contact(Surface to Surface)
Thermal Conductivity(k)
Density(ρ)
9 Conduction
열전달 Specific Heat(Cp)
9 Convection
(Heat Transfer) Heat Transfer Coefficient(h)
9 Radiation
Emissivity(ε)
Stefan Boltzmann Constant((σ)
표[1.2] 각 요소 특성 및 자유도
NKTP Element Static Dynamic Buckling NL/NT DOF
1 2-D Plane Stress √ √ √ √ 2
2 2-D Plane Strain √ √ √ √ 2
3 Axisymmetric Solid √ √ √ √ 2
4 3-D Solid √ √ √ √ 3
5 3-D Thick Shell √ √ √ √ 3
7 3-D Composite Solid √ √ - - 3
9 3-D Hybrid Solid √ √ √ - 3
10 Membrane √ √ √ -
11 3-D Tapered Beam √ √ √ - 6
12 3-D Beam √ √ √ - 6
13 2-D Beam √ √ √ - 3
14 3-D Spar √ √ √ √ 3
15 2-D Spar √ √ √ - 2
17 3-D Translational Spring √ √ √ √ 3
18 2-D Translational Spring √ √ √ - 2
20 3-D General Shell √ √ √ √ 6
21 3-D Torsional Spring √ √ √ - 3
22 2-D Torsional Spring √ √ √ - 2
25 2-D Point Mass √ √ √ √ 2
26 3-D Point Mass √ √ √ √ 3
27 2-D General Point Mass √ √ √ √ 2
28 3-D General Point Mass √ √ √ √ 3
29 2-D Point Mass with Inertia √ √ √ √ 3
30 3-D Point Mass with Inertia √ √ √ √ 3
3-D Laminated Composite
32 √ √ √ √ 3
General Shell
33 3-D Sandwich Shell √ √ √ - 6
34 Axisym. Solid with Loading √ - - - 3
36 Axisymmetric Shell √ √ √ - 3
37 Axisym. Shell with Loading √ - - - 6
38 3-D General Spring √ √ √ √ 6
39 3-D General Beam √ √ √ √ 6
40 3-D Thin Shell √ √ √ √ 6
3-D Mixed Interpolation
41 √ √ √ - 6
General Shell
42 2-D or Axisymmetric Gap √ - - - 2
43 3-D Gap √ - - - 3
45 3-D Cable Element - - - √ 3
46 3-D Straight Pipe Element √ √ - - 6
47 3-D Elbow Element √ √ - - 6
48 3-D Dmaper - √ - √ 3
49 2-D Gap/Friction Element - - - - 2
50 3-D Gap/Friction Element - - - √ 3
NKTP : Element Type in NISA FEM S/W, NL : Nonlinear Static, NT : Nonlinear Direct Transient
DOF : Degree of Freedom
2단계, 요소에 대한 수식을 정립한다(Develop the Element Equation).
망(Mesh)을 이루고 있는 각각의 전형적인 요소에 대해서 수식을 유도한다. 이 과정에서

는 대체로 변분법(Variation Methods) 또는 그에 준하는 법칙을 적용하여 강성(Stiffness)행
렬 또는 유연성(Flexibility)행렬 및 절점(Node)의 상당 부하(Equivalent Loads)를 결정하게
되는데, 평형방정식(Equilibrium Equations)과 적합조건(Compatibility Conditions)이 만족되
도록 한다.
수식의 유도를 위해서는 접근 방법에 따라 변위법, 하중법 및 혼합법의 세가지 방법이 이
용되는데 이를 좀 더 구체적으로 설명하면 다음과 같다.
• 변위법(Displacement Method)
최소 위치 에너지 법(Minimum Potential Energy Theorem)에 기초한 이 방법은 변위를
미지수로 간주하고 적합조건이 만족되도록 요소에 대한 변위를 가정하여 최소 위치 에너
지 법칙을 적용, 각 요소에 대한 수식을 유도하는 방법이다.
이를 변위법 또는 강성법(Stiffness Method)이라고 하는데 결과적으로 얻게 되는 요소
에 대한 식은 근사 평형방정식(Approximate Equilibrium Equation)이 된다.
• 하중법(Force Method)
최소 보충 에너지 법칙(Minimum Complementary Energy Theorem)에 기초한 이 방법
은 힘을 미지수로 간주하고, 평형방정식이 만족되도록 요소에 대한 응력을 가정하여 최소
보충 에너지 법칙을 적용, 각 요소에 대한 수식을 유도하는 방법이다.
이를 하중법 또는 유연성법(Flexibility Method)이라고 하며 결과적으로 얻게 되는 식은
근사 적합방정식(Approximate Compatibility)이 된다.
• 혼합법(Hybrid Method 또는 Mixed Method)

특별한 에너지 법칙에 기초한 이 방범은 라이스너의 원리(Reissner’s Principal)와 같은
특별한 에너지 법칙을 이용하는 것으로서, 변위 와 응력을 부분적으로 가정하여 적용하므
로 혼합법이라고 한다.
이 상의 세 가지 방법은 문제의 성질에 따라서 한 가지 방법이 다른 방법에 비해서 보다

더 적합한 경우가 있으나, 대부분의 문제에 있어서는 변위법이 가장 적용하기가 간단하다.
따라서 구조해석용으로 개발된 대부분의 범용(General) 유한요소해석 프로그램은 변위법
을 사용하고 있다.
그래서 좀더 변위법에 대하여 설명하면 다음과 같은 단계가로 이루어진다.
1) 각 요소내의 변위를 국부 좌표(Local Coordinates)를 이용하여 특정의 연속함수로 가

정한다. 이를 위해서는 주로 다항식(Polynomials)이 사용된다.
2) 가정한 변위 함수를 절점의 값으로 표시한다.

이를 수식으로 표시하면,
{u} = [N ] {q} (1-1)
{u} : 변위 벡터(Displacement Vector)

[N] : 형상 함수 행렬(Shape Function Matrix)
{q} : 절점 변위 벡터(Nodal Displacement Vector)
3) 변형률과 변위(Strain-Displacement) 및 응력과 변형률(Stress-Strain)관계식이 만족

되도록 한다.
{ε} = [B] {q}, (1-2)

{σ} = [D] {ε}, (1-3)
{ε} : 변형률 벡터
{σ} : 응력 벡터
[B] : 형상 함수를 미분하여 유도되는 행렬
[D] : 탄성계수(Elastic Coefficient, Constants, Modulus)로 이루어진 행렬
열역학 법칙(Thermodynamics Origin)에 따라 대칭 행렬(Symmetric Matrix)
이 된다.
4) 요소에 대한 강성 행렬 및 상당 절점 부하를 일 또는 에너지 법칙을 이용하여 결정한

다. 여기서는 최소 위치 에너지 법칙을 이용하기도 한다.
전체 위치 에너지(Total Potential Energy)는 내부 변형 에너지(Strain Energy)와 외부
의 힘에 의한 위치 에너지(Potential Energy)의 합으로 정의된다.
(1) 내부의 변형 에너지(U)
⎡ ε T ⎤
U = ∫vol ⎢⎣∫0 {σ} {dε}⎥⎦ dV (1-4)
로 정의된다. 여기서 상첨자(Superscript) T는 전치(Transpose) 행렬을 뜻하며, 적분

은 부피(Voulme)에 대해서 행해진다. 특히, 혹의 법칙(Hook’s Law)을 따르는 선형
(Linear) 물질애 대해서는
U =
1
∫ {σ}T {ε} dV (1-5a)
2 vol
로 주어지는데 3차원 문제에 대한 변형에너지를 성분(Component)별로 나타내 보면,

직교 좌표계(Cartesian Coordinate System)에서는 다음과 같으며,
1
{σ x ε x + σ y ε y + σ z ε z + τ yz γ yz + τ xz γ xz + τ xy γ xy }dV
2 ∫vol
U =
직교 좌표계(Cartesian Coordinate System)에서는,
1
{σ r ε r + σ θ ε θ + σ z ε z + τ θz γ θz + τ rz γ rz + τ rθ γ rθ }dV
2 ∫vol
U =
와 같이 표시된다.
따라서 식(2-5a)는
U =
1
∫ {σ} {ε}T dV (1-5b)
2 vol
로 나타낼 수 있다. 여기서 직교 좌표계인 경우,
{σ}T = {σ x , σ y , σ z , τ yz , τ xz , τ xy },
{ε}T = {ε x , ε y , ε z , γ yz , γ xz , γ xy 이다.
}
식(2.3)을 식(2.5a)와 (2.5b)에 각각 대입하면,
{ε} [D ] {ε} dV = {ε} [D ]{ε} dV

1 1
∫ ∫
U =
T T T (1-6)
2 vol 2 vol
따라서 [D]T=[D]가 된다. 이것이 식(2.3)의 탄성계수 행렬[D]가 대칭 행렬임에 대한
수학적 증명이다. 여기에서는 열팽창 계수를 포함시키지 않았으나, 열 팽창에 의한 효과
를 고려하여도 [D]의 대칭성은 변하지 않는다.
식(1.2)를 식(1.6)의 오른편 식에 대입하면,
{q} [B ] [D ][B ] {q}dV

1
U = ∫
T T
(1-7)
2 vol
U = {q} ⋅ ∫vol [B ]T [D ][B ] dV ⋅ {q}

1 T (1-8)
2
(2) 외부의 힘에 의한 위치 에너지(W)
W = cons tan t − ∑ f i qi (1-9)

i
= cons tan t − { f }T {q} = cons tan t
여기서 {f}는 절점 하중 벡터(Nodal Load Vector)이다. 특히 하중이 절점에 가해지지

않고 요소 내의 어떤 위치에 작용할 때에는 이를 상당 절점 하중(Equivalent Nodal
Loads)으로 환산해서 더해 주어야 한다.
(3) 최소 위치 에너지 법칙 또는 변분 이론을 적용하면
∂
(U + W ) = 0 ; i = 1, 2, ⋅ ⋅ ⋅ (1-10)
∂qi
식 (2.8)과 (2.9)를 식(1.10)에 대입하여 정리하면 다음 식을 얻을 수 있다.
[k ] {q} = { f } (1-11)
여기서,
(1-12)
[k ] = ∫vol [B] [D][B] dV
T
이며, [D]가 대칭 행렬이므로 [k]도 대칭 행렬이 된다. 한 가지 기억할 것은 상당 절

점 하중이 존재할 경우 실제로 절점에 작용하는 힘(Actual Nodal Forces)은 식(1.11)의
{f}로부터 상당 절점 하중을 뺀 값이 되다는 점이다. 즉, 상당 절점 하중을 {f}0 라 하면
실제 절점 하중 {f}0 는
{ f }a = { f }− { f }0 (1-13)
따라서 변형에너지 U는 식(1.8)로부터,
U = {q} [k ][q ]
1 T
(1-14)
2
5) 필요하다면 국부 좌표계에서 표현된 각 요소에 대한 강성 행렬[k] 및 절점 행렬{f} 등
을 대표 좌표계(Global Coordinates)에 대한 것으로 변환한다.
{q} = [T ] {Q} (1-15)

여기서 [T]는 변환(Transformation) 행렬이다. 이 식을 변형 에너지 식(1.14)에 대입하

면,
U =
1
{Q}T [T ]T [k ][T ] {Q} = 1 {Q}T [k ] {Q}
2 2
즉, 요소의 대표 좌표 상의 강성 행렬 [K]는,
[K ] = [T ]T [K ][T ] (1-16)
식(2.11)의 양변의 앞에 [T]T를 곱한 뒤 식(1.15)에 대입하여, 식(1.16)을 이용하여 정

리하면,
[K ] {Q} = [F ] (1-17)
여기서
{F } = [T ]T { f } (1-18)
덧붙여 식(1.13)을 대표 좌표로 표시하면
{F }a = {F } − {F }0 (1-19)
여시서, {F}a는 실제로 절점에 작용하는 대표 좌표상의 힘(Actual Global Nodal Force
Vector)이고, {F}0는 대표 좌표 상의 상당 하중 벡터(Equivalent Global Nodal Load
Vector)이다. 또한 앞에서 언급한 여러 값, 이를 테면, {u}, {ε}, 등도 식(1.14)를 대입하
여 대표 좌표계에서의 값으로 쉽게 변환할 수 있다.
2단계에 포함된 식을 정리하면,
{u} = [N ] {q} (1-1)

{ε} = [B] {q} (1-2)
{σ} = [D] {ε} (1-3)
[k ] {q} = { f } ; [k ] = ∫vol [B] [D][B]
T
(1-11,12)
{q} = {T }{Q}
[K ] {Q} = {F }; [K ] = [T ]T [k ][T ], {F } = [T ]T { f } (1-16,17,18)
덧붙여, 앞으로 요소의 번호를 구분해야 하 필요가 있을 때는 첨자를 행렬 기호와 바깥

에 붙이도록 한다. 예를 들면, [K]i 는 요소 i에 대한 강성 행렬을 의미한다.
3단계, 요소에 대한 수식을 결합하여 전체 계에 대한 식을 얻고 경계조건(B/C)을 부여한다

(Assemble the Element Equation to Obtain the Global Equation and Impose the B/C).
3단계에서는 요소에 대한 수식을 결합하여 요소망(Element Mesh)으로 모형화된 전체 계

(Overall Systems)즉, 해석하고자 하는 전체 물체 또는 영역에 대한 평형 방정식을 유도하
는 단계이다.
이렇게 결합하여 얻어지는 전체 계에 대한 변위와 힘의 관계식(Displacement-Force
Relations은 각각의 요소에 대한 수식과 같은 형태로 표시되며 다만 모든 결점을 다 포함하
는 만큼 행렬식의 크기가 아주 커지게 된다.
[K ]g {Q} = [F ] (1-20)
여기서, [K]g는 대표 강성 행렬(Global or Total Stiffness Matrix), {Q}는 모든 절점 변위로

이루어진 벡터(The Vector of Global Nodal Displacement) 그리고, {F}는 모든 절점에서의
힘으로 이루어진 벡터(The Vector of Global Nodal forces)이다. 주의할 점은 강성 행렬을
결합할 때 연결되는 절점의 위치에 맞게 합산 해야 한다는 것이다.
컴퓨터상의 결합은 알맞은 알고리즘에 의하여 체계적으로 이루어지게 된다.
다음은 경계조건 또는 제한조건(Constraints)을 부여하게 되는데, 전체 계의 자유도

(Degree of Freedom)는 식(1.20)에서 경계조건을 제외하고 남는 벡터{Q}의 크기가 된다.
대표 강성 행렬 [K]g는 다음과 같은 성질을 가지고 있다.
1) [K]g는 대칭 행렬이다. [K]g의 대칭성은 식(1.12) 및 (1.16)으로부터 분명해진다. 그

런데 식(1.12)는 변형 에너지가 식(1.16)과 같은 2차형(Quadratic Form)으로 주어지
는 일반적인 경우에 대해서 유도된 식인데, 특수한 경우에는 변형 에너지 또는 그
에 준하는 범함수(Functional)가 식(1.16)과 같은 형태로 되지 않는 경우가 있다. 이
때의 [K]g는 비대칭이 되며, 이 경우 외에는 항상 대칭이 된다.
2) 경계조건(Boundary Conditions)이 부여되지 않는 상태에서의 [K]g는 특이 행렬

(Singular Matrix)이 된다. 이는 물리적으로 전체 계의 강체 운동(Rigid Body Motion)
때문인데 경계조건에 의하여 이러한 강체운동은 제거된다.
3) [K]g의 주 대각선(Main Diagonal) 값은 향상 양(Positive)이 된다. 그렇지 않을 경우

어떤 절점에서의 힘이 양일 때, 그 절점에서의 변위가 음이 되는, 현실에서는 있을
수 없는 결과가 초래되게 된다.
4) 변위법에 의해서 유도된 강성은 실제와 같거나 크게0 된다. 바꾸어 말하면, 변형 에
너지가 실제와 같거나 크게 된다. 이는 변위가, 가정한 변수 함수 식(1.1)로 표시할
수 있는 형상으로만 일어날 수 있으므로 이것이 하나의 제한 조건이 되어 강성이
강하게 되는 효과를 가져오기 때문이다.
따라서, 유한요소법으로 구한 변위는 하한계(Lower Bound)가 된다. 고유 진동수는
상한계(Upper Bound)가 된다.
그런데, 한 가지 기억할 것은 실지로 유한 요소법으로 값을 구해 보면 절점의 수가
많은 경우 즉 자유도가 큰 경우에는 꼭 그렇게 되지 만은 않는다는 사실이다. 이는
주로 큰 행렬의 처리 과정에서 발생하는 수치 상의 오차(Round-off Error) 때문이다.
4단계, 평형 방정식의 해를 구한다(Solver the Equilibrium Equations).
소거법(Elimination Method)이나 반복법(Iterative Method)등을 이용한 행렬식 처리를 통해

연립방장식(1.20)으로부터 절점의 변위를 구하고 필요할 경우 절점에서의 힘을 구한다.
5단계, 기타 필요한 값을 구한다

(Calculate Additional Parameters if Desired, Data Recovery).
필요할 경우 요소 내부의 어떤 점에서의 값을 계산한다. 즉 식(1.1)을 이용하여 요소의

내부의 변위를, 나아가 변형률과 변위 관계 (1.2) 및 응력과 변형률 관계 식(1.3) 등을 애용
하여 변형률, 응력 등의 필요한 값을 계산한다
6단계, 결과를 분석한다(Interpret the Results).
필요한 결과를 알아보기 쉽게 표 또는 그림으로 나타낸다. 이를 후처리(Post Processor)

라고 한다. 후처리를 통해서 결과를 검증하고 나아가 이를 설계 변경 등의 필요한 목적에
적절하게 반영한다.
3. 요소의 형상 함수
3.1 형상 함수(Shape Function) 특성
유한요소법(FEM)을 적용하는 과정에 있어서 아주 중요한 작업중의 하나가 형상 함수를

결정하는 작업이다. 유한요소 해석 결과의 정확성이나 수렴 속도 등은, 요소 수를 같이 했
을 때 거의 전적으로 형상 함수에 의하여 좌우된다.
대부분의 경우에 있어서는 다항식이 사용되고 있어 이를 간략하게 언급하고자 한다. 형상
함수의 결정은 때로는 매우 번거롭고 난해한 작업이 될 수 있다.
간단한 요소에 대한 형상 함수를 결정함에 있어서는 그다지 어려움이 없으나 자유도가 증
가 함에 따라서 복잡성은 기하급수적으로 증가하게 된다.
여기서는 형상 함수의 특성을 고찰해 함으로써 경우에 따라서는 이를 이용하여 번거로운

과정을 거치지 않고 보다 쉽게 형상함수를 결정할 수 있다고 본다.
한 예로 축방향 요소에 대하여 설명하면 다음과 같다.
x1 x2
q1 q2
1 2
그림[3.1] 축 방향 보 요소
(a) 전통적인 방법(Convectional Way)
두 절점사이의 변위는 다음과 같이 가정할 수 있다.
⎧C ⎫
u = C1 + C 2 x = [1, x ] ⎨ 1 ⎬ (3-1)
⎩C 2 ⎭
이 식에서 절점에서의 값을 대입하여 행렬식으로 정리하면,
⎧ q1 ⎫ ⎡1 x1 ⎤ ⎧ C1 ⎫
⎨ ⎬ = ⎢ ⎥⎨ ⎬
⎩q2 ⎭ ⎣1 x2 ⎦ ⎩C 2 ⎭
따라서 미지수 Ci의 해는,

⎧ C1 ⎫ 1 ⎡ x2 − x1 ⎤ ⎧ q1 ⎫
⎨ ⎬ = ⎢ − 1 1 ⎥ ⎨q ⎬
⎩C 2 ⎭ x2 − x1 ⎣ ⎦ ⎩ 2⎭
이 해를 식(3.1)에 대입하면
x − x1 ⎤ ⎧ q1 ⎫ ⎧q ⎫
u =
1
[1 x ] ⎡⎢ 2 ⎥⎨ ⎬ = [N1 N2 ]⎨ 1 ⎬
x2 − x1 ⎣ − 1 1 ⎦ ⎩q2 ⎭ ⎩q2 ⎭
여기서
x2 − x − x2 + x
N1 = , N2 = (3-2)
x 2 − x1 x 2 − x1
(b) 형상 함수의 특성을 이용하는 방법(using the characteristic of shape function)
형상함수 Ni 는 다음과 같이 가정할 수 있다.
N i = a1 + a 2 x (3-3)
( )
그런데 N i x j = δ ij 이므로, 먼저 N1에 대해서 이를 적용하면,
⎡1 x1 ⎤ ⎧ a1 ⎫ ⎧0 ⎫
⎢1 x ⎥ ⎨a ⎬ = ⎨1 ⎬
⎣ 2⎦ ⎩ 2⎭ ⎩ ⎭
따라서
⎧ a1 ⎫ 1 ⎡ x2 − x1 ⎤ ⎧1 ⎫ 1 ⎧ x2 ⎫
⎨ ⎬ = ⎢− 1 1 ⎥ 0 ⎨ ⎬ = ⎨ ⎬
⎩a 2 ⎭ x2 − x1 ⎣ ⎦⎩ ⎭ x2 − x1 ⎩− 1⎭
이 식을 식(3.3)에 대입하면,
x2 − x
N1 =
x 2 − x1
N2는 같은 방법으로 구하거나, 첨자 1과 2를 서로 바꾸면,
x1 − x
N2 =
x1 − x 2
위의 방법을 전통적인 방법과 비교해 보면, 이 경우는 수식의 전개 과정이 대체로 비슷

하게 되어 형상함수의 특성을 이용하는 방법이 특별히 유리한 점이 없다. 그러나, 형상함수
가 식(3.3)와 같이 1차원 선형함수인 경우에는 그림[3.2]와 같이 형상함수의 특성을 도형으
로 나타내면, 비례관계로부터 간단히 형상함수가 구해진다. 즉,
N1 1 N2 1
= , = (3-4)
x2 − x x2 − x1 x − x1 x2 − x1
따라서
x2 − x x − x1
N1 = , N2 =
x 2 − x1 x 2 − x1
N1 N2
1 1
x1 x2 x1 x2
x x
그림[3.2] 1차원 선형 형상 함수
3.2 내분함수(Interpolation Functions)
향상함수의 특성을 이용함에 있어서 대체로 두 가지 개념이 사용된다. 하나는 내분함수

개념이고 다른 하나는 자연 좌표 개념이다.
내분함수란 용어는 종종 형상함수(Shape Functions)와 같은 의미로도 쓰인다. 먼저 선형요
소(Linear Element)에 대해서 내분함수의 개념을 이용하는 방법을 보이고, 나아가 유한요소
법에서 자주 쓰이는 라그란지 다항식(Lagrange Polynomials)과 허마이트(Hermite
Polynomials)에 대해서 언급하고자 한다.
3.2.1 선형요소에 대한 내분 함수
식 (3.3) 또는 (3.4)를 검토해 보면 Ni는 (i는 1 또는 2) x로부터 다른 절점까지 길이의 요

소의 길이에 대한 비이다. 요소의 길이를 그림[3.3]에서와 같이 ″l″ 이라 두면,
N1 =
l1
=
1
(x2 − x )
l l
N1
l
= 2 =
1
(x − x1 )
l l
l2 l1 절점2
절점1
x1 x x2
그림[3.3] 1차원 내분 함수
이러한 접근 방법을 그림[3.4]의 삼각형 요소에 대해서 적용하면,
1 ⎡1
N1 =
A1
= ⎢ (x3 − x )( y3 + y ) + 1 (x − x2 )( y + y 2 ) − 1 (x3 − x2 )( y3 + y 2 )⎤⎥
A A ⎣2 2 2 ⎦
=
1
[(x2 y3 − x3 y 2 ) + ( y 2 − y3 )x + (x3 − x2 ) y ] (3-5)
2A
여기서 A는 절점을 꼭지점으로 하는 삼각형의 면적이고 다음의 식으로 계산할 수 있다.
1 x1 y1
1
A = 1 x2 y2
2
1 x3 y3
N2 및 N3는 1, 2, 3을 순환적으로 치환(Cyclically Permuting)함으로써 구해진다. 주의할
것은 1, 2, 3의 순서를 그림[3.4]에서와 같이 시계 반대 방향으로 지켜야만 한다.

만약, 순서를 바꾸게 되면 부호가 변하여 형상함수의 특성을 잃어 버리게 된다.
y 1(x1,y1)
A2
3(x3,y3)
A3
A1
2(x2,y2)
그림[3.4] 2차원 선형 요소
내분 함수 개념은 또한 그림[3.5]와 같은 사면체 요소에 대해서도 그대로 적용된다. 즉,

사면체 내부의 점을 (x,y,z)라 하면 N1은 이 점과 꼭지점 2, 3, 4로 이루어지는 사면체 부피
의 요소의 부피에 대한 비가 된다. 요소의 부피는,
1 x1 y1 z1
1 1 x2 y2 z2
V = (3-6)
6 1 x3 y3 z3
1 x4 y4 z4
이 식을 이용하여 N1을 구하면,
1 x y z
V1 1 1 x2 y2 z2
N1 = = (3-7)
V 6V 1 x 3 y3 z3
1 x4 y4 z4
=
1
(a1 + b1 x + c1 y + d 1 z )
6V
여기서
x2 y2 z2 1 y2 z2
a1 = x3 y3 z3 , −1 y3 z3
x4 y4 z4 1 y4 z4
1 x2 z2 1 x2 y2
c1 = 1 x3 z3 , − 1 x3 y3
1 x4 z4 1 x4 y4
N2, N3 및 N4도 같은 방법으로 구해진다.
z
(x, y, z)
y 2 3
4
그림[3.5] 3차원 선형 요소
3.2.2 라그란지 다항식(Lagrange Polynomial)
라그란지 다항식은 다음과 같이 정의된다.
Lk ( x ) =
n
x − xi (x − x1 )(x − x2 )Λ (x − xk −1 )(x − xk +1 )(x − xn )
∏x =
(xk − x1 )(xk − x2 )Λ (xk − xk −1 )(xk − xk +1 )(xk − xn )
(3-8)
i =1 k − xi
k =i
여기서 (n-1)은 다항식의 차수가 된다. xi를 절점 좌표라하면, x=xk일 때는 분모와 분자가
서로 같게 되므로 Lk는 1이 되고, x=xi 이고 i ≠ k이면 Lk는 0이되어, Lk는 형상함수의 특성을
지니고 있음을 알 수 있다. 따라서 Lk는 각 절점 당 자유도가 1인 요소의 형상함수가 된다.
즉, N k ( x ) = Lk ( x ) 이며, 이때 n은 절점의 수가 된다. 앞의 식(3.2)은 그러니까 n=2인 경
우이고, 예를 들어 n이 3인 경우에는 형상함수는 그림 3.6(b)와 같은 모양이 되며, 식(3.7)
을 이용하면 다음과 같은 형상 함수를 바로 얻을 수 있다.
N1 =
(x − x 2 )(x − x3 ) , N = (x − x1 )(x − x3 ) , N = (x − x1 )(x − x 2 ) (3-9)
(x1 − x 2 )(x1 − x3 ) 2 (x2 − x1 )(x2 − x3 ) 3 (x3 − x1 )(x3 − x 2 )
1
1 i
2 n
1
N1
x1 x2 x3
N2
1
x1 x2 x3
N3
1
x1 x2 x3
그림[3.6] Lagrange Interpolation Polynomials
3.2.3 허마이트 다항식(Hermite Polynomial)
라그란지 다항식이 각 절점 당의 자유도가 1인 경우에 n 절점 요소에 대한 형상함수를

결정하는데 사용하는데 비해서, 허마이트 당항식은 2장에서의 굽힘 보 요소와 같이 절점에
서의 미분의 계수(Order)가 증가할 경우에 사용된다. n계 허마이트 다항식은 변수가 x일때
Hn(x)와 같이 표시되며, 다항식의 차수는 (2n+1)이 된다. 여기서는 미분의 계수 m과 절점번
호 i를 나타내는 두개의 첨자를 덧붙여서 H min ( x ) 와 같이 표시하기로 한다. 그러면,
d l H min (x j )
l
= δ lmδ ij (3-10)
dx
이 허마이트 다항식 역시 형상함수의 특성을 가지고 있고 따라서 형상함수로서 사용할 수

있다.
예를 들어 굽힘 보 요소에 대한 형상함수를 허마이트 다항식으로 나타내보면 다음과 같다.
여기서 절점 좌표를 x1, x2로 표시하면,
N 1 = H 01
1
(ξ ) = (1 + 2ξ )(1 − ξ )2
N 2 = H 111 (ξ ) = (x2 − x1 )(ξ − 1)2 (3-11)
N3 = H 1
02 (ξ ) = (3 − 2ξ )ξ 2
N4 = H 1
12 (ξ ) = (x2 − x1 )(ξ − 1)ξ 2
여기서,
x − x1 ∂ 1 ∂
ξ = , 0 ≤ ξ ≤ 1, =
x2 − x1 ∂x x2 − x1 ∂ξ
3.3 자연 좌표계(Natural Coordinate System)
요소 내의 어떤 위치의 값을 0과 1사이의 무차원 수로 표시하는 국부 좌표계를 자연 좌

표계라 부른다. 자연 좌표계는 보통 해당 절점에서는 값이 1이고 다른 절점에서는 값이 0이
되도록 하며, 두 절점 사이의 값은 선형적으로 변한다.
자연 좌표계는 수학적으로는 임의의 차원의 요소에 대해서 적용할 수 있으나, 일반적으로
요소는 3차원을 넘지 않으므로 여기서도 3차원을 넘는 초공간(Hyperspace)에 대해서는 취
급하지 않기로 한다. 자연 좌표계를 사용하게 되면 형상함수를 쉽게 결정할 수 있을 뿐만
아니라, 적분값을 간단한 수식으로 나타낼 수 있는 이점이 있다.
3.3.1 선형 자연 좌표
그림[3.7]에 보이는 1차원 요소에 대해서, N1 및 N2를 자연좌표라 하면 요소 내의 어떤

점의 좌표 x는
x = N 1 x1 + N 2 x2 (3-12)
로 나타낼 수 있다.
여기서 x가 x1일 때는 N1 = 1, N2 = 0가 되고, x가 x2일 때는 N1 = 0, N2 = 1이 되므로 형상함
수의 특성을 갖고 있음을 알 수 있다. 나아가 N1과 N2는 그림에서와 같이 각각 0에서 1까지
변하는데, 둘을 합하면 항상 1이 된다.
N1 + N 2 = 1 (3-13)
1 2
x1 x x2
N1
x1 x x2
1 0.75 0.5 0.25 0
N2
x1 x x2
0 0.25 0.5 0.75 1
그림[3.7] 1차원 자연 좌표
식(3-12)와 (3-13)을 행렬식으로 정리하면,
⎡1 1 ⎤ ⎧ N1 ⎫ ⎧1 ⎫
⎢x ⎥ ⎨ ⎬ = ⎨ ⎬
⎣ x2 ⎦ ⎩ N 2 ⎭ ⎩x⎭
따라서 {N}의 해는,
⎧ N1 ⎫ 1 ⎡ x2 − 1⎤ ⎧1 ⎫ 1 ⎧ x2 − x ⎫
⎨ ⎬ = ⎢− x ⎥ ⎨ ⎬ = ⎨ ⎬ (3-14)
⎩N 2 ⎭ x2 − x1 ⎣ 1 1 ⎦⎩ x⎭ x2 − x1 ⎩− x1 + x ⎭
이 자연 좌표개념을 이용하여 구한 {N}은 식(3.2)의 1차원 선형요소에 대한 형상함수와

정확하게 일치함을 알 수 있다. 그런데 N1 및 N2는 길이의 비이므로 종종 길이 좌표
(Length Coordinate)라도 불리기도 한다. 덧붙여, 여기서는 내분함수 개념과 자연좌표 개념
을 분리해서 설영을 했으나, 이는 어디까지나 형상 함수를 구하는데 있어서의 접근 방법의
차이를 설명하기 위한 것일 뿐 함수자체(같은함수)에 대한 분류를 뜻하는 것은 아니다.
1차원 자연좌표의 길이에 대한 적분은 다음과 같은 간편한 수식으로 주어진다.
(x2 − x1 ) α ! β !
x2
∫x1
N1α N 2β dx =
(α + β + 1)!
(3-15)
여기서 α와 β는 자연수이며 α!은 α의 계승(Factorial)이다.

x에 대한 미분은 연쇄 법칙(Chain Rule)으로부터
d 2
∂N i ∂ ∂N 1 ∂ ∂N 2 ∂ 1 ⎛ ∂ ∂ ⎞
dx
= ∑
i =1 ∂ x ∂N i
=
∂ x ∂N 1
+
∂x ∂N 2
= ⎜⎜ − ⎟ (3-16)
x2 − x1 ⎝ ∂N 2 ∂N 1 ⎟⎠
1 3
0.75 0
0.5
0 0.25
1 0.25
1
0.75 0.5 N2
y
0.5
0.75
x 0.25
1 그림[3.8] 2차원 자연 좌표
0 2
1차원 선형요소에 대한 위의 접근 방법을 그림[3.8]의 삼각형 요소에 적용하면,
⎡1 1 1 ⎤ ⎧ N1 ⎫ ⎧1 ⎫
⎢x ⎥ ⎪ ⎪ ⎪ ⎪
⎢ 1 x2 x3 ⎥ ⎨ N 2 ⎬ = ⎨ x ⎬ (3-17)
⎣⎢ y1 y2 y3 ⎦⎥ ⎪⎩ N 3 ⎪⎭ ⎪ y⎪
⎩ ⎭
따라서 {N}의 해는,
⎧ N1 ⎫ ⎡ a1 b1 c ⎤ ⎧1 ⎫
⎪ ⎪ 1⎢ ⎪ ⎪
⎨N2 ⎬ = ⎢a2 b2 c2 ⎥⎥ ⎨ x ⎬ (3-18a)
⎪N ⎪ Δ
⎩ 3⎭ ⎢⎣ a3 b3 c3 ⎥⎦ ⎪⎩ y ⎪⎭
여기서
ai = x j y k − xk y j , bi = y j − yk , ci = x k − x j
I, j, k는 1, 2, 3의 순환 순열(Cyclic Permutation)
1 1 1
Δ = x1 x2 x3 = 2A = (a i + a j + ak )
y1 y2 y3
식(3.18a)를 간편히 다시 쓰면,
Ni =
1
(ai + bi x + ci y ) (3-18b)
2A
삼각형 요소에 대한 자연좌표는 넓이와 관계 되므로 종종 넓이 좌표(Area Coordinate)라

고 불리기도 한다. 넓이 좌표의 삼각형 요소의 넓이에 대한 적분은 다음과 같은 간편한 수
식으로 얻어진다.
α ! β !γ !
∫A
N1α N 2β N 3γ = 2A
(α + β + γ + 2 )! (3-19)
여기서 α, β, γ는 자연수이며, α!은 α의 계승(Factorial)이다.

x 및 y에 대한 미분은,
∂ 1 ⎛ ∂ ∂ ∂ ⎞
= ⎜⎜ b1 + b2 + b3 ⎟
∂x 2 A ⎝ ∂N 1 ∂N 2 ∂N 3 ⎟⎠
(3-20)
∂ 1 ⎛ ∂ ∂ ∂ ⎞
= ⎜⎜ c1 + c2 + c3 ⎟
∂x 2 A ⎝ ∂N 1 ∂N 2 ∂N 3 ⎟⎠
그림[3.5]의 사면체 요소에 대해서 자연좌표 개념을 도입하면,
⎡1 1 1 1 ⎤ ⎧ N1 ⎫ ⎧1 ⎫
⎢x ⎥ ⎪ ⎪
x4 ⎥ ⎪ N 2 ⎪ ⎪x⎪
⎢ 1 x2 x3 ⎪ ⎪
⎨ ⎬ = ⎨ ⎬ (3-21)
⎢ y1 y2 y3 y4 ⎥ ⎪ N 3 ⎪ ⎪ y⎪
⎢ ⎥⎪ ⎪ ⎪⎩ z ⎪⎭
⎣ z1 z2 z3 z4 ⎦ ⎩ N 4 ⎭
따라서
⎧ N1 ⎫ ⎡ a1 b1 c1 d1 ⎤ ⎧ 1 ⎫
⎪N ⎪ ⎢ d 2 ⎥⎥ ⎪⎪ x ⎪⎪
⎪ 2⎪ 1 ⎢a2 b2 c2
⎨ ⎬ = ⎨ ⎬ (3-22)
⎪N3 ⎪ Δ ⎢ a3 b3 c3 d3 ⎥⎪ y⎪
⎪⎩ N 4 ⎪⎭ ⎢ ⎥
⎣a4 b4 c4 d 4 ⎦ ⎪⎩ z ⎪⎭
여기서
x2 x3 x4 1 1 1
a1 = y2 y3 y4 b1 = − y2 y3 y4
z2 z3 z4 z2 z3 z4
(3-23)
1 1 1 1 1 1
c1 = x2 x3 x4 d1 = − x2 x3 x4
z2 z3 z4 y2 y3 y4
1 1 1 1
x1 x2 x3 x4
Δ = = 6V = (a1 + a2 + a3 + a 4 ) (3-24)
y1 y2 y3 y4
z1 z2 z3 z4
그리고 기타 상수 a2, b2,… c4, d4 는 식(3.23)와 비슷하게 주어진다. 그런데 한 가지 주의할

것은 첨자 1, 2, 3, 4는 순환 순열이 아니라는 점이다. 식(3.22)의 상수는 식(3.24)의 디터미
넌트(Determinant)의 코펙터(Cofactor)이므로, 순환 순열로 취급했을 때는 부호가 바뀔 경
우가 발생한다. 식(3.23)은 좌표축이 오른 손 법칙을 따를 때, 절점 1, 2, 3이 절점 4에서
볼 때 시계 반대 방향인 경우에 대한 것이다.
사면체에 대한 자연 좌표는, 종종 부피 좌표(Volume Coodinate)라고 불린다. 부피 좌표에
대한 적분 및 미분은 다음 식으로 주어진다.
α ! β !γ !δ !
∫N
α
N 2β N 3λ N 4δ dV = 6V (3-25)
v
1
(α + β + λ + δ + 3)!
∂ 1 4
∂
∂x
=
6V
∑b
i =1
i
∂N i
∂ 1 4
∂
∂y
=
6V
∑c
i =1
i
∂N i
(3-26)
∂ 1 4
∂
∂z
=
6V
∑d
i =1
i
∂N i
3.3.2 중(重) 자연 좌표
그림[3.9a]와 같은 임의의 사각형 요소에 대해서는 그림[3.9b]와 같은 정사각형으로 좌

표변환하여 취급하게 된다. 정사각형에 요소에 대한 자연 좌표는 각각의 축 방향에 대한 1
차원 자연 좌표를 곱하여 얻을 수 있다. 그림에서와 같이 정사각형의 중심을 원점으로 택하
여 변의 좌표를 ξ = ±1, η = ±1 로 하면, 절점 1과 2사이의 자연 좌표Li(ξ)는 식(3.14)로부터
Li (ξ ) = (1 − ξ ), L2 (ξ ) = (1 + ξ )
1 1
2 2
η 방향에 대해서도 같은 방법으로 주어진다. 따라서 정사각형 요소에 대한 형상함수는
N 1 = L1 (ξ )L1 (η ) = (1 − ξ )(1 − η )
1
4
N 2 = L2 (ξ )L1 (η ) (1 + ξ )(1 − η )
1
=
4 (3-27a)
N3 =
1
(1 + ξ )(1 + η )
4
N1 =
1
(1 − ξ )(1 + η )
4
이를 다시 쓰면,
Ni =
1
(1 − ξ iξ )(1 + η iη ) (3-27b)
4
여기서 ξi , ηi 는 절점의 좌표를 뜻한다. 이러한 함수는 2중 선형(Bilinear) 함수라 한다. 그
런데 여기서 한 가지 주목할 것은 국부 좌표 ξ , η 와 대표 좌표 x, y 사이의 좌표 변환관계
식이 지역 좌표에 대한 형상함수와 같은 함수로 주어진다는 점이다. 즉,
x = N 1 x1 + N 2 x 2 + N 3 x3 + N 4 x4
(3-28)
y = N 1 y1 + N 2 y 2 + N 3 y3 + N 4 y 4
이와 같은 좌표변환을 위한 함수와 형상함수가 같을 때 이러한 요소를 등매개

(Isoparametric) 요소라 한다.
η
4(x4,y4)
4(-1,1) 3(1,1)
3(x3,y3)
ξ
1(x1,y1)
2(x2,y2) 1(-1,-1) 2(1,-1)
(a) 대표좌표 (b) 자연좌표
그림[3.9] 사각형 요소에 대한 자연좌표

이러한 자연좌표 개념은 또한 그림[3.10a]와 같은 육면제 요소에 대해서도 적용된다.
그림[3.10b]와 같이 지역 좌표 ξ, η로 좌표 변환하면,
8 8 8
x = ∑N x,
i =1
i i y = ∑N y,
i =1
i i z = ∑N z
i =1
i i
여기서
Ni =
1
(1 + ξ iξ )(1 + η iη )(1 + ς i ς )
8
Ni는 또한 그림[3.10b]의 정육면제 요소에 대한 형상함수이기도 하다. Ni는 삼중 선형 자
연 좌표라기도 한다.
x
(a) 대표좌표
ζ
(-1,-1,1) (-1,1,1)
(1,1,1)
(1,-1,1)
η
(-1,-1,-1)
(-1,1,-1)
(1,-1,-1)
ξ (1,1,-1)
(b) 자연좌표
그림[3.10] 육면체 요소에 대한 자연좌표

4. 요소 형상 및 특징
일반적으로 각각의 다양한 구조 특징 및 형상에 적용될 수 있게 많은 요소의 형식

(Element Type, NKTP)이 제공되고 있으나 크게 구조물을 2차원 및 3차원 구조로 보고 그
구조에 적용될 수 있는 요소의 선택이 요구되어진다.
다음은 다양한 요소의 형식으로 일반적으로 주로 사용하고 있는 요소는 선형요소(Linear
Element, 1차 요소), 2차형 요소(Quadratic, 2차 요소) 및 3차형 요소(Cubic, 3차 요소)이다.
특히 특이한 형상의 주도가 아니면 선형요소를 주로 사용한다고 보아도 무리는 아니다.
4.1 1차원요소(One-Dimensional Elements)
선분요소(Line Element)로 철골구조물, 철탑, 강교, 트러스 구조물 등과 같이 단면에 비해

길이가 긴 구조에 적용된다고 보면 된다. 요소형상의 특징은 선분으로 이루어져 있기 때문
에 어떤 면적 및 입체의 형상을 띄우지 않는 요소라 보면 된다. 단 내부적으로 그 선분에
대한 각각의 특성(단면적, 관성모우멘트, 단면계수 등)이 입력이 된다
아래의 2차형 요소는 선형요소(1차형 요소)의 중간에 새로운 절점이 1개가 존재하는 요소
이며, 3차형 요소는 선형요소를 3등분하는 부분에 새로운 절점 2개가 존재하는 요소이다.
예) Beam, Truss, Bar, Spring, Cable,…
1차형(Linear)요소 2차형(Quadratic)요소 3차형(Cubic)요소
그림[4.1] 1차원 요소의 요소 형상
Line Element
Line Element
3D Shell Element
(b) Steel Box with Arch Cable Bridge
(a) 강교(Beam및 Truss로 이루어진 교량) (c) 주차타워(Beam및 Truss로 이루어진 구조물)
그림[4.2] 1차원 요소의 구조물

4.2 2차원요소(Two-Dimensional Elements)
일정한 면을 형성하는 평면(판)요소로 실 구조물은 두께가 존재하는 입체형상의 Solid 형

태이지만 구조물의 특징상 Shell요소로 전환할 수 있는 구조물에 적용이 되며 내부적으로
각각 구조물의 특징 부위별로 두께가 입력이 된다.
4.2.1 삼각형 요소(Triangular Elements)
3 3 3
6
5
5 41
7 4
1 2 1 2 1 2
6 9
8
선형(Linear)요소 2차형(Quadratic)요소 3차형(Cubic)요소
3 3
8 7
7 6
9 6
8 5
10 5
9 4
11 4
1 2 2
1
10 11 12 12 13 14 15
4차형(Quartic)요소 5차형(Quintic)요소
그림[4.3] 삼각형 2차원 요소의 형상
4.2.2 사각형 요소(Quadrilateral Elements)
4차형(Quartic)요소 5차형(Quintic)요소
그림[4.4] 사각형 2차원 요소의 형상
[Purge Tank]
[Gas Reboiler]
[Skirt Support Junction]
[Gas Cleaning Plant]
Stress Analysis [Cyclone Reactor]

of a Ventury
[Cyclone Reactor]
[Artificial Teeth Implant]

그림[4.5] 구조 부분의 2차형 요소의 2D/3D구조물
Water
Convection
Convection
Air
Radiation
Convection
1200°C
Steel Plate
Refractory
Gas Velocity(1200°C) Cooling Tube Lines Temperature(°C)
Water
Air
Fluid Flow in Chuteway with Aerator
Over a Weir
Over a Wall
Water Gate & Spillway Flow
Hydraulic Jump Flow
그림[4.6] 비구조 부분의 2차형 요소의 2D/3D 열/유체 유동

4.3 3차원요소(Three-Dimensional Elements)
일정한 여러 개의 면으로 이루어진 하나의 체적을 형성하는 입체(Solid)요소로 실 구조물

의 모든 특징 및 형상을 그대로 반영하는 요소를 말한다. 다음은 각각 4차 및 5차형의 요소
형상은 생략이 되었다.
4.3.1 삼각뿔요소(Tetrahedronal Elements)
그림[4.7] 삼각뿔 3차원 요소의 형상
4.3.2 육면체(Hexahedronal Elements)
그림[4.8] 육면체 3차원 요소의 형상
4.3.3 삼각기둥(Prismatical Elements)
그림[4.9] 삼각기둥 3차원 요소의 형상

Planetary Gear
Automesh using surface of Pro/E
G-Gas Tank
Heavy Equipment
그림[4.10] 구조 부분의 3차형 요소의 3D구조물

Heavy Equipment Knuckle Parts
Automobile Aluminum Wheel
Oil Screw Conveyer
Regulator Stopper Wheel
Automobile ET Condenser
Milling Machine
Train Rail Road Base Structures
그림[4.11] 구조 부분의 3차형 요소의 3D구조물

Fluid Flow of Cooling Guide

Fluid Flow of Dust Absorber Plant
Fluid Flow of Automobile Evaporator
Natural Cooling Problem of Rail(900°C)
External Fluid Flow of G7 Car
그림[4.12] 비구조 부분의 2/3차형 요소의 2D/3D 열/유체 유동

5. 요소의 생성 기법 및 결과 분석 방법
Theoretical Course 1 5. 요소의 생성기법 및 결과분석
5.1 구조해석의 단계별 순서
다음의 단계별 해석적인 접근 순서는 구조해석 프로그램마다 상이하지만 거의 유사한 접

근 방법으로 이루어지기 때문에 이를 숙지하는 것이 필요하리라 본다.
단, 처음 구조해석 프로그램(CAE Tool)을 접하는 학생이나 엔지니어는 다소 이해가 되지
않는 부분이 있지만 차후에 예제를 통한 해석적인 접근을 해 보고 이 부분을 다시 숙지한다
면 쉽게 이해하리라 본다. (2), (3), (4)의 순서는 무관하다(편리한 순서 선택).
다음의 해석적인 접근은 구조 모델(비구조 제외)에 대한 기준으로 설명이 되었다.
(1) 모델의 차원(Model Dimension) 결정
구조물의 형상 및 해석의 주안점에 따라 유한요소모델 접근을 2차원으로 할 것인지 아니

면 3차원의 모델로 접근 할 것인지 결정한다.
예)
구조형상이 축에 대해서 아니면 어떤 좌표를 기준으로 대칭이면 그 기준면을 중심으로 1/2 또는 1/4
유한요소모델만 구축한다.
(2) 경계조건 부여(Boundary Condition)
전체 구조 모델이든 아니면 대칭형 모델이든 구조물에 경계조건을 부여한다.

예)
모든 구조물에는 그 구조물을 고정하고 있는 부분(구속부분)이 꼭 있기 때문에 이 부분이 유한요소모
델로 구현이 되었다면 이 부분을 그 부분의 특성에 따라 경계조건을 부여한다. 단, 전체 모델이 아닌
대칭형의 모델로 구현이 되었다면 그 대칭면에 대칭모델에 적합한 경계조건을 부여한다. 경계조건의
부여는 Displacement 6개의 자유도(D.O.F) 경계조건이 존재하는데(Ux, Uy, Uz, Rx, Ry, Rz, 모델의 형
상에 따라 자유도에 따라 상이) 0.0 입력 시에는 각 방향 별 완전 구속이며, 0.0 이 아닌 숫자가 입력
이 되면 그 부분이 강제 변위 하중(Forced Displacement Loading)으로 입력이 된다는 점에 주의해야
한다.
(3) 하중조건 부여(Loading Conditions)
경계조건과 같이 아주 중요한 부분으로 하중의 종류는 다음과 같다.

예)
-. 자중(Shelf Weight) : 구조물의 무게(밀도, 중력가속도, 중력 방향 입력)
-. 집중하중(Concentrated Nodal Force) : 일정한 지점에 집중적인 하중 작용
축방향 하중(Fx, Fy, Fz), 모멘트(Mx, My, Mz)
-. 온도하중(Nodal Temperature Loading) : 열적인 조건(열팽창계수 입력)
-. 균일분포하중(Distributed Pressure Loading) : 일정한 영역에 균일하중 작용
Uniform Pressure(등 분포하중), Non-Uniform Pressure(불 균일 분포하중)
Hydrostatic Pressure(수 두압), Wind Load(바람 하중), Seismic(지진,가진하중)…
(4) 기계적 물성치 부여(Material Properties)
일반적인 Steel의 경우에는 재질을 균질체(Homogeneous)로 본다면 3축방향(X,Y,Z)에 대

한 탄성계수(Elastic Modulus, Young’s Modulues) 및 포이송비(Poisson’s Ratio)가 같기 때
문에 등방성(Isotropic)이고 외력에 의한 정적해석(Linear Static)을 하기 위해서는 단순히 탄
성계수(Ex) 및 포이송비(νxy) 두개의 특성만 입력이 되면 된다.
예)
등방성 재질(Isotropic Material) : 탄성계수(Ex), 포이송비(νxy)
정방성 재질(Orthotropic Material) : 탄성계수(Ex, Ey, Ez), 포이송비(νxy, νyz, νxz)
이방성 재질(Anisotropic Material) : Laminated / Solid / Sandwich Composite
-. 세로탄성계수(Ex, Ey, Ez), 가로탄성계수(Gx, Gy, Gz), 포이송비(νxy, νyz, νxz)

-. Tensile/Compressive Failure Stress(x/y/z-direction)
-. Shear Failure Stress(Inplane/out-of-plane)
-. Coupling Coefficient(Tsai-Wu Theory)
(5) 모델특성 입력(Model Properties)
일반적으로 처음 접하는 학생이나 엔지니어는 어떤 모델에 어떤 특성을 입력해야 하는지

그리 쉽지는 않지만 각각 요소형식이 같고 있는 특성을 이해하고 생각하면서 입력하면 쉬우
리라 생각된다.
예)
Beam, Truss, Bar, Spring, Cable,… : 단면적, 단면계수, 관성모멘트,…
Solid Element(NKTP 4): 입력없음(실형상 입력)
Shell Element(NKYP 20) : 두께입력(두께모델 없이 종이와 같은 얇은 경우)
(6) 입력화일 작성
NISA Software를 이용하여 상기의 과정 (1) ~ (5)가 이상 없이 수행이 이루어 졌다면 다

음과 같은 다양한 포멧으로 출력할 수 있다. 이는 곳 외부의 프로그램에 제공되기도 한다.
예)
NISA Solver를 수행하기 위한 입력 파일 : ***.nis(절점,요소,재질,하중,경계 등)
현재 까지의 작업 내용 수록 : ***.ses(모든 단계 별 명령 수행 저장)
죄종 작업까지의 내용 저장 : ***..dbs(Database 파일, ***.nis + 기하학적 형상정보)
….
(7) Solver 수행
상기의 ″(6)″에서 생성된 ***.nis file을 입력하고 실행한다.

예)
-. ***26.dat, ***27.dat 파일 생성 : 그래픽으로 모든 응력, 변형, 모멘트, 전단력 등을 확인할 수 있
는 정보 파일 생성(Binary File)
-. ***.out 파일 생성 : 모든 결과 정보를 Text로 제공
(8) 결과분석(Result Evaluation)
″(7)″의 수행에서 생성된 ***26.dat, ***27.dat를 DisplayIII/IV Post/Processor에서

Loading 한 후 결과 분석 절차에 따라 분석한다.
예)
File ⇒ Read ⇒ Post Files ⇒ ***.26.dat ⇒
Post ⇒ Contours ⇒ Pick Result ⇒ Displacement(X,Y,Z, Resultant), Stress(Sxx, Syy, Szz, Sxy, Syz,
Sxz, Von-Mises, Principal 1/2/3, Maximum Shear Stress,…)
Post ⇒ Deformation ⇒ Pick Results,…
(9) 결과그림 출력(Results Report/Output)
결과 그림의 출력 및 각 종 에니메이션으로 보고서 작성이나 프래젠테이션

(Presentations)용으로 편집할 수 있다.
예)
그림 파일(Figures) : BMP, JPG, TIF,…
에니메이션화일 : AVI,…
5.2 유한요소 모델 기법
다음의 그림[5.1]에서 보여지듯이 해석으로 바로 들어가지 전에 가능한 역학적인 지식을

동원하여 수 계산적인 공학 계산을 해 본다든지 하는 전 단계 절차가 필요하다. 공학해석
결과의 모든 것을 바로 설계에 반영하여 제작 및 생산에 적용하기 전에 많은 부분에서 검토
가 이루어져야 하기 때문에 다음의 사항의 단계를 고려하여 설계에 반영 하도록 하는 것이
필요하리라 본다.
(1) 역학적인 개념을 도입하여 수 계산적인 방법이 가능한 모델이면 분석한다.

(2) 완벽한 해석 결과를 도출하였다 하여도 그 모델이 생산성과 제작이 가능한 모델인지
분석한다. 그리고 아주 튼튼하고 문제가 없는 제품의 모델를 구현 했다 하더라고 원가
분석 측면에서 경제성 Cost 면에서 분석을 해야 한다.
(3) 상기의 (1)(2)를 고려할 수 있는 엔지니어가 동시에 설계 및 해석 수행이 가능하면 이
상적인 모델 및 단 시간 내에 제품을 구현할 수 있지만 현재 국내의 현장 여건은 설계
및 해석의 범주가 분리되어 운영되는 경우가 대부분이기 때문에 많은 문제를 내포하고
있다고 보여진다.
그래서 구조해석 및 설계가 분리된 현장에서는 해석을 전담하는 연구소 와 설계부서 간
에 감정적인 대립이 발생하곤 한다. 해석적인 차원에서는 이상적인 모델이 설계 및 디
자인을 전담하는 부분에서는 수용할 수 없는 부분이 있기 때문이다.
Answer(%)
100%
50%
More Detailed
F.E. Analysis
Concept
F.E. Analysis
Analysis Effort(%)
Hand Calculations 100%
σ = (My I ) + (P A)
Information(%)
Need for
Analysis Results
100%
Availability
of Good FE Results
Time
Start of Design Start of Poduction
그림[5.1] 단계 별 해석 접근 및 제품 설계
5.3 유한요소 모델의 주의 사항
앞서 설명한 여러 가지의 특이성을 고려하여 유한요소모델로 구조물을 구현 했다 해도

100% 실 구조물을 구현하기란 아주 매우 힘든 일이다. 아주 단순한 모델은 100% 구현이
가능하지만 굴곡, 유선형적인 원형이 존재하는 부분에 대해서는 실 구조물에 비해 100%를
초과 할 수도 있고 100% 이하의 모델을 구현할 수 있기 때문에 다음을 고려하여 유한요소
모델로 구현된 완성된 모델을 다시 확인하는 과정이 요구되어 진다.
아니면 다음의 사항을 숙지하여 유한요소모델로 구현 하기 전에 미리 숙지하여 반영하면

서 모델을 구현하면 아주 이상적인 방법이라 생각이 된다.
주) 여기서 100%의 구현이라는 의미는 모델의 면적, 체적이라 생각하면 간단하다. 결국은 이러한
면적 및 체적에 밀도(Density)가 부여되어 무게가 되기 때문이다. 따라서 다음의 여러 가지 모델을 참
고하여 유한요소모델을 구현할 때문 100%를 기준으로 본 다면 110% ~ 90%의 모델에 준하는 접근
모델 기법이 필요하다. 이는 곧 어느 정도 정확한 모델(FE Model)를 구현했고 결과에 대한 신뢰성과
깊은 관련이 있는 아주 중용한 인자임을 인식해야 한다.
(1) 실 구조물 보다 추가되는 경우(100% 이상) 모델
구조물의 내부 Hole(φd)과 같은 부분에 대한 유한요소모델의 구현은 아래 그림[5.2]와 같

이 실 구조물에는 존재하지 않는 허공의 일 부분을 구조의 유한요소모델로 포함하여 분석
하는데 이는 절점(Node)과 절점(Node)의 사이는 직선으로 연결되는 요소의 특성상 피할 수
없는 모델의 구현 방법이다.
따라서 실 형상의 면적이나 체적을 기준으로 10%를 초과하지 않는 부분에서 유한요소모
델(Finite Element Model)를 구현해야 한다.
E6 E5
E7 Internal Hole
E4
φd
E8 E3 E3(요소 ″3″)
E1 E2
면적이 추가된 부분 모델
그림[5.2] 이하)
(2) 실 구조물 보다 삭감된 경우(100% 실 형상보다
모델 추가된 모델
구조물의 외부 원형(φd)과 같은 부분에 대한 유한요소모델의 구현은 아래 그림[4.3]과 같

이 실 구조물의 내부 일 부분을 구조의 유한요소모델로 포함하여 분석하는데 이는 절점
(Node)과 절점(Node)의 사이는 직선으로 연결되는 요소의 특성상 피할 수 없는 모델의 구
현 방법이다.
따라서 실 형상의 면적이나 체적을 기준으로 10% 이상 삭감 되지 않는 부분에서 유한요
소모델(Finite Element Model)를 구현해야 한다.
면적이 삭감된 부분 모델
E6 E5 φd
E7 E12 E11 E4
E9 E10 E3
E8 E3(요소 ″3″)
E1 E2 Shaft Outside Diameter(φd)
그림[5.3] 실 형상보다 삭감된 모델
(3) 원형 부분의 모델은 30°이내로 모델
그림[5.4]와 같은 부분의 모델은 기하학적 모델[Line(선), Patch(면), Hyperpatch(입체)]

부터 30° 이내로 모델를 구현하고 요소를 생성할 때에도 최소 30° 이내에서 요소가 생성이
되도록 주의한다.
E6 E5
45°>30°
E7 E4 22.5° < 30°
E8 E3
E1 E2
(a) Bad Model (b) Good Model
그림[5.4] 원형 부분의 요소 생성 모델
(4) 원형 축과 같은 부분의 모델 기법
그림[5.5]와 같은 부분의 모델은 생성은 삼각형 요소(Triangular Element) 및 삼각뿔 요소

(Tetrahedron Element)를 피하고 사각형(Quadrilateral) 및 직육면체 요소(Hexahedron
Element)를 생성하도록 주의한다. 응력 분석 부분(Evaluation Point)이 아니면 삼각형으로
생성해도 무관하나 되도록 이면 사각형으로 요소를 생성한다.
이러한 부분의 모델에 있어서 삼각형 모델의 구현은 쉽지만 사각형 요소의 구현에는 어려
움이 있지만 노력이 필요하다. 차 후에 설명이 되겠지만 일반적으로 삼각형이 사각형 보다
실제의 응력에 비유하면 10% Over Stress가 발생하기 때문에 주의 해야 한다. 특히 응력이
집중이 되거나 평가 부분에서의 삼각형 요소 모델은 생성하지 않도록 주의해야 한다.
삼각형 요소 사각형 요소
그림[5.5] 원형 축과 같은 부분의 요소 생성 모델
(5) 삼각형 요소보다 사각형 요소 사용
그림[5.6]와 같은 구조 모델에서는 응력의 기울기가 심하게 일어나는 부분이기 때문에 되

도록 이면 삼각형의 요소보다 사각형의 요소 형상을 사용해야 한다. 가능한 한 앞서 언급하
였듯이 삼각형을 배제하고 사각형의 요소 생성에 집중할 필요가 있다.
응력집중 부위 Fine Mesh
삼각형 요소 사각형 요소
그림[5.6] 응력 집중 부위 모델
(6) 자유도가 상이 구조의 결합
일반적으로 자유도(Degree of Freedom)가 상이한 구조 와 구조의 결합문제는 많은 이론

적인 접근 경험에 의한 모델링 접근 방법이 제시되고 있지만 현장 적용차원에서의 제시 방
안은 다음과 같다.
z MPC(Multiple Point Constraint) 사용
예를 들어 입체요소(Solid Element)와 판요소(Shell Element)의 결합에는 서로 접속이 되

는 부분에서는 자유도가 상이(Solid : 3 DOF, Shell : 6 DOF)하기 때문에 그 접속 부분에 하
나의 중간 매개 역할을 하는 구조가 삽입이 되어 서로 상이한 자유도의 결합이 요구
(Singular Matrix)되어지는 부분에 적용할 수 있다.
상이한 자유도 접속부분

6 DOF(Ux, Uy, Uz, Rx, Ry, Rz)
3 DOF(Ux, Uy, Uz)
되도록 거리 짧게
일반절점
Rigid Link(무한 강서의 강체)
Real Shell Thickness
주절점(Master Node)
종속절점(Slave Node) 그림[5.7] MPC 사용 모델
MPC에서 주절점(Master Node)과 종속절점(Slave Node)은 각각 6개의 자유도(DOF)를 소

유하고 있고, 종속절점은 자유도를 제어할 수 있다. 예를 들어 상기와 같은 부분에서 종속
절점 6개의 DOF 중 영향을 받는 부분의 3개의 자유도(Ux, Uy, Uz)이기 때문에 ″Slave DOF
Code″ 항에 ″1/2/3″을 입력하면 된다. 이는 곧 나머지 Rx, Ry,Rz는 주절점에 영향을 받지
않는다는 의미 이며 MPC Equation항에서 이항이 무시되어 같은 자유도로 인식이 된다.
종속관계⇒Dependent : DOF(Ux,Uy,Uz), 비 종속관계⇒Independent : DOF(Rx, Ry, Rz)
z 쉘요소(Shell Element)의 입체요소 속에 삽입
아래 그림과 같이 쉘요소를 입체의 요소 속에 통과 시켜 새로운 쉘요소를 만든다. 단 입

체 속에 포함되는 쉘요소의 두께는 Thickness Property에 입력한다.
그리고 입체요소 속에 포함이 되는 쉘요소 두께 부분에 반듯이 입체요소의 경계면이 존재
해야 한다.
상이한 자유도 접속부분

6 DOF(Ux, Uy, Uz, Rx, Ry, Rz)
3 DOF(Ux, Uy, Uz)
입체요소의 경계면
Shell Element 삽입
Real Shell Thickness
그림[5.8] 입체 속에 삽입 쉘 모델
(7) 쉘 요소(Shell Element)의 법선 방향 확인
일반적인 CAE 소프트웨어에서 Shell요소를 생성할 때 주의를 하면서 기하학적 모델링 및

요소의 생성 순서를 시계 반대방향(Counterclockwise)으로 형성 하면서 구조물을 완성해야
한다. 그 이유는 쉘 요소는 형상적으로 종이와 같은 얇은 면으로 형성이 되어 있어 시각적
으로 두께를 고려하여 입체화 한 모델링으로 구현하는 방식이 아닌 내부 적으로 두께를 입
력하고 모델은 면(Shell)으로 구성하는 모델링 기법이다.
이러한 모델링 방법에 의해서 쉘의 면(Face)이 상, 하 두개 및 응력의 결과는 3가지의 면

에서 산출할 수 있다.
-. 중간면(Middle Plane) : (ζ = 0)
-. 윗면(Top Plane) : Face 2(ζ = +1)
-. 아랫면(Bottom Plane) : Face 1(ζ = -1)
따라서 모델링을 구현할 때는 이웃 하는(같은 Plane) 같은 구조의 면 형성은 같은 방향으

로 요소의 방향이 형성이 되어야 한다.
η
ζ = + 1(Top Surface, Face 2)
ξ
ζ = 0(Middle Surface)
ζ = -1(Bottom Surface, Face 1)
그림[5.9] 쉘 요소의 법선 방향
다음의 그림의 요소 모델에서 요소의 법선 방향을 확인하면 다음과 같은 방향으로 표현이

되었다면 요소 1/4 와 2/3의 경계면에서는 같은 Plane 끼리 형성이 되지않고 다른 면으로
형성이 되어있기 때문에 평균 응력에서 문제가 발생한다.
특히 Graphic Post Processor에서 응력 성분을 Contours로 확인하면 Error가 발생한 이

상한 응력의 분포를 발견할 수 있다. 이러한 이유는 결과의 응력선도를 표시할 때는 Nodal
Average Stress를 사용하기 때문이다. 그리고 절점의 응력의 값은 그 절점에 연결된 요소의
응력의 평균값을 사용한다.
각 요소의 응력 분포는 Membrane Stress, Bending Stress를 합한 분포가 된다. 실 구조
물에서는 연속적이지만, Post Processor에서 Top Surface Stress Contours를 그리면 실제와
는 달리 응력의 값을 계산할 때는 사용자가 지정한 대로 좌측의 Top Surface에서의 응력은
σm+ σb, 우측의 Top Surface에서의 응력은 σm- σb이다.
그래서 절점에서의 응력은 (σm+ σb+ σm- σb)/2 = σm 이 된다. 이렇게 절점의 연결 순서를
반대로 하면 결과 검토 시 Error가 발생하는데 주의 하여야 한다.
쉘 요소에서 Membrane Stress가 주요성분이 되고 Bending Stress의 성분이 미약할 때는
Error의 크기가 미약하여 문제가 없지만 Bending Stress가 큰 경우에는 결과의 평가에 큰
영향을 미치는데 주의 해야 한다.
Face 1 Face 1
E4 E3
평균응력
E1 E2 [(σ m + σ b ) + (σ m − σ b )] 2
=σm
σm σb σm σb
Top Surface Bottom Surface
E1 E2
Bottom Surface Top Surface
그림[5.10] 쉘 요소의 절점에서의 응력 계산법
그러나 상기의 문제는 현재 NISA 프로그램에서는 시계반대방향으로 모델를 고려하지 않

고 엔지니어의 편리한 방법으로 모델를 구현 한 다음 꼭 방향성에 대하여 통일을 하라는 명
령을 수행하면 모든 같은 평면 내에서는 같은 방향으로 통일을 수행하여 주기 때문에 아주
편리하다고 볼 수 있다.
예)
Element ⇒ Modify ⇒ NISA Parameter ⇒ Element Face ⇒ Element Reference를 선택하여 요소
모델링 중의 임의 요소를 선택하여 모든 요소가 선택한 요소의 방향성에 대하여 통일 하라는 작업을
수행하면 상이한 면 들이 색상이 변하면서 일체의 방향성을 통일 시켜준다.

(8) 요소의 형상 Error 및 추천 방안
상기의 여러 부분에서 삼각형 요소(Triangular/Tetrahedron Element)보다 사각형요소

(Quadrilateral/Hexahedron Element)를 추천하고 있는데 사각형요소를 생성하였다 하여도
다음의 여러 경우를 참조하여 유한요소 모델를 구현하는데 인지 하여야 한다. 다소 프로그
램에 따라 차이는 있지만 대체로 같은 범주에 존재한다고 본다.
z 추천 요소의 순서
다음은 유한요소모델 구현하는데 높은 정도의 요소에서 낮은 정도 순서로 요소의 형상을

나열하였다. 정사각형이 가장 정도가 높다고는 하지만 구조물의 형상에 따라 불 가능한 부
분이 많기 때문에 대부분 직사각형 및 경사진 마름모 꼴의 형상을 많이 사용하고 있다.
b b
θ θ
a a
(a) 정사각형 (b) 직사각형 (d) 정삼각형 (e) 직 삼각형

[a : b= 1:1] [a : b = 1 : 3(3:1)] [모든 내각θ = 60°] [내각 120°>θ> 60°]
c
가로(세로) : 세로(가로) 어느 내각(기울기)이든
길이의 비가 b (120°>θ>60°)
응력의 평가에서는 1:3(3:1)
변형의 평가에서는 1:7(7:1) a
(c) 마름모꼴
[a : b= 1 : 3(3:1), a → b(c >a/2)
그림[5.11] 정도에 따른 요소 생성
z 세장비(Aspect Ratio) : 1:3(3:1), 1:7(7:1)
직사각형에서의 요소의 가로(세로) 대 세로(가로)의 비가 응력 분석을 원할 때에는 1 : 3

(3 : 1)를 기준으로 요소를 형성하고, 변형을 주로 분석을 원할 경우에는 응력 분석 기준 보
다 다소 세장비가 커도 결과에는 크게 영향을 주지 않기 때문에 1 : 7(7 : 1) 이 내에서 요소
를 형성해야 한다.
z 삼각형의 내각 크기 : 120° > θ > 60°
직삼각형에서의 요소 내의 어떤 내각의 크기가 최대 120° 보다는 작고 최소 60° 보다는

커야 한다.
z Skew Index(S.I) : 30°
그림에서
Skew Index = 90° - δ
7 9
5
⎡ v26 × v84 ⎤ 8 δ
δ = cos −1 ⎢ ⎥,
⎣⎢ v26 v84 ⎦⎥ 4
1 2 3
그림[5.12] 요소의 Skew Index
1-3-5-7 : Original quad

1-3-9-7 : Parallelogram enclosing the original quad
1,3,5,7 : Corner Nodes of original quad
2,4,6,8 : Mid-side edge points original quad
V26 : The vector joining the two mid-side points 2 and 6
V84 : The vector joining the two mid-side points 8 and 4
z Warping Index(W.I) : 0.2
5′
W .I = h
5
A
θ = sin (h / l )
−1
6
4
3
7
A 3′
7′
8 2
1′
θ
h
l 1
그림[5.13] 요소의 Warping Index
θ : θ = 0 for flat element

1-3-5-7 : Original quad
1′-3′-5′-7′ : Mean plane of the original quad with area ″A″
h : is the maximum value of the distance 11′,33′,55′ and 77′
z Twist Index(T.I) : 60°
그림에서
Twist Index = Max(θ1, θ2)
7 9
θ2
6
5
8
θ1
4
1 2 3
그림[5.14] 요소의 Twist Index
Maximum Angle : Maximum angle subtended by any two sides of the quad
Minimum Angle : Minimum angle subtended by any two sides of the quad
z Taper Index(Ta.I)
그림에서
Taper Index = (α/A)×4
ta = 1 ta → 0
A3
A4
A2
A1
그림[5.15] 요소의 Tape Index

[Taper]
α : Smallest triangular area(min of A1,A2,A3,A4)

A : = A1+A2+A3+A4
5.4 CAD 프로그램과의 연관관계
현재까지의 CAE 프로그램의 개발 동향의 흐름을 보면 프로그램의 특성을 갖기 위해서 3

단계의 흐름으로 크게 볼 수 있는데 전체적인 개발의 흐름 틀로 구분한 다면 다음과 같다.
(1) 1단계(전처리 단계)
앞의 장에서 언급 하였듯이 기하학적모델링을 통한 요소 및 절점 생성, 경계조건, 기계적

물성치, 모델특성 등을 부여하여 하나의 유한요소해석을 위한 전처리 단계를 주로 개발의
목표하는 일명 Modeler라 한다.
예) Hypermesh, I-DEAS, GridGen,…등
(2) 2단계(Solver 단계)
외부의 프로그램에서 만들어진 입력 데이터를 처리하기 위한 Solver 만을 제공하는 프로

그램으로 특히 구조부분에서 비선형 부분에 많이 이용하고 있다.
예) ABAQUS, ADINA, MARC(NASTRAN에 흡수),…등
(3) 3단계(후 처리 단계)
상기의 두 부분이 없이 단순히 외부에서 상기의 1, 2 단계를 수행한 결과를 좀더 강력한

그래픽 및 기타 결과의 시각성을 부여한 Post Processor 전용 프로그램이다. 특히 비 구조
부분, 유체 부분에서 많이 사용하나 구조 부분에서는 사용 빈도가 많지는 않다.
예) FIELDVIEW, Ensight,…
(4) 상기의 3단계 통합
사용상 편리 및 데이터의 호환성을 갖기 위해서 상기의 3단계(전처리, Solver, 후처리)를

통합한 프로그램이 일반적으로 현장에서 많이 사용하고 있다. 이러한 프로그램은 사용하기
쉽고 전체 프로그램의 흐름을 파악하기는 쉽지만 때로는 비선형부분(Nonlinear Static)에 약
한 부분이 있어 그러한 특수 목적에서는 앞에서 언급한 ABAQUS등의 프로그램과 연계하여
사용하기도 한다.
예) NISA, ANSYS, NASTRAN(PATRAN, MARC), COSMOS/M, RUSAS,…등
(5) CAD 테이타 이용
피상적으로 보면 CAD 프로그램의 데이터를 바로 CAE와 연계하여 사용하면 아주 쉽다고

는 볼 수 있지만 근본적인 모델링의 기법 등 근본 바탕이 상이하기 때문에 쉽게 이용하기에
는 아직 문제점이 많이 있다.
CAD 프로그램의 근본은 도면을 작성하기 위한 Drawing에 초점을 두는 것이고, CAE 프로
램은 그러한 점, 선, 면, 입체를 이용하여 최종적으로 요소 및 절점을 생성하기 위한 체계적
인 사각형 및 직육면체를 근본으로 하는 모델링 기법이다.
어떤 정량적인 결과 보다는 경향적인 결과를 원한다면 외부의 CAD 데이타를 이용하여
바로 Automesh 기능으로 요소를 생성하고 해석을 추구할 수 는 있다. 단, 경향적인 결과
분석에 주안점을 두어야 한다.
다음의 단계는 외부의 CAD 프로그램과의 연계관계 이다.
Autocad CAD System

Pro/Engineer
I-Deas DXF
Solid-Works IGES
- STL
CAD Data File
- Neutral
- Solid Surface Data
-
-
Translator
2D Plane Mesh
(Triangular, Quadrilateral)
NISA Automesh 3D Surface Mesh
Execution (Triangular, Quadrilateral)
3D Solid Mesh
(Tetrahedron, Hexahedron)
Mesh/Node ″Analysis″
Generation Boundary Condition
Material Data
-
-
-
그림[5.16] CAD와의 Interface
예) AUTOCAD DXF FILE Interface


Theoretical Course ″B″

1-1 : 선형 용수철 요소(Linear Spring Element)
1-2 : 축 방향 보 요소(Axial Bar Element)
1-3 : 굽힘 보 요소(Bending Beam Element)
1-4 : 평면 구조 요소(Plane Frame Element)
1-5 : 평면 격자 요소(Plane Lattice Element)
1-6 : 평면 요소(Plane Element)

1. 각 요소에 대한 평형 방정식
1-1 선형 용수철 요소(Linear Spring Element)
요소에 대한 수식을 유도함에 있어 앞서 언급한 바와 같이 대체로 4가지 접근 방법이 있

다. 그 중에서 직접법은 절점(Node)에서의 힘과 변위 사이의 관계식을 이용하여 직접 강성
행렬과 요소에 대한 평형 방정식을 얻는 방법이다.
이 방법은 간단한 대신 사용할 수 있는 범위가 상대적으로 간단한 문제에 국한되어 있다.
여기서는 선형 용수철에 대해서 직접법으로 수식을 유도하고 유한요소법(Finite Element
Method)을 적용하는 과정을 단계별로 보이고자 한다.
이를 통해서 유한 요소법의 기본적인 개념을 보다 쉽게 파악할 수 있을 것이다.
q2
Q2Y
x
f2
Node 2 (b) θ
Y 2
Q2X
q2
q1 f2
Q1Y F2Y
θ θ
1 2
Q1X F2X
Node 1 f1
F1Y
f1 q1 (c)
X θ
1
F1X
(a)
그림[1.1] 선형 용수철 요소
그림[1.1(a)]에 강성 계수(Stiffness Coefficient) 또는 용수철 상수(Spring Constant)가 k

인 용수철 요소의 양단에 국부 좌표계 x를 따라서 힘 f1과 f2가 각각 절점(Node) 1과 2에 작
용하는 경우를 보이고 있다.
이 때 절점에서의 x축상의 변위를 각각 q1과 q2라고 하면, 용수철에 있어서 힘과 변형 사
이의 관계식은 혹의 법칙에 따라 다음과 같다.
f = k⋅q
여기서 f는 힘, q는 변형이다. 절점 1 및 절점2에 대해서 각각 이 식을 적용하면,
f1 = k (q1 − q 2 ),
f 2 = k (q 2 − q1 )
이를 행렬식으로 표시하면,
[k ] {q} = {f } (1-1)
여기서 {q}T = {q1 , q 2 }, { f }T = { f1 , f 2 }, 그리고,
⎡ 1 − 1⎤
[k ] = k⎢ ⎥ (1-2)
⎣− 1 1 ⎦
가 된다.
변위 q1, q2를 대표 죄표계 X, Y 방향의 성분으로 나타내면,
그림[2.1(b)]에서,
⎧ Q1 X ⎫
⎧ q1 ⎫ ⎡C S 0 0 ⎤ ⎪⎪ Q1Y ⎪⎪ (1-3)
⎨ ⎬ = ⎢ ⎨ ⎬
⎩q 2 ⎭ ⎣0 0 C S ⎥⎦ ⎪Q2 X ⎪
⎪⎩ Q2Y ⎪⎭
가 되어 {q} = [T ] {Q}의 형태가 된다. 여기서 C = cosθ, S = sinθ, 그리고 θ 는 x 축과 X축

사이의 각도이다.
그러니까 여기서 변환 행렬 [T]는,
⎡C S 0 0⎤
{T } = ⎢
⎣0 0 C S ⎥⎦
대표 좌표계 방향의 힘의 성분을 f1 및 f2로서 나타내면, 그림[2.1(c)]에서
⎧ F1 X ⎫ ⎡C 0⎤
⎪F ⎪ ⎢S 0 ⎥⎥ ⎧ f 1 ⎫
⎪ 1Y ⎪
⎨ ⎬ = ⎢ ⎨ ⎬ (1-4)
F
⎪ 2X ⎪ ⎢0 C⎥ ⎩ f2 ⎭
⎪⎩ F2Y ⎪⎭ ⎢ ⎥
⎣0 S⎦
가 되어 식 {F } = [T ]T { f } 와 같은 형태가 된다. (1-5)
식(1.1)의 양변의 앞에 각각 [T]T를 곱하고 식(1.3)을 대입하여 정리하면, 다음과 같은

형태의 대표좌표계에 대한 평형방정식 {K }{Q} = {F } 가 구해지고, 다음이 유도된다.
⎡C 2
CS −C2 − CS ⎤
⎢ ⎥
S2 − CS − S2 ⎥
[K ] = [T ]T [k ][T ] = k⎢
⎢ (1-6)
C2 CS ⎥
⎢ ⎥
⎢⎣ S 2 ⎥⎦
식(1.2)에 주어진 강성 행렬은 축 방향 및 비틀림 강성 계수가 그림[1.2]에서와 같이 주

어짐을 이용하면, 그대로 구조물(Framed Structures)에 대한 축 방향 보요소(Axial Bar
Element) 및 비틀림 보요소(Torsional Bar Element)에 대한 강성 행렬이 된다.
k = EA , E = young ' s mod ulus

(a) 축 방향 강성계수 l
A = cross − sec tion area
l = length
(b) 비틀림 강성계수 k = GJ , GJ = torsional rigidity

l
l = length
그림[1.2] 축 방향 및 비틀림 강성계수
Beam은 굽힘강성(Bending Stiffness)이 있는 경우에 대해서, Bar는 굽힘 강성이 없는 경

우에 대해서 쓰인다.
축 방향 강성 계수 k = EA / l 이므로, 축 방향 요소에 대한 강성 행렬은
⎡ 1 − 1⎤
[k ] =
EA
⎢− 1 1 ⎥ (1-7)
l ⎣ ⎦
여기서 E는 영률(Young’s Modulus), A는 단면적, l은 요소의 길이를 나타낸다.
비틀림 강성 계수는 k = GJ / l 이므로, 비틀림 요소에 대한 강성 행렬은
GJ ⎡ 1 − 1⎤
[k ] = (1-8)
l ⎢⎣ − 1 1 ⎥⎦
여기서 G는 전단계수(Shear Modulus), J는 극 관성 모맨트(Polas Moment of Inertia)이다.

평면 구조물에 대한 대표 좌표계에서의 강성 행렬은 식(1.6)으로 주어진다.

1-2 축 방향 보 요소(Axial Bar Element)
앞 절에서는 직접법으로 축 방향 요소에 대한 강성 행렬을 식(1.7)과 같이 구했다. 여기

서는 같은 강성 행렬을 변위법을 사용하여 유도하는 과정은 다음과 같다.
Y
x
Q2Y
q2
Q2X
2
l
y Q1Y
q1
Q1X
1
그림[2.1] 2차원 축 방향 보 요소
그림[2.1]에 보이는 축 방향 보요소(Axial Bar Element)의 탄성계수(Modulus of Elasticity,

Young’s Modulus)를 E, 초기의 길이를 l 이라 하고, 단면적 A가 일정한 것으로 간주한다.
보요소의 양단을 절점 1 및 2라 두고, 절점에서의 요소의 길이 방향 변위를 그림에서와
같이 각각 q1 및 q2 라고 하기도 한다.
(1) 요소 내의 변위를 국부 좌표를 사용하여 연속 함수로 가정한다.
절점 1로부터의 길이 방향 좌표축 x를 따라 변위 u가 선형적으로 일어나는 것으로 가정하

면,
u = C1 + C 2 x, (2-1a)
여기서 미지 지수 계수 Ci의 개수는 항상 요소의 자유도와 같아야 한다.

식(2.1a)를 행렬식으로 나타내면
⎧C ⎫
u = [1, x] ⎨ 1 ⎬ (2-1b)
⎩C 2 ⎭
2) 가정한 변위 함수를 절점(Node)의 값으로 표시한다.
미지 계수 Ci는 양단에서의 조건에 의하여 결정된다. 즉 x = 0 에서 u = q1,

x = l 에서 u = q2 이므로 식(2.9a)로 부터
q1 = C1 , q 2 = C1 + C 2 l
이를 행렬식으로 나타내면,
⎧ q1 ⎫ ⎡1 0 ⎤ ⎧ C1 ⎫
⎨ ⎬ = ⎢ ⎥⎨ ⎬
⎩q2 ⎭ ⎣1 l ⎦ ⎩C 2 ⎭
따라서 계수 Ci의 해는,
⎧ C1 ⎫ ⎡ 1 0 ⎤ ⎧ q1 ⎫
(2-2)
⎨ ⎬ = ⎢− 1 1⎥ ⎨ ⎬
⎩C 2 ⎭ ⎢⎣ l l ⎥⎦ ⎩q 2 ⎭
식(2.2)을 (2.1b)에 대입하면,
⎡ 1 0⎤ ⎧ q ⎫ ⎧q ⎫
u = [1, x ] ⎢ 1 1 ⎥ ⎨ 1 ⎬ = [N 1 , N 2 ] ⎨ 1 ⎬ = [N ] {q} (2-3)
⎢⎣ − l l ⎥⎦ ⎩q 2 ⎭ ⎩q 2 ⎭
여기서 N1 = 1 – x / l, N2 = x / l 과 같은 형태의 식을 얻을 수 있다.
N1 N2
1 1
x
1 2
그림[2.2] 1차원 형상 함수
그림[2.2]에 형상 함수의 모양을 보여지고 있다. 한 가지 주목할 것은 형상 함수 N1은 절

점 1에서는 값이 1이고 점점 2에서는 0, 반면에 N2는 점점 2에서는 1이고 절점 1에서는 0
이라는 사실이다.
이를 수식으로 표시하면,
N i (x j ) = δ ij (2-4)
여기서 δij 는 Kronecker Delta 이다.
(3) 변형률과 변위 및 응력과 변형률 관계식이 만족되도록 한다.
변형률과 변위 관계식(Strain-Displacement Relationship)은
∂u ⎡ ∂N ∂N 2 ⎤ ⎧ q1 ⎫ ⎧q ⎫
εx = = ⎢ 1, ⎥⎨ ⎬ = [B1 , B2 ] ⎨ 1 ⎬ = [B ] {q} (2-5)
∂x ⎣ ∂x ∂x ⎦ ⎩ q 2 ⎭ ⎩q 2 ⎭
여기서 B1 = ∂N 1 / ∂x = −1 / l , B2 = ∂N 2 / ∂x = 1 / l 와 같은 형태의 식을 얻을 수 있다.
응력과 변형률 관계식(Stress-Strain Relationship)은
σ x = Eε x =
E
(q 2 − q1 ) (2-6)
l
(4) 요소에 대한 강성 행렬 및 절점 부하를 일 또는 에너지 법칙을 이용하여 결정한다.
요소의 변형에너지는 훅의 법칙을 따르는 선형 물질에 대해서는
∫ {q} [B ] E [B ]{q}dV
1 1 1
2 ∫vol ∫
U = σ x ε x dV = ε xT Eε x dV =
T T
2 vol 2 vol
=
1
{q}T [k ]{q} (2-7)
2
여기서
[k ] ∫ [B ] E [B ] dV ∫ [B ] E [B ]A dx
l
= =
T T
vol 0
이므로, [B] = [− 1 / l , 1 / l ] 을 대입하면
⎡ 1 − 1⎤
[k ] =
EA
⎢− 1 1 ⎥ (2-8)
l ⎣ ⎦
이 식은 (1)절에서 직접법으로 구한 식(1.7)과 꼭 같다.

절점 부하{f}라 하면, 이 힘에 의한 일은,
W = ∑f i qi = {q}T { f } (2-9)
i
그런데 힘이 x축 방향으로 절점에 작용할 때는 그 값이 그대로 절점 부하가 된다.

그러나 힘이 요소의 길이를 따라서 분포하중으로 주어질 때는 절점에서의 상당 부하로 환
산해 주어야 한다. 예를 들어, 단위 길이 당 힘 Px가 x축을 따라서 선형적으로 변하는 분포
하중의 경우에 대해서 생각해 보기로 한다. 절점 1,2에서의 각각 P1, P2라 하면,
Px = P1 + (P2 − P1 ) x / l (2-10)
이 힘에 의한 일은,
∫ P {q} [N ] {q}T { f }0
l l
= ∫ Px u dx = dx =
T T
W x
0 0
여기서 상당 절점 부하 {f}0 는,
⎧ x ⎫
∫ {P + (P − P ) x l }⎨ x
⎧ f1 ⎫ ⎪1 − l ⎪ dx
∫ Px [N ] dx =
l l
⎨ ⎬ = 1 2 1 ⎬
⎩ f 2 ⎭0 0 0
⎪⎩ l ⎪⎭
l ⎧2 P1 + P2 ⎫ (2-11)
= ⎨ ⎬
6 ⎩ P1 + 2 P2 ⎭
덧붙여, 상당 절점 하중이 존재할 경우, 실제에 작용하는 힘은 상당 절점 하중을 뺀 값,

즉 {f}-{f}0 이 된다.
최소 에너지 법칙을 적용하면,
⎡ ∂ (U + W ) ⎤
⎢ ⎥ = 0 , i = 1, 2
⎣ ∂q i ⎦
여기서 식(2.7)와 (2.9)을 대입하면 다음과 같은 형태의
[k ] {q} = {f } (2-12)
를 얻는다.
(5) 국부 좌표계로 나타낸 강성 행렬 및 절점 하중을 대표 좌표에 대한 것으로 변환한다.
평면 트러스(Plane Truss)와 같은 2차원 문제에 대한 대표 강성 행렬 및 대표 절점 하중

등은 앞 절에서와 꼭 같이 변환 된다. 이를 다시 정리하면,
{q} = [T ] {Q}, {F } = [T ]T { f }
⎡C 2 CS −C2 − CS ⎤
⎢ ⎥
S2 − CS − S2 ⎥
[K ] =
EA ⎢ (2-13)
l ⎢ C2 CS ⎥
⎢ ⎥
⎢⎣ S 2 ⎥⎦
[K ] {Q} = {F }
여기서
{Q}T = {Q1 X , Q1Y , Q2 X , Q2Y }, {F }T = {F1 X , F1Y , F2 X , F2Y }
⎡C S 0 0⎤
[T ] = ⎢ (2.14)
⎣0 0 C S ⎥⎦
이고, C=cosθ, S=sinθ, 그리고 θ는 국부 좌표축의 대표 좌표축에 대한 각도이며, {Q}({F})

의 방향은 그림[2.3]에 보이는 바와 같다.
z′ x
y θZ
θY
y′
θX z
x′
Q1X
그림[2.3] 3차원 축 방향 보 요소
X
이제 3차원 공간에서의 임의의 방향으로 기울어져 있는 요소에 대한 변환에 대해서 알아

보자.
그림[2.3]에 보이는 바와 같이, 대표 좌표축 X, Y, Z에 대한 절점 1의 변위를 Q1X, Q1Y, Q1Z,
절점 2의 변위를 Q2X, Q2Y, Q2Z라 하면, 변환 행렬 [T]는,
⎡C CY CZ 0 0 0⎤
[T ] = ⎢ X (2.15)
⎣ 0 0 0 CX CY C Z ⎥⎦
여기서 CX=cosθX, CY=cosθY, CZ=cosθZ, 이고, θX, θY, θZ 는 요소와 대표 좌표축 사이의 각도
이다.
따라서 강성 행렬 [K]는 다음과 같다.
[K ] = [T ]T [k ][T ]
⎡C X2 C X CY C X CZ − C X2 − C X CY − C X CZ ⎤
⎢ ⎥
⎢ CY2 CY C Z − C X CY − CY2 − CY C Z ⎥
EA ⎢ C Z2 − C X CZ − CY C Z − C Z2 ⎥ (2.15a)
= ⎢ ⎥
l ⎢ C X2 C X CY C X CZ ⎥
⎢ CY2 CY C Z ⎥
⎢ ⎥
⎢⎣ C Z2 ⎥⎦
덧붙여, 식(2.15a)를 검토해 보면, 절점과 연관하여 같은 값은 값이 반복됨을 알 수 있다.

이런 경우에는 부분 행렬(Submatrices)로 나누어 표시하면 간편하게 된다. 그러면 식
(2.15a)는 다음과 같이 나타낼 수 있다.
⎡[K ], [K12 ]⎤
[K ] = ⎢ 11 (2.15b)
⎣[K 21 ], [K 22 ]⎥⎦
따라서
⎡C X2 C X CY C X CZ ⎤
EA ⎢ ⎥
[K 11 ] = [K 22 ] = − [K 12 ] = − [K 21 ] = C 2
CY C Z ⎥
l ⎢
Y
⎢ C Z2 ⎥⎦
⎣
1-3 굽힘 보요소(Bending Beam Element or Flexural Beam Element)
보(Beam)는 보통 앞절에서 취급한 축 방향 변형보다 훨씬 많은 변형을 일으키고, 실제에

있어서 자주 횡 방향 하중에 의한 굽힘 변형이 축 방향 변형보다 더 큰 문제가 된다. 이 절
에서는 단순 보(Simple Beam or Slender Beam) 요소에 횡 방향 하중이나 모멘트가 주어졌
을 경우에 대해서 설명하기로 한다.
q1
ν q3
2
1 x
q2 q4
x
z
그림[3.1] 굽힘 보 요소
그림[2.6]에 보이는 바와 같이 길이가 l 인 직선 보 요소의 양단을 절점 1및 2라 하고,

보의 축 방향 지역 좌표를 x, 횡 방향 국부 좌표를 y라 하자. 그리고, 강성 행렬에 대한 수
식을 유도하기 전에 먼저 부호 규약을 다음과 같이 정하기로 한다.
-. 회전과 모멘트는 시계 반대 방향일 때 양이다.

-. 변위와 힘은 양의 y방향일 때 양이다.
(1) 요소 내의 변위를 지역 좌표를 사용하여 연속 함수로 가정한다.
요소의 각 절점의 자유도를 횡 방향 변위 ν 와 회전 φ (2자유도)로 가정한다. 그러면 절

점 변위 벡터 {q}는 다음과 같다.
{q}T = {q1 , q 2 , q3 , q 4 } = {ν 1 , φ1 , ν 2 , φ 2 } (3.1)
여기서 νi 는 절점 i 에서의 횡 방향 변위이며, 회전 φi 는
φ i = dν i dx
(3.2)
로 주어진다.
요소의 자유도가 4이므로, 두 절점 사이의 변위 ν 는 3차 다항식(Cubic Polynominal)으로
가정할 수 있다.
ν ( x ) = C1 + C 2 x + C 3 x 2 + C 4 x 3 = [S ] {C } (3.3)
여기서 S [ ] [ ]
= 1, x, x 2 , x 3 이고, {C}
T
= {C1 , C 2 , C 3 , C 4 } 이다.
(2) 가정한 변위 함수를 절점의 값으로 표시한다.
절점 조건은 ν 1 = ν (0 ), φ1 = (∂ν ∂x )x =0 , ν 2 = ν (l ), φ 2 = (∂ν ∂x )x =l 이므로,

이 들의 식(3.3)에 대입하여 행렬식으로 정리하면,
⎧ν 1 ⎫ ⎡1 0 0 0 ⎤ ⎧ C1 ⎫
⎪φ ⎪ ⎢0 1 0 0 ⎥⎥ ⎪⎪C 2 ⎪⎪
⎪ 1⎪ ⎢
⎨ ⎬ = ⎨ ⎬ (3.4)
ν
⎪ 2⎪ ⎢1 l l 2 l 3 ⎥ ⎪C 3 ⎪
⎪⎩φ 2 ⎪⎭ ⎢ ⎥
⎣0 1 2l 3l 2 ⎦ ⎪⎩C 4 ⎪⎭
{C}의 해를 구하면,
⎧ C1 ⎫ ⎡ l3 0 0 ⎤ ⎧ν 1 ⎫
0
⎪C ⎪ ⎢ ⎥⎪ ⎪
⎪ 2⎪ 1 ⎢ 0 l3 0 ⎥ ⎪φ 1 ⎪
0
⎨ ⎬ = 3 ⎢ ⎨ ⎬
⎪C 3 ⎪ l − 3l − 2l 2 3l − l 2 ⎥ ⎪ν 2 ⎪ (3.5)
⎪⎩C 4 ⎪⎭ ⎢ ⎥
⎢⎣ 2 l −2 l ⎥⎦ ⎪⎩φ 2 ⎪⎭
이 식을 식(3.3)에 대입하면,
ν = [N ] {q} = [N1 N 2 N3 N4 ] {q} (3.6)
여기서
2 3
N1 = 1 − 3x 2 + 2x ,
l l3
N2 = lx (l − 2x
2
l 2 + x
3
l3
),
2 3
N3 = 3x 2 − 2x ,
l l3
N4 = l−x(2
l2)+ x3
l3
형상 함수 Ni 는 그림[3.2]에 보이는 바와 같이, 해당되는 절점에 해당 단위 변위(회전) 만

주어지고 다른 절점 또는 다른 변위는 0일때의, 요소의 길이를 따라 변위 ν가 변하는 모양
을 나타낸다.
따라서, N i (qi ) = δ ij
N1
ㅣㅣ
(a)
N2
(b)
N3
ㅣ
(c) ㅣ
N4
(e)
그림[3.2] 굽힘 보 요소의 형상 함수
(3) 변형률과 변위 및 응력과 변형률 관계식이 만족되도록 한다.
변형률과 변위 관계식은 단순 보 이론(Simple Beam Theory)의 가정 중 평면 단면은

변형 후에도 그대로 평면이라고 하는 가정으로부터 주어진다. 그림[3.3]에 보이는 것과 같
이 단면 위의 임의의 점에서의 x방향의 변위 u는 중립면으로부터의 거리를 y라하면,
∂ν
u = −y (3.7)
∂x
이므로
∂u ∂ 2ν
εx = = −y 2 = − y [N ′′] {q} (3.8)
∂x ∂x
여기서
″ ″ ″ ″
[B ] = − y ⎡ N1 , N 2 , N 3 , N 4 ⎤
⎢⎣ ⎥⎦
″
N1 = (− 6 + 12 x l ) / l 2 ,
″
N2 = (− 4 + 6 x l ) / l ,
″
N3 = (6 − 12 x l ) / l 2 ,
″
N4 = (− 2 + 6 x l ) / l
y,ν
dA
y x,u
z dv/dx
dx
그림[3.3] 굽힘 변형
응력과 변형률 관계식은 다음과 같다.
σ x = E ε x = E [B ] {q} (3.9)
이 된다.
(3) 요소에 대한 강성 행렬 및 절점 부하를 일 또는 에너지 법칙을 이용하여 결정한다.
∫ {q} [B ] E [B ] E [B ]{q}dV
1 1 1
U = ∫ σ x ε x dV = ∫ ε xT Eε x dV = T T T
2 vol 2 vol 2 vol
=
1
{q}T [k ]{q} (3.10)
2
여기서
[k ] ∫ ∫ [B ] E [B ] dAdx y 2 dA ∫ [N ′′] [N ′′] dx

l A A l
= = E ∫
T T
0 0 0 0
⎡12 6l − 12 6l ⎤
⎢ 4l 2 − 6l 2l 2 ⎥⎥ (3.11)
EI ⎢
= 3
l ⎢ 12 − 6l ⎥
⎢ ⎥
⎣ 4l 2 ⎦
윗 식에서 I는 단면의 중립축(Neutral Axis)에 대한 관성 모멘트이며, 단면적은 길이를 다

라 일정한 것으로 간주했다.
평형 방정식
[k ] {q} = {f } (3.12)
여기서 { f }T = { f1 , m1 , f 2 , m2 } 이며 f1, f2는 절점에서의 전단력, m1, m2는 절점에서의 굽힘
모멘트이다. 그런데 {f}가 절점에서의 값으로 주어질 때는 그것이 그대로 절점 부하가 되지

만, 요소가 길이를 따라 변하는 하중을 받고 있는 경우에는 이를 절점에서의 상당 부하로
환산해서 더해 주어야 한다. 따라서 이 경우에 실제로 절점에 작용하는 하중은 상당 절점
부하를 뺀 값이 된다.
여기서는 선형적으로 변하는 분포 하중의 경우에 대해서 생각해 보기로 한다.
절점 1 및 2에서의 단위 길이 당 힘을 p1 및 p2라 하면 분포하중은 px는 다음과 같다.
px = p1 + ( p 2 − p1 ) x (3.13)
l
이 힘에 의한 일은,
p x {q} [N ] dx = {q} { f }0
l l
W = ∫ p x νdx = ∫
T T T
(3.14)
0 0
{f } = ∫ p x [N ] dx 이며,
l T
따라서 상당 절점 부하 {f}0 는
0
⎧ f1 ⎫ ⎧ 21 p1 + 9 p 2 ⎫
⎪m ⎪ ⎪ ⎪
⎪ ⎪ l ⎪ (3 p1 + 2 p 2 )l ⎪
{ f }0 = ⎨ 1⎬ = ⎨ ⎬
(3.15)
⎪ f2 ⎪ 60 ⎪ 9 p1 + 21 p 2 ⎪
⎪⎩m2 ⎪⎭ ⎪⎩(− 2 p1 − 3 p 2 )l ⎪⎭
요소가 대표 좌표축에 대해서 기울어져 있을 때, 대표 좌표계에 대한 강성 행렬은 굽힘

요소자체 만으로는 실제적인 의미가 없으므로 생략한다.

1-4 평면 구조 요소(Plane Frame Element)
그림[4.1]에 보이는 평면 구조 요소는 축 방향과 횡 방향 변형이 모두 일어날 수 있다.

따라서 각 절점의 자유도는 3(UX, UY, RM)이 된다. 그런데 국부 좌표계에서는 축 방향 변형
과 횡 방향 변형이 서로 독립이므로, 이러한 경우에는 축 방향 변위와 횡 방향 변위를 연립
하여 취급하게 된다. 바꾸어 말하면, 앞에서 언급한 축 방향 보 요소에 대한 수식과 굽힘
요소에 대한 수식을 단순히 결합함으로써 평면구조요소의 국부 좌표계에 대한 수식을 얻을
수 있다.
그러므로 변위는 식(2.3)과 (3.6)으로부터 다음과 같이 정리될 수 있다.
⎧u ⎫ ⎡ N1 0 0 N4 0 0⎤
⎨ ⎬ = ⎢ {q} (4.1)
⎩v ⎭ ⎣0 N2 N3 0 N5 N 6 ⎥⎦
여기서 {q}T = {q1x q1 y φ1 q 2 x q 2 y φ 2 },
x x
N1 = 1 − , N4 = ; 축 방향 요소
l l
2 3
N2 = 1− 3 x 2 +2x ,
l l3
(
2
N3 = l x − 2 x 2
l l
+x
3
l3
),
2 3
굽힘 요소
N5 = 3 x −2x 3,
l2 l
N6 = l−x (
2
l2
+x
3
l3
),
이며, u는 축 방향 변위, v는 횡 방향 변위를 뜻한다.
하중 벡터를 { f }T = {f1x f 1 y m1 f 2 x f 2 y m2 }라 하면,
구조 요소(Frame Element)에 대한 평형방정식은
(4.2)
[k ] {q} = {f }
Y
x
Q2Y
q2x
q2y
Q2X
l m2, M2
y Q1Y
q1y q1x
Q1X
1
m1, M1
X
그림[4.1] 평면 구조 요소(Plane Frame Element)
여기서 강성 행렬은 축 방향 요소에 대한 강성 행렬, 식(2.8)과 굽힘 요소에 대한 강성 행

렬, 식(3.11)를 결합함으로써 다음과 같디 얻어지게 된다.
⎡ EA EA ⎤
⎢ l 0 0 − 0 0⎥
l
⎢ 12 EI 6 EI 12 EI 6 EI ⎥
⎢ 0 − ⎥
⎢ l3 l2 l3 l2 ⎥
⎢ 4 EI
0
6 EI
− 2
2 EI ⎥
⎢ l ⎥
[k ] = ⎢ l
EA
l
⎥ (4.3)
⎢ 0 0 ⎥
⎢ l ⎥
⎢ 12 EI 6 EI ⎥
− 2
⎢ l3 l ⎥
⎢ 대 칭 4 EI ⎥
⎢ ⎥
⎣ l ⎦
이제 평면 구조 요소의 대표 좌표계에 있어서의 수식을 얻기 위하여 변위{q}를 {Q}의 성

분으로 나타내면 그림[2.9]에서
[q ] = [T ] {Q} (4.4)
여기서 {Q} = {Q1 X Q1Y Φ 1 Q2 X Q2Y Φ 2 }

T
⎡C S 0 0 0 0⎤
⎢− S C 0 0 0 0⎥⎥
⎢ (4.5)
⎢ 0 0 1 0 0 0⎥
[T ] = ⎢ ⎥
⎢ 0 0 0 C S 0⎥
⎢ 0 0 0 −S C 0⎥
⎢ ⎥
⎣⎢ 0 0 0 0 0 1⎦⎥
이고 C=cosθ, S=sinθ를 뜻하며, θ는 국부 좌표축 x와 대표 좌표축 X사이의 각도이다.

대표 좌표계에서의 평형방정식은
[k ] {Q} = {F } (4.6)
여기서 {F }T = {F1 X F1 y M 1 F2 X F2 y m2 }이며
{F } = [T ]T { f } (4.7)
⎡[K ] [K12 ]⎤
[K ] = [T ]T [k ][T ] = ⎢ 11
⎣[K 21 ] [K 22 ]⎥⎦ (4.8)
⎡ 12 I 2 ⎛ 12 I ⎞ 6I ⎤
⎢ AC + l 2 S , ⎜ A − 2 ⎟CS , −
2
l ⎥⎥
S
⎢ ⎝ l ⎠
[K11 ] =
E ⎢ 12 I
AS 2 + 2 C 2 ,
6I ⎥
C
l ⎢ l l ⎥
⎢ 4I ⎥
⎢ 대 칭 ⎥
⎣ ⎦
⎡ ⎛ 12 I 2 ⎞ ⎛ 12 I ⎞ 6I ⎤
⎢ − ⎜ AC + l 2 S ⎟, − ⎜ A − 2 ⎟CS , −
2
l ⎥⎥
S
⎢ ⎝ ⎠ ⎝ l ⎠
E ⎢ ⎛ 12 I ⎞ ⎛ ⎞ 6I ⎥
[K12 ] = [K 21 ]T = − ⎜ A − 2 ⎟CS ,
12 I
− ⎜ AS 2 + 2 C 2 ⎟, C
l ⎢ ⎝ l ⎠ ⎝ l ⎠ l ⎥
⎢ 6I 6I ⎥
⎢ S, − C 2I ⎥
⎣⎢ l l ⎥⎦
⎡ 12 I 2 ⎛ 12 I ⎞ 6I ⎤
⎢ AC + l 2 S , ⎜ A − 2 ⎟CS ,
2
S ⎥
⎢ ⎝ l ⎠ l ⎥
[K 22 ] =
E ⎢ 12 I 2
AS + 2 C , − C
2 6I ⎥
l ⎢ l l ⎥
⎢ 4I ⎥
⎢ 대 칭 ⎥
⎣ ⎦

1-5 평면 격자 요소(Plane Lattice Element)
앞 절의 평면 구조 요소는 평면 구조물(Plane Frame)에 힘이 평면 내에서 작용할 경우에

대한 요소이다. 본 요소는 같은 종류의 평면 구조물에 그림[5.1]와 같이 힘이 평면에 수직
인 방향으로 작용할 경우에 대해서 고찰하면 다음과 같다.
Z
그림[5.1] 평면 격자 구조(Plane Lattice Structure)
이 경우의 변형은 주로 비틀림 변형과 힘 방향 변형이 된다. 그런데 국부 좌표계에서 각

요소에 대한 두 가지 변형은 서로 독립적이므로, 평면 구조 요소에서와 마찬가지로 두 변형
에 대한 수식을 단순히 연립하여 취급하면 다음과 같다.
변위는
⎧φ x ⎫ ⎡0 N1 0 0 N4 0⎤
⎨ ⎬ = ⎢ {q} (5.1)
⎩ν ⎭ ⎣N2 0 N3 N5 0 N 6 ⎥⎦
여기서 {q}T = {q1 y , φ1x , φ1z , q 2 y , φ 2 x , φ 2 z },
x x
N1 = 1 − , N4 = ; 비틀림 요소
l l
N2 = 1 − 3 x3 l 2 + 2 x3 l 3 ,
(
N3 = l x l − 2x2 l 2 + x3 l 3 )
굽힘 요소
N 5 = 3x 2 l 2 − 2 x3 l 3
(
N 6 = l − x2 l 2 + x3 l 3 )
이며, φx는 회전, ν는 평면에 수직인 방향으로의 변위를 뜻한다.

하중 벡터 { f }T = {f1 y , m1x , m1z , f 2 y , m2 x , m2 z }라 두면, 격자 요소에 대한 평형 방정식
은 다음과 같다.
[k ] {q} = {f } (5.2)
여기서 강성 행렬은 비틀림 강성 행렬, 식(1.8)과 같은 강성 행렬, 식(3.11)를 결합함으로

써 얻어지게 된다. 즉,
⎡12 EI 6 EI 12 EI 6 EI ⎤
⎢ l3 0 − 0
l2 l3 l2 ⎥
⎢ GJ GJ ⎥
⎢ 0 0 0 ⎥
⎢ l l ⎥
⎢ 4 EI
−
6 EI
0
2 EI ⎥
⎢ l ⎥
[k ] = ⎢ l l2
12 EI 6 EI ⎥
(5.3)
⎢ 0 ⎥
⎢ l3 l2 ⎥
⎢ GJ
0 ⎥
⎢ l ⎥
⎢ 대 칭 4 EI ⎥
⎢ ⎥
⎣ l ⎦
여기서 한 가지 상기할 것은 전단 계수(G)와 단성계수(E) 사이의 관계식은 포이송 비

(Poisson’s Ratio)를 ν라 나타내면 G = E/2(1+ν)이다.
{Q}T = {Q1Y , Φ1 X , Φ 1Z , Q2Y , Φ 2 X , Φ 2 Z , } 로 두면 {q} = [T ]{Q} 에서의 변환 행렬[T]는
그림[5.2]에서 다음과 같이 주어진다.
⎡1 0 0 0 0 0⎤
⎢0 C S 0 0 0 ⎥⎥
⎢
⎢0 − S C 0 0 0⎥
[T ] = ⎢ ⎥ (5.4)
⎢0 0 0 1 0 0⎥
⎢0 0 0 C S⎥
⎢ ⎥
⎢⎣0 0 0 0 −S C ⎥⎦
여기서 C=conθ, S=sinθ 를 뜻하며 θ 는 국부 좌표계의 x축과 대표 좌표계의 X축 사이의

각도이다.
φ1x
q1y Q1Y
Y
Φ1X
θ
φ1z
Φ1Z
X
Z
q2y Q2Y
Φ2X
φ2z φ2x
Φ2Z
그림[5.2] 평면 격자 요소(Plane Lattice Structure)

(2중 화살표는 회전, 방향은 오른손 법칙)
대표 좌표계에서의 강성 행렬은 다음과 같다.
⎡[K ] [K12 ]⎤
[K ] = [T ]T [k ][T ] = ⎢ 11
⎣[K 21 ] [K 22 ]⎥⎦ (5.5)
여기서
⎡12 EI 6 EI 6 EI ⎤
⎢ l , − S, C ⎥
1⎢ l l
[K11 ] = ⎢ GJC + 4 EIS , (GJ − 4 EI )CS ⎥⎥
2 2
l
⎢ GJS 2 + 4 EIC 2 ⎥
⎢⎣ 대 칭 ⎥⎦
⎡ 12 EI 6 EI ⎤ 6 EI
⎢− l , −
l
S, ⎥ l
C
1 ⎢ 6 EI ⎥
[K12 ] = [K 21 ] = ⎢ S GJC 2 + 4 EIS 2 , − (GJ − 2 EI )CS ⎥
l⎢ l ⎥
⎢ − 6 EI C − (GJ + 2 EI )CS , − GJS + 2 EIC ⎥
2 2⎥
⎢⎣ l ⎦
⎡12 EI 6 EI 6 EI ⎤
⎢ l2 , S, − C ⎥
1⎢ l l
[K 22 ] = ⎢ GJC 2
+ 4 EIS 2
, (GJ − 4 EI )CS ⎥⎥
l
⎢ GJS 2 + 4 EIC 2 ⎥
⎢⎣ 대 칭 ⎥⎦

1-6 평면 요소(Plane Element)
(1) 일반화된 훅의 법칙(Generalized Hook’s Law)
많은 물질에 있어서 일축(Uniaxial) 상태에서의 변형은 힘에 비례한다. 이러한 성질을 발

견자의 이름을 따서 훅의 법칙(Hook’s Law)이라고 한다. 이러한 힘과 변형 사이의 관계는
그대로 식(2.6)와 같은 응력과 변형률 사이의 비례관계로 표시될 수 있다.
3차원 물체의 경우에는 응력 과 변형률은 각각 6개의 성분으로 이루어진다. 훅의 법칙을 3
차원 물체로 연장하면 변형률의 각 성분은 응력 성분들의 선형 함수로 표시하고, 또한 응력
의 각 성분은 변형률 성분의 선형 함수로 표시될 수 있다. 이를 일반화된 훅의 법칙이라고
한다.
이를 행렬식으로 표시하면,
{ε } = [C ] {σ } (6.1)
{σ } = [D ] {ε } (6.2)
여기서 변형률 벡터 {ε }은 {ε }T = {ε x , ε y , ε z , γ yz , γ xz , γ xy },
응력 벡터 {σ }는 {σ }T = {σ x , σ y , σ z , τ yz , τ xz , τ xy },
[D ] = [C ]−1
여기서 [C], [D]는 각각 6×6 행렬이다.

선형 등방성 균질 물질(Linear Isotropic Uniform Material)에 대해서 [C] 및 [D]는 영률
(Young’s Modulus)과 포이송비(Poisson’s Ratio)만으로 표시할 수 있다. 변형률과 응력 관
계식(Strain-Stress Relationship)은 다음과 같다.
⎧εx ⎫ ⎡1 −ν −ν 0 0 0 ⎤ ⎧σ x ⎫
⎪ε ⎪ ⎢− ν ⎪ ⎪
⎪ y⎪ ⎢ 1 −ν 0 0 0 ⎥⎥ ⎪σ y ⎪
⎪⎪ ε z ⎪⎪ 1 ⎢− ν −ν 1 0 0 0 ⎥ ⎪⎪σ z ⎪⎪
⎨ ⎬ = ⎢ ⎥⎨ ⎬ (6.3)
⎪γ yz ⎪ E⎢ 0 0 0 2(1 + ν ) 0 0 ⎥ ⎪τ yz ⎪
⎪γ xz ⎪ ⎢ 0 0 0 0 2(1 + ν ) 0 ⎥ ⎪τ xz ⎪
⎪ ⎪ ⎢ ⎥⎪ ⎪
⎪⎩γ xy ⎪⎭ ⎣⎢ 0 0 0 0 0 2(1 + ν )⎦⎥ ⎪⎩τ xy ⎪⎭
따라서 응력 과 변형률 관계식(Stress-Strain Relationship)은

⎡1 − ν ν ν 0 0 0 ⎤
⎧σ x ⎫ ⎢ ν 1 −ν ν 0 0 0 ⎥ ⎧εx ⎫
⎪σ ⎪ ⎢ ⎥⎪ ⎪
⎪ ⎪ y ⎢ ν ν 1 −ν 0 0 0 ⎥ ⎪ε y ⎪
⎪⎪σ z ⎪⎪ 1 ⎢ 1 − 2ν ⎥⎪ ⎪
⎨ ⎬ = ⎢ 0 0 0 0 ⎥ ⎪⎨ ε z ⎪⎬ (6.4)
⎪τ yz ⎪ ⎥ ⎪γ yz ⎪
E⎢ 2
1 − 2ν
⎪τ xz ⎪ ⎢ 0 0 0 0 0 ⎥ ⎪γ xz ⎪
⎪ ⎪ ⎢ 2 ⎥⎪ ⎪
⎪⎩τ xy ⎪⎭ ⎢ 1 − 2ν ⎥ ⎪⎩γ xy ⎪⎭
⎢⎣ 0 0 0 0 0
2 ⎥⎦
여기서 전단계수G=E/2(1+ν)를 이용하면, 전단 응력과 전단 변형률(Shear Strain)사이의

관계식은 τyz = Gγyz 와 같은 간편한 형태로 쓰일 수 있다.
(2) 평면응력(Plane Stress, NKTP=1)
평면 탄성(Plane Elasticity)문제에 있어서 만약 평면에 수직인 방향으로 응력의 분포가 변

하지 않는다는 가정의 평면구조 이다.
따라서 두께가 일정한 얇은 평판이 x-y 평면 내에 놓여 있다고 하면, 그림[6.1]에 보이는

것과 같이 평판의 모서리에 평면 내에서(In-Plane) 힘이 작용한다면 힘이 작용하는 모서리
근처에서는 힘의 전달로 말미암아 z축 방향 응력 성분들이 평면 내에 작용하는 응력 성분들
에 비해서 극히 작게 된다. 이러한 경우를 평면 응력 구조라 한다.
그림[6.1] 평면 응력 구조 모델
그러므로 평면 응력 상태에 있어서는
σ z = τ xz = τ yz = 0 (6.5)
로 가정하여 이를 식(2.55)에 대입하면,
⎧εx ⎫ ⎡ 1 −ν 0 ⎤ ⎧σ x ⎫
⎪ ⎪ 1 ⎢ ⎪ ⎪
⎨ε y ⎬ = −ν 1 0 ⎥⎥ ⎨σ y ⎬ (6.6)
⎪γ ⎪ E ⎢
⎩ xy ⎭ ⎢⎣ 0 0 2(1 + ν )⎥⎦ ⎪⎩τ xy ⎪⎭
덧붙여 나머지 다른 변형률 성분은
ν ν
εx = − (σ +σ y ) = − (ε +εy) (6.7)
1 −ν
x x
E
γ yz = γ xz = 0 (6.8)
εz는 식(6.7)에서 보는 바와 같이 εx 와 εy의 선형 함수이기 때문에 식(6.6)에는 포함시키

지 않았다.
식(6.6)으로부터 응력과 변형률 관계식은
⎧σ x ⎫ ⎡ ⎤
⎪ ⎪ E ⎢1 ν 0 ⎥ ⎧εx ⎫
⎪ ⎪
⎨σ y ⎬ = ⎢ν 1 0 ⎥ ⎨ε y ⎬ (6.9)
⎪σ ⎪ 1 −ν 2 ⎢ 1 −ν ⎥ ⎪ ⎪
⎩ xy ⎭ ⎢0 0 ⎥ γ xy
⎣ 2 ⎦⎩ ⎭
(3) 평면변형률(Plane Strain, NKTP=2)
그림[6.2]에 보이는 것과 같이 길이 방향으로 단면의 형상이나 하중 상태가 거의 변하지

않는 문제를 평면 변형률 구조라 하며, 이런 구조의 양 끝에서는 어느 정도 떨어진 단면에
서는 길이 방향의 변형률 성분들이 단면 내의 변형률 성분들에 비해서 매우 작게 된다.
x
z
그림[6.2] 평면 변형률 구조 모델
그러므로 평면 변형률 상태에 있어서는 길이 방향을 z 방향으로 택하면,
εz = γ yz = γ xz = 0 (6.10)
로 가정하여 이를 식(6.4)에 대입하면,
⎧σ x ⎫ ⎡ ⎤⎧ ⎫
⎪ ⎪ E ⎢1 − ν ν 0 ⎥ ⎪εx ⎪ (6.11)
⎨σ y ⎬ = ⎢ ν 1 −ν 0 ⎥ ⎨ε y ⎬
⎪τ ⎪ (1 + ν )(1 − 2ν ) ⎢ 1 − 2ν ⎥ ⎪ ⎪
⎩ xy ⎭ ⎢ 0 ⎥ γ xy
2 ⎦⎩ ⎭
0
⎣
나아가 다른 응력 성분은
Eν
σz = (ε + ε ) = ν (σ x + σ y ) (6.12)
(1 + ν )(1 − 2ν ) x y
τ z = τ xz = 0 (6.13)
σz 는 식(6.12)와 같이 σz 와 σz의 선형 함수 이기 때문에 독립 변수가 아닌 관계로 식

(6.11)에는 포함시키지 않았다.
변형률과 응력 관계식을 정리하면, 식(6.11)으로부터
⎧εx ⎫ ⎡1 − ν −ν 0⎤ ⎧σ x ⎫
⎪ ⎪ 1 +ν ⎢ −ν ⎪ ⎪
⎨ε y ⎬ = ⎢ 1 −ν 0 ⎥⎥ ⎨σ y ⎬ (6.14)
⎪γ ⎪ E ⎪τ ⎪
⎩ ⎭
xy
⎢⎣ 0 0 2 ⎥⎦ ⎩ xy ⎭
(4) 평면 요소(Plane Element)
두께가 t인 평면 삼각형 요소가 그림[6.3]에 보이는 바와 같이 x-y 평면 위에 놓여 있다

고 가정하자. 각 절점(Node) 당 자유도가 각각의 좌표축 방향에 대해서 1일 때 강성 행렬
및 절점 하중 등을 유도하면 다음과 같은 과정이 도출 된다.
v3
Theoretical Course 2 1-6 평면 요소

3 u3
v1
1 u1
y
v2
x
u2
2
그림[6.3] 상수 변형률 삼각형 요소
• 요소내의 변위를 국부 좌표계를 사용하여 연속 함수로 가정한다.

• 가정한 변위 함수를 절점의 값으로 표시한다.
x 방향의 변위를 u, y 방향의 변위를 v라 하면,
{u} = [N ] {q} (6.14)
여기서
{u}T = {u , v}, {q} = {u1 , v1 , u 2 , v 2 , u 3 , v}
⎡N 0 N2 0 N3 0 ⎤
[N ] = ⎢ 1
⎣0 N1 0 N2 0 N 3 ⎥⎦
형성함수 Ni는 내분 함수 개념 및 자연좌표 개념을 이용하면,
Ni =
1
(ai + bi + C iy ) (6.15)
2A
여기서
ai = x j y k − x k y j , bi = y i − y k , C i = x k − x j ,
2 A = ai + a j + a k = a1 + a 2 + a
I, j, k는 1, 2, 3의 순환 순열(Cyclic Permutation)이다.
• 변형률과 변위 및 응력과 변형률 관계식이 만족되도록 한다.
변형률과 변위 관계식은 식(6.3)으로부터
∂u ∂v ∂u ∂v
εx = , εy = , γ xy = + (6.16)
∂x ∂y ∂y ∂x
이므로
{ε } = [B ] {q} (6.17)
여기서
{ε }T = {ε x , ε y , γ xy },
⎡ N1, x 0 N 2, x 0 N 3, x 0 ⎤
⎢ ⎥
[B ] = ⎢ N1, y 0 N 2, y 0 N 3, y ⎥
⎢ N1, x N1, x N 2, y N 2, x N 3, y N 3, x ⎥⎦
⎣
⎡b1 b2 b3 ⎤
1 ⎢
= c1 c2 c3 ⎥⎥ (6.18)
2A ⎢
⎢⎣c1 b1 c2 b2 c3 b3 ⎥⎦
위의 행렬에서 Ni,x는 ∂N 1 ∂x 를 의미한다. 행렬[B]는 좌표 x 또는 y를 포함하지 않으

므로, 요소내의 변형률은 식(6.17)에서 상수가 됨을 알 수 있다.{q} 는 절점의 변위로서
상수, 이런 이유로 이 요소는 상수 변형률 삼각형 요소가 불린다.
응력과 변형률 관계식은 식(6.9)과 (6.11)로부터
{σ } = [D ] {ε } (6.19)
여기서 {σ }T = {σ x σ y τ xy } 이고,
⎡e1 ν 0⎤
[D ] =
E ⎢ν e1 0 ⎥⎥ (6.20)
(1 + ν )e 2 ⎢
식(6.20)의 e1, e2, e3는 각각,
⎢⎣ 0 0 e3 ⎥⎦
-. 평면 응력 문제일 때는 e1 = 1, e2 = 1-ν, e3 = (1-ν)/2
-. 평면 변형률 문제일 때는 e1 = 1-ν, e2 = 1-2ν, e3 = (1-2ν)/2
• 요소에 대한 강성 행렬 및 절점 부하를 일 또는 에너지 법칙을 이용하여 결정한다.
요소 내의 변형에너지는
1
[ ] { } 1
[ ] [D ]{ε } dV
T
2 ∫v 2 ∫v
U = B
T
ε dV = ε
=
1
∫ [q ]T [B ]T [D ][B ] {q}dV = 1 {q}T [k ] {q} (6.21)
2 v 2
여기서
[k ] = ∫ [B ] [D ][B ] dV
v
T
= t ∫ [B ] [D ][B ] dA
A
T
= [k ] + [k ]
n s
(6.22)
식(6.18)의 [B] 와 (6.20)의 [D]를 대입하여 강성 행렬을 구하면, 강성 행렬은 식

(6.22)의 오른쪽과 같이 두 부분으로 나눌 수 있다. [kn]는 수직 변형률(Normal Strain)에
의한 값이며 [ks]는 전단 변형률(Shear Strain)에 의한 값이다. 계산 결과는 다음과 같이
주어진다.
⎡e1b12 νb1 c1 e1b1 b2 νeb1 c2 e1b1 b3 νb1c3 ⎤

⎢ ⎥
⎢ e1c12 νb2 c1 e1c1b2 νb3 c1 e1c1c3 ⎥
⎢ νb2 c2 νb2 c3 ⎥
[k ]
n
=
Et
2 ⎢
4 A(1 + ν )e ⎢
e1b22
e1c22
e1b2b3
νb3 c2
⎥
e1c2 c3 ⎥
(6.23)
⎢ e1b32 νb3 c3 ⎥
⎢ 대칭 ⎥
⎣⎢ e1c32 ⎦⎥
⎡ c12 b1c1 c1c 2 b2 c1 c1c 3 b3 c1 ⎤

⎢ ⎥
⎢ b12 b1c 2 b1b2 b1c 3 b1b3 ⎥
⎢ b3 c 2 ⎥
[k ]
2
Et c b2 c 2 c 2 c3 (6.24)
s
= ⎢
2
⎥
8 A (1 + ν ) ⎢ b 2
2 b2 c 3 b2 b3 ⎥
⎢ c32 b3 c 3 ⎥
대칭
⎢ ⎥
⎢⎣ b32 ⎥⎦
3 3 3
1
l12N2
1 1
y pg px
l12
x 2 2 2
(a) 자중에 의한 힘 (b) 일정한 분포 하중 (c) 자연 좌표
그림[6.4] 상수 변형률 요소의 상당 절점 하중
이제 절점 하중에 대해서 생각해 보면, 그림[6.4]에 보이는 자중에 의한 힘과 두 절점

사이의 변에 x방향으로 일정한 분포하중이 작용할 때의 두 가지의 경우에 대해서 상당 절
점 하중에 대해서 고찰하면 다음과 같다.
a) 자중에 의한 힘
단위 부피당 무게를 Pg라 두면, Pg는 –y 방향으로 작용한다.
이 중력에 의한 일은,
⎧ 0 ⎫
∫v{p}{u} dV ∫v{u} {p} dV ∫A {u , v} ⎨
T
W = = = t
T
⎬ dA
⎩ − p g ⎭
⎧ 0 ⎫
= t ∫ {q} [N ]
T T
⎨ ⎬ dA
A
⎩ − p g ⎭
따라서
⎧ 0 ⎫
{f } = t ∫A [N ]
T
⎨ ⎬ dA
⎩− p g ⎭
넓이 좌표의 적분식을 이용하여 계산하면 다음과 같다.
{ f }T = {f 1x f1 y f 2 x f 2 y f 3 x f 3 y } = ⎜ − p g tA ⎟ [0 1 0 1 0 1]
⎛ 1 ⎞
(6.25)
⎝ 3 ⎠
b) 절점 사이의 분포 하중
절점 1과 2사이의 변에 x방향으로 단위 길이당 힘이 px인 분포 하중이 작용할 때, 이
힘에 의한 일은 다음과 같다.
⎧p ⎫
W = ∫ {u, v} ⎨⎩ 0 ⎬⎭dl
x
l12
여기서 l은 길이, l12는 절점 1과 2사이의 길이를 말한다.

a)에서와 같은 과정을 거치면 상당 절점 하중은
T ⎧p ⎫
{f } = ∫ [ N ] ⎨ x ⎬dl
l12
⎩0⎭
[N]를 대입하여 정리하면,
{ f }T = px ∫ [N
l12
1 0 N 2 0 N 3 0 ] dl
그런데 그림[6.4](c)에서 l=l12N2이므로 dl= l12dN2이다. 그리고 절점 1과 2사이의 변에

서는 N3=0이므로 N1+N2+N3=1에서 N1=1-N2이다.
따라서,
1
⎛ N2⎞
∫0 (1 − N 2 )l12 dN 2
1 l12
∫l12
N 1 dl = = l12 ⎜⎜ N 2 − 2 ⎟⎟
⎝ 2 ⎠0
=
2
1
1 N 22 l12
∫l12
N 2 dl = ∫N
0
l dN 2 =
2 12 l12
2
=
2
0
그러므로 상당 절점 하중은,
⎛1 ⎞
{ f }T = ⎜ p x l12 ⎟[1 0 1 0 0 0] (6.26)
⎝2 ⎠
• 국부 좌표로 나타낸 강성행렬 및 절점 하중 등을 대표좌표에 대한 것으로 전환한다.
만약 요소의 국부 좌표축이 대표 좌표축과 평행이 아니면 국부 좌표계의 강성 행렬이나

절점 변위, 절점 하중 등을 대표좌표에 대한 것으로 변환해 주어야 한다. 변환 행렬은 앞장
에서 설명한 과정으로 유도될 수 있다.
그러나 모든 요소가 한 평면 안에 있을 때는 국부 좌표축을 택할 때 대표 좌표축과 평행
하게 택할 수 있으므로, 이 때는 국부 좌표계에 대한 수식이 그대로 대표 좌표계에 대한 수
식이 된다.
Theoretical Course C
예제 1 : 스프링 계로 이루어진 구조물의 유한요소해석
예제 2 : 트러스(Truss) 구조물의 유한요소해석
예제 3 : 빔(Beam) 구조물의 유한요소해석

스프링 계로 이루어진 구조의 유한요소해석
Grid Create Method

Line Create Method
Line Element Generation
Boundary Condition Create Method
Loading Condition Create Method
Element Type Create Method
Element Property Data Create Method
Solver Create Method
Graphic Post-Processing Method

[예제 1] 스프링 계로 이루어진 구조의 유한요소해석(FEA)
◆ 문제의 정의
아래 그림[1]과 같이 용수철(Spring) 상수가 K1, K2, K3인 3개의 용수철로 이루어진 계

(System)의 형상 구조이다. 그리고 세 용수철이 만나는 ″3″ 부분에 집중하중이 작용할 때
각 지점에서의 반력과 ″3″ 지점의 최대 거동 변형량을 수치 계산 과 FEA 결과와 비교해
보자.
: 절점에서의 힘과 변위 사이의 관계식인 직접법을 이용한 직접 강성행렬과 요소에

대한 평형방정식으로 유도하는데 간편한 유도과정이 주어지지만 이는 간단한 문제에 국한된
다.
283.24 N
4
Y
K1 = 150 N/mm
K2 = 200 N/mm
80mm
K3 ③ K3 = 100 N/mm
X
3
0.5KN 500 N
40mm
① ②
K1 -1000 N
1KN K2
1 2 18.78 N
518.78 N
30mm 30mm 25.05 N 691.71 N
그림[1] Spring System 형상 및 힘의 분력 도
◆ 이론 및 FEM NISA와 결과 비교
표[1] 결과 비교
″3″ 지점의 최대 거동
각 지점 Rx(N) Ry(N)
변형량(mm)
반력(N)
이론 NISA 이론 NISA 이 론 NISA
1 18.78 1.8783E+01 25.05 2.5044E+01
X=3.429 X=3.4287
2 -518.78 -5.1878E+02 691.71 6.9171E+02
Y=-2.832 Y=-2.83244
4 0.0 0.0 283.24 2.8324E+02
1. 강성 행렬을 이용한 유도과정

283.24 N
☞ 연속체 또는 계를 분할한다. 4
Y
그림과 같이 각각의 용수철을 하나의 K1 = 150 N/mm

요소로 간주하여 세 개의 요소로 이루어 K2 = 200 N/mm K3 ③
X
K3 = 100 N/mm
진 계로 간주하고, 각 요소와 절점의 번호
및 강성을 고려한다. 3
500 N
① ②
• ① : 절점 1과 3을 연결하는 요소 1KN K2
K1
• ② : 절점 2와 3을 연결하는 요소
18.78 N 1 2
• ③ : 절점 3과 4를 연결하는 요소 518.78 N
• K1, K2, K3 : 각 요소 ①,②,③의 강성 25.05 N 691.71 N
☞ 요소에 대한 수식을 정립한다.
각 요소에 대한 강성행렬을 각각의 요소 ①,② 와 ③에 적용한다.
¨ 요소 ①
• 절점 : 1과 3으로 이루어져 있다.
• C= Cosθ = 3/5, S = Sinθ = 4/5
• 강성(K1) =150 N/mm 150 × 3 ( 5 ) = 54
2
⎡C 2 CS −C2 − CS ⎤ ⎡54 72 − 54 − 72 ⎤
⎢ 2⎥ ⎢
S2 − CS −S ⎥ 96 − 72 − 96 ⎥⎥
[K ]1 = k1 ⎢
⎢
= ⎢
C2 CS ⎥ ⎢ 54 72 ⎥
⎢ 대칭 2 ⎥ ⎢ 대칭 ⎥
⎣⎢ S ⎦⎥ ⎣ 96 ⎦
¨ 요소 ②
• C= Cosθ = -3/5, S = Sinθ = 4/5
• 강성(K1) =200 N/mm
⎡C 2 CS −C2 − CS ⎤ ⎡72 − 96 − 72 96 ⎤
⎢ 2⎥ ⎢
S2 − CS −S ⎥ 128 96 − 128 ⎥⎥
[K ]2 = k 2 ⎢⎢ = ⎢
C2 CS ⎥ ⎢ 72 − 96 ⎥
⎢ 대칭 2 ⎥ ⎢ 대칭 ⎥
¨ 요소 ③ ⎣⎢ S ⎦⎥ ⎣ 128 ⎦

• C= Cosθ = 0, S = Sinθ = 1
• 강성(K1) =100 N/mm
⎡C 2 CS −C2 − CS ⎤ ⎡0 0 0 0 ⎤
⎢ 2⎥ ⎢ 100 0 100 ⎥
S2 − CS −S ⎥
[K ]3 = k3 ⎢
⎢
= ⎢ ⎥
C2 CS ⎥ ⎢ 0 0 ⎥
⎢ 대칭 ⎥ ⎢ 대칭 ⎥
⎣⎢ S 2 ⎦⎥ ⎣ 100 ⎦
☞ 각 요소에 대한 식을 결합하여 전체 계에 대한 식을 유도하고 경계조건을 부여한다.
¨ 각 절점에서의 변위
• 절점 1의 변위 : Q1X, Q1Y
• 절점 2의 변위 : Q2X, Q2Y
¨ 힘과 변위를 강성행렬에 결합하여 전체 계에 대한 평형방정식을 유도하면 다음과 같다.
1 2 3 4
⎡54 72 0 0 − 54 − 72 0 0 ⎤ ⎧ Q1 X ⎫ ⎧ F1 X ⎫
1⎢
96 0 0 − 72 − 96 0 ⎥ ⎪
0 ⎥ ⎪ Q1Y ⎪ ⎪ ⎪F ⎪
⎢ ⎪ 1Y ⎪
⎢ 72 − 96 − 72 96 0 0 ⎥ ⎪Q2 X ⎪ ⎪ F2 X ⎪
2⎢ ⎥⎪ ⎪ ⎪ ⎪
⎢ 128 96 − 128 0 0 ⎥ ⎪ Q 2Y ⎪ ⎪ F2Y ⎪
⎨ ⎬ = ⎨ ⎬
⎢ 126 − 24 0 0 ⎥ ⎪Q3 X ⎪ F3 X ⎪
3⎢ ⎥ ⎪
⎢ 324 0 100 ⎥ ⎪ Q3Y ⎪ ⎪ F3Y ⎪
⎢ ⎪ ⎪ ⎪ ⎪
대칭 0 0 ⎥ ⎪Q4 X ⎪ ⎪ F4 X ⎪
4⎢ ⎥
⎣⎢ 100 ⎦⎥ ⎪⎩ Q 4Y ⎪⎭ ⎪F ⎪
⎩ 4Y ⎭
☞ 주어진 힘과 경계조건을 부여하여 2자유도인 수식으로 정리한다.
¨ 외력 힘의 분포 및 경계 조건
• F3X = 500 N
• F3Y = -1000 N
• Q1X = Q1Y = Q2X = Q2Y = Q4X = Q4Y = 0
⎡ 126 − 24 ⎤ ⎧Q3 X ⎫ ⎧ 500 ⎫

⎢ − 24 324 ⎥ ⎨ Q ⎬ = ⎨− 1000 ⎬
⎣ ⎦ ⎩ 3Y ⎭ ⎩ ⎭
☞ 절점 ″3″에서 각 방향에 대한 변위를 구하여 보자.
¨ 각 방향에 대한 변위는 Cramer’s Rule이나 Gauss Elimination Method를 이용하여 Q3x,

Q3의 변위를 구하면 다음과 같다.
⎧Q3 X ⎫ ⎧ 3.429 ⎫
⎨ ⎬= ⎨ ⎬ mm
⎩ Q3Y ⎭ ⎩− 2.832 ⎭
☞ 각 절점에서의 반력을 강성 행렬을 이용하여 구한다.
¨ 앞서 언급하였듯이 요소에 대한 평형 방정식을 이용하여 각 절점에서의 반력을 구할 수

도 있지만 또한 상기에서 구한 변위량을 이용하여 전체 계에 대한 평형방정식으로 다음과
같이 유도할 수 있다.
1 2 3 4
⎡54 72 0 0 − 54 − 72 0 0 ⎤⎧ 0 ⎫ ⎧ F1 X ⎫
1⎢
96 0 0 − 72 − 96 0 ⎥ ⎪
0 ⎥⎪ 0 ⎪ ⎪ ⎪ F ⎪
⎢ ⎪ 1Y ⎪
⎢ 72 − 96 − 72 96 0 0 ⎥⎪ 0 ⎪ ⎪ F2 X ⎪
2⎢ ⎥⎪ ⎪ ⎪ ⎪
⎢ 128 96 − 128 0 0 ⎥⎪ 0 ⎪ ⎪ F2Y ⎪
⎢ ⎨ ⎬ = ⎨ ⎬ Newtons
126 − 24 0 0 ⎥ ⎪ 3.429 ⎪ ⎪ 500 ⎪
3⎢ ⎥
⎢ 324 0 100 ⎥ ⎪− 2.832 ⎪ ⎪− 1000 ⎪
⎢ 대칭 ⎪ ⎪ ⎪ ⎪
0 0 ⎥⎪ 0 ⎪ ⎪ 0 ⎪
4⎢ ⎥
⎢⎣ 100 ⎥⎦ ⎪⎩ 0 ⎪⎭ ⎪ F ⎪
⎩ 4Y ⎭
⎧ F1 X ⎫ ⎧ 18.78 ⎫
⎪F ⎪ ⎪ 25.05 ⎪
⎪ 1Y ⎪ ⎪ ⎪
⎪ F2 X ⎪ ⎪− 518.78⎪
⎪ ⎪ ⎪ ⎪
⎪ F2Y ⎪ ⎪ 691.71 ⎪
⎨ ⎬ = ⎨ ⎬ Newtons
⎪ F3 X ⎪ ⎪ 500.00 ⎪
⎪ F3Y ⎪ ⎪ − 1000 .0 ⎪
⎪ ⎪ ⎪ ⎪
⎪ F4 X ⎪ ⎪ 0 .0 ⎪
⎪F ⎪ ⎪ 283.24 ⎪
⎩ 4Y ⎭ ⎩ ⎭
2. NISA Software를 이용한 해석 과정
☞ 기하학적형상의 구현단계는 다음과 같은 단계를 따른다.
사용자의 편리성에 따라 Menu 방식이나 Command 방식을 선택할 수 있으나 다음의

단계는 Menu 과정으로 접근 하였다(대/소문자 구분 없음).
Ä 좌표축을 윗면(Top)으로 설정한다.

• ⇒ : 우축하단부에 좌표 표시
• X-Y(Top)가 기본 좌표계이나 상이한 경우에는 이를 정의한다.
Ä Grid를 생성한다(G1 to G4)
Menu를 이용하여 Gird를 생성한다.

• GEOMETRY > GRID > CREATE > XYZ Coord
► XYZ-Data[0/0/0] > OK ---- G1 생성
► XYZ-Data[60/0/0] > OK ---- G2 생성
► XYZ-Data[30/40/0] > OK ---- G3 생성
► XYZ-Data[30/120/0] > OK ---- G4 생성
► Cancel ---- Grid 생성이 완료되면 Grid 생성화면 지움
GRIDS
Ä Line을 생성한다(L1 to L3)

• GEOMETRY > LINE > CREATE > 2 GRIDS
► Grid IDs[1,3] > OK ---- L1 생성
► Grid IDs[2,3] > OK ---- L2 생성
► Grid IDs[3,4] > OK ---- L3 생성
► Cancel
LINES
상기와 다른 Line 생성 방법은 다음과 같다(참조).

► Grid IDs[1/2/3,3/3/4] > OK -- L1to L3 생성
► Cancel
Grid ″1″ 은 상대 Grid ″3″으로 연결되어 선이 생성
Grid ″2″ 은 상대 Grid ″3″으로 연결되어 선이 생성
Grid ″3″ 은 상대 Grid ″4″으로 연결되어 선이 생성 [1 / 2 / 3, 3 / 3 / 4]
• GEOMETRY > LINE > CREATE > 2 GRIDS [대응되는 GRID]

► 마우스로 두 Grid를 선택하여 선을 연결한다.
⇒ 좌측마우스로 Grid ″1″ ⇒ Grid ″3″ -- L1

⇒ Grid ″2″ ⇒ Grid ″3″ -- L2
⇒ Grid ″3″ ⇒ Grid ″4″ -- L3 L3
⇒ Cancel하기 위해 우측마우스 버턴 선택
L1 L2
Ä Spring 요소를 생성한다(E1 to E3)
Command를 이용하여 요소를 생성한다.

• FEG,BAR,1,,,1,18/1/1/1 > Return ⇒ E1
• FEG,BAR,2,,,1,18/1/1/2 > Return ⇒ E2
• FEG,BAR,3,,,1,18/1/1/3 > Return ⇒ E3
FEG : Finite Element Generation(요소생성 명령어)

BAR : Spring, Beam, Pipe, Cable등의 Line 요소
″1″ : Line Number
″,,,″ : 자동으로 생성되는 요소 및 절점 번호
″1″ : 앞에서 지정한 Line에 요소 1개 생성
″18/1/1/1″ : 2차원 스프링(NKTP=18)/요소차원
/재질번호/단면특성(스프링상수)
ELEMENT
NODES
Ä 구조물 일체화(절점 Merge) 한다.
Command를 이용하여 절점을 일체화 한다.

• NODE,MER,ALL > Return > Y
Menu 방식을 이용하여 절점을 일체화 한다(참조)

• FEG > NODES > NODE MANAGEMENT > MERGE NODES > NODE Ids > METHOD(ALL)
> OK > OK(MAIN MENU > YES ALL> CANCEL : 중복되어 있는 절점들은 가장 낮은 절점만 남는다.
NODE
MERGE
Ä 경계조건을 정의한다.
메뉴를 이용하여 생성하는데 계속 좌측 마우스 버턴을 선택한다.

• FE GENERATION > STRUCTURAL > DISPLACEMENT > ADD AN NODES
> BC Data > Click ″Ux″ > Click ″Uy” > OK
> Node Ids > Input Ids[1/3/6] > OK > OK > Cancel
• BC DATA > Ux = Uy = 0(Constrain,고정)
구조물의 좌,우측 및 상단 부분을 완전고정을

하기위해 Ux,Uy 부분을 선택한다.
2차원 Spring요소는 자유도(DOF)가 2 방향성만
있기 때문에 Uz,Rx,Ry,Rz는 선택하지 않는다.
여기서 ″0″이라는 숫자가 나타나는데 ″0″은 곳
구조물 선택 부분이 변형이 없다는 의미이고 인
위적으로 숫자를 입력하면 강제변위의 하중으로
정의된다는데 주의한다.
[경계조건 부분의 기호 설명]

Ux,Uy,Uz : 변형에 대한 구속 DOF
Ry
Rx,Ry,Rz : Rotation에 대한 구속 DOF
흑색으로 표시된 DOF가 구속이 되었다는 의미 Rz Rx
이고 나머지 흰색 부분은 구속(고정)이 없다는
의미 이다. 단, 각각의 요소마다 정해진 DOF가 Uz
Ux
있기 때문에 본 스프링 요소에서는 Ux, Uy를 제
외한 DOF는 의미가 없다. 예, 다른 DOF를 구 Uy
속하여도 의미는 없고 경고문이 주어진다.
Ä 외력(하중)조건을 정의한다.
• FE GENERATION > STRUCTURAL > FORCE.MOMENT> ADD AN NODES
• BC DATA > Ux[500.0] > Uy[-1000] > OK
Fx,Fy,Fz : 각각 X,Y,Z 축 방향의 집중하중

Mx,My,Mz : 각각 X,Y,Z축 방향의 모멘트
Fx : 좌표상에서 +X축 방향 ► 500
Fy : 좌표상에서 –Y축 방향 ► -1000
• NODE IDS > Input Ids[2] > OK > OK

> Cancel
집중하중(외력)조건은 절점에 적용되기 때문에

우측의 그림에서 보여주듯이 선택의 방법은 직접
절점의 정보를 입력하여 선택도 가능하고 여러가
지 마우스의 선택에 의한 방법도 가능하다.
Ä 단면(재질)특성을 정의한다.
메뉴를 이용하여 생성한다.

• ANALYSIS/STATIC > PROPERTY
> ADD(화살표 부분선택)
> Spar,Spring and Cable(2D trans. S)
> PROP(선택)
> Axial Spring Constant(150.0)
> OK ► ″1번요소 스프링상수″
► PROPERTY ID가 ″2″로 바뀜
> PROP(선택)
> PROP(선택)
> OK
스프링은 재질(단면)특성에는 스프링

상수만 입력이 된다.
입력 방법은 앞의 요소 생성에서 단면특성
번호를 요소 별 ″1/2/3″으로 정의 하였기
때문에 각 스프링 상수를 달리 입력한다.
반복 수행
Ä 해석 및 결과 화일을 정의한다.

• ANALYSIS/STATIC > EXCUTIVE
기본이 정적해석이기 때문에 달리 해석범

위를 설정하지 않고 바로 실행의 항으로
접근한다.
단, 동, 열, 유체 해석 등은 달리 설정해
야 한다.
Save/Access NISA > File 26, 27 부분선택

Enter Filename > spring1
> 확인
Ä 입력 NISA INPUT 화일을 생성한다.
• FILE > SAVE AS > 파일 이름[spring1]

> 파일 형식(*.nis) 선택 > 저장(S)
Ä Solver를 실행한다.
각각의 정보들이 Tree 방식으로 저장이 되는

좌측 메뉴 화면의 아래 부분에 있는 화살표를
선택하여 Executables이 선택이 되도록 설정한
다.
⇒
설정이 되면 이 부분을 마우
스로 선택한다.
그러면 우측의 그림의 화면과 같이 여러 가지의 Solver의 기능

이 선택이 된다.
> 마우스로 NISA2를 두 번 연속 선택하여 실행한다.
> Type in NISA data file name > spring1.nis > Enter
Solver 실행을 위한 아래와 같은 창이 선택이 되면 앞서
선정한 파일명을 입력(spring1.nis)하고 Enter를 한다.
> Type in output file name > spring1.out > Enter [NISA SOLVER]
각종 결과를 Text로 출력 할 수 있는 화일명을 입력한다.
enter
enter
Ä NISA Input 화일 확인(spring1.nis, 참조)
ANALYSIS = STATIC
SOLV = FRON
FILE = spring1 ** File Name
SAVE = 26,27
*ELTYPE
1, 18, 1
*RCTABLE
1, 5,1, 0 ** Property ID(1), 입력 항(5)
1.5000000E+02, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 스프링상수(150)
0.0000000E+00, ** Property ID(2)
2, 5,1, 0 스프링상수(200)
2.0000000E+02, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00, ** Property ID(3)
0.0000000E+00,
스프링상수(100)
3, 5,1, 0
1.0000000E+02, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00,
0.0000000E+00,
*NODES ** 각 절점의 좌표값
1,,,, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0
2,,,, 3.00000E+01, 4.00000E+01, 0.00000E+00, 0
3,,,, 6.00000E+01, 0.00000E+00, 0.00000E+00, 0
6,,,, 3.00000E+01, 1.20000E+02, 0.00000E+00, 0
*ELEMENTS ** 각 요소의 절점 연결성
1, 1, 1, 1, 0
1, 2,
2, 1, 1, 2, 0
3, 2,
3, 1, 1, 3, 0
2, 6,
*LDCASE, ID= 1 ** 하중 제한 조건
0,1,1,0,-1,2,0
*SPDISP ** 경계조건
** SPDISP, SET = 1 각 절점의 경계조건 구속
1,UX , 0.00000E+00,,,,,,,, 0
1,UY , 0.00000E+00,,,,,,,, 0
3,UX , 0.00000E+00,,,,,,,, 0
3,UY , 0.00000E+00,,,,,,,, 0
6,UX , 0.00000E+00,,,,,,,, 0
6,UY , 0.00000E+00,,,,,,,, 0
** 집중하중 조건
*CFORCE
** CFORCE, SET = 1
2,FX, 5.00000E+02,,, 0, 0
2,FY,-1.00000E+03,,, 0, 0
*ENDDATA
Ä 각 절점에서 반력 및 변형량을 확인한다(spring1.out, 참조).
***** REACTION FORCES AND MOMENTS AT NODES *****
LOAD CASE ID NO. 1
NODE FX FY FZ MX MY MZ
1 1.87835E+01 2.50447E+01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00

3 -5.18784E+02 6.91711E+02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
6 0.00000E+00 2.83244E+02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
SUMMATION OF REACTION FORCES IN GLOBAL DIRECTIONS

-------------------------------------------------------------------------------------
FX FY FZ
-5.000000E+02 1.000000E+03 0.000000E+00
*** EMRC NISA *** -- Version 10.0 (07/09/00-128MEG) LOAD CASE ID NO. 1 FEB/27/2001
12:45:19
****** DISPLACEMENT SOLUTION ******
LOAD CASE ID NO. 1
NODE UX UY UZ ROTX ROTY ROTZ
1 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00

2 3.42874E+00 -2.83244E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
3 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
6 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
► 각 절점에서의 반력(Reaction Force) 및 변형량(Displacement)을 확인하기 전에 주의해

할 부분은 구조물에 제하한 Total 하중이 경계(구속)조건이 주어진 부분의 반력의 합과 동
일한지는 우선 확이 해야 한다. 그 이유는 하중의 조건이 정확하게 주어졌는지 확인의 절
차이기 때문이다(Fx=-500N, Fy=1000.0N). 상반된 부호는 제하된 하중의 반력 부분이기
때문에 반대적으로 반력이 발생한다.
► 반력 및 변형 결과 분석
절점 번호(Nodes) 반 력(Reaction Force, N) 변형량(mm)

(FE MODEL에서) Rx Ry X Y
1 1.87835E+01 2.50447E+01 0.0 0.0
2 -5.18784E+02 6.91711E+02 3.4287 -2.83244
6 0.0 2.83244E+02 0.0 0.0
Ä 결과 Post File를 Loading한다.
• FILE > NEW(새로운 결과 File Loading)

• FILE > OPEN > 파일 형식(*.dat)
> spring126 > 열기 > spring127
> 열기
Ä 그래픽으로 변형을 확인한다
• POST-PROCESSING/STATIC > DEFORMATIONS > PICK DEFORMATIONS

트러스(Truss) 구조물의 유한요소 해석
Grid Create Method

Line Create Method

[예제 2] 트러스(Truss) 구조물의 유한요소 해석
아래 그림[1]과 같이 각각의 단면적이 상이한 네 개의 Truss구조가 절점″5″에서 외력을

받고 있는 형상으로, ″5″ 지점의 최대 거동 변형량 및 각 지점에서의 반력을 수치 계산에
의한 결과와 비교해 보자.
: 각각의 Truss는 자유도가 3개인 요소로 모멘트의 항이 없다.
5
5N A1 = 49mm2
2
10 N A2 = 24.5mm
2
A3 = 98mm
A4 = 49mm2
④ 7m
Z 15 N 탄성계수 : 210000 MPa
①
4
6m
② 1 Y
③
Z
4m
6m
Y
2
3
X
X
그림[1] Truss 구조물의 형상
각 지점 Rx(N) Ry(N) Rz(N)
의 반력(N) 이론 NISA 이론 NISA 이론 NISA
1 -0.926 -9.2592E-01 -1.389 -1.38889E+00 -2.778 -2.77778E+00
2 # 1.85185E-02 # -1.80041E-03 # -2.70062E-03
3 # -5.7407E+00 # -8.61111E+00 # 1.72222E+01
4 # -1.5740E+00 # 2.36111E+00 # -4.72222E+00
변형량(mm) X 방향 Y 방향 Z 방향
절점″5″ 0.01852 1.8518E-02 -0.0018 -1.8004E-03 -0.0027 -2.70062E-03
**부재 ″1″에서의 절대값의 축력(N) : NISA=3.24074, 이론=3.241
☞ 연속체 또는 계를 분할한다.
앞의 그림에서와 같이 각각의 봉은 하나의 요소와 2개의 절점으로 이루어진 계로 간주한

다. 그리고 각 요소의 길이는 7m이다.
• ① : 절점 1과 5을 연결하는 요소
• ② : 절점 2와 5을 연결하는 요소
• ③ : 절점 3과 5를 연결하는 요소
• ④ : 절점 4와 5를 연결하는 요소
• A1, A2, A3, A4 : 각 요소 ①, ②, ③, ④의 단면적
각 요소에 대한 강성행렬을 각각의 요소 ①, ②,③, ④에 적용한다.
¨ 요소 ①
• CX= CosθX = 2/7, CY= CosθY = 3/7, CZ= CosθZ = 6/7
⎡[K ] [K 15 ]⎤
[K ]1 = ⎢ 11
⎣[K 51 ] [K 55 ]⎥⎦
여기서
⎡ 4 6 12 ⎤
= (1470 ) ⎢ 9 18 ⎥⎥
1 ⎢
[K 11 ] = [K 55 ] = − [K 15 ] = − [K 51 ]
49 대칭
⎣⎢ 36 ⎥⎦
¨ 요소 ②
• 절점 : 2와 5으로 이루어져 있다.
• CX= CosθX = -2/7, CY= CosθY = 3/7, CZ= CosθZ = 6/7
⎡[K ] [K 25 ]⎤
[K ]2 = ⎢ 22
⎣[K 52 ] [K 55 ]⎥⎦
여기서
⎡4 − 6 − 12 ⎤
1 ⎢
[K 22 ] = [K 55 ] = − [K 25 ] = − [K 52 ] = (735 ) ⎢ 9 18 ⎥⎥
49 대칭
⎣⎢ 36 ⎦⎥
¨ 요소 ③
• CX= CosθX = -2/7, CY= CosθY = -3/7, CZ= CosθZ = 6/7
⎡[K ] [K 35 ]⎤
[K ]3 = ⎢ 33
⎣[K 53 ] [K 55 ]⎥⎦
여기서
⎡ 4 6 − 12 ⎤
[K 33 ] = [K 55 ] = − [K 35 ] = − [K 53 ] = (2940 ) ⎢ 9 − 18⎥⎥
1 ⎢
49 대칭
⎣⎢ 36 ⎥⎦
¨ 요소 ④
• 절점 4와 5으로 이루어져 있다.
• CX= CosθX = 2/7, CY= CosθY = -3/7, CZ= CosθZ = 6/7
⎡[K ] [K 45 ]⎤
[K ]4 = ⎢ 44
⎣[K 54 ] [K 55 ]⎥⎦
여기서
⎡4 − 6 12 ⎤
[K 44 ] = [K 55 ] = − [K 45 ] = − [K 54 ] = (1470 ) ⎢ 9 − 18⎥⎥
1 ⎢
49 대칭
⎢⎣ 36 ⎥⎦
☞ 각 요소에 대한 식을 결합하여 전체 계에 대한 식을 유도하고 경계조건을 부여

한다.
{Qi }T = {QiX , QiY , QiZ }

i = 1, 2 , 3, 4 , 5
{Fi }T = {FiX , FiY , FiZ }
라 정의 되면,
⎡[K 11 ] 0 0 0 [K 15 ] ⎤ ⎧ 0 ⎫ ⎧[F1 ]⎫
⎢ ⎥
⎢ [K 22 ] 0 0 [K 25 ] ⎥ ⎪ 0 ⎪ ⎪F ⎪
⎪⎪ ⎪⎪ ⎪⎪ 2 ⎪⎪
⎢ [K 33 ] 0 [K 35 ] ⎥ ⎨ 0 ⎬ = ⎨ F3 ⎬
⎢
⎢ 대칭
[K 44 ] [K 45 ] ⎥⎥ ⎪ 0 ⎪ ⎪F ⎪
⎪ ⎪ ⎪ 4⎪
⎢
⎣ [K 55 ]g ⎥⎦ ⎪⎩{Q5 }⎪⎭ ⎪⎩ F5 ⎪⎭
여기서
{F5 }T = {10, 5, − 15}, [K 55 ]1 + [K 55 ]2 + [K 55 ]3 + [K 55 ]4
결국, 3자유도인 수식이 얻어진다.
⎡540 270 − 180 ⎤ ⎧Q5 X ⎫ ⎧ 10 ⎫

⎢ ⎪ ⎪ ⎪ ⎪
⎢ 1215 − 810⎥⎥ ⎨ Q5Y ⎬ = ⎨ 5 ⎬
⎢⎣ 4860 ⎥⎦ ⎪⎩Q5 Z ⎪⎭ ⎪⎩− 15⎪⎭
⎧Q5 X ⎫ ⎧ 18.52 ⎫
⎪ ⎪ ⎪ ⎪ −3
⎨ Q 5Y ⎬ = ⎨ − 1.80 ⎬ × 10 mm
⎪Q ⎪ ⎪− 2.70 ⎪
⎩ 5Z ⎭ ⎩ ⎭
☞ 절점 ″1″에 작용하는 반력을 구해보자
{F1} = [K15 ]{Q5 } 이므로,
⎧ F1 X ⎫ ⎡120 180 360 ⎤ ⎧ 18.52 ⎫ ⎧0.926 ⎫

⎪ ⎪ ⎢ ⎥ ⎪ ⎪ −3 ⎪ ⎪
⎨ F1Y ⎬ = ⎢ 270 540 ⎥ ⎨ − 1.80 ⎬ × 10 = ⎨1.389 ⎬ Newtons
⎪F ⎪ ⎢⎣ 1080 ⎥⎦ ⎪⎩− 2.70 ⎪⎭ ⎪2.778 ⎪
⎩ 1Z ⎭ ⎩ ⎭
요소 ①에 대해서 구해보면, 변환행렬 [T]는,
1 ⎡ 2 3 6 0 0 0⎤
[T ]1 =
7 ⎢⎣0 0 0 2 3 6⎥⎦
요소 ①에 대한 대표 좌표계에서의 변위는,
{Q}1T = {0, }
0, 0, 18.52 × 10 −3 , − 1.80 × 10 −3 , − 2.70 × 10 −3 mm
이므로,
⎧q ⎫ ⎧ 0 ⎫
{q}1 = ⎨ 1⎬ = [T ]1 {Q1 } = ⎨ −3
⎬ × 10 mm
⎩q 5 ⎭ ⎩2.205 ⎭
따라서,
⎧ f1 ⎫ EA ⎡ 1 − 1⎤ ⎧ 0 ⎫ ⎧− 3.241⎫
⎨ ⎬ = [k ]1 {q}1 = ⎢ ⎥ ⎨ −3 ⎬
= ⎨ ⎬ Newtons
⎩ f5 ⎭ l ⎣ − 1 1 ⎦ ⎩2.205 × 10 ⎭ ⎩ 3.241 ⎭

Ä 좌표축을 윗면(Isotropic)으로 설정한다.


► XYZ-Data[0/0/0] > OK ---- G1 생성
► XYZ-Data[4000/0/0] > OK ---- G2 생성
► XYZ-Data[0/6000/0] > OK ---- G3 생성
► XYZ-Data[4000/6000/0] > OK ---- G4 생성
► XYZ-Data[2000/3000/6000] > OK ---- G5 생성
► CANCEL ---- Grid 생성이 완료되면 Grid 생성화면 지움
GRIDS

► Grid IDs[1/2/3/4,5/5/5/5] > OK ---- L1 to L4 생성
► Cancel
LINES
Ä 3D Bar 요소를 생성한다(E1 to E4)
• FEM > MESH > FEG OPTION > BAR/BEAM/SPRING

► No.of Elements[1]
► > 3D Spar > OK
► Line Ids > Input Ids[1] > OK > OK ⇒ E1 생성
► > 3D Spar > Property ID[2] > OK
► Line Ids > Input Ids[3] > OK > OK > Cancel
⇒ E4 생성
ELEMENT
NODES
[3D Bar 요소생성] [3D Bar Property 별 색상 표현]


Menu 방식을 이용하여 절점을 일체화 한다(참조).
> OK > OK(MAIN MENU > YES ALL > CANCEL
NODE
MERGE

> BC Data > Click ″Ux″ > Click ″Uy” > Click ″Uz″ > OK
> Node Ids > Input Ids[1/3/5/7] > OK > OK > Cancel
• BC DATA > Ux=Uy=Uz=0(Constrain,고정)
구조물의 하단부분을 완전고정을 하기 위해

Ux,Uy,Uz 부분을 선택한다.
3차원 Spar요소는 자유도(DOF)가 3 방향성만
있기 때문에 Rx,Rx,Rz는 선택하지 않는다. 여기
서 선택시 ″0″이라는 숫자가 나타나는데 ″0″은
곳 구조물 선택 부분이 변형이 없다는 의미이고
인위적으로 숫자를 입력하면 강제변위의 하중으
로 정의된다는데 주의한다.
경계조건 부분의 기호 설명
Ry
Ux
있기 때문에 본 Spar 요소에서는 Ux,Uy,Uz를
제외한 DOF는 의미가 없다. 예, 다른 DOF를 Uy
구속하여도 의미는 없고 경고문이 주어진다.
Ä 외력(하중)조건을 정의한다.
• FE GENERATION > STRUCTURAL > FORCE. MOMENT> ADD AN NODES
• BC DATA > Ux[10] > Uy[5] > Uz[-15] OK

Fx : 좌표상에서 +X축 방향 ► 10
Fy : 좌표상에서 +Y축 방향 ► 5
Fz : 좌표상에서 -Z축 방향 ► -15

> Cancel

Ä 단면특성을 정의한다.

• ANALYSIS/STATIC > PROPERTY
> Spar,Spring and Cable(3D Spar)
> PROP(선택)
> First Node(49.0)
> Second Node(49.0)
> OK ► ″1번요소 단면특성″
> PROP(선택)
> First Node(24.5)
> Second Node(24.5)
> PROP(선택)
> First Node(98.0)
> Second Node(98.0)
> PROP(선택)
> First Node(49.0)
> Second Node(49.0)
> OK
반복 수행
Ä 재질특성을 정의한다.

• ANALYSIS/STATIC > MATERIAL
> Ex = 210000
> Nuxy = 0.3
> OK
재질특성은 기본이 등방성(Isotropic)이기

때문에 Ex, Nuxy 두개만 입력을 하면 자
동적으로 나머지 부분들은 입력이 된다.
단, 정방성(Orthotropic) 및 이방성
(Anisopropic)의 재질은 입력 방법이 상
이함에 유의해야 한다.
Ä Print Control를 정의한다.

각 부재(Spring, Bar, Beam, Cable등 Line요소)에 대한 축력(Axial Force)를 Text Output으
로 확인하기 위해서는 기본이 ″Off″로 되어 있는 부분을 ″On″으로 정의한다.
• ANALYSIS/STATIC > PRINT CONTROL

> ADD > SUMMARY OF LINE ELMENT FORCES(SLFO) > PRINT ALL > OK
Ä 해석 및 결과화일을 정의한다.


접근한다.
야 한다.

Enter Filename > 3dspar1
> 확인
• FILE > SAVE AS > 파일 이름[3dspar1]


다.
⇒

> Type in NISA data file name > 3dspar1.nis > Enter
선정한 파일명을 입력(3dspar1.nis)하고 Enter를 한다.
> Type in output file name > 3dspar1.out > Enter
각종 결과를 Text로 출력 할 수 있는 파일명을 입력한다.
enter
enter
Ä NISA Input 파일 확인(3dspar1.nis)
**EXECUTIVE data deck

ANALYSIS = STATIC
SOLV = FRON
FILE = 3dspar1
SAVE = 26,27,29
*ELTYPE
1, 14, 1 ** 14 : Element Type(3D Spar)
*RCTABLE
1, 3,1, 0
4.9000000E+01, 4.9000000E+01, 0.0000000E+00, ** 각 3D Spar Element에 대한 단면
2, 3,1, 0 특성
2.4500000E+01, 2.4500000E+01, 0.0000000E+00,
3, 3,1, 0
9.8000000E+01, 9.8000000E+01, 0.0000000E+00,
4, 3,1, 0
4.9000000E+01, 4.9000000E+01, 0.0000000E+00,
*NODES
** 각 절점에서의 3차원 좌표값
1,,,, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0
2,,,, 2.00000E+03, 3.00000E+03, 6.00000E+03, 0
3,,,, 4.00000E+03, 0.00000E+00, 0.00000E+00, 0
5,,,, 4.00000E+03, 6.00000E+03, 0.00000E+00, 0
7,,,, 0.00000E+00, 6.00000E+03, 0.00000E+00, 0
*ELEMENTS ** 각 요소의 절점과 연결성
1, 1, 1, 1, 0
1, 2,
2, 1, 1, 2, 0
3, 2,
3, 1, 1, 3, 0
5, 2,
4, 1, 1, 4, 0
7, 2,
*MATERIAL ** 재질특성
EX , 1,0, 2.10000E+05,(탄성계수)
NUXY, 1,0, 3.00000E-01,(포이송비) ** 하중제하 조건
*LDCASE, ID= 1
0,1,1,0,-1,2,0
*SPDISP
** SPDISP, SET = 1 ** 경계조건
1,UX , 0.00000E+00,,,,,,,, 0
1,UY , 0.00000E+00,,,,,,,, 0
1,UZ , 0.00000E+00,,,,,,,, 0
3,UX , 0.00000E+00,,,,,,,, 0
3,UY , 0.00000E+00,,,,,,,, 0
3,UZ , 0.00000E+00,,,,,,,, 0
5,UX , 0.00000E+00,,,,,,,, 0
5,UY , 0.00000E+00,,,,,,,, 0
5,UZ , 0.00000E+00,,,,,,,, 0
7,UX , 0.00000E+00,,,,,,,, 0
7,UY , 0.00000E+00,,,,,,,, 0
7,UZ , 0.00000E+00,,,,,,,, 0
*CFORCE ** 집중하중 조건
** CFORCE, SET = 1
2,FX, 1.00000E+01,,, 0, 0
2,FY, 5.00000E+00,,, 0, 0
2,FZ,-1.50000E+01,,, 0, 0
*PRINTCNTL ** Print Control 조건
SLFO,0 Line Element Axial Force Printing
*ENDDATA
Ä 각 절점에서 반력 및 변형량을 확인한다(3dspar1.out).
LOAD CASE ID NO. 1
1 -9.25926E-01 -1.38889E+00 -2.77778E+00 0.00000E+00 0.00000E+00 0.00000E+00

3 -1.75926E+00 2.63889E+00 5.27778E+00 0.00000E+00 0.00000E+00 0.00000E+00
5 -5.74074E+00 -8.61111E+00 1.72222E+01 0.00000E+00 0.00000E+00 0.00000E+00
7 -1.57407E+00 2.36111E+00 -4.72222E+00 0.00000E+00 0.00000E+00 0.00000E+00

----------------------------------------------------------------------------------------
FX FY FZ
-1.000000E+01 -5.000000E+00 1.500000E+01
LOAD CASE ID NO. 1
1 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00

2 1.85185E-02 -1.80041E-03 -2.70062E-03 0.00000E+00 0.00000E+00 0.00000E+00
3 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
5 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
7 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
STRESS RESULTANTS FOR LINE ELEMENTS - LOAD CASE ID NO. 1
ELE ELE NODE FORCE SHEAR SHEAR TORQUE MOMENT MOMENT

NUMBER NKTP NUMBER AXIAL LOCAL-Y LOCAL-Z AXIAL LOCAL-Y LOCAL-Z
1 14 1 -3.24074E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00

2 3.24074E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
2 14 3 6.15741E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00

2 -6.15741E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
3 14 5 2.00926E+01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00

2 -2.00926E+01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
4 14 7 -5.50926E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00

2 5.50926E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
[3D Spar Element Type ID]
[요소(부재) 번호]

일한지는 우선 확이 해야 한다. 그 이유는 하중의 조건이 정확하게 주어졌는지 확인의 절
차이기 때문이다(Fx=-10N, Fy=-5N, Fz=15N). 상반된 부호는 제하된 하중의 반력 부분이
기 때문에 반대적으로 반력이 발생한다.
► 각 절점에서의 반력 및 변형 결과 분석
절점번호 반 력(Reaction Force, N) 변형량(mm)
(Nodes) Rx Ry Rz X Y Z
1 -9.25926E-01 -1.38889E+00 -2.77778E+00 # # #
3 1.85185E-02 -1.80041E-03 -2.70062E-03 # # #
5 -5.74074E+00 -8.61111E+00 1.72222E+01 # # #
7 -1.57407E+00 2.36111E+00 -4.72222E+00 # # #
2 # # # 1.85185E-02 -1.80041E-03 -2.70062E-03
► 각 부재(3D Spar)에서의 축력(Axial Force)

부재번호
[방향] 축 력(Axial Force, N)
(Elements)
1 1 -3.24074E+00 2 3.24074E+00 인장부재
2 3 6.15741E+00 2 -6.15741E+00 압축부재
3 5 2.00926E+01 2 -2.00926E+01 압축부재
4 7 -5.50926E+00 2 5.50926E+00 인장부재
Z
④
Z 15 N
①
Y
1 7 Y
X ③
②
[인장/압축 부재의 판단 기준]
부재의 방향은 상기의 절점 순서
이다.(예, 부재 ″1″)은 절점 ″1″ ⇒
3 5 ″2″로 가는데 절점″1″의 방향이 ″-
″이고(부재방향 반대) 절점″2″는
X ″+″(부재방향)이므로 그림과 같이
양쪽으로 인장이 되는 인장 부재
[해석에 적용된Truss 구조물의 형상] 로 정의된다.

> 3dspar126 > 열기 > 3dspar127
> 열기

빔(Beam) 구조물의 유한요소 해석
Grid Create Method

Line Create Method

[예제 3] 빔(Beam) 구조물의 유한요소해석
아래 그림[1]과 같이 각각의 단면적이 일정한 세 개의 Beam구조가 절점 ″3″에 집중하중

이 작용하고, 요소 ③에 균일한 분포하중 작용하는 형상으로, 각 지점에서의 반력 및 축력
을 수치 계산에 의한 결과와 비교해 보자.
: 2차원 Beam의 요소는 3 자유도(DOF) UX, UY 및 MZ를 갖는다.
5KN
1.5 KN/m
4
3
봉의 반지름 : 13mm ③
탄성계수 : 208000 MPa
① ② 2m
4.8m
Y
1 Z 2
X
2m 2m
그림[1] Beam 구조물의 형상
각 지점 Rx Ry Mz
반력(N) 이론 NISA 이론 NISA 이론 NISA
1 1736.48 1.73655E+03 3713.19 3.71322E+03 -302.88 -3.02993E+02
2 -1473.61 -1.47365E+03 3991.20 3.99129E+03 -303.05 -3.03036E+02
4 -262.83 -2.62905E+02 3295.58 3.29549E+03 -2394.24 -2.39402E+03
☞ 연속체 또는 계를 분할한다.
앞의 그림에서와 같이 각각의 Beam은 하나의 요소와 2개의 절점으로 이루어진 계로 간

주한다. 그리고 각 요소의 특성은 다음과 같다.
• ① : 절점 1과 3을 연결하는 요소(길이=5.2m, 반경=13mm)

• ② : 절점 2와 3을 연결하는 요소(길이=5.2m, 반경=13mm)
• ③ : 절점 3과 4를 연결하는 요소(길이=4.0m, 반경=13mm)
• A1, A2, A3, : 각 요소 ①, ②, ③의 단면적(πr2 = 169π mm2)
• Ι : 관성모우멘트(πr4/4 = 7140.25 π mm2)
각 요소에 대한 강성행렬을 각각의 요소 ①, ②, ③, ④에 적용한다.
¨ 요소 ①
• C = Cosθ = 2/5.2, S = Sinθ = 4.8/5.2
⎡[K ] [K 13 ]⎤
[K ]1 = ⎢ 11
⎣[K 31 ] [K 33 ]⎥⎦
여기서
⎡ ⎛ 12 I 2 ⎞ ⎛ 12 I ⎞ 6I ⎤
⎢ − ⎜ AC + l 2 S ⎟, − ⎜ A − 2 ⎟CS , −
2
l ⎥⎥
S
⎢ ⎝ ⎠ ⎝ l ⎠
E⎢ ⎛ 12 I ⎞ ⎛ ⎞ 6I ⎥
[K13 ] = − ⎜ A − 2 ⎟CS ,
12 I
− ⎜ AC 2 + 2 C 2 ⎟, C
l ⎢ ⎝ l ⎠ ⎝ l ⎠ l ⎥
⎢ 6I 6I ⎥
⎢ S, C, 2I ⎥
⎣⎢ l l ⎦⎥
⎡ − 1000 .108, − 2399 .95500 , − 304 .2 ⎤

[K 13 ] ⎢
= π ⎢ − 2399 .95500 , − 5760 .01875 , 126 .75 ⎥⎥
⎢⎣ 304 .2, − 126 .75, 571220 .0⎥⎦
208000 ⎡ ⎤ 12 × 7140 .25 ⎛ 4.8 ⎞ 2

2
⎛ 2 ⎞
⎢169 × ⎜ ⎟ ⎥+ ×⎜ ⎟
5200 ⎢⎣ ⎝ 5 .2 ⎠ ⎥⎦ 5200 2 ⎝ 5 .2 ⎠
⎡1000 .108, 2399 .95500 , 304 .2 ⎤

[K 33 ] ⎢
= π ⎢ 5760 .01875 , − 126 .75 ⎥⎥
⎢⎣ 대칭 1142440 .0⎥⎦
¨ 요소 ②
• 절점 : 2와 3으로 이루어져 있다.
• C = Cosθ = -2/5.2, S = Sinθ = 4.8/5.2
⎡[K ] [K 23 ]⎤
[K ]2 = ⎢ 22
⎣[K 32 ] [K 33 ]⎥⎦
여기서
⎡ − 1000 .108, 2399 .95500 , − 304 .2 ⎤

[K 23 ] ⎢
= π ⎢ 2399 .95500 , − 5760 .01875 , − 126 .75 ⎥⎥
⎢⎣ 304 .2, 126 .75, 571220 .0⎥⎦
⎡1000 .108, − 2399 .95500 , 304 .2 ⎤

[K 33 ] ⎢
= π ⎢ 5760 .01875 , 126 .75 ⎥⎥
⎢⎣ 대칭 1142440 .0⎥⎦
¨ 요소 ③
• C = Cosθ = 1, S = 0
• 요소 ③은 국부 좌표계와 대표 좌표계가 같으므로, 국부 좌표계에 대한 수식(강성 행렬,
절점 하중, 변위 등)이 바로 대표 좌표계에 대한 수식이 된다.
⎡[K ] [K 34 ]⎤
[K ]3 = ⎢ 33
⎣[K 42 ] [K 44 ]⎥⎦
여기서
⎡8788, 0, 0 ⎤
[K 33 ] ⎢
= π ⎢ 0.27847 , 556 .9395 ⎥⎥ = [k 33 ]
⎢⎣ 대칭 1485172 .0⎥⎦
⎡ − 888, 0, 0 ⎤
[K 33 ] ⎢
= π ⎢ 0, − 0.27847 , − 556 .9395 ⎥⎥ = [k 43 ]
⎢⎣ 0, 556 .9395, 742586 .0 ⎥⎦
그리고 절점 부하 {f }는 횡 방향 하중만 존재하므로 f 3 x = 0, f 4 x = 0 이며, 분포하중에

의한 상당 절점 부하는 다음과 같다.
{f 3x , f 3 y , m 3 , f 4 x , f 4 y , m 4 }0 = {0, − 3 N , − 2000 N ⋅ mm , 0, − 3 N , 2000 N ⋅ mm }× 10 3

4
따라서 절점 부하는,
{F3 X , F3Y , M 3 , F4 X , F4Y , M 4 }T = {f 3x , f 3 y , m3 , f 4 x , f 4 y , m 4 }

T
= {0, − 8 N , − 2000 N ⋅ mm , 0, − 3 N , 2000 N ⋅ mm }× 10 3
☞ 각 요소에 대한 식을 결합하여 전체계에 대한 식을 유도하고 경계조건을 부여한다.
⎡[K11 ] 0 [K13 ] 0 ⎤ ⎧ 0 ⎫ ⎧ F1 ⎫
⎢ 0 [K 22 ] [K 23 ] 0 ⎥⎥ ⎪⎪ 0 ⎪⎪ ⎪F ⎪
⎢ ⎪ ⎪
= ⎨ 2⎬
⎢[K 31 ] [K 32 ] [K 33 ]g [K 34 ]⎥ ⎨⎪Q3 ⎬⎪ ⎪ F3 ⎪
⎢
⎣ 0 0 [K 43 ] [K 44 ]⎥⎦ ⎪⎩ F4 ⎪⎭ ⎪⎩ F4 ⎪⎭
여기서 {Qi }T = {QiX , QiY , Φ i }, {Fi }T = {FiX , FiY , M i } ,
[K 33 ]g = [K 33 ]1 + [K 33 ]2 + [K 33 ]3
⎡10788 .216, 0, 608 .4 ⎤

= π ⎢⎢ 11520 .31597 , 556 .9395 ⎥⎥
⎢⎣ 대칭 3770052 .0 ⎥⎦
[k 33 ]g {Q3 } = {F3 }
로부터
⎡Q3 X ⎤ ⎧ 0.00952 mm ⎫
⎢ Q ⎥ = ⎪− 0.021288 mm ⎪
⎢ 3Y ⎥ ⎨ ⎬
⎢⎣ Φ 3 ⎥⎦ ⎪ ⎪
⎩ − 0.16883 rad ⎭
☞ 각 절점에서의 반력을 구하여 보자.
⎧ F1 X ⎫ ⎧ 1736 .48 N ⎫
⎪ ⎪ ⎪ ⎪
⎨ F1Y ⎬ = [K13 ] {Q3 } = ⎨ 3713 .19 N ⎬
⎪M ⎪ ⎪− 300 .88 N ⋅ m ⎪
⎩ 1⎭ ⎩ ⎭
⎧ F2 X ⎫ ⎧ − 1473 .61N ⎫
⎪ ⎪ ⎪ ⎪
⎨ F1Y ⎬ = [K 23 ] {Q3 } = ⎨ 3991 .20 N ⎬
⎪M ⎪ ⎪− 303 .05 N ⋅ m ⎪
⎩ 2⎭ ⎩ ⎭
⎧ F4 X ⎫ ⎧ − 262 .83 N ⎫
⎪ ⎪ ⎪ ⎪
⎨ F4Y ⎬ = [K 43 ] {Q3 } = ⎨ 295 .58 N ⎬
⎪M ⎪ ⎪− 394 .24 N ⋅ m ⎪
⎩ 4⎭ ⎩ ⎭
그런데, 절점 4에서는 요소 ③의 분포하중이 상당 절점 하중으로 더해졌으므로 실제로 절

점 ″4″에서 작용하는 하중은 이를 뺀 값이 된다.
즉 다음과 같다.
⎧ F4 X ⎫ ⎧− 262 .83⎫ ⎧ 0 ⎫ ⎧ − 262 .83 N ⎫

⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪
⎨ F4Y ⎬ = ⎨ 295 .58 ⎬ − ⎨− 3000 ⎬ = ⎨ 3295 .58 N ⎬
⎪M ⎪ ⎪− 394 .24 ⎪ ⎪ 2000 ⎪ ⎪− 2394 .24 N ⋅ m ⎪
⎩ 4 ⎭a ⎩ ⎭ ⎩ ⎭ ⎩ ⎭
참고로 평형 관계를 확인하면,
∑F X = 1736 .48 − 173.61 − 262.83 = 0.04 ≈ 0
∑F Y = 3713 .19 + 3991 .20 − 5000 − (1500 )( 4) + 3295 .58 = − 0.03 ≈ 0
∑ M 1 = − 302.88 − 303.05 + (3991.20)(4) − (5000)(2) − ∫ (1500)(2 + x)dx

0
− 2394 .24 + ( 262.83)( 4.8) + (3295 .58)(6) = − 0.306 ≈ 0
여기서 ΣM1은 절점1에 대한 모멘트의 합을 의미하며 다른 절점 또는 임의의 위치에 대하

여 계산해도 결과는 마찬가지가 된다. 그런데, 힘과 모멘트의 합이 각각 정확하게 0이 되지
않는 것은 수치상의 오차(Round-off Error) 때문이다.
☞ 각 요소에 대한 국부 좌표계에서의 절점 하중을 계산해 보자.
각 요소에서의 변환 행렬은 다음과 같다.
⎡ 5 12 0 0 0 0⎤
⎢ − 12 5 0 0 0 0 ⎥⎥
⎢
1 ⎢ 0 0 13 0 0 0⎥
[T ]1 = ⎢ ⎥
13 ⎢ 0 0 0 5 12 0 ⎥
⎢ 0 0 0 − 12 5 0 ⎥
⎢ ⎥
⎢⎣ 0 0 0 0 0 13⎥⎦
⎡ − 5 12 0 0 0 0⎤
⎢− 12 − 5 0 0 0 0 ⎥⎥
⎢
1 ⎢ 0 0 13 0 0 0⎥
[T ]2 = ⎢ ⎥
13 ⎢ 0 0 0 − 5 12 0 ⎥
⎢ 0 0 0 − 12 − 5 0 ⎥
⎢ ⎥
⎢⎣ 0 0 0 0 0 13⎥⎦
따라서 각 요소에 대한 절점 변위 {q}는
{Q}1T = {0, Q3 }, {Q}T2 = {0, Q3 }
로 다음의 식에 대입하면,
{q} = [T ] {Q}
⎧ q1 x ⎫ ⎧ 0 ⎫ ⎧ q2 x ⎫ ⎧ 0 ⎫
⎪q ⎪ ⎪ 0 ⎪ ⎪q ⎪ ⎪ 0 ⎪
⎪ 1y ⎪ ⎪ ⎪ ⎪ 2y ⎪ ⎪ ⎪
⎪⎪ φ1 ⎪⎪ ⎪⎪ 0 ⎪⎪ ⎪⎪ φ 2 ⎪⎪ ⎪⎪ 0 ⎪⎪
⎨ ⎬ = ⎨ ⎬ , ⎨ ⎬ = ⎨ ⎬
⎪ q3 x ⎪ ⎪ − 0.19284 mm ⎪ ⎪ q3 x ⎪ ⎪− 0.20017 mm ⎪
⎪ q3 y ⎪ ⎪− 0.09067 mm ⎪ ⎪ q3 y ⎪ ⎪− 0.07309 mm ⎪
⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪
⎪⎩ φ 3 ⎪⎭1 ⎪⎩ − 0.16883 mm ⎪⎭ ⎪⎩ φ 3 ⎪⎭ 2 ⎪⎩ − 0.16883 mm ⎪⎭
강성 행렬(국부 좌표계)을 구하면,
⎡[K ] [K 13 ]⎤ ⎡[K ] [K 23 ]⎤
[K ]1 = ⎢ 11 [K ]2 = ⎢ 22
[K 33 ]⎥⎦ [K 33 ]⎥⎦
,
⎣[K 13 ] ⎣[K 32 ]
여기서
⎡ − 6760, 0, 0 ⎤
[k13 ] = [k 23 ] ⎢
= π ⎢ 0, − 0.12675, 329 .55 ⎥⎥
⎢⎣ 0, − 329 .55, 571220 .0⎥⎦
⎡6760, 0, 0, ⎤
[k33 ]1 = [k33 ]2 ⎢
= π ⎢ 0.12675 − 329 .55 ⎥⎥
⎢⎣ 대칭 1142440 .0⎥⎦
따라서 각 절점에서의 반력을 구하면 다음과 같다.
⎧ f1x ⎫ ⎧ 4095 .38 N ⎫

⎪f ⎪ ⎪ − 174 .76 N ⎪
⎪ 1y ⎪ ⎪ ⎪
⎪⎪ m1 ⎪⎪ ⎪⎪ − 302 .88 N ⋅ m ⎪⎪
⎨ ⎬ = ⎨ ⎬
⎪ f 3x ⎪ ⎪ − 4095 .38 N ⎪
⎪ f3y ⎪ ⎪ 174 .76 N ⎪
⎪ ⎪ ⎪ ⎪
⎪⎩ m 3 ⎪⎭1 ⎪⎩− 605 .85 N ⋅ m ⎪⎭
⎧ f 2x ⎫ ⎧ 4251 .04 N ⎫
⎪f ⎪ ⎪ − 174 .76 N ⎪
⎪ 2y ⎪ ⎪ ⎪
⎪⎪ m 2 ⎪⎪ ⎪⎪ − 302 .90 N ⋅ m ⎪⎪
⎨ ⎬ = ⎨ ⎬
⎪ f 3x ⎪ ⎪ − 4251 .04 N ⎪
⎪ f3y ⎪ ⎪ 174 .76 N ⎪
⎪ ⎪ ⎪ ⎪
⎪⎩ m 3 ⎪⎭ 2 ⎪⎩− 605 .87 N ⋅ m ⎪⎭
⎧ f 3x ⎫ ⎧ 262 .83 ⎫ ⎧ 0 ⎫ ⎧ 262 .83 N ⎫
⎪f ⎪ ⎪− 295 .58 ⎪ ⎪ − 3000 ⎪ ⎪ 2704 .42 N ⎪
⎪ 3y ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪
⎪⎪ m 3 ⎪⎪ ⎡[K 33 ]3 [K 34 ]⎤ ⎧q 3 ⎫ ⎪⎪− 788 .10 ⎪⎪ ⎪⎪− 2000 ⎪⎪ ⎪⎪ 1211 .90 N ⋅ m ⎪⎪
⎨ ⎬= − { f }0 =⎨ ⎬−⎨ ⎬=⎨
⎪ f 4x ⎪
⎢ [K ]
⎣ 43 [K 44 ]⎥⎦ ⎨⎩ 0 ⎬⎭ ⎪− 262 .83⎪ ⎪ 0 ⎪ ⎪ − 262 .83 N ⎪
⎬
⎪ f4y ⎪ ⎪ 295 .58 ⎪ ⎪ − 3000 ⎪ ⎪ 3295 .58 N ⎪
⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪
⎪⎩ m 4 ⎪⎭ ⎪⎩− 394 .24 ⎪⎭ ⎪⎩ 2000 ⎪⎭ ⎪⎩− 2394 .24 N ⋅ m ⎪⎭
각 요소에 대한 자유 물체도를 그려 보면 다음 그림[3]과 같이 된다. 각 요소에 대한 평

형관계 역시 만족됨을 알 수 있다.
예를 들어 요소 ①에 대한 모멘트의 평형을 분석하면 다음과 같다.
∑m 1 = − 302 .88 − 605 .85 + (174 .76)(5.2) = − 0.022 ≈ 0
그런데, 여기에 각 절점에 대한 힘의 합이 대표 좌표계에서의 힘과 일치 해야 하는데 다

소 다른 경우를 보게된다.
간단한 예로, 절점 ″2″ 에서의 모멘트는 국부 좌표게에서는 302.90, 대표 좌표계에서는
303.05를 얻었는데 이는 수치 계산상의 오차(Round-off Error) 때문이다.

Ä 좌표축을 윗면(Top)으로 설정한다.

• X-Y(Top)가 기본 좌표계이나 상이한 경우에는 이를 정의한다.
Menu를 이용하여 Gird를 정의하는데 각 좌표를 ″m″ 단위로 생성한다.

► XYZ-Data[0/0/0] > OK ---- G1 생성
► XYZ-Data[4/0/0] > OK ---- G2 생성
► XYZ-Data[2/4.8/0] > OK ---- G3 생성
► XYZ-Data[6/4.8/0] > OK ---- G4 생성
► Cancel ---- Grid 생성이 완료되면 Grid 생성화면 지움
GRIDS

► Grid IDs[1/2/3,3/3/4] ⇒ OK ---- L1 to L3 생성
► Cancel
LINES
Ä 2D Beam 요소를 생성한다(E1 to E3)
• FEM > MESH > FEG OPTION > BAR/BEAM/SPRING

► No.of Elements[1]
► > 2D Beam > OK
► Line Ids > Input Ids[1/2/3] > OK > OK > Cancel
⇒ E1 to E3 생성
ELEMENT
NODES
[2D Beam 요소생성] [절점을 지운 상태 요소 표현]


Menu 방식을 이용하여 절점을 일체화 한다(참조).

> OK > OK(MAIN MENU > YES ALL > CANCEL
NODE
MERGE

> BC Data > Click ″Ux″ > Click ″Uy” > Click ″Rz″ > OK
> Node Ids > Input Ids[1/3/6] > OK > OK > Cancel
• BC DATA > Ux=Uy=Rz=0(Constrain,고정)
구조물의 하단부분을 완전고정을 하기 위해

Ux,Uy,Rz 부분을 선택한다.
2차원 Beam요소는 자유도(DOF)가 3 방향성만
있기 때문에 Uz,Rx,Ry는 선택하지 않는다. 여기
서 선택시 ″0″이라는 숫자가 나타나는데 ″0″은
곳 구조물 선택 부분이 변형이 없다는 의미이고
인위적으로 숫자를 입력하면 강제변위의 하중으
로 정의된다는데 주의한다.
경계조건 부분의 기호 설명
Ry
Ux
있기 때문에 본 Beam 요소에서는 Ux,Uy,Rz를
제외한 DOF는 의미가 없다. 예, 다른 DOF를 Uy
구속하여도 의미는 없고 경고문이 주어진다.
Ä 외력 집중(하중)조건을 정의한다.
• FE GENERATION > STRUCTURAL > FORCE.MOMENT> ADD AN NODES
• BC DATA > Uy[5000] > OK

Fy : 좌표상에서 -축 방향 ► -5000

> Cancel

집중하중
Ä 외력 균일분포(하중)조건을 정의한다.
• FE GENERATION > STRUCTURAL > BEAM LOADS > ADD AN ELEMENTS
• BC DATA
► Loc. of Load(from
⇒ In Terms of Beam Length
► Dir. of Load
⇒ Local XY Plane
► Distributed Load
⇒ Start Loc : 0, End Loc : 1
⇒ Load at Start : -1500
Load at End : -1500
► OK
• ELEMENT IDS > Input Ids[3] > OK > OK

> Cancel
Beam의 균일분포 하중은 상기에서 정의한 방법 집중하중

은 Beam의 전체 길이를 Unit로 보고 각 양단의
절점에 작용하는 Value를 입력하면 된다.
상기의 경우는 Beam ″3″부분에 양끝단의 분포
하중이 같은(1500 to 1500) 균일 분포 하중의 경
우이다.
균일분포하중
Ä 단면특성을 정의한다.

• FEM > FE MODEL DATA
> BEAM SECTIONS
> GENERAL SECTION
¾ 선택한다.
> DIAMETER[0.013] : 환봉의 직경(0.013m)

> OK
지름이 13mm인 환봉의 단면특성 확인
> 직경을 입력한 후 > VIEW

Ä 재질특성을 정의한다.

• ANALYSIS/STATIC > MATERIAL
> Ex = 2.08E11 : (2.08e11 N/m2)
> Nuxy = 0.3
> OK
재질특성은 기본이 등방성(Isotropic)이기

때문에 Ex, Nuxy 두개만 입력을 하면 자
동적으로 나머지 부분들은 입력이 된다.
단, 정방성(Orthotropic) 및 이방성
(Anisopropic)의 재질은 입력 방법이 상
이함에 유의해야 한다.
Ä Print Control를 정의한다.

각 부재(Spring, Bar, Beam, Cable등 Line요소)에 대한 축력(Axial Force)를 Text Output으
로 확인하기 위해서는 기본이 ″Off″로 되어 있는 부분을 ″On″으로 정의한다.
• ANALYSIS/STATIC > PRINT CONTROL

> ADD > SUMMARY OF LINE ELMENT FORCES(SLFO) > PRINT ALL > OK
Ä 해석 및 결과화일을 정의한다.


접근한다.
야 한다.

Enter Filename > 2dbeam1
> 확인
• FILE > SAVE AS > 파일 이름[2dbeam1]

> Yes All :Beam의 단면 형상이 주어질
경우 결과 그림에서 실 3차원 형상
구현이 가능하다.
> 2dbeam1 : Shapes내용 수정
> 저장(S)

다.
⇒

> Type in NISA data file name > 2dbeam1.nis > Enter
선정한 파일명을 입력(2dbeam1.nis)하고 Enter를 한다.
> Type in output file name > 2dbeam.out > Enter
각종 결과를 Text로 출력 할 수 있는 파일명을 입력한다.
2-13
enter
enter
Ä NISA Input 파일 확인(2dbeam1.nis)
**EXECUTIVE data deck

ANALYSIS = STATIC
SOLV = FRON
FILE = 2dbeam1
SAVE = 26,27,35 ** File Name
*ELTYPE
1, 13, 1 ** NKTP : 13
*RCTABLE
1, 43,1, 0 ** Property ID(1), 입력 항(4)
1.3273230E-04, 1.4019850E-09, 1.4019850E-09, 2.8039699E-09, 단면적,관성모멘트,Transverse,..
0.0000000E+00, 0.0000000E+00, 1.3000000E-02, 1.3000000E-02,
0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00,
0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00,
0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00,
0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00,
0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00,
0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00,
0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00,
0.0000000E+00, 0.0000000E+00, 0.0000000E+00, 0.0000000E+00,
0.0000000E+00, 0.0000000E+00, 0.0000000E+00,
*NODES
** 각 절점의 좌표값
1,,,, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0
2,,,, 2.00000E+00, 4.80000E+00, 0.00000E+00, 0
3,,,, 4.00000E+00, 0.00000E+00, 0.00000E+00, 0
6,,,, 6.00000E+00, 4.80000E+00, 0.00000E+00, 0
*ELEMENTS ** 각 요소의 절점 연결성
1, 1, 1, 1, 0
1, 2,
2, 1, 1, 1, 0
3, 2,
3, 1, 1, 1, 0
2, 6,
*MATERIAL ** 재질특성
EX , 1,0, 2.08000E+11,
NUXY, 1,0, 3.00000E-01, ** 하중 제한 조건
*LDCASE, ID= 1
0,1,1,0,-1,2,0 ** 경계조건
*SPDISP 각 절점의 경계조건 구속
** SPDISP, SET = 1
1,UX , 0.00000E+00,,,,,,,, 0
1,UY , 0.00000E+00,,,,,,,, 0
1,ROTZ, 0.00000E+00,,,,,,,, 0
3,UX , 0.00000E+00,,,,,,,, 0
3,UY , 0.00000E+00,,,,,,,, 0
3,ROTZ, 0.00000E+00,,,,,,,, 0
6,UX , 0.00000E+00,,,,,,,, 0
6,UY , 0.00000E+00,,,,,,,, 0
6,ROTZ, 0.00000E+00,,,,,,,, 0
*BEAMLOAD ** Beam Loading
**_DISP3_: BEAMLOAD, SET = 1
3, 0,0,0,,, 0.000E+00, 1.000E+00,-1.500E+03,-1.500E+03,0
*CFORCE ** 집중하중 조건
** CFORCE, SET = 1
2,FY,-5.00000E+03,,, 0, 0
*PRINTCNTL
SLFO,0
*ENDDATA
Ä 각 절점에서 반력,축력 및 변형량을 확인한다(spring1.out).
LOAD CASE ID NO. 1
1 1.73655E+03 3.71322E+03 0.00000E+00 0.00000E+00 0.00000E+00 -3.02993E+02

3 -1.47365E+03 3.99129E+03 0.00000E+00 0.00000E+00 0.00000E+00 -3.03036E+02
6 -2.62905E+02 3.29549E+03 0.00000E+00 0.00000E+00 0.00000E+00 -2.39402E+03

--------------------------------------------------------------------------------------
FX FY FZ
2.273737E-13 1.100000E+04 0.000000E+00
LOAD CASE ID NO. 1
1 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00

2 3.80908E-05 -8.51535E-04 0.00000E+00 0.00000E+00 0.00000E+00 -2.70168E+00
3 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
6 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
STRESS RESULTANTS FOR LINE ELEMENTS - LOAD CASE ID NO. 1
ELE ELE NODE FORCE SHEAR SHEAR TORQUE MOMENT MOMENT

NUMBER NKTP NUMBER AXIAL LOCAL-Y LOCAL-Z AXIAL LOCAL-Y LOCAL-Z
1 13 1 4.09549E+03 -1.74808E+02 0.00000E+00 0.00000E+00 0.00000E+00 -3.02993E+02

2 -4.09549E+03 1.74808E+02 0.00000E+00 0.00000E+00 0.00000E+00 -6.06010E+02
2 13 3 4.25106E+03 -1.74825E+02 0.00000E+00 0.00000E+00 0.00000E+00 -3.03036E+02

2 -4.25106E+03 1.74825E+02 0.00000E+00 0.00000E+00 0.00000E+00 -6.06053E+02
3 13 2 2.62905E+02 2.70451E+03 0.00000E+00 0.00000E+00 0.00000E+00 1.21206E+03

6 -2.62905E+02 3.29549E+03 0.00000E+00 0.00000E+00 0.00000E+00 -2.39402E+03
[3D Spar Element Type ID]
[요소(부재) 번호]

일한지는 우선 확이 해야 한다. 그 이유는 하중의 조건이 정확하게 주어졌는지 확인의 절차
이기 때문이다(Fy=-11000N). 상반된 부호는 제하된 하중의 반력 부분이기 때문에 반대적으
로 반력이 발생한다.
► 각 절점에서의 반력 및 변형 결과 분석
절점번호 반 력(Reaction Force, N) 변형량(mm)
(Nodes) Rx Ry Mz X Y Rz
1 1.73655E+03 3.71322E+03 -3.02993E+02 # # #
2 # # # 3.80908E-05 -8.51535E-04 -2.70168E+00
3 -1.47365E+03 3.99129E+03 -3.03036E+02 # # #
6 -2.62905E+02 3.29549E+03 -2.39402E+03 # # #
► 각 부재(2D Beam)에서의 축력(Axial Force)

부재번호 축 력(Axial Force, N)
[방향]
(Elements)
1 1 4.09549E+03 2 -4.09549E+03 압축부재
2 3 4.25106E+03 2 -4.25106E+03 압축부재
3 2 2.62905E+02 6 -2.62905E+02 압축부재
6
2
③
[인장/압축 부재의 판단 기준]

① ② 모든 부재의 축력 방향이 부재의
중심으로 향하고 있어 본 구조의
모든 부재는 압축을 받고있다.
Y
1 3
X
[해석에 적용된Truss 구조물의 형상]


> 2dbeam126 > 열기 > 2dbeam127
> 열기


KEC plant pressure vessel training manual 4 of 5 권 PDF

Uploaded by

Copyright:

Available Formats

KEC plant pressure vessel training manual 4 of 5 권 PDF

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

KEC plant pressure vessel training manual 4 of 5 권 PDF

Uploaded by

Copyright:

Available Formats

What is the document about?

What is the document about?

What topics are covered in the document?

What topics are covered in the document?

KOREA ENGINEERING CONSULTING, INC.

PRESSURE VESSEL BASIC CONCEPT

Pressure Vessel Basic Concept

Korea Engineering Consulting, Inc.

==== Contents ===

[1권Book] : Theoretical Strength Calculation & FEA

[[2권Book] : Theoretical Strength Calculation & FEA

[3권Book] : Theoretical Strength Calculation & FEA

[4권Book] : FEA(Finite Element Analysis) Theoretical Course

z Theoretical Course ″B″

z Theoretical Course ″C″

[5권Book] : Crack Propagation & Facture Mechanics

1. Crack Propagation ..................................................................................................................... 7

″Finite Element Theoretical & Analysis by NISA″

Theoretical Course ″A″

This page internally left blank

1. 유한요소법(Finite Element Method)의 개념

유한요소법(Finite Element Method, 有限要素法)은 1930년대에 개발된 행렬(Matrix)의

유한요소법은 초기에는 공학 및 물리학에서 많이 사용되었으나 현재에는 공학의 여러 분

유한요소법은 초기에는 고전적인 이론식으로 해석할 수 없는 복잡한 형상의 구조물의 응

기술의 발전과 더불어 더욱 복잡한 형상이나 하중 상태 및 물성치 들을 포함하는 문제를

좀더 유한요소법을 설명하기 전에 처음 입문자를 위하여 전산기를 이용해 현재 공학에 적

단) 현재 CAD, CAM 및 CAE가 복합적으로 통합화 된 프로그램도 있고, CAD 및 CAE만

[CAD] : ″Computer Aided Design″

[CAM] : ″Computer Aided Manufacturing″

[CAE] : ″Computer Aided Engineering″

이러한 알고리즘의 적용범위는 결과의 정도 와 결과를 도출하기 까지의 시간에 영향을

[FEM] : 유한요소법(Finite Element Method)

단 유체의 유동 문제는 많은 요소 및 내부 알고리즘의 문제로 많은 시간이 소요되기 때문

앞서 언급하였듯이 유한요소법은 거의 모든 분야에 적용이 가능하기 때문에 적용분야 별

그리고 아래의 소프트웨어는 대부분 외국에서 개발된 프로그램이며 국내에서 개발된 프로

[FVM] : 유한체적법(Finite Volume Method)

-. 내부 및 외부 유동, 사출, 주물, 화학반응 등의 압축 및 비압축성 유체

[FDM] : 유한차분법(Finite Differential Method)

-. Moldflow, Flow-3D, C-Mold, MAGMA(주물),… 등

z 여러 가지 복합재료로 이루어진 연속체에도 사용이 가능하다.

z 재료의 물성과 구조물 거동의 비선형에도 적용이 가능하다.

z 응력, 좌굴, 진동, 열전달, 유동/열유동, 전자기장, 운동학 및 피로 및 해석이 가능하다.

z 국부응력해석에 적용할 시에는 많은 경험이 필요하다.

z 실구조물을 완벽하게 유한요소로 치환하는 것은 어려우므로 구조물을 단순화 한다.

z 대형구조물인 경우 높은 성능의 전산기가 필요하다.

위에서 설명한 내용은 반드시 읽어보아야 하고 기억하여야 한다. 아무리 많은 경험자일지

This page internally left blank

유한요소법(FEM)은 응력, 진동, 피로 및 파괴 해석 등의 구조 분야 와 열, 유체 및 전자

1단계, 기하학적 형상 및 연속체를 정립한다(전처리 단계, Pre/Processor)

응력, 변형 등을 분석하고자 하는 구조형상에서 기하학적 형상을 기반으로 해석 영역 부

외부의 CAD 데이터를 이용하

그림[2.1] 유한요소해석의 단계별 기능

1) 기하학적 형상 생성(Geometric Modeling)

모든 형상을 생성하는 명령어는 소, 대문자 구분 없이 2차 및 3차원의 공간을 사용한다.

그림[2.5] Hyperpatch 생성하기

앞의 기하학적 형상(Grid, Line, Patch, Hyperpatch)을 이용하여 연속체로 분할하는데 절

(1) 절점(Node)을 이용한 요소생성 : Mass(Point) Element

(3) 면(Patch)을 이용한 요소생성 :

(4) 입체(Hyperpatch)를 이용한 요소생성 : Solid, Composite,