KR101028978B1 - 해제된 자유도를 갖는 제한적으로 변형가능한 시스템을위한 솔버 - Google Patents

해제된 자유도를 갖는 제한적으로 변형가능한 시스템을위한 솔버 Download PDF

Info

Publication number
KR101028978B1
KR101028978B1 KR1020060031941A KR20060031941A KR101028978B1 KR 101028978 B1 KR101028978 B1 KR 101028978B1 KR 1020060031941 A KR1020060031941 A KR 1020060031941A KR 20060031941 A KR20060031941 A KR 20060031941A KR 101028978 B1 KR101028978 B1 KR 101028978B1
Authority
KR
South Korea
Prior art keywords
node
frame
finite element
freedom
displacement
Prior art date
Application number
KR1020060031941A
Other languages
English (en)
Other versions
KR20060107407A (ko
Inventor
일레랭 샤를-앙드르 드
세바스띠앙 보겔르
Original Assignee
다솔 시스템므
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from EP05290790.4A external-priority patent/EP1710721B1/en
Application filed by 다솔 시스템므 filed Critical 다솔 시스템므
Publication of KR20060107407A publication Critical patent/KR20060107407A/ko
Application granted granted Critical
Publication of KR101028978B1 publication Critical patent/KR101028978B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F5/00—Orthopaedic methods or devices for non-surgical treatment of bones or joints; Nursing devices; Anti-rape devices
    • A61F5/37—Restraining devices for the body or for body parts, e.g. slings; Restraining shirts
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F5/00—Orthopaedic methods or devices for non-surgical treatment of bones or joints; Nursing devices; Anti-rape devices
    • A61F5/0003—Apparatus for the treatment of obesity; Anti-eating devices
    • A61F5/0013—Implantable devices or invasive measures
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61H—PHYSICAL THERAPY APPARATUS, e.g. DEVICES FOR LOCATING OR STIMULATING REFLEX POINTS IN THE BODY; ARTIFICIAL RESPIRATION; MASSAGE; BATHING DEVICES FOR SPECIAL THERAPEUTIC OR HYGIENIC PURPOSES OR SPECIFIC PARTS OF THE BODY
    • A61H19/00—Massage for the genitals; Devices for improving sexual intercourse
    • A61H19/40—Devices insertable in the genitals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61H—PHYSICAL THERAPY APPARATUS, e.g. DEVICES FOR LOCATING OR STIMULATING REFLEX POINTS IN THE BODY; ARTIFICIAL RESPIRATION; MASSAGE; BATHING DEVICES FOR SPECIAL THERAPEUTIC OR HYGIENIC PURPOSES OR SPECIFIC PARTS OF THE BODY
    • A61H2205/00—Devices for specific parts of the body
    • A61H2205/08—Trunk
    • A61H2205/087—Genitals
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/16—Cables, cable trees or wire harnesses
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S706/00—Data processing: artificial intelligence
    • Y10S706/902—Application using ai with detail of the ai system
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S706/00—Data processing: artificial intelligence
    • Y10S706/902—Application using ai with detail of the ai system
    • Y10S706/919—Designing, planning, programming, CAD, CASE
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S706/00—Data processing: artificial intelligence
    • Y10S706/902—Application using ai with detail of the ai system
    • Y10S706/919—Designing, planning, programming, CAD, CASE
    • Y10S706/92—Simulation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S715/00—Data processing: presentation processing of document, operator interface processing, and screen saver display processing
    • Y10S715/961—Operator interface with visual structure or function dictated by intended use
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S715/00—Data processing: presentation processing of document, operator interface processing, and screen saver display processing
    • Y10S715/961—Operator interface with visual structure or function dictated by intended use
    • Y10S715/964—CAD or CAM, e.g. interactive design tools

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Evolutionary Computation (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Nursing (AREA)
  • General Physics & Mathematics (AREA)
  • Vascular Medicine (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Computer Hardware Design (AREA)
  • Reproductive Health (AREA)
  • Rehabilitation Therapy (AREA)
  • Physical Education & Sports Medicine (AREA)
  • Pain & Pain Management (AREA)
  • Epidemiology (AREA)
  • Obesity (AREA)
  • Child & Adolescent Psychology (AREA)
  • Complex Calculations (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

시스템의 하나의 노드에서 국소 좌표계 대해 자유도가 해제된 소정의 부하 및 변위 경계 조건을 따르는 시스템의 거동을 결정하기 위하여, 본 방법은 유한 요소 해석을 사용한다.
기하학적으로 비선형인 솔루션 방식에 대하여, 본 방법은 복수의 증분 단계를 가진다. 각 증분 단계에서는, 유한 요소의 노드에서 시스템의 유한 요소에 대한 변위 벡터를 계산한다. 시스템이 소정의 부하와, 소정 변위의 일부를 따를 경우에, 시스템의 에너지를 최소화함으로써, 변위 벡터의 이러한 계산이 수행된다.
기하학적으로 선형인 솔루션 방식에 대하여, 본 방법은 단일 솔빙 (solving) 단계에 의해 수행된다. 이 단계에서는, 유한 요소의 노드에서 시스템의 유한 요소에 대한 변위 벡터를 계산한다. 시스템이 소정의 부하와 소정의 변위를 따를 경우에, 시스템의 에너지를 최소화함으로써 상기 계산이 수행된다.
국소 기준 프레임에서는, 해제된 자유도를 갖는 노드에 대해 최소화 단계가 수행된다. 노드에서의 해제된 자유도는 국소 기준 프레임의 하나의 축을 따라 표현된다.
이러한 하이브리드 기준 프레임을 사용하는 것은, 제약조건을 인자화 하지 않고도, 해제되지 않은 자유도에 대한 제약조건을 전체 공식에서 간단하게 표현 및 참조할 수 있도록 한다.
솔버, 부품 설계, 자유도, 제약조건, 거동

Description

해제된 자유도를 갖는 제한적으로 변형가능한 시스템을 위한 솔버{SOLVER FOR A RESTRAINED DEFORMABLE SYSTEM WITH RELEASED DEGREES OF FREEDOM}
도 1 은 시작 컨피규레이션과 최종 컨피규레이션에서 다수의 중간 컨피규레이션을 도시하는 플렉시블 시스템의 개략도.
도 2 는 시작 컨피규레이션과 제 1 잠정 컨피규레이션에서, 플렉시블 시스템의 유한 요소 모델을 부분적으로 표현한 도면.
도 3 은 도 2 의 유한 요소의 확대도.
도 4 및 도 5 는 각종 기준 프레임을 갖는 유한 요소의 도면.
도 6 은 시스템의 노드가 자유도를 가질 경우에 사용된 프레임을 갖는 두 개의 유한 요소의 도면.
도 7 내지 도 22 는 본 방법의 각종 예를 도시한 도면.
본 발명은 컴퓨터 프로그램 및 시스템 분야에 관한 것으로, 더욱 상세하게는 부품 설계 프로그램 및 이러한 제약조건-기반 시스템 (constraints-based system) 에 관한 것이다.
부품의 설계 또는 부품의 조립을 위한 시장에서는, CATIA 라는 상표로 출원인에 의해 제공되는 것과 같은 다수의 시스템 및 프로그램이 제공되고 있다. 이러한 소위 캐드 (CAD : Computer-aided design) 시스템은 사용자가 부품 또는 부품 어셈블리의 복잡한 3 차원 (3D) 모델을 구성하고 조작할 수 있도록 한다. 이러한 시스템 및 프로그램은 모델을 정의하기 위한 다양한 제약조건을 사용한다. 제약조건의 세트 (set) 를 구하기 위해 사용되는 프로그램 또는 시스템은 일반적으로 "솔버 (solver)" 라고 불린다. 이러한 솔버는 CAD/CAM/CAE 시스템이나, 보다 일반적으로는, 임의의 종류의 대상물 (object) 을 정의하기 위한 제약조건을 사용하는 시스템에 사용된다. 설계중인 대상물과, 이들 대상물에 적용되는 제약조건의 종류에 대해 솔버가 적합하게 된다. CATIA 에서 사용되는 것과 같은 솔버는 고체 대상물의 설계에 대해 적합하게 되고, 제약조건은 치수상의 제약조건, 응력, 대상물간의 접촉 등을 포함한다.
소정의 부하 및 변위 경계 조건하의 변형가능한 고체 물리 시스템의 변형된 형상에 대한 시뮬레이션을 솔버에 의해 가능하게 할 필요성이 있다.
CAD/CAM/CAE 시스템에서 통상 설계되는 일부 플렉시블 (flexible) 시스템 및 대상물은, 별개의 지점 (point) 이 부하의 작용하에서 임의의 큰 회전을 거칠 수 있다는 의미에서, 그 높은 변형 가능성 ("기하학적 비선형성 (geometrical non-linearity)") 에 특징이 있다. 바꾸어 말하면, 이러한 플렉시블 시스템에 대하여, 시스템의 변형된 형상은 시스템의 변형되지 않은 형상과 실질적으로 상이할 수도 있다. 다른 한편으로, 현존하는 CAD/CAM/CAE 시스템에서 기하학적으로 선형 인 솔버는 대상물의 변형된 형상이 변형되지 않은 형상에 근접해 있다는 가정하에서 작동한다.
보다 일반적으로, 큰 회전을 거치게 되어, 라디언 (radian) 단위의 회전 각도 (대표적으로 5 도이지만, 이 수치는 용인되는 모델링 에러량에 따라 변동될 수 있음) 가 접선에 의해 근사화될 수 없는 시스템은 "기하학적으로 비선형인 시스템"으로 통상 간주된다. 이와 달리, 아래의 Green-Lagrange 계측:
Figure 112006024464303-pat00001
과 같은 "큰 변형률 (large strain)" 정의식이, Green-Lagrange 계측의 곱셈항(product term)
Figure 112006024464303-pat00002
을 무시할 수 없을 경우의 변형률 텐서 성분 (strain tensor component) :
Figure 112006024464303-pat00003
의 "작은 변형률 (small strain)" 정의식에 의해 대체될 수 없을 때에, 시스템이 비선형 (non-linear) 이라고 말할 수 있다.
일반적으로, 시스템의 일부 제한된 지점에서 자유도의 부분적인 해제를 허용하는 것이 솔버에서는 필요하다. 예를 들어, 하나의 방향을 따라 슬라이딩하는 것이 허용되면서 시스템 변위가 두 방향으로 제한되어 있고, 슬리이딩 방향 주위로 회전하는 것이 가능한 슬라이딩 지점을 제공할 수 있으며, 이 예에서는, 시스템이 병진 (translation) 시에 하나의 자유도를 가지고, 그 지점의 회전 시에 하나의 자 유도를 가진다. 아래에 상세하게 논의되는 바와 같이, 종래 기술 (일례 : Use of the FEM for the Design of Flexible Parts, C.A. de Hillerin, Proceedings of NAFEMS Word Congress 1999, pp. 345-356) 은 이러한 상황을 설명하기 위한 제약조건을 인자화 (factorizing) 하는 것에 대해 교시하고 있다.
비선형 시스템인 슬렌더 바디 (slender body) 의 플렉시블 시스템에 존재하는 이러한 필요성은, 시스템의 변형된 형상이 변형되지 않은 형상에 근접한 것으로 간주되는 선형 시스템에 대해서도 발생한다.
슬렌더 바디의 플렉시블 시스템의 또 다른 특징은, 유한 요소 해석에서의 시스템의 요소가 병진 및 회전 성분을 가지는 변위를 하기 쉽고, 그 요소는 소위 "구조적 (structural)" 요소이며, 구조적 요소의 다른 예는 빔 (beam) 과 쉘 (shell) 을 포함한다. 제한된 지점에서 자유도의 부분적인 해제를 허용할 필요성은 시스템의 요소에 대한 변위 벡터가 병진 성분만을 갖는 소위 "연속 (continuous)" 요소를 이용하여 모델링된 시스템에도 존재한다. 연속 시스템의 예는 막대 (bar), 막 (membrane) 및 고체 (solid) 를 포함하며, 구조적 및 연속 시스템간의 구별은 Finite Elements Procedures in Engineering Analysis, K-J. Bathe, Prentice Hall 1982 에서 논의되어 있다. 물론, 복잡한 시스템은 동일 모델에서 양쪽 타입의 요소를 포함할 수도 있다.
제약조건을 인자화 하는 종래 기술의 솔루션은, 특히, 해제된 자유도가 전체 기준 프레임의 축과 대응하지 않을 경우에 길고 복잡한 계산을 포함할 수도 있다. 따라서, 솔버는, 계산을 단순화하면서, 시스템의 일부 제한된 지점에서 자유도의 부분적인 해제를 허용할 필요가 있다.
본 발명의 일 실시 형태에 따르면, 이 문제에 대한 솔루션 (solution) 이 제공된다. 이 솔루션은 구조적이고 연속적인 요소 수식 양쪽과, 비선형 및 선형 해결 방식 양쪽에 적용될 수 있다.
연속 시스템에 대하여, 본 발명의 제 1 실시 형태는, 정적 유한 요소 솔루션의 구조 내에서, 부분적으로 제한된 노드의 전체적인 병진 자유도를 포함하는 제약조건 관계의 수가 1 또는 2 인 선형 제약조건 수식의 시스템을 처리할 수 있는 방법에 관한 것이다. 이러한 제약조건 관계는 전체 좌표축과 정렬되지 않은 하나 또는 두 개의 경사진 방향을 따르는 특정 노드의 병진 변위의 부분적인 제한을 나타낸 것이다.
이러한 상황에 대해, 이 제 1 실시 형태는 표준 << 제약조건 인자화 >> 절차를 새로운 절차로 대체하는 것을 제시하고 있다.
표준 "제약조건 인자화" 절차는, 전체 변위 성분의 항으로 기록된 쌍을 이룬 제약조건 수식의 세트를 생성하고, 모든 종속적인 자유도가 제거되어 있는 등가의 형태로 변환하며, 이 형태를 (어셈블리 이전에 전체 기준 프레임에서 회전하게 되는 기초 강성도 (elementary stiffness) 행렬과 부하 벡터로 구성되는 것으로 간주된) 평형 수식 (equilibrium equation) 의 전체 시스템에서 치환하며, 마지막으로 쌍을 이룬 제약조건 수식 세트의 영 공간 (nullspace) 상에 이 시스템을 투영함으 로써, 모든 종속적인 자유도와 미지의 라그랑주 승수 (Lagrange multiplier) 를 제거하고 평형 수식의 전체 시스템을 더욱 간결하고 가역성인 수식의 시스템으로 변환하는 것으로 구성되어 있다.
제 1 실시 형태에 따른 새로운 절차에서는, 우선, 각각의 부분적으로 제한된 노드의 병진 변위의 전체 성분이 국소 경사 제한 방향으로부터 구성된 정규직교 (orthonormal) 기준 프레임에 대한 국소 성분으로 변환된다. 국소 성분에서는, 이러한 부분적으로 제한된 노드의 병진 변위를 포함하는 제약조건 관계가 보다 간단하고 쌍을 이루지 않는 형태로 표현된다. 평형 수식의 대응하는 시스템의 구성은 어셈블리 이전에 "하이브리드 (hybrid)" 기초 강성도 행렬 및 부하 벡터의 계산을 필요로 한다.
"하이브리드" 행렬 및 벡터 회전에 의해, 즉, 상이한 기준 프레임에서의 상이한 노드에 대응하는 강성도 행렬의 대각선 블록을 선택적으로 회전시키고, 부하 벡터에 대해 유사하게 처리하여 비대각선 블록에 그 효과를 적절하게 전파시킴으로써, "하이브리드" 기초 강성도 행렬 및 부하 벡터가 얻어진다.
일단 평형 수식의 적절한 시스템이 조립되면, 쌍을 이루지 않은 제약조건 수식은 조립된 강성도 행렬 및 부하 벡터의 직접적인 검사 후의 간단한 로우 및 컬럼 제거 동작에 의해 참조된다. 결과적으로 얻어지는 수식의 시스템은, 그 미지의 부분이 유일한 기준 프레임에 대한 변위 성분이 아니기 때문에, "하이브리드"로 불린다. 그럼에도 불구하고, 모든 종속적인 자유도와 미지의 라그랑주 승수는 제거되었고, 평형 수식의 변환된 시스템은 다시 간결하고 가역성으로 된다.
그 방법은, 증분 단계의 부재가 제한이 오직 한번만 고려되어야 함을 의미하는 비선형 솔루션 뿐만 아니라 선형 솔루션에도 적용된다.
(아래의 예에서 설명되는 "케이블" 요소와 함께, 빔, 평판 및 쉘 요소를 포함하는) 구조적인 요소에 대하여, 연속 요소에서의 주요 차이점은, 구조적인 요소가 병진 노드의 자유도와 함께 회전 노드의 자유도를 나타낸다는 사실에 있다. 구조적인 요소에 대하여, 본 발명의 제 2 실시 형태는, 정적 유한 요소 솔루션의 구조 내에서, 부분적으로 제한된 노드의 전체적인 병진 자유도를 포함하는 제약조건 관계의 수가 1 또는 2, 및/또는 부분적으로 제한된 노드의 전체적인 회전 자유도를 포함하는 제약조건 관계의 수가 1 또는 2인 선형 제약조건 수식의 시스템을 처리할 수 있는 방법에 관한 것이다. 이러한 제약조건 관계는 전체적인 좌표축과 정렬되지 않은 경사진 방향을 따르는 특정된 노드의 병진 및/또는 회전 변위의 부분적인 제한을 나타낸 것이다.
이러한 상황에 대하여, 본 발명의 제 2 실시 형태는 표준 << 제약조건 인자화 >> 절차를 앞서 설명한 것과 유사한 새로운 절차로 대체하는 것을 제시하고 있다. (선형화된 회전 운동학을 갖는 구조적인 시스템이 비선형 증분 솔로션인 경우에는 업데이트를 필요로 한다는 사실과 관련된) 이러한 특수한 차이를 제외하고, 그 밖의 모든 것은 부분적인 이완을 갖는 경사 제약조건의 처리에 대한 하이브리드화 방법에 관하여 동일하다.
구체적으로, 본 발명의 제 1 실시 형태에서는, 유한 요소 해석을 이용하여, 소정의 부하 및 변위 경계 조건을 따르는 시스템의 거동 (behaviour) 을 결정하기 위한 방법이 제공되며, 상기 소정의 변위 경계 조건은, 상기 방법이 복수의 증분 단계를 포함하는 상기 시스템의 하나의 노드에서의 하나 이상의 자유도 해제를 포함하며, 각 증분 단계는, 소정의 변위의 일부를 따를 경우에 시스템의 에너지를 최소화함으로써, 유한 요소의 노드에서 시스템의 유한 요소에 대하여 변위 벡터를 계산하는 단계를 포함하며, 상기 최소화 단계는 국소 기준 프레임에서 해제된 자유도를 갖는 노드에 대해 수행되며, 상기 노드에서의 해제된 자유도는 국소 기준 프레임의 하나의 축을 따라 표현된다.
본 발명의 제 2 실시 형태에서는, 유한 요소 해석을 이용하여, 소정의 부하 및 변위 경계 조건을 따르는 시스템의 거동을 결정하기 위한 방법이 제공되며, 상기 소정의 변위 경계 조건은 상기 시스템의 하나의 노드에서의 하나 이상의 자유도 해제를 포함하며, 상기 방법은, 상기 시스템이 소정의 변위를 따를 경우에 시스템의 에너지를 최소화함으로써, 유한 요소의 노드에서 시스템의 유한 요소에 대한 변위 벡터를 계산하는 단계를 포함하며, 상기 최소화 단계는 국소 기준 프레임에서 해제된 자유도를 갖는 노드에 대해 수행되며, 상기 노드에서의 해제된 자유도는 국소 기준 프레임의 하나의 축을 따라 표현된다.
제 1 및 제 2 실시 형태에서는,
- 변위 벡터가 병진 성분만을 가지는 연속 요소에 의해 시스템이 모델링되고,
- 변위 벡터가 병진 및 회전 성분을 가지는 구조적인 요소에 의해 시스템이 모델링되며,
- 유한 요소의 노드에는 노드 기준 프레임이 제공되고, 노드에서의 변위 벡터의 회전 성분에 의해 노드 기준 프레임을 회전시킴으로써, 유한 요소의 변위된 노드에 대하여 업데이트된 노드 기준 프레임을 계산하는 것을 더 포함하며,
- 상기 방법은, 유한 요소에 대한 노드의 위치와 노드 기준 프레임에 의거하여, 시스템에 대해 평탄화된 형상을 삽입 (interpolation) 하는 단계를 더 포함하며,
- 상기 방법은 삽입된 평탄화 형상을 표시하는 단계를 더 포함하며,
- 시스템의 최소화된 에너지는, 노드에서의 노드 기준 프레임과, 노드를 포함하는 유한 요소의 기초 기준 프레임 사이의 차이의 함수이며,
- 상기 방법은, 계산 단계 중에, 유한 요소의 노드에서, 기초 기준 프레임을 노드 기준 프레임으로 변환하는 회전을 계산하고, 계산된 회전에 따라 에너지를 최소화하는 것을 더 포함하며,
- 시스템의 최소화된 에너지는, 노드에서의 노드 기준 프레임과, 노드를 포함하는 유한 요소의 기초 기준 프레임 사이의 차이에 의한 스케일링 인자의 곱의 함수이고, 상기 방법은 하나의 증분 단계로부터 또 다른 증분 단계까지 스케일링 인자를 변경하는 것을 더 포함하며, 상기 스케일링 인자는 수동으로 변경되거나, 증분 단계에서 에너지를 최소화하기 위한 수개의 반복과정에 따라 변경될 수 있으며, 그 밖에, 스케일링 인자는, 유한 요소에 대한 기초 기준 프레임과, 유한 요소에 대한 노드 기준 프레임 사이의 차이에 따라 변경될 수 있으며,
- 에너지 최소화 단계는 반복적으로 수행되고, 나머지 (residual) 와 변위 변동의 스칼라 곱을 계산하는 단계와, 스칼라 곱의 값에 따라 반복과정을 중단하는 단계를 포함하는,
이러한 추가적인 특징 중에서 하나 이상을 제공할 수 있다.
마지막으로, 본 발명은, 컴퓨터 판독가능 매체 상에 상주하여, 이러한 방법의 모든 단계를 컴퓨터상에서 실행하도록 구성된 컴퓨터 프로그램 코드를 포함하는 컴퓨터 프로그램 제품을 제공한다.
이하, 본 발명을 실시하는 시스템에 대하여, 비 한정적인 예에 의하여 첨부한 도면을 참조하여 설명한다.
도 1-5 를 참조한 슬렌더 바디의 플렉시블 시스템에 대한 설명은, 한정된 자유도의 문제를 구체적으로 처리하는 것이 아니라 종래 기술의 전체 증분 라그랑주 수식에서 동시 회전의 업데이트에 관한 것이지만, 솔버의 전반적인 동작을 이해하기에 유용하다. 명료하게 하기 위하여, NAFEMS 세계 대회 논문 (이하, NAFEMS 논문) 에서 논의된 방법을 여기서 다시 논의한다. 본 설명의 첫 번째 부분에서는, 플렉시블 시스템의 양 단부에서 소정의 기준 변위가 완전히 결정된 경우에 대해 참조한다.
본 방법의 목적은, 플렉시블 시스템이 소정의 변위 경계 조건을 따를 경우에 플렉시블 시스템의 거동을 결정하기 위한 것이다. 도 1 의 참조번호 2 는 플렉시블 시스템의 시작 컨피규레이션 (configuration) 을 나타낸다. 시스템의 양 단부에는, 기준 프레임 (4 및 6) 이 표시되어 있다. 이들 기준 프레임은 국소 기준 프레임이고, 양 단부에서 플렉시블 시스템의 방위 (orientation) 를 표시하고 있다. 따라서, 용어 "변위" 는 플렉시블 시스템의 단부의 위치를 포함할 뿐만 아니라, 플렉시블 시스템의 방위도 포함하는 것으로 이해할 수 있다. 이에 따라, 한 지점의 변위는 6개의 자유도를 포함한다.
도 1 의 참조번호 12 는 플렉시블 시스템의 최종 컨피규레이션을 나타내며, 참조번호 14 및 16 은 플렉시블 시스템의 양쪽 단부에서의 기준 프레임을 나타낸다. 제 1 단부 (기준 프레임 6 및 16) 에서, 변위는 위치의 변위를 포함하고, 기준 프레임 (6 및 16) 의 축은 평행하다. 제 2 단부 (기준 프레임 4 및 14) 에서는, 변위는 위치 및 방위에서의 변위를 포함하고, 기준 프레임 (14) 의 축은 기준 프레임 (4) 의 축에 비해 회전되어 있다.
또한, 도 1 은 플렉시블 시스템의 잠정 컨피규레이션 (20, 22 및 24) 을 도시한 것이다. 도 1 의 잠정 컨피규레이션의 수는 오로지 설명을 위해 선택된 것이다. 하나의 컨피규레이션으로부터 다음 컨피규레이션으로 통과 (이것은 하나의 증분 단계를 나타냄) 하기 위하여, 소정 변위의 일부가 플렉시블 시스템에 인가된다. 도 1 의 예에서는, 3 개의 잠정 컨피규레이션과 이로 인한 4 개의 증분 단계에 의하여, 소정 변위의 1/4 을 플렉시블 시스템에 적용함으로써 하나의 컨피규레이션으로부터 다음 컨피규레인션으로 통과하게 된다. 예를 들어, 제 1 증분 단계에서는, 병진 및 회전 상태인 소정 변위의 1/4 을 플렉시블 시스템의 단부에 적용함으로써, 시작 컨피규레이션 (2) 으로부터 제 1 잠정 컨피규레이션 (20) 으로 통과하게 된다.
플렉시블 시스템의 각 잠정 컨피규레이션과 최종 컨피규레이션은 기하학적으로 비선형인 유한 요소 해석을 이용하여 계산되며, 구체적으로, 플렉시블 시스템은 일련의 유한 요소로서 표현되어 있다. 설명을 위하여, 제 1 증분 단계가 참조되며, 이 제 1 증분 단계에서는, 시작 컨피규레이션 (2) 으로부터 제 1 잠정 컨피규레이션 (20) 으로 통과하게 된다. 도 2 는 기준 프레임 (4) 과 일부 요소를 갖는 시작 컨피규레이션 (2) 의 일부를 도시한 것이다. 수개의 요소 (30, 35, 40, 45) 가 도시되어 있고, 각 요소는 두 개의 노드 (31, 36, 41, 46 및 51) 사이에 연장되어 있다. 노드 (31) 는 실제로 플렉시블 시스템의 단부이다. 각 요소 (30, 35, 40, 45) 는 도 2 에서 32, 37, 42 및 47 로 참조부호가 붙여진 기초 기준 프레임과 연관되어 있다. 또한, 도 2 는 제 1 잠정 컨피규레이션 (20) 과 각 노드를 부분적으로 도시한 것이고, 요소와 기준 프레임은 30 으로 증가된 참조번호를 가진다.
각 증분 단계에서는, 플렉시블 시스템이 소정의 변위 부분을 따를 경우에 플렉시블 시스템의 전체 위치 에너지를 최소화한다. 플렉시블 시스템의 전체 위치 에너지는 플렉시블 시스템을 구성하는 각각의 유한 요소의 전체 위치 에너지의 합으로서 계산된다. 이것은 NAFEMS 자료에 개시된 반복과정 처리에 의해 수행된다. 이 반복과정 처리는 이하에 논의되는 바와 같이 수행될 수도 있다.
각 증분 단계에서는, 동시 회전 업데이트를 갖는 전체 라그랑주 공식 (total Lagrangian formalism) 을 사용한다. NAFEMS 자료의 섹션 4.3.2.1에서 논의된 바와 같이, 소정의 증분 단계 내의 변위 (또는 NAFEMS 자료에서 사용된 용어 중에서 "palier")는 이전 증분 단계에서 얻어진 솔루션에 대응하는 컨피규레이션에 대 해 측정된다. 바꾸어 말하면, 하나의 증분 단계에서는, 증분 단계의 초기에 플렉시블 시스템의 컨피규레이션의 요소에 연관된 기초 기준 프레임에서, 각 요소에 대하여 플렉시블 시스템의 전체 위치 에너지의 최소화가 수행된다. 도 2 의 예에서는, 증분 단계의 초기의 컨피규레이션이 시작 컨피규레이션 (2) 이다. 그러므로, 도 2 의 기초 기준 프레임 (32, 37, 42 및 47) 인 시작 컨피규레이션에 대응하는 기준 프레임에서 제 1 잠정 컨피규레이션을 얻기 위하여 플렉시블 시스템의 전체 위치 에너지를 반복적으로 최소화하게 된다.
그러나, 다음 증분 단계에 대하여, 제 1 잠정 컨피규레이션에 대응하는 업데이트된 기준 프레임을 참조하게 될 것이다. 바꾸어 말하면, 그 이후의 증분 단계는 제 1 잠정 컨피규레이션 (20) 으로부터 제 2 잠정 컨피규레이션 (22) 으로 통과하게 된다. 제 1 잠정 컨피규레이션과 연관된 기초 기준 프레임 (62, 67, 72 및 77) 에서는 제 2 증분 단계의 계산이 수행된다. 이것은 다음 증분 단계 이전에 각 증분 단계의 단부에서 새로운 기준 프레임을 계산하는 것을 필요로 한다. NAFEMS 자료는 각 요소에 대한 업데이트된 기준 프레임의 계산에 대해 더 상세한 설명을 제공하지 않는다. 바꾸어 말하면, 동시 회전 업데이트가 제시되더라도, 이 동시 회전 업데이트의 수행에 관해서는 설명이 제공되지 않는다.
본 발명은, 변위 벡터의 휨 회전 성분을 참조하지 않고도, 유한 요소의 노드에서 계산된 변위 벡터의 성분에 의거하여, 하나의 주어진 유한 요소에 대한 동시 회전 업데이트를 진행하는 것을 제시하고 있다.
본 발명의 솔루션은, 업데이트된 기초 기준 프레임에서, 하나의 축이 유한 요소의 노드와 결합하는 라인과 평행하다는 점을 보증하는 장점을 가진다. 이 특성은 각 증분 단계에서 유지되며, 모든 계산을 유사하게 하여, 보다 간단하고 용이하게 한다.
도 3 은 도 2 의 유한 요소 (40 및 70)의 확대도이다. 시작 컨피규레이션에는 노드 (41 및 46) 가 도시되어 있고, 유한 요소 (40) 는 노드 (41 및 46) 사이에 연장되어 있다. 유한 요소 (40) 의 기초 기준 프레임 (42) 은 축 x, y 및 z를 가진다. x-축은 노드 (41 및 46) 를 결합하는 라인과 평행하다.
플렉시블 시스템의 전체 위치 에너지를 최소화하는 것은 플렉시블 시스템의 각 노드에 변위 벡터를 제공한다. 변위 벡터는 병진 성분 (ux, uy, uz) 과 회전 성분 (θx, θy, θz) 을 포함한다. 병진 성분 (ux, uy, uz) 은 각 노드의 위치에서의 변위를 나타내고, 병진 성분은 용이하게 알 수 있으며, 도 3의 예에서는, 노드 (41) 내지 노드 (71) 또는 노드 (46) 내지 노드 (76) 의 위치에서의 변위에 대응한다.
회전 성분 (θx, θy, θz) 은 각 노드의 방위에서의 변위를 나타내고, 이 회전 성분은, 시스템이 플렉시블 시스템이어서 주어진 노드에서의 시스템의 방위가 하나의 증분 단계에서 변경된다는 사실을 나타내고 있다.
회전 성분은, 증분 단계의 단부에서 얻어진 컨피규레이션의 주어진 노드에서 플렉시블 시스템에 인가되는 휨 (bending) 을 나타내고 있는 휨 회전 성분 (bending rotational component) 을 포함한다.
또한, 회전 성분은 주어진 노드에서 플렉시블 시스템에 인가되는 비틀림 (torsion) 을 나타내고 있는 비틀림 회전 성분을 포함한다.
노드 (41) 에서의 회전 성분 (θx, θy, θz) 이 요소 (40) 의 기준 프레임 (42) 에서 표현되어 있으면, x-축은 노드 (41) 로부터 노드 (46)로 연장되어 있는 라인과 평행하며, 비틀림 성분은 축 x 를 따라 계측된 성분 θx 이다. 휨 회전 성분은 성분 θy 및 θz 이다. 바꾸어 말하면, 휨 회전 성분은, 유한 요소의 노드와 결합하는 라인에 대해 평행한 기초 기준 프레임의 x-축을 가로지르는 축을 따르는 노드의 회전을 나타내고 있다.
비틀림 회전 성분은 θx 이고, 기초 기준 프레임의 x-축을 따르는 노드에 인가되는 비틀림 (torsion) 을 나타내고 있다.
제 1 실시예에서는, 오로지 유한 요소의 노드에서의 변위 벡터의 병진 성분에 의거하여 업데이트된 기초 기준 프레임이 계산된다. 이 실시예는 2 차원으로만 변형하는 플렉시블 시스템에 적합하게 되어 있으며, 이러한 경우, 플렉시블 시스템은 임의의 비틀림 회전 제약조건에 종속되지 않고, 2D 평면에서만 휘어진다. 이러한 평면 변형의 예는 양 단부에서 클램프(clamp) 된 수평 빔에 의해 제공되고, 중력에 의해 작용을 받아서, 일 단부에 부과된 수평 압축 변위를 따르게 된다.
이 경우, 업데이트된 기초 기준 프레임은 아래와 같이 계산될 수도 있다. 증분 단계의 초기의 컨피규레이션에서는, 두 개의 노드 P1 및 P2 를 결합하는 유한 요소 P1P2 를 참조한다. 이 유한 요소에 대한 기초 프레임은 3개의 축 x, y 및 z 를 가진다. x-축은 두 개의 노드를 결합하는 라인 (P1P2)과 평행하다. y-축은 2D 평면의 x-축에 대해 수직이다. z-축은 2D-평면에 대해 수직이다.
(ux1, uy1, uz1) 을, 유한 요소의 제 1 노드 P1 에서 계산되어 유한 요소의 기초 기준 프레임에서 계측된 변위 벡터의 병진 성분이라고 한다. (θx1, θy1, θz1) 을, 유한 요소의 제 1 노드 P1 에서 계산된 변위 벡터의 회전 성분이라고 한다. 이와 유사하게, (ux2, uy2, uz2) 와 (θx2, θy2, θz2) 를, 유한 요소의 기초 기준 프레임에서 계측된, 제 2 노드 P2 에 대한 변위 벡터의 회전 및 병진 성분이라고 한다. P1' 및 P2' 는 증분 단계에서 계산된 변위된 노드이다.
2D 인 경우의 플렉시블 시스템은 어떠한 비틀림 상수를 가지지 않으므로, θx1 = θx2 = 0 이다. 제 1 및 제 2 노드에서의 비틀림 회전 제약조건은 제로 (zero) 이다. 플렉시블 시스템은 2D 평면 내에서 이동하여 uz1 = uz2 = 0 이 된다.
따라서, 제 1 및 제 2 노드에서의 변위 벡터의 휨 회전 성분은 θz1 및 θz2 에 의해 표현된다.
이 예에서는, 업데이트된 기초 기준 프레임이 아래와 같이 계산된다. 기초 기준 프레임의 x-축은 변위된 노드 P1' 및 P2' 를 결합하는 라인 (P1'P2') 에 대하여 평행하게 설정된다. x-축은 P1' 으로부터 P2' 의 방향으로 방위가 정해져 있다. 플렉시블 시스템이 이동하고 있는 2D 평면에서는, y-축이 x-축에 대해 수직이다. z-축은 다시 2D 평면에 대해 수직이다. 업데이트된 기초 기준 프레임을 계산하기 위하여, 변위 벡터의 회전 성분이 참조되는 것은 아니다.
제 2 실시예에서는, 유한 요소의 노드에서의 변위 벡터의 병진 성분 및 비틀림 회전 성분에 의거하여, 업데이트된 기초 기준 프레임이 계산된다. 이 실시예는 3 차원에서 변형하는 플렉시블 시스템에 적합하게 되어 있다. 이것은 예를 들어, 케이블에 적용되며, 케이블에 인가된 비틀림 제약조건은 케이블의 3D 변형으로 귀결되며, 보다 일반적으로는, 직선 중립 라인이 비-평면 (non-planar) 의 곡선으로 변형하는 임의의 상황으로 귀결된다.
노드 P1 및 P2 에 대한 변위 벡터는 제 1 실시예에서와 같이 참조되며, 이 예에서는, θx1 = θx2 = 0 또는 uz1 = uz2 = 0 이라고 가정할 이유가 없다. 업데이트된 기초 기준 프레임은 아래와 같이 계산된다. 먼저, 노드 P1 및 P2 의 비틀림 성분 θx1 및 θx2 를 참조한다. 업데이트된 기초 기준 프레임은,
- 기초 기준 프레임에 대하여, x-축을 따르는 회전, 즉, 노드 P1 및 P2 의 평균 비틀림 회전 성분 (θx1 + θx2)/2 와 동일한 각도를 갖는 라인으로서 노드 P1 및 P2 를 결합하는 라인 (P1 P2) 을 따르는 회전을 인가하고,
- 회전된 프레임에 대하여, 노드 P1 으로부터 노드 P2 로 향하는 단위 벡터를 노드 P1' 으로부터 노드 P2' 으로 향하는 단위 벡터로 변환하는 최소 회전과 동일한 회전을 인가하며, 이 회전은 변위 벡터의 병진 성분 (ux1, uy1, uz1) 및 (ux2, uy2, uz2) 에 의거하여 계산되며, 변위 벡터의 회전 성분 (θx1, θy1, θz1) 및 (θx2, θy2, θz2) 를 참조하지 않는,
위의 두 단계에 의해 계산된다.
바꾸어 말하면, 증분 단계 도중에 사용된 기초 기준 프레임은 노드에 인가되는 평균 비틀림 회전 성분에 대응하는 양만큼 일단 회전하게 되고, 다음으로, 그 회전된 프레임이 더 회전하게 되어, 업데이트된 기초 기준 프레임의 x-축이 변형된 노드 P1' 및 P2' 를 결합하는 라인 (P1'P2') 과 평행하게 된다.
제 1 실시예에서와 같이, 업데이트된 기준 프레임의 x-축은 변위된 노드 P1' 및 P2' 를 결합하는 라인과 평행하지만, 제 1 실시예와 반대로, y-축 및 z-축은 노드에 인가된 비틀림 변위를 참조하여 회전하게 되므로, xy-평면은 일정하게 유지할 이유가 없다.
두 실시예에서는, 유한 요소의 노드에서 변위 벡터의 휨 회전 성분을 참조하지 않고도 업데이트된 기초 기준 프레임이 계산된다. 또한, 두 실시예에서는, 기초 기준 프레임의 x-축이 유한 요소의 노드를 결합하는 라인에 대해 평행하게 유지된다.
위에서 논의된 방법은 각 유한 요소에 대하여 업데이트된 기준 프레임을 제공한다. 이 업데이트된 기준 프레임은 다음 증분 단계에서의 반복적인 계산을 위해 사용될 수도 있다. 기초 기준 프레임의 업데이트는 유한 요소 해석에서 사용된 근사법이 유효함을 보증한다. 구체적으로, NAFEMS 자료의 섹션 4.1.1.1 에 개시된 케이블 운동학
u(x, y, z) = u0(x) - y.θz(x) + z.θy(x)
v(x, y, z) = v0(x) - z.θy(x)
w(x, y, z) = w0(x) + y.θz(x) 는,
회전 θy 및 θz 가 충분히 작은 동안에는 유효하게 유지되므로, sinθ
Figure 112006024464303-pat00004
θ 이며, 이 가정은 기초 기준 프레임이 위에서 논의된 바와 같이 업데이트 될 경우에 유효하다.
기초 기준 프레임의 동시 회전 업데이트는 각 노드에서의 변위 벡터의 휨 성분을 참조하지 않는다. 변위 벡터의 휨 성분은 노드 기준 프레임의 업데이트 처리를 위해 참조된다.
또한, 증분 단계는 시스템의 경사-불연속 (slope-discontinuous) 표시인 플 렉시블 시스템의 유한 요소 표시를 제공한다. 플렉시블 시스템의 평탄한 그리고 경사-연속적인 표시를 얻기 위하여, 본 발명은 노드 기준 프레임의 이용을 제시하고 있다. 노드 기준 프레임은 플렉시블 시스템의 각 노드와 연관되어 있다. 노드 프레임은, 변위 벡터의 휨 회전 성분을 참조함으로써, 기초 기준 프레임과 별개로 업데이트 된다.
또한, 도 3은 노드 기준 프레임을 도시한 것이다. 시작 컨피규레이션 (2) 에서는, 노드 기준 프레임이 아래와 같이 정의된다:
- 각 노드 기준 프레임의 x-축은 플렉시블 시스템에 대해 접한다.
- 각 노드 기준 프레임의 y-축은 중립 라인을 따르는 케이블 시스템의 비틀림을 나타내고, 바꾸어 말하면, 플렉시블 시스템의 일 단부, 예를 들어, 제 1 노드 (31) 에서 y-축의 소정의 임의의 방향이 선택된다. 다음 노드에 대하여, 플렉시블 시스템에 인가되는 비틀림 제약조건을 참조하여, y-축은 y-축에 대하여 수직인 것으로 결정된다. 따라서, 플렉시블 시스템이 평면 내에서 임의의 비틀림을 가지지 않으면, 각종 노드 기준 프레임에 대한 y-축은 모두 평면에 놓여 있거나, 모두 평면에 대해 직교한다.
- 각 노드 기준 프레임의 z-축은 프레임의 모든 축이 수직이 되도록 되어 있다.
각 증분 단계에서는, 노드 기준 프레임이 아래와 같이 업데이트 된다. 노드 P1 에서의 노드 기준 프레임은 그 노드에서의 변위 벡터의 회전 성분 (θx1, θy1, θz1) 에 의한 회전에 의하여 업데이트 된다. 업데이트된 노드 기준 프레임의 원점은 변위된 노드 P1' 에 있다. 이것은, 각각의 노드 (71 및 76)에 대한 노드 기준 프레임 (73 및 78) 뿐만 아니라, 각각의 노드 (41 및 46) 에 대한 노드 기준 프레임 (43 및 48) 과 함께, 도 3 에 표시되어 있다. 또한, 도 3 은 두 개의 점을 갖는 단절된 라인에서, 플렉시블 시스템의 평탄화된 컨피규레이션을 도시한 것이다. 노드 기준 프레임의 업데이트는 변위 벡터의 휨 및 비틀림 회전 성분을 참조한다.
노드 기준 프레임을 제공하는 것은, 각 증분 단계에서, 플렉시블 시스템의 형상을 계산할 수 있도록 한다. 플렉시블 시스템의 형상은 노드의 위치와, 노드 기준 프레임에 의거하여 계산된다. 구체적으로, 플렉시블 시스템의 형상은 각종 노드에서의 노드의 위치와 노드 기준 프레임을 이용하여 삽입된다. 아래의 삽입 규칙이 사용될 수도 있다:
- 플렉시블 시스템은 각 노드를 통과해야 한다.
- 각 노드에서는, 플렉시블 시스템이 노드 기준 프레임의 x-축에 대해 접해야 한다.
- 각 노드에서는, y-축 및 z-축이 케이블에 인가되는 비틀림 제약조건을 나타내고 있다.
플렉시블 시스템의 형상을 삽입하기 위하여, 예를 들어, 큐빅 스플라인 삽입 (cubic spline interpolation) 과 같은 삽화 삽입 툴의 상태를 이용할 수도 있다.
삽입 단계는 플렉시블 시스템의 변형된 형상을 사용자에게 표시하기 이전에 수행되는 것이 바람직하다. 따라서, 사용자에게는 변형된 플렉시블 시스템의 평탄화된 도면이 제공될 수도 있고, 플렉시블 시스템의 유한 요소 도면이 제공되지 않을 수도 있다.
각 증분 단계에서는, 노드 기준 프레임으로 인해, 단계의 시작 시점에 플렉시블 시스템에 인가되는 초기 응력 (initial stress) 이 참조될 수도 있다. 위에서 설명한 바와 같이, 유한 요소의 노드에 인가되는 변위 벡터의 휨 회전 성분을 참조하지 않고도, 기초 기준 프레임이 업데이트 된다. 다른 한편으로, 노드에 인가되는 변위 벡터의 모든 회전 성분에 의거하여, 노드 기준 프레임이 업데이트 된다. 기초 기준 프레임과 노드 기준 프레임을 업데이트 한 후, 다른 노드의 한 쪽에 있는 기초 기준 프레임과 노드 기준 프레임 사이에 있는 하나의 노드에서의 각도 차이는 휨 회전 성분의 양을 나타내고 있다. 이 차이는 플렉시블 시스템에 인가되는 초기 응력의 양을 또한 나타내고 있다.
또한, 본 발명은 증분 단계의 시작 시점에 플렉시블 시스템에 인가되는 초기 응력을 계측하는 수단으로서 노드 기준 프레임을 이용하는 것에 대해 제시하고 있다. 이것은, 플렉시블 시스템의 각 유한 요소의 각 노드에서, 노드 기준 프레임을 유한 요소의 기초 기준 프레임과의 비교에 의해 수행된다. 상기 비교 (또는 노드 기준 프레임과 기초 기준 프레임 사이의 차이) 는 플렉시블 시스템에 인가되는 초기 응력을 나타내고 있다. 하나의 가능한 비교는 기초 기준 프레임을 노드 기준 프레임으로 변환하는 회전을 계산하는 것에 있다. 이 회전의 각도는 노드의 초기 응력을 나타내고 있으며, 증분 단계 도중에 최소화되는 에너지에 대응하는 기여도를 계산하기 위해 사용될 수도 있다. 노드를 포함하는 두 개의 유한 요소에 대한 노드 기준 프레임과 기초 기준 프레임 사이의 어느 하나의 주어진 노드에서의 차이를 참조할 수도 있다. 이와 달리, 주어진 요소에 대한 기초 기준 프레임과, 유한 요소를 구성하는 노드의 노드 기준 프레임 사이의 차이를 참조할 수도 있다. 두 개의 예는, 그 차이가 전체 플렉시블 시스템을 넘어서 합산될 경우, 동일한 결과가 되거나 비례하는 결과가 된다.
예를 들어, x-y 평면에서의 평면 운동에 대해, 기초 변형률 에너지로서 아래의 수식을 사용할 수도 있다:
Figure 112006024464303-pat00005
이 수식에서, 윗 첨자가 표시된 항목은, 증분 단계의 시작 시점의 노드 기준 프레임과 이에 대응하는 기초 기준 프레임간의 요소 노드에서의 각도 차이를 나타내고 있다. 윗 첨자가 표시되지 않은 항목은 현재의 반복과정 노드 변위를 나타내고 있다.
최소화되는 에너지에서, 노드 기준 프레임과 기초 기준 프레임 사이의 차이를 나타내고 있는 항목을 사용하는 것은 플렉시블 시스템의 거동을 계산함에 있어서 불안정을 발생시킬 수도 있으며, 이러한 불안정은 하나의 증분 단계로부터 다음 단계까지 플렉시블 시스템의 형상의 폭 넓은 변경으로서 나타난다. 이러한 항목을 사용하는 것은 증분 단계 도중에 수렴된 솔루션을 발견하는 것을 더욱 어렵게 할 수도 있다.
이들 문제를 회피 또는 제한하기 위하여, 대응하는 항목 또는 항목들은 스케일링 인자 μIS 로 곱해질 수도 있으며, 스케일링 인자는 0 및 1 사이의 값으로 이루어진다. 기초 변동률 에너지에 대한 예시된 수식은 아래와 같이 된다.
Figure 112006024464303-pat00006
스케일링 인자 μIS 가 1 이면, 초기 응력이 참조되며, 스케일링 인자 μIS 가 0 이면, 초기 응력이 참조되지 않는다. 스케일링 인자 μIS 는 그 결과에 따라 증분 단계 상에서 변경되는 것이 바람직하다. 수렴된 솔루션이 증분 단계 도중에 발견되거나, 플렉시블 시스템이 불안정한 것으로 판명되면, 스케일링 인자를 감소시키는 것이 전형적이다. 이것은 사용자에 의해 수행될 수도 있으며, 이와 달리, 예를 들어, 증분 단계에서의 현재의 반복과정 횟수에 따라 스케일링 인자가 자동으로 적합하게 될 수도 있다. 또 다른 예에서, 스케일링 인자는 요소의 평균 노드 기준 프레임 (그 단부 노드에서의 노드 기준 프레임의 평균) 과 그 기초 기준 프레임 사이의 차이 값에 의존하며, 이 예는, 작은 차이가 증분 단계에서의 반복과정 처리의 결과에 대한 수렴도 (convergence) 를 또한 나타내고 있다는 사실을 이용하고 있으며, 바꾸어 말하면, 차이가 작을수록, 유한 요소는 플렉시블 시스템의 거동에 대한 시뮬레이션을 더욱 잘 행한다.
도 4 및 5 는 증분 단계의 단부에서의 유한 요소에 대한 예를 도시한 것이다. 두 도면은 설명을 위하여 2D 이며, 평균 노드 프레임 (80) 뿐만 아니라, 유한 요소 (70), 노드 프레임 (73 및 78), 기초 프레임 (72) 을 도시한 것이다. 위에서 설명한 바와 같이, 평균 노드 프레임은 노드 프레임의 평균으로 계산되고, 그 원점은 노드의 중앙이며, 이에 따라, 기초 프레임과 동일한 원점이 된다. 평균 노드 프레임의 각 벡터는 노드 프레임의 대응하는 벡터의 정규화된 합이다. 도 4에서, 기초 프레임 (72) 과 평균 노드 프레임 (80) 사이의 각도 차이는 작으며, 본 예에서는 몇 도 (degree) 이다. 이것은 양호한 수렴도와, 스케일링 인자 μIS 를 올바르게 선택하였음을 나타내고 있다. 도 5에서는, 기초 프레임 (72) 과 평균 노드 프레임 (80) 사이의 각도 차이가 중요하며, 45°근처의 값은 스케일링 인자 μIS 가 너무 높은 값임을 나타내고 있다. 3 차원 예에서, 평균 노드 프레임과 기초 기준 프레임 사이의 차이는 기초 기준 프레임을 평균 노드 프레임으로 변환하는 회전 각도로서 계산될 수도 있다. 10°미만인 각도의 값은 스케일링 인자의 값이 올바르고 적절한 수렴도임을 나타내고 있다.
위에서 논의한 바와 같이, 각 증분 단계에서는, 평형 솔루션에 대한 반복 탐색을 진행한다. 각 반복과정은 2 단계 처리이며, 처리과정의 첫 번째 단계는 에너지 표면상에서의 하강 방향을 구하는 것이고, 처리과정의 두 번째 단계는 최적의 단계를 구하기 위한 하강 방향의 라인-탐색 (line-searching) 을 행하는 것이다. 최적의 단계를 구하는 것은 정지 기준을 필요로 하며, NAFEMS 논문의 섹션 4.3.2.4 는 각 반복과정에서 나머지 벡터와 전체 에너지를 평가하는 것을 제시하고 있다. 다음으로, 나머지의 놈 (norm) 과 전체 에너지 변동의 모노토니 (monotony) 에 대하여 각각의 수렴도 검사가 수행된다. 나머지의 놈과 에너지 변동의 모노토니 양자에 기초한 이 정지 기준은 유효하지만, 덜 엄격한 정지 기준을 제공하는 것이 유용할 수도 있으며, 플렉시블 시스템이 균형을 찾기 전에 불안정한 상태를 통과할 경우에 특히 유용하다.
본 발명은, 나머지와 변위 변동간의 스칼라 곱으로 정의된 나머지 에너지 놈을 정지 기준으로서 사용하는 것에 대해 제시하고 있다. 이러한 정지 기준을 사용하는 것은 시스템이 불안정한 상태를 보다 용이하게 통과할 수 있도록 한다. 이것은 시스템이 분기 (bifurcation) 를 겪게 되어 하나의 소정 컨피규레이션으로부터 매우 상이한 컨피규레이션으로 냉정하게 통과하는 경우에 특히 뚜렷하다. 제안된 정지 기준을 사용하는 것은 각 증분 단계에서 수행되는 에너지의 반복적인 최소화 과정 도중에 솔루션을 구하는 것을 보다 용이하게 하며, 따라서, 계산이 더 짧다. 이 정지 기준을 사용하는 것은 위에서 개시된 동시 회전 업데이트와 무관하다.
지금까지의 처리과정에 대한 설명에서는, 플렉시블 시스템의 양 단부에서의 소정의 기준 변위가 완전히 결정되며, 바꾸어 말하면, 플렉시블 시스템의 위치 및 방위가 완전히 제약되어 있다. 본 발명에서는, 플렉시블 시스템이 그 단부에서도 얼마간의 자유도를 갖는 경우에 대해 참조한다. 플렉시블 시스템의 하나의 예는, 플레시블 시스템의 하나의 지점, 주로 그 하나의 단부가 위치에 있어서 고정 되지만 회전은 자유로운 경우이며, 이것은 그 지점이 병진에 있어서는 고정되지만 회전에 있어서는 3 개의 자유도를 가진다는 것을 의미한다. 또 다른 예는 슬라이딩 시스템에 대한 슬라이딩 지점이며, 시스템은 하나의 방향을 따라 슬라이딩 하는 것이 자유롭고, 슬라이딩 방향 주위로 회전하는 것이 가능하며, 이 예에서는, 플렉시블 시스템이 병진에 있어서 하나의 자유도를 가지고 회전에 있어서 하나의 자유도를 가진다.
이 문제는 비선형이며 구조적인 시스템인 슬렌더 바디의 플렉시블 시스템에 존재한다. 이 문제는 상대적으로 강성 (rigid) 이거나 플렉시블 하지 않은 빔 또는 쉘과 같은 선형 구조 시스템에도 존재한다.
이 문제는 선형 또는 비선형인 연속 시스템에도 존재한다. 선형 연속 시스템은, 통상의 CAD 로 모델링되는 것, 예를 들어, 대형 비행기 설비, 엔진 캡, 실린더 블록, 크랭크 샤프트, 또는 다른 복잡한 형상의 대형 기계 부품 등과 같은 복잡한 고체 형상을 포함한다. 비선형 연속 시스템은 자동차 산업에 사용되는 고무 사일런트 블록 (rubber silent block) 또는 윈드서핑 마스트에 사용되는 조인트와 같은 대형의 변형가능한 대상물을 포함한다. 거시적인 대상물의 모델링된 변형 형상은 전체적으로 회전을 포함하는 것이 가능하지만, 각 유한 요소에 대한 변위 벡터는 병진 성분만을 가진다.
선형 시스템의 솔루션은 종래 기술에 존재한다. 이러한 솔루션은
Figure 112007078066222-pat00045
, J.F. Imbert, CEPADUES Editions, 1979-1984, page 276 또는 MSC NASTRAN Version 64, Handbook for linear Analysis, 1985, The MacNeal-Schwendel Corp., page 4.4.1, Chapter 4.4, Multipoint Constraint Operations 에 기재되어 있다. 이러한 문헌은 제약조건을 인자화 하는 것, 즉, 제약조건 또는 자유도의 일부를 다른 자유도의 함수로서 표현하는 것에 대해 제시하고 있다. 구체적으로, 소정의 선형 시스템에서의 자유도의 컬럼 벡터를 X 라고 한다. 풀어야 할 일반적인 문제는 수식 (1)과 같이 표현될 수 있다.
K.X + AT.λ = Fapp 및 A.X = b (1)
여기서, K 는 강성도 행렬, A 는 변위 제약조건의 야코비 행렬 (jacobian matrix) 이며, λ 는 라그랑주 승수의 컬럼 벡터이며, b 는 제약조건 우측의 컬럼 벡터이며, Fapp 는 시스템에 가해지는 힘의 컬럼 벡터이다 (AT 는 행렬 A 의 전치행렬을 나타냄). 종래 기술의 솔루션은 제약조건을 인자화 하는 것, 즉, 아래의 수식 (2)와 같이, 행렬 N 및 M과 함께, 자유도의 독립적인 서브세트 Xa 를 구하는 것이다.
X = N.Xa + M (2)
이것은 수식 (1)의 X 값을 대체하는 것을 가능하게 한다. A.N = 0 이라고 가정하고, NT 를 미리 곱하면, 수식 (1)은 아래의 수식과 같이 되며,
NT.K.N.Xa = NT.(Fapp - K.M)
이 수식은 종래 기술의 솔루션을 이용하여 독립적인 자유도 Xa 에서 구해질 수도 있다. (이완된 자유도의 세트를 포함하는) Xa 가 일단 구해지면, 전체 변위는 수식 (2) 에 의해 복구될 수도 있다.
자유도의 세트를 독립 및 종속 서브세트로 자동 분할하는 것이 선형의 문제일 경우에 값비싼 처리 과정이라는 점을 제외하고도, 이 솔루션은 비선형의 문제일 경우에 몇몇 추가적인 문제를 발생시킨다. 우선, 비선형 시스템에서는, 기준 프레임이 위에서 설명한 바와 같이 업데이트 될 필요가 있다. 이것은, 기준 프레임이 업데이트 될 때마다 제약조건이 다시 인자화 되어야 함을 의미하며, 바꾸어 말하면, 수식 (2)의 행렬 N 의 값은 각 증분 단계를 지난 후에 다시 계산되어야 할 필요가 있다. 다음으로, 숫자가 제한된 자릿수로 기억되는 컴퓨터 시스템에서는, 제약조건의 인자화를 위해 요구되는 행렬의 반전이 계산의 난이도를 낮출 수도 있다. 마지막으로, 강성도 행렬 K 는 각 증분 단계에서 변동하며, 이것은 NT.K.N 을 반복적으로 계산하는 것을 필요로 하며, 이것은 비용이 많이 드는 것으로 판명되어 있다.
선형 시스템에서는, 동시 회전 업데이트가 필요하지 않으며, 본 발명은 또한 제약조건을 인자화할 필요성을 회피하고 있으며, 기준 프레임의 간단한 변경만을 필요로 한다.
따라서, 기준 프레임의 업데이트와 함께, 비선형 시스템에서 제약조건이 이완되도록 한다는 요건이 있다. 본 발명은, 일부 이완되고 일부 부과된 자유도로 이루어진 제약조건을 시스템의 이러한 노드에 대한 국소 기준 프레임에 의해, 에너지를 최소화하기 위한 하이브리드 기준 프레임을 이용하는 것에 대해 제시하고 있다.
해제된 자유도가 국소 기준 프레임의 축을 따라 표현되도록 국소 기준 프레임이 선택된다. 다음으로, 제약되지 않은 자유도 Xa 의 서브세트를 구하는 것은, 이완된 (부과되지 않은) 자유도에 대응하는 라인을 수식 (1) 에서 선택함으로써 (바꾸어 말하면, 강성도 행렬로부터, 부과된 자유도, 즉, 이완되지 않은 자유도에 대응하는 라인 및 컬럼을 제거함으로써), 간단히 수행될 수도 있다.
이것은 위에서 논의한 문제에 대한 솔루션을 제공한다. 본 발명의 솔루션은 각 증분 단계에서 기준 프레임의 구체적으로 선택된 변경을 수행하는 것만을 필요로 하지만, 에너지를 최소화하기 위하여, 각 단계에서는 여하튼 기준 프레임의 변경이 필요하며, 이하 논의하는 바와 같이, 하이브리드 기준 프레임의 세트에 대한 변경은 계산의 관점에서 더 복잡한 것은 아니다.
각각의 노드 (36, 41) 및 (41, 46) 을 갖는 두 개의 유한 요소 (35 및 40) 를 도시하고 있는 도 6 을 참조하여 상기 솔루션이 예시되어 있다. 또한, 도 6은 제 1 및 제 2 유한 요소 (35 및 40) 에 대한 기초 기준 프레임 (37 및 42) 을 도시하고 있으며, 프레임 (37 및 42) 은 간단히 R1 및 R2 로 표기되어 있다. 지점 (36) 에서의 자유도의 컬럼 벡터를 X1 이라고 하고, 지점 (41 및 46) 에서의 자유도의 컬럼 벡터를 각각 X2 및 X3 이라고 한다.
자유도의 이완이 없을 경우에는, 우선, 유한 요소에 대한 기초 기준 프레임 또는 임의의 프레임에서, 각 유한 요소 (35 및 40)에 대해 강성도 행렬이 계산될 것이다. 대응하는 기초 프레임에서 각 요소에 대한 강성도 행렬이 일단 계산되면, 각종 행렬이 전체 기준 프레임에서 표현되며, 그 다음으로, 이하의 수식 (3) 내지 (6) 에서 설명된 바와 같이 합산된다.
구체적으로, 국소 프레임 R1 의 제 1 요소 (35) 에 대하여, 평형 수식 K.X = F 는 아래와 같다.
Figure 112006024464303-pat00007
(3)
이와 유사하게, 제 2 요소에 대한 대응하는 국소 프레임 R2 에서는, 제 2 유한 요소에 대해 아래의 수식을 얻으며,
Figure 112006024464303-pat00008
(4)
여기서, 각 Xk 가 6 개의 자유도를 가지면, 각 Kij 는 실제로 6 x 6 행렬이다. 수식 (3) 및 (4) 에서, 행렬에 대한 아래 첨자 1 및 2 는 행렬이 국소 기준 프레임 R1 및 R2 에서 표현되어 있음을 나타내고 있다.
제약조건이 이완되지 않으면, 기준 프레임의 변경을 수식 (3) 및 수식 (4) 에 간단하게 적용하여, 도 6 의 참조부호 82로 표시된 전체 기준 프레임 RG 에서 상기 제약조건을 모두 읽어낸다. 국소 기준 프레임 R1 로부터 전체 기준 프레임 RG 로 통과시키기 위한 행렬을 P1 로 표기하는 것으로 가정하면, 아래의 수식이 얻어지며,
Figure 112006024464303-pat00009
여기서, 아래 첨자 G 는 자유도가 전체 기준 프레임 RG 에서 표현되어 있음을 나타낸다. 따라서, 수식 (3) 은 아래와 같이 변경되며,
Figure 112006024464303-pat00010
여기서, P1 은 국소 기준 프레임 R1 로부터 전체 기준 프레임 RG 로 통과하기 위한 기준 프레임 행렬의 치환행렬이다. 따라서, 수식 (3) 은 프레임 RG 에서 아래와 같이 표현되고,
Figure 112006024464303-pat00011
(5)
수식 (4) 는 프레임 RG 에서 아래와 같이 표현되며,
Figure 112006024464303-pat00012
(6)
P2 는 국소 기준 프레임 R2 로부터 전체 기준 프레임 RG 로 통과하기 위한 기준 프레임 행렬의 치환행렬이다. 모든 유한 요소에 대해 수식 (5) 및 (6) 을 합산함으로써, 전체 평형 수식이 얻어진다. 이것은 두 개의 유한 요소만을 갖 는 도 6의 예에서는 아래의 수식을 제공한다.
Figure 112006024464303-pat00013
(6)
이러한 수식에서는 약간의 모호함이 있으며, 예를 들어, F2 는 R1 으로 표현된 제 2 노드 (41) 에 가해지는 힘과, 이와 동시에, R2 로 표현된 제 2 노드 (41) 에 가해지는 힘을 나타낸다. 이러한 모호함은 행렬 P1 및 P2 의 사용으로 인해 야기되며, 따라서, 간단하게 하기 위하여 모호함을 회피하기 위한 추가적인 아래 첨자를 사용하지 않는다. 수식 (6) 은, 전체 기준 프레임에서 계산된 바와 같이, 강성도 행렬의 형태를 도시한 것이다. 도 6 의 예는, 각각의 국소 기준 프레임에서 계산된 강성도 행렬의 합산에 의해 적어도 둘 이상의 유한 요소로 일반화된다. 기준 프레임의 필수적인 동시 회전 업데이트로 인해, 플렉시블 시스템의 각 요소에 대한 프레임의 치환행렬 Pi 에 대한 계산은 각 증분 단계에서 수행된다.
본 발명은 에너지를 최소화하기 위하여 하이브리드 기준 프레임을 사용하는 것에 대해 제시하고 있고, 일부 부과되고 일부 해제된 자유도를 갖는 노드에 대하여, 전체 프레임을 사용하는 대신에 국소 프레임을 사용한다. 자유도가 국소 프레임의 축을 따라 표현되도록 국소 프레임이 선택된다.
예를 들어, 적어도 자유도가 노드 (41) 에서 해제되어 있고, 그렇지 않으면, 고정되어 있다고 가정한다. 해제된 자유도가 하나의 축을 따라 표현되어 있는 기준 프레임을 RL 이라고 하고, 이 프레임은 도 6 에서 참조부호 83 으로 표기되어 있다. 예를 들어, 이완된 제약조건 또는 자유도가 접선이어서, 시스템이 자유롭게 하나의 방향으로 슬라이딩 한다면, 접선 노드 프레임을 국소 기준 프레임 RL 으로 사용할 수 있다. 시스템이 하나의 방향과 동일한 방향을 중심으로 한 회전을 따라 자유롭게 슬라이딩 하게 된다면, 동일한 노드 프레임이 사용될 수 있다. 이러한 경우, 행렬 X1 에 대해, 프레임의 치환행렬 P1 을 동일한 전체 기준 프레임 RG 에 사용하며, 이와 유사하게, 행렬 X3 에 대해, 프레임의 치환행렬 P2 를 동일한 전체 기준 프레임 RG 에 사용하게 된다. 그러나, 행렬 X2 에 대해서는, 프레임이 국소 프레임 RL 로 변경된다. 상기 수식 (3) 은 아래와 같이 변경되며,
Figure 112006024464303-pat00014
여기서, 적용 가능하다면, 마지막 아래 첨자는 X1 및 X2 를 표현하기 위해 사용되는 기준 프레임을 나타내고,
- P1L 은 국소 기준 프레임 R1 로부터 국소 기준 프레임 RL 로 변경하기 위한 행렬을 나타내고,
- P2L 은 국소 기준 프레임 R2 로부터 국소 기준 프레임 RL 로 변경하기 위한 행렬을 나타낸다.
국소 프레임 RL 이 국소 프레임 R1 이면, P1L 이 간단하게 단위 행렬이 될 수 있음을 알 수 있다. 이것은 프레임 R1 이 하나의 축을 따라 자유도를 표현하도록 되어 있는 경우일 수도 있다.
따라서, 하이브리드 기준 프레임에서는, 수식 (3) 및 (4)가 각각 아래와 같이 된다.
Figure 112006024464303-pat00015
(7)
및
Figure 112006024464303-pat00016
(8)
이것은 두 개의 유한 요소에 대해 아래의 강성도 행렬로 유도된다.
Figure 112006024464303-pat00017
(8)
수식 (8) 은, 제약조건이 이완되지 않은 경우에 사용되는 것과 거의 동일한 계산에 의해, 강성도 행렬과, 시스템에 인가되는 힘이 표현될 수 있음을 증명하고 있다. 유일한 차이는, 국소 기준 프레임 RL 을 구해야 하고, 제약조건이 이완되어 있는 부분적으로 제약된 노드에서 국소 기준 프레임 R1 및 R2 로부터 새로운 국소 기준 프레임으로 통과하기 위한 두 개의 추가적인 행렬 P1L 및 P2L 을 계산해야 한다는 점이다.
그러나, 일단 수식 (8) 이 계산되면, 에너지를 최소화하는 것은 위에서 주어진 수식 (2) 의 행렬 N 및 M 에 대한 계산과 같은 것을 필요로 하지 않는다. 그래서, 이완된 자유도가 기준 프레임의 하나의 축을 따라 표현되도록 국소 기준 프레임 RL 이 선택된다. 이것은, 이완되어 있는 모든 제약조건을 선택함으로써, 서브세트 Xa 의 선택이 간단하게 수행될 수도 있음을 의미한다. 바꾸어 말하면, 강성도 행렬의 라인 및 컬럼과, 제약조건이 이완되지 않은 국소 기준 프레임의 축 및 축들에 대응하는 부하 벡터의 라인을 수식 (8) 에서 간단하게 제거한다는 것이다.
구체적으로, 도 6 의 예에서는, 각 노드에 대하여 6 개의 자유도, 즉, 병진 (u, v, w) 에 있어서 3개와 회전 (θx, θy, θz) 에 있어서 3개를 참조하고 있다. 각각의 Xi 는 아래와 같다.
Figure 112006024464303-pat00018
이와 같이, 하이브리드 기준 프레임에서 표현되어 있는 수식 (8) 의 강성도 행렬은 24 x 24 행렬이다. 노드 (41) 에서 하나의 자유도가 있다고 가정하면, 그 노드는 축 u 주위로 자유롭게 회전할 수 있으며, 다시 말하면, θx2 가 이완되어 있다. 이것은, u2, v2, w2, θy2 및 θz2 가 고정되어 있고, θx2 와, X1 및 X3 의 모든 변수가 계산될 필요가 있음을 의미한다. 이러한 효과를 위하여, 도 2 및 3 을 참조하여 위에서 논의한 바와 같이, 단부 노드로부터 떨어진 노드에서의 변위의 계산을 행하도록 진행된다.
보다 구체적으로, 우선, 수식 8 은 아래와 같이 상세하게 표현된다.
Figure 112006024464303-pat00019
이 수식에서, kij (1 ≤ i ≤ 18 이고 1 ≤ j ≤ 18) 는 위에서 논의된 바와 같이 계산되는 강성도 행렬의 각종 계수이다. u, v, w, θx, θy 및 θz 에 대한 첫 번째 아래 첨자는 노드를 나타내고, u, v, w, θx, θy 및 θz 에 대한 두 번째 아래 첨자는 변위가 표현되어 있는 프레임을 나타내며, G 는 전체 기준 프레임 RG 를 대표하고, L 은 국소 기준 프레임 RL 을 대표한다.
인가된 부하와 관련하여, Fx app, Fy app 및 Fz app 는 노드에 인가되는 힘의 성분을 나타내는 한편, Mx app, My app 및 Mz app 는 노드에 인가되는 토크의 성분을 나타내고, 첫 번째 첨자 (아라비아 숫자로 된 첨자) 는 부하가 인가되는 노드를 나타내며, 두 번째 첨자는 변위가 표현되어 있는 프레임을 나타내며, G 는 전체 기준 프레임 RG 를 대표하고, L 은 국소 기준 프레임 RL 을 대표한다.
부하 좌표 시스템 RL 에 표현되어 있는 아래와 같은 5 개의 제약조건과 연관지어 참조하면,
수식은 아래와 같이 된다.
Figure 112006024464303-pat00021
일단 에너지가 최소화되면 (이에 따라, 중간 노드에서 국소 u 축을 주위의 이완된 회전을 포함하는 자유도의 Xa 세트, 즉, 다른 두 개의 노드의 모든 전체 자유도와 함께, 국소 자유도 θx2 를 얻으면), 자유도의 완전한 벡터 X 는 아래와 같 이 복구될 수도 있다.
Figure 112006024464303-pat00022
이 복구 수식에서는, 변위 벡터 성분이 국소 기준 프레임, 즉, 중간 노드에서 부분적으로 표현되어 있고, 전체 기준 프레임, 즉, 다른 두 개의 노드에서 부분적으로 표현되어 있음을 유의해야 하며, 물론, 전체 기준 프레임에서 변위 벡터를 완전히 표현하는 것이 이제는 가능하다.
Xa 벡터의 모든 성분이 계산되면, 필요에 따라, 노드 기준 프레임의 업데이트 뿐만 아니라 동시 회전 업데이트가 위에서 논의된 바와 같이 수행될 수도 있다. 유일한 차이점은, 자유도를 갖는 노드에 대해, 계산된 성분과 제약된 성분의 조합이 업데이트에서 사용된다는 것이다. 도 6 의 예에서는, u2, v2, w2, θy2 및 θz2 가 고정되어 있지만, θx2 는 X1 및 X2 의 성분과 동시에 계산된다. 기준 프레임과 노드 프레임의 업데이트를 위하여, 이러한 성분이
- 소정의 변위 경계 조건의 일부분과,
- 증분 단계의 도중에 에너지를 최소화하는 것으로부터 얻어진다는 점을 참조하지 않고, 변위 성분 (u2, v2, w2, θx2, θy2, θz2) 을 이용하여 위에서 개시된 바와 같이 진행할 수도 있다. 예를 들어, 선형 시스템이나 연속 시스템의 경우에는, 업데이트가 필요하지 않으면, 직접 다음 단계로 진행할 수도 있다.
본 발명은 시스템의 유한 요소 표시의 제약된 노드에서 일부 자유도를 이완하는 것을 가능하게 한다. 에너지가 최소화되는 프레임을 신중하게 선택하는 것은 제약조건의 인자화 회피를 가능하게 한다는 사실을 이용하고 있다.
본 방법의 예는 도 7 내지 도 22 를 참조하여 아래와 같이 주어진다. 이러한 4 개의 예는 Dassault Systemes, Suresnes France 의 CATIA V5 에 의해 통합된 비선형 솔버 코어 (FORTRAN code) 를 이용하고 있다. 각 예에 대하여, 시뮬레이션은 75 개의 증분 단계 (시작부터 초기 컨피규레이션까지인 45 개의 단계, 초기부터 최종 컨피규레이션까지인 30 개의 단계) 로 구성된다. 이들 단계 중에서 몇 개만 도면에 표시되어 있다. 모든 실행은 IBM Thinkpad T42 (1.8 GHz Pentium IV processor, 1.0 GB 의 RAM) 를 통해 수행되었다.
제 1 실시예는 도 7 내지 10 에 도시되어 있다. 시스템은 단일 세그먼트이고, 양 단부 지점은 완전히 부과되어 있다. 시스템은 12 개의 요소 모델로서 분해되어 있다. 전체 실행은 CPU 시간으로 0.28 sec 내에 수행되며, 초기-최종 단계의 각 증분 단계에서는, 평균적으로 6 개의 반복과정에 의해 수렴도가 달성된다.
노드 위치와 노드 기준 프레임이 모든 도면에 표시되어 있다. 도 7 은 시스템의 형상이 사용자에 의해 완전히 특정되어 있는 시작 컨피규레이션을 도시한 것이다.
도 8 은 초기 컨피규레이션을 도시한 것이고, 단부 지점 위치 및 방위는 사용자에 의해 특정되며, 도면에는 광 기준 프레임으로서 표시되어 있다. 시스템의 변형된 형상은 프로그램에 의해 계산된다.
도 9 는 최종 컨피규레이션을 도시한 것이고, 단부 지점 위치 및 방위는 사용자에 의해 특정되고, 두꺼운 기준 프레임으로서 표시되어 있다. 변형된 형상은 프로그램에 의해 계산된다. 도 8 및 도 9 의 단부 지점 위치 및 방위간의 차이는 시스템이 따르는 소정의 변위 경계 조건을 나타낸다.
도 10 은 3 개의 중간 컨피규레이션을 도시한 것이고, 단부 지점 위치 및 방위가 삽입되어 있으며 소정의 변위 경계 조건의 일부를 나타낸다. 변형된 형상은 프로그램에 의해 계산된다.
제 2 실시예는 도 11 내지 도 14 에 나타내어져 있다. (변형 평면에 직교하는 국소 축 주위의) 회전 자유도가 초기 및 최종 컨피규레이션 양쪽의 좌측 단부 지점에서 해제되어 있다는 점을 제외하고는, 제 1 실시예와 유사하다.
전체 실행은 CPU 시간으로 0.27 sec 내에 수행되며, 제 1 실시예에서와 같이, 평균적으로 증분 당 6 개의 반복과정에 의해 수렴도가 달성된다.
도 11-14 는 도 7-10 에 각각 대응하며, 다시 설명하지 않는다. 제 1 실시예에 관하여, 모든 도면에서 노드 위치와 노드 기준 프레임이 표현되어 있다.
제 3 예는 도 15 내지 도 18 에 나타내어져 있다. 이 예에서는, 시스템이 상이한 특징의 3 개의 세그먼트를 포함하며, 하나의 단부 지점과 하나의 중간 지점은 완전히 부과되어 있다. 이것은 도면에 나타나 있으며, 시스템의 슬렌더 부분은 제약되어 있지 않고 중력의 작용하에 아래로 휘어 있음을 명확하게 알 수 있다.
유한 요소 해석은 3 x 12 요소 모델에 의해 수행된다.
전체 실행은 CPU 시간으로 6.45 sec 내에 수행되고, 평균적으로 증분 당 40 개의 반복과정에 의해 수렴도가 달성된다.
도 15-18 은 도 7-10 에 각각 대응하며, 다시 설명하지 않는다. 노드 위치와 노드 기준 프레임이 모든 도면에 표시되어 있는 것은 아니다.
제 4 예는 도 19 내지 도 22 에 표시되어 있고, 3개의 회전 자유도가 모두 중간 지점에서 해제되어 있다는 점을 제외하고는, 이 예는 제 3 예와 동일하다.
전체 실행은 CPU 시간으로 5.2 sec 내에 수행되고, 각 증분 단계에서는, 평균적으로 증분 당 25 개의 반복과정에 의해 수렴도가 달성된다.
도 19-22 는 도 15-18 에 각각 대응하며, 다시 설명하지 않는다. 시스템에 인가되는 상이한 제약조건으로 인해 최종 컨피규레이션에서 차이가 있음에 유의해야 한다.
본 방법은 위에서 제공된 예에 한정되지 않는다. 예를 들어, 도 1 의 예 에서는, 소정의 변위가 플렉시블 시스템의 양 단부에 존재하는 것으로 간주된다. 시스템의 두 지점 이상에서 (완전히 제약된 또는 그렇지 않은) 소정의 변위가 존재한다고 간주할 수도 있다.
도 1 내지 도 3 을 참조하여 개시된 방법에서는, 각 증분 단계 이후 또는 다음 증분 단계의 시작전에 동시 회전 업데이트가 수행된다. 각 단계에서 동시 회전 업데이트를 수행하지 않는 것도 가능하며, 구체적으로, 필요한 경우, 즉, 각 기초 기준 프레임과, 이에 대응하는 업데이트 된 기초 기준 프레임 사이의 각도 변화가 소정의 값보다 높을 경우에는, 동시 회전 업데이트만 수행할 수도 있다. 실제로, 이 소정의 값은, 의사-시스템 업데이트 (quasi-systematic updating) 에 대응하는 낮은 값, 예를 들어, 1 sec 를 갖는 0 과 30 도 사이 일 수도 있다.
증분 단계 도중에 수렴도를 제어하기 위한 스케일링 인자가 위에 개시되어 있고, 이러한 스케일링 인자는 플렉시블 시스템의 거동을 제어할 목적으로 사용될 수도 있다. 예를 들어, 하나의 소정의 증분 단계에서 제로 값을 스케일링 인자에 대해 적용하는 것은 점탄성 (visco-elastic) 플렉시블 시스템을 시뮬레이션 하는 결과가 되며, 모든 초기 응력은 이완되어 있다.
상기 개시된 각종 실시 형태는 플렉시블 시스템을 더욱 효율적으로 시뮬레이션 하기 위하여 결합될 수도 있다. 각종 실시 형태를 독립적으로 수행할 수도 있다. 예를 들어, 정지 기준은 노드 프레임에 무관하게 사용될 수도 있다. 상기 논의된 정지 기준은 도 3 및 도 4 를 참조하여 논의된 동시 회전 업데이트에 무관하게 사용될 수도 있고, 정지 기준은 비구조적 (즉, 연속적) 요소에도 사용될 수 있으며, 이러한 비구조적 또는 연속적 요소는 실제로 "고체" 요소이며, 노드 회전 자유도를 가지지 않는다.
상기 논의한 바와 같이, 도 6 을 참조하여 개시된 하이브리드 기준 프레임의 사용은 이전에 논의된 동시 회전 업데이트와 무관하게 이루어질 수도 있다. 이러한 하이브리드 기준 프레임을 선형 솔버에 사용할 수도 있다.
이상과 같이 설명한 본 발명에 의하면, 계산을 단순화하면서, 시스템의 일부 제한된 지점에서 자유도의 부분적인 해제를 허용하는 솔버를 제공할 수 있다.

Claims (15)

  1. 유한 요소 해석을 이용하여, 시작 컨피규레이션에서부터 최종 컨피규레이션까지 시뮬레이션되는 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법으로서, 상기 시스템은 소정의 부하 및 변위 경계 조건을 따르고, 상기 소정의 변위 경계 조건이 상기 시스템의 하나의 노드에서 하나 이상의 자유도 해제를 포함하며,
    상기 방법은 상기 시작 컨피규레이션에서부터 상기 최종 컨피규레이션까지의 복수의 증분 단계를 포함하고,
    상기 각 증분 단계는, 상기 시스템이 소정 변위의 일부를 따를 경우에 상기 시스템의 에너지를 최소화함으로써, 유한 요소의 노드에서, 상기 시스템의 상기 유한 요소에 대한 변위 벡터를 계산하는 단계를 포함하며,
    상기 최소화하는 단계는, 국소 기준 프레임에서, 해제된 자유도를 갖는 노드에 대해 수행되며,
    상기 노드에서의 해제된 자유도는 상기 국소 기준 프레임의 하나의 축을 따라 표현되어 있고,
    상기 단계들은 컴퓨터 프로세서에 의해 실행되는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  2. 유한 요소 해석을 이용하여, 시작 컨피규레이션에서부터 최종 컨피규레이션까지 시뮬레이션되는 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법으로서, 상기 시스템은 소정의 부하 및 변위 경계 조건을 따르고, 상기 소정의 변위 경계 조건이 상기 시스템의 하나의 노드에서 하나 이상의 자유도 해제를 포함하며,
    상기 방법은 상기 시작 컨피규레이션에서부터 상기 최종 컨피규레이션까지 증분되기 위하여, 상기 시스템이 소정 변위의 일부를 따를 경우에 상기 시스템의 에너지를 최소화함으로써, 유한 요소의 노드에서, 상기 시스템의 상기 유한 요소에 대한 변위 벡터를 계산하는 단계를 포함하고,
    상기 최소화하는 단계는, 국소 기준 프레임에서, 해제된 자유도를 갖는 노드에 대해 수행되며,
    상기 노드에서의 해제된 자유도는 상기 국소 기준 프레임의 하나의 축을 따라 표현되어 있고,
    상기 단계들은 컴퓨터 프로세서에 의해 실행되는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  3. 제 1 항 또는 제 2 항에 있어서,
    상기 시스템은 상기 변위 벡터가 병진 성분만을 가지고 있는 연속 요소에 의해 모델링되는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  4. 제 1 항 또는 제 2 항에 있어서,
    상기 시스템은 상기 변위 벡터가 병진 및 회전 성분을 갖는 구조적 요소에 의해 모델링되는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  5. 제 1 항 또는 제 2 항 있어서,
    유한 요소 (40) 의 노드 (41, 46) 에는 노드 기준 프레임 (43, 48) 이 제공되고,
    상기 계산 단계는, 상기 노드 (41, 46) 에서의 변위 벡터의 회전 성분에 의해 상기 노드 기준 프레임 (43, 48) 을 회전함으로써, 상기 유한 요소의 변위된 노드 (71, 76) 에 대하여, 업데이트된 노드 기준 프레임 (73, 78) 을 계산하는 것을 더 포함하는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  6. 제 5 항에 있어서,
    상기 유한 요소에 대한 상기 노드의 위치와, 상기 노드 기준 프레임을 기준으로, 상기 시스템에 대한 평탄화된 형상을 삽입하는 단계를 더 포함하는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  7. 제 6 항에 있어서,
    상기 삽입된 평탄화 형상을 표시하는 단계를 더 포함하는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  8. 제 1 항 또는 제 2 항 있어서,
    상기 시스템의 최소화된 에너지는, 노드에서의 노드 기준 프레임과 상기 노드를 포함하는 유한 요소의 기초 기준 프레임간의 차이 함수인, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  9. 제 8 항에 있어서,
    계산 단계 도중에, 유한 요소의 노드에서, 상기 기초 기준 프레임을 상기 노드 기준 프레임으로 변환하는 회전을 계산하는 것을 더 포함하고,
    상기 최소화된 에너지는 상기 계산된 회전에 의존하는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  10. 제 8 항에 있어서,
    상기 시스템의 최소화된 에너지는, 노드에서의 노드 기준 프레임과 상기 노드를 포함하는 유한 요소의 기초 기준 프레임간의 차이와, 스케일링 인자의 곱의 함수이고, 상기 방법은 하나의 증분 단계로부터 또 다른 증분 단계로 상기 스케일링 인자를 변경시키는 단계를 더 포함하는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  11. 제 10 항에 있어서,
    상기 스케일링 인자는 수동으로 변경되는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  12. 제 10 항에 있어서,
    상기 스케일링 인자는 증분 단계에서 에너지를 최소화하는 수개의 반복과정에 따라 변경되는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  13. 제 10 항에 있어서,
    상기 스케일링 인자는 유한 요소에 대한 기초 기준 프레임과 상기 유한 요소에 대한 노드 기준 프레임간의 차이에 따라 변경되는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  14. 제 1 항 또는 제 2 항 있어서,
    상기 에너지를 최소화하는 단계는
    반복적으로 수행되고,
    나머지와 변위 변동의 스칼라 곱을 계산하는 단계와,
    상기 스칼라 곱의 값에 따라 반복과정을 정지하는 단계를 포함하는, 시스템의 거동을 결정하는 컴퓨터로 실행되는 방법.
  15. 제 1 항 또는 제 2 항에 기재된 컴퓨터로 실행되는 방법의 모든 단계를 컴퓨터상에서 실행하도록 구성된 컴퓨터 프로그램 코드 수단을 포함하는, 컴퓨터 프로그램 제품이 상주하는 컴퓨터 판독가능 매체.
KR1020060031941A 2005-04-08 2006-04-07 해제된 자유도를 갖는 제한적으로 변형가능한 시스템을위한 솔버 KR101028978B1 (ko)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
EP05290790.4 2005-04-08
EP05290790.4A EP1710721B1 (en) 2005-04-08 2005-04-08 Non-linear solver for the simulation of a flexible slender-body
EP06290450.3A EP1710723A3 (en) 2005-04-08 2006-03-21 Solver for a restrained deformable system with released degrees of freedom
EP06290450.3 2006-03-21

Publications (2)

Publication Number Publication Date
KR20060107407A KR20060107407A (ko) 2006-10-13
KR101028978B1 true KR101028978B1 (ko) 2011-04-12

Family

ID=37101477

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020060031941A KR101028978B1 (ko) 2005-04-08 2006-04-07 해제된 자유도를 갖는 제한적으로 변형가능한 시스템을위한 솔버

Country Status (4)

Country Link
US (1) US8222581B2 (ko)
JP (1) JP2006294038A (ko)
KR (1) KR101028978B1 (ko)
CA (1) CA2541948C (ko)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8050897B2 (en) * 2008-12-09 2011-11-01 Livermore Software Technology Corporation Solid finite elements suitable for simulating large deformations and/or rotations of a structure
AU2011253779A1 (en) * 2011-12-01 2013-06-20 Canon Kabushiki Kaisha Estimation of shift and small image distortion
US20170140072A1 (en) * 2014-06-17 2017-05-18 Masashi Kubota Method and system for determining a configuration of a model having a collection of entities and satisfying a set of constraints
CN105260499A (zh) * 2015-08-20 2016-01-20 南京理工大学 一种三维柔性梁系统力学的运动仿真方法及系统
US10599879B2 (en) * 2016-06-17 2020-03-24 Dassault Systemes Simulia Corp. Optimal pressure-projection method for incompressible transient and steady-state navier-stokes equations
EP3483755B1 (en) * 2017-11-09 2022-07-13 Dassault Systèmes Additive manufacturing of a 3d part
CN116090113B (zh) * 2022-10-08 2023-09-08 盐城工业职业技术学院 一种块体结构可变节点悬挂单元网格有限元模型

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5443569A (en) 1993-03-05 1995-08-22 Daiwa Can Company Foam dispensing pump container
KR20030004299A (ko) * 1999-10-15 2003-01-14 몰드플로우 피티와이. 리미티드. 구조해석 장치 및 방법
KR20060107420A (ko) * 2005-04-08 2006-10-13 다솔 시스템므 플렉시블 슬렌더-바디 시스템에 대한 비선형 솔버

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9302450D0 (en) * 1993-02-08 1993-03-24 Ibm Cumputer aided design system
US5581489A (en) * 1994-01-05 1996-12-03 Texas Instruments Incorporated Model generator for constructing and method of generating a model of an object for finite element analysis
US5901072A (en) * 1996-10-23 1999-05-04 Nelson Metal Products Corporation Method for incorporating boundary conditions into finite element analysis
US5956500A (en) * 1996-10-23 1999-09-21 Nelson Metal Products Corporation Method for incorporating boundary conditions into finite element analysis
JP3643504B2 (ja) * 1998-08-25 2005-04-27 株式会社東芝 アセンブリモデル作成方法およびアセンブリモデル作成処理プログラムを記録した記録媒体
EP1077431A1 (en) * 1999-08-16 2001-02-21 Mtu Motoren- Und Turbinen-Union MàœNchen Gmbh Generating a three-dimensional mesh for finite element analysis
US6876956B1 (en) * 1999-08-31 2005-04-05 California Institute Of Technology Method and system for thin-shell finite-element analysis
JP2001126087A (ja) * 1999-10-27 2001-05-11 Toyota Central Res & Dev Lab Inc 構造物の設計方法及び記録媒体
US6904395B1 (en) * 2000-08-16 2005-06-07 Ford Global Technologies, Llc System and method of generating a finite element mesh for a threaded fastener and joining structure assembly
JP2002231074A (ja) * 2001-01-26 2002-08-16 Furukawa Electric Co Ltd:The ワイヤーハーネスの設計方法およびこの方法をコンピュータに実行させるプログラム
US7286972B2 (en) * 2001-04-17 2007-10-23 Livermore Software Technology Corporation Implicit-explicit switching for finite element analysis
US7167816B1 (en) * 2002-02-22 2007-01-23 Livermore Software Technology Corporation Eulerian-Lagrangian mapping for finite element analysis
US7027048B2 (en) * 2002-05-31 2006-04-11 Ugs Corp. Computerized deformation analyzer
US20040054509A1 (en) * 2002-09-12 2004-03-18 Breit Stephen R. System and method for preparing a solid model for meshing
US20040122630A1 (en) * 2002-12-19 2004-06-24 Wallace Ronald Hugh Fife Method and apparatus for linking finite element models to computer-aided design models

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5443569A (en) 1993-03-05 1995-08-22 Daiwa Can Company Foam dispensing pump container
KR20030004299A (ko) * 1999-10-15 2003-01-14 몰드플로우 피티와이. 리미티드. 구조해석 장치 및 방법
KR20060107420A (ko) * 2005-04-08 2006-10-13 다솔 시스템므 플렉시블 슬렌더-바디 시스템에 대한 비선형 솔버

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
물리탐사 7권4호 2004년 한국지구물리탐사학회, 양동우 외1명, "고성능 병렬 유한요소 솔버를 이용한 3차원 주시와 진폭계산" pp.234-244(2004.)

Also Published As

Publication number Publication date
US8222581B2 (en) 2012-07-17
KR20060107407A (ko) 2006-10-13
CA2541948C (en) 2014-09-09
JP2006294038A (ja) 2006-10-26
US20070007382A1 (en) 2007-01-11
CA2541948A1 (en) 2006-10-08

Similar Documents

Publication Publication Date Title
KR100954307B1 (ko) 플렉시블 슬렌더-바디 시스템에 대한 비선형 솔버
Vu-Bac et al. A NURBS-based inverse analysis for reconstruction of nonlinear deformations of thin shell structures
KR101028978B1 (ko) 해제된 자유도를 갖는 제한적으로 변형가능한 시스템을위한 솔버
Yu et al. Variational asymptotic beam sectional analysis–an updated version
Zhang et al. Analysis of three-dimensional curved beams using isogeometric approach
Eugster et al. Director‐based beam finite elements relying on the geometrically exact beam theory formulated in skew coordinates
US10311180B2 (en) System and method of recovering Lagrange multipliers in modal dynamic analysis
Ferradi et al. A model reduction technique for beam analysis with the asymptotic expansion method
Choi et al. Higher order analysis of thin-walled beams with axially varying quadrilateral cross sections
Arredondo-Soto et al. The compliance matrix method for the kinetostatic analysis of flexure-based compliant parallel mechanisms: Conventions and general force–displacement cases
Borković et al. Geometrically exact static isogeometric analysis of an arbitrarily curved spatial Bernoulli–Euler beam
Zupan et al. The three-dimensional beam theory: finite element formulation based on curvature
Han Configurational forces and geometrically exact formulation of sliding beams in non-material domains
Pourboghrat et al. A hybrid membrane/shell method for calculating springback of anisotropic sheet metals undergoing axisymmetric loading
US6684168B1 (en) Body displacement and vibration analysis method
Kawecki et al. Numerical results quality in dependence on Abaqus plane stress elements type in big displacements compression test
Rosen et al. A transfer matrix model of large deformations of curved rods
Sopanen et al. Studies on the stiffness properties of the absolute nodal coordinate formulation for three-dimensional beams
Oñate Introduction to the finite element method for structural analysis
Dwarshuis et al. Kinematically started efficient position analysis of deformed compliant mechanisms utilizing data of standard joints
Pirmanšek et al. Material softening and strain localization in spatial geometrically exact beam finite element method with embedded discontinuity
Rungamornrat et al. Adaptive, Small‐Rotation‐Based, Corotational Technique for Analysis of 2D Nonlinear Elastic Frames
Hartloper et al. Best-fit constraint equations for coupling mixed-dimension simulation models with wide flange cross sections
dos Santos Tool descriptions and contact strategies used in the static explicit fem code ITAS3D for simulation of 3-D sheet metal forming processes
Ljukovac et al. Geometrically exact 3D beam theory with embedded strong discontinuities for modeling of localized failure in bending

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
FPAY Annual fee payment
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20150302

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20160328

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20170324

Year of fee payment: 7