• 전산구조공학
• /
• 제3권3호
• /
• pp.113-123
• /
• 1990

본 연구에서는 개선된 degenerated 쉘 유한요소의 탄소성 및 기하학적 비선형 해석에의 적용성을 고찰하였다. 본 연구의 개선된 degenerated 쉘요소는 shear locking 해결에 우수한 결과를 보인 가정된 전단변형도를 대치사용하고, membrance locking 현상을 제거하기위해 평면내 변형도의 구성시 감차적분을 행하며, 쉘요소 자체의 거동을 보완하기위해 비적합변위형을 선택적으로 추가하였다. 본 요소는 shear/membrance locking이 발생하지 않으며, 전달가능한 거짓 영에너지모드도 나타나지않는다. 소성변형 정형화에서는 적층모델을 사용하며, 재료는 von Mises항복조건을 따른다고 가정한다. 유한변형을 고려한 기하학적 비선형 방정식을 total lagrangian수식화를 사용하여 정형화 하였고, 비선형 방정식은 하중제어 및 변위제어법을 사용한 Newton-Raphson 반복법으로 반복 계산한다. 여러 예제해석을 통하여 본 개선된 degenerated 쉘 유한요소의 정확도를 고찰하였다.

• 대한조선학회지
• /
• 제27권1호
• /
• pp.45-56
• /
• 1990

본 해석에서는 압력용기를 설계하기 위하여 복합적층된 원통형 구각의 비선형 해석을 유한요소법으로 수행하였다. 적층순서의 변화에 따라 최소변위 또는 최대압력을 갖는 최적의 적층구조를 얻기 위하여 8절점 Isoparametric 격하요소를 사용하며 구조요소의 비선형거동은 Total Lagrangian 수식과 하중증분법을 적용하여 해석하며 평형반복수렴은 Newton-Raphson Method를 이용하였다. 선형해석의 경우에 9가지 적층구조를 선정하여 하중조건이 내압일때 최소변위를 나타내는 적층구조를 조사한 결과 $[50^{\circ}/-50^{\circ}]$의 최적구조를 구하였고 적층순서를 $[{\theta}^{\circ}/{-\theta}^{\circ}]$로 하여 비선형해석과 동시에 Quadratic Failure Criteria를 적용하였을 때 하중조건이 외압만을 고려하는 상태에서도 $\theta=50^{\circ}$가 최소 변위비를 나타내었고 외압과 축하중을 동시에 고려한 상태에서는 $\theta=90^{\circ}$였다.

• 품질경영학회지
• /
• 제9권2호
• /
• pp.31-36
• /
• 1981

점증되고 있는 고신뢰성 제품의 설계에 있어 주어진 선형제약조건 내(內)에서 체계의 신뢰성을 최대화하는 방법을 소개하고 이를 해결하는 비선형계획법을 반복단계로 하여 Computer Programming을 하였다. 단, 본 논문에서 다루는 체계는 병렬중복구조를 갖는 직렬다단계 구조이다. 타당성 검토를 위한 예제를 해결하였으며 Computer Programming은 지면관계로 생략하였다.

• 대한기계학회논문집
• /
• 제17권12호
• /
• pp.3007-3019
• /
• 1993

For the same configuration, the stiffness of 6-node two dimensional isoparametric is stiffer than that of 8-node two dimensional isoparametric element. This phenomenon may be called 'Relative Stiffness Stiffening Phenomenon.' In this paper, the relative stiffness stiffening phenomenon was studied, and could be corrected by modifying the position of Gauss integral points used in the numerical integration of the stiffness matrix. For the same deformation (bending) energy of 6-node and 8-node two dimensional isoparametric elements, Gauss integral points of 6-node element have to move closer, in comparison with those of 8-node element, in the case of numerical integration along the thickness direction.

• 대한조선학회논문집
• /
• 제39권1호
• /
• pp.28-37
• /
• 2002

본 연구에서는 저장력 예인케이블의 비선형 동적거동을 수치적으로 해석하였다. 고장력 케이블해석에서는 흔히 무시되는 굽힘강성의 효과가 저장력 케이블에서는 중요한 역할을 하므로 본 연구에서는 이를 고려하였다. 또한 저장력 케이블에서는 대변위가 발생하기 쉬우므로 기하학적인 비선형 및 유체 비선형 효과가 크므로 이를 고려하였다. 저장력 예인케이블에 대한 3차원 비선형 운동방정식을 수립하고 유한차분법을 적용하여 이산화 시켰다. 시간적분에 있어서 안정적인 해를 얻을 수 있는 음해법(implicit method)을 적용하였으며 비선형 해를 구하기 위하여 Newton-Raphson 반복법을 사용하였다. 케이블과 같이 양단경계조건을 갖고 대각선 주변 성분만 있는 행렬식을 계산하는 경우에는 Gauss-Jordan 방법 등과 같이 일반적인 방법 보다 블록삼중대각행렬 풀이법이 계산시간을 상당히 줄일 수 있음을 알 수 있었다. 몇 가지 예제해석을 수행하였으며 실해역 실험결과에 의해 이미 검증되어 있는 케이블 해석프로그램인 WHOI Cable 프로그램의 해석결과와 비교 검토한 결과 서로 잘 일치함을 알 수 있었다.

• 대한기계학회논문집
• /
• 제15권2호
• /
• pp.523-536
• /
• 1991

This dissertation analyzes the running mechanism of flexible and thin tape above rotating head through the numerical simulation and the experiment. The scope of analysis is confined to the phenomena of two dimensional elasto hydrodynamic lubrication between the protruded bump on a rotating cylinder and the running tape. This model is based on the elastic deformation equation of plate and shell and Reynolds equation. Finite difference method is employed as a numerical technique to calculate (1) the distribution of pressure between the running tape and rotating bump and (2) the vertical deformation of elastic thin tape over he rotating bump under hydrodynamic pressure. In numerical analyses, the effects of bump size on flying characteristics of the tape were evaluated and examined considering the influence of tension and stiffness of tape.

• 한국생물환경조절학회:학술대회논문집
• /
• 한국생물환경조절학회 2001년도 봄 학술발표논문집
• /
• pp.31-33
• /
• 2001

잠열축열 그린하우스 난방 시스템의 난방특성을 분석하기 위하여 이에 대한 열평형 이론을 정립하고 수치해석에 의하여 컴퓨터 시뮬레이션 모델을 개발하고자 잠열축열 그린하우스 난방 시스템의 열저항 회로망을 구성하였다. 그리고 그린하우스의 피복재, 내부 공기, 토양표면, 잠열 축열재와의 열평형 방정식을 구성하였으며, Newton-Raphson반복법을 이용하여 수치해석을 하였고, 실험 분석을 통하여 수치해의 타당성을 검증하였다. 시뮬레이션 모델을 위하여 C언어를 사용하였으며, 겨울철 (11월-2월)의 기후 조건이 유사한 여러 날을 선정하여 온도, 태양강도, 상대습도, 토양 수분함량 등을 자료로 하여 모델링을 하였다. 여기에 사용된 토양 조건은 사양토로 건조한 상태를 유지하였다. 이상과 같은 분석에 의하여 그린하우스내 경시적 공기온도 변화와 열전달 현상의 실험치와 이론분석 결과가 잘 일치하고 있음을 알 수 있었다.

• 대한토목학회논문집
• /
• 제32권3B호
• /
• pp.175-183
• /
• 2012

본 연구에서는 천수방정식에 Galerkin법과 Petrov-Galerkin 기법의 일종인 SU/PG 기법을 적용하여 유한요소 천수흐름 해석 모형을 개발하고, 다양한 실험수로에서 이송이 지배적인 흐름을 수치 모의하였다. 수로 내부에 얇은 판 형태의 구조물이 존재하는 경우 Fr 수와 Re 수가 매우 낮은 흐름에서는 Galerkin 기법과 SU/PG 기법이 동일한 결과를 나타냈으나, Fr=1.58인 경우 Galerkin법은 발산하여 해를 얻을 수 없었다. 이 경우 SU/PG법은 Newton-Raphson법에 의한 5회의 반복에 의해 수렴된 유속결과를 구할 수 있었다. 사류와 상류가 혼재하여 천이류가 나타나는 단면확대 수로 모의에서 SU/PG 기법을 이용한 본 연구의 경우 상류단 수심조건이 잘 유지되며, 도수가 발생하는 지점 및 도수에 의한 수심 경사, 도수 후의 수심이 모두 Khalifa(1980)의 실험결과와 매우 근사하였다. 이송이 지배적인 사류(Fr=2.74)에 의한 사각도수 모의의 경우에도 Galerkin 기법은 최초 모의시간의 첫 번째 반복 이후 발산하였으나, SU/PG 기법은 도수 경계면을 수치진동 없이 잘 포착하였으며, 해석해와 비교한 수심 및 유속의 최대 오차는 0.2% 이내로 나타나 기존 연구(Levin 등, 2006; Ricchiuto 등, 2007)에 비해 더욱 정확한 결과를 도출하였다.

• 한국안전학회지
• /
• 제11권4호
• /
• pp.143-150
• /
• 1996

본 연구는 선형, 비선형 hygrothermal 응력 문제를 위한 explicit-Implicit 유한요소 해석 모델 개발에 관한 것이다. 부가적으로 moilsture 확산 방정식, J-적분 평가를 위한 균열 요소 및 가상 균열 진전법이 도입된다. 시간 변화에 따른 균열 추진력을 계산하기 위하여 선형 탄성 파괴 역학(LEFM)이론이 고려되며 재료의 기공은 실온에서 액체 상태의 습기로 포화되어 있으며 온도가 상승함에 따라 증기화된다는 가정하에서 균열 추진력과 증기 효과의 관계가 연구된다. 이상 기체방정식은 각 시간 단계에서 증기에 의한 열역학적 압력을 계산하기 위하여 이용된다. 다공질 재료의 시간 종속 응답을 지배하는 방정식들은 혼합이론에 기초하며 다공질 재료의 유체 흐름을 위한 Darcy의 법칙과 Von-Mises 항복 기준을 포함하고 있는 Perzyna의 점소성 모델이 첨가된다. 또한 Green-Naghdi 응력률이 중첩된 강체 운동하에서 응력 텐서 invariant로 사용되며, 모델링을 위하여 사각요소가 이용되고 비선형 지배 방정식을 풀기 위하여 full Newton-Raphson법에 의한 반복법이 사용된다. 본 연구를 통하여 얻은 결과는 다음과 같다. 1) 본 유한요소 프로그램은 복합재의 hygrothermal 파괴 해석에 매우 유용하게 적용될 수 있다. 2) 습기의 온도에 의한 영향을 가지는 재료의 J-적분을 정확히 예측하기 위하여는 증기 효과를 고려하여야 한다. 왜냐하면 초기단계에 균열 전파력이 가속되기 때문이다. 3) 본 해석을 위해 Uncoupled scheme에 의한 결과도 Coupled scheme에 결과에 비해 아주 타당하므로 CPU 측면에서 매우 경제적인 Uncoupled scheme이 추천된다.

• 대한토목학회논문집
• /
• 제36권6호
• /
• pp.947-954
• /
• 2016

고정하중을 받는 사장교의 최적 초기형상 결정을 위한 간략해석법을 제시한다. 이 방법은 비선형 유한요소 해석을 수행하지 않고 최적의 초기형상을 구현할 수 있으며, 최소의 휨모멘트를 발생시키는 무응력길이를 제공한다. 초기형상 해석 후 여타의 하중조합 또는 시공단계해석을 수행하기 위해 무응력길이를 변화시키지 않고 Newton-Raphson 반복법을 수행하는 무응력길이법을 이용한다. 수치예제는 제작캠버를 포함한 인천대교 2차원 모델이며, 제시된 해석법과 상용 유한요소해석 프로그램 MIDAS의 초기형상 해석결과를 비교 분석한다. 주로 케이블의 장력과 주부재의 휨모멘트 및 변위를 관찰하였으며, 두 해석결과는 축방향 변위를 제외하고 모든 결과 값이 서로 잘 일치함을 확인할 수 있었다.