• Title/Summary/Keyword: Stiffness matrix

Search Result 929, Processing Time 0.023 seconds

A Study on Composite Blade Analysis Library Development through Dimension Reduction/Recovery and Calculating Energy Release Rate (단면의 차원축소/복원해석과 에너지 해방률 계산을 위한 복합재 블레이드 해석 라이브러리 개발에 대한 연구)

  • Jang, Jun Hwan;Ahn, Sang Ho
    • Composites Research
    • /
    • v.29 no.1
    • /
    • pp.16-23
    • /
    • 2016
  • In this paper, numerical results of sectional analysis, stress recovery and energy release rate were compared with the results of VABS, 3-D FEM through the blade analysis library. The result of stress recovery analysis for one-dimensional model including the stiffness matrix is compared with stress results of three-dimensional FEM. We discuss the configuration of the blade analysis library and compare verifications of numerical analysis results of VABS. Blade analysis library through dimensional reduction and stress recovery is intended to be utilized in conjunction with pre- and post-processing of the analysis program of the composite blade, high-altitude uav's wing, wind blades and tilt rotor blade.

Performance validation and application of a mixed force-displacement loading strategy for bi-directional hybrid simulation

  • Wang, Zhen;Tan, Qiyang;Shi, Pengfei;Yang, Ge;Zhu, Siyu;Xu, Guoshan;Wu, Bin;Sun, Jianyun
    • Smart Structures and Systems
    • /
    • v.26 no.3
    • /
    • pp.373-390
    • /
    • 2020
  • Hybrid simulation (HS) is a versatile tool for structural performance evaluation under dynamic loads. Although real structural responses are often multiple-directional owing to an eccentric mass/stiffness of the structure and/or excitations not along structural major axes, few HS in this field takes into account structural responses in multiple directions. Multi-directional loading is more challenging than uni-directional loading as there is a nonlinear transformation between actuator and specimen coordinate systems, increasing the difficulty of suppressing loading error. Moreover, redundant actuators may exist in multi-directional hybrid simulations of large-scale structures, which requires the loading strategy to contain ineffective loading of multiple actuators. To address these issues, lately a new strategy was conceived for accurate reproduction of desired displacements in bi-directional hybrid simulations (BHS), which is characterized in two features, i.e., iterative displacement command updating based on the Jacobian matrix considering nonlinear geometric relationships, and force-based control for compensating ineffective forces of redundant actuators. This paper performs performance validation and application of this new mixed loading strategy. In particular, virtual BHS considering linear and nonlinear specimen models, and the diversity of actuator properties were carried out. A validation test was implemented with a steel frame specimen. A real application of this strategy to BHS on a full-scale 2-story frame specimen was performed. Studies showed that this strategy exhibited excellent tracking performance for the measured displacements of the control point and remarkable compensation for ineffective forces of the redundant actuator. This strategy was demonstrated to be capable of accurately and effectively reproducing the desired displacements in large-scale BHS.

Buckling Analysis of Axisymmetric Shells by Incremental Finite Element Mothod (증분형(增分形) 유한요소법(有限要素法)에 의한 축대칭(軸對稱) Shell구조(構造)의 좌굴해석(挫屈解析))

  • J.B.,Kim;C.Y.,Kim
    • Bulletin of the Society of Naval Architects of Korea
    • /
    • v.22 no.1
    • /
    • pp.21-30
    • /
    • 1985
  • This paper deals whth the buckling as well as postbuckling analysis of axisymmertric shells taking the initial deflection effects into account. Incremental equilibrium equations, based on the principle of virtual work, were derived by the finite element method, the successive step-by-step Newton-Raphson iterative technique was adopted. To define the transition pattern of postbuckling behavior from the prebuckling state more accurately, a simple solution method was developed, i.e. the critical load was calculated by the load extrapolation method with the determinant of tangent stiffness matrix and the equilibrium configuration in the immediate postbuckling stage was obtained by perturbation scheme and eigenvalue analysis. Degenerated isoparametric shell elements were used to analyse the axisymmetric shell of revolution. And by the method developed in this paper, the computer program applicable to the nonlinear analysis of both thin and moderately thick shells was constructed. To verify the capabilities and accuracies of the present solution method, the computed results were compared with the results of analytical solutions. These results coincided fairly well in both the small deflection and large deflection ranges. Various numerical analyses were done to show the effect of initial deflection and shape of shells on buckling load and postbuckling behavior. Futhermore, corrected directions of applied loads at every increment steps were used to determine the actual effects of large deflection in non-conservative load systems such as hydrostatic pressure load. The following conclusions can be obtained. (1) The method described in this paper was found to be both economic and effective in calculating buckling load and postbuckling behavior of shell structure. (2) Buckling and postbuckling behavior of spherical caps is critically dependent upon their geometric configuration, i.e. the shape of spherical cap and quantities of the initial deflection. (3) In the analysis of large deflection problems of shells by the incremental method, corrections of the applied load directions are needed at every incremental step to compensate the follower force effects.

  • PDF

A Four-node General Shell Element with Drilling DOFs (면내회전자유도를 갖는 4절점 곡면 쉘요소)

  • Chung, Keun-Young;Kim, Jae-Min;Lee, Eun-Haeng
    • Journal of the Earthquake Engineering Society of Korea
    • /
    • v.16 no.4
    • /
    • pp.37-52
    • /
    • 2012
  • In this study, a new 4-node general shell element with 6 DOFs per node is presented. Drilling rotational degrees of freedom are introduced by the variational principle with an independent rotation field. In formulation of the element, substitute transverse shear strain fields are used to avoid shear locking, while four nonconforming modes are applied in the in-plane displacement fields as a remedy for membrane locking. In addition, a direct modification method for nonconforming modes is employed in the numerical implementation of nonconforming modes to represent constant strain states. A 9-points integration rule is adopted for volume integration in the computation of the element stiffness matrix. With the combined use of these techniques, the developed shell element has no spurious zero energy modes, and can represent a constant strain state. Several numerical tests are carried out to evaluate the performance of the new element developed. The test results show that the behavior of the elements is satisfactory.

Parallel Nonlinear Analysis of Prestressed Concrete Frame on Cluster System (클러스터 시스템에서 프리스트레스트 콘크리트 프레임의 병렬 비선형해석)

  • 이재석;최규천
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.14 no.3
    • /
    • pp.287-298
    • /
    • 2001
  • Analysis of nonlinear behavior of prestressed concrete frame structures on PC is a time-consuming computing job if the problem size increase to a certain degree. Cluster system has emerged as one of promising computing environments due to its good extendibility, portability, and cost-effectiveness, comparing it with high-end work-stations or servers. In this paper, a parallel nonlinear analysis procedure of prestressed concrete frame structure is presented using cluster computing. Cluster system is configured with readily available pentium III class PCs under Win98 or Linux and fast ethernet. Parallel computing algorithms on element-wise processing parts including the calculation of stiffness matrix, element stresses and determination of material states, check of material failure and calculation of unbalanced loads are developed using MPL. Validity of the method is discussed through typical numerical examples. For the case of 4 node system, maximum speedup is 3.15 and 3.74 for Win98 and Linux, respectively. Important issues for the efficient use of cluster computing system based un PCs and ethernet are addressed.

  • PDF

On the local stability condition in the planar beam finite element

  • Planinc, Igor;Saje, Miran;Cas, Bojan
    • Structural Engineering and Mechanics
    • /
    • v.12 no.5
    • /
    • pp.507-526
    • /
    • 2001
  • In standard finite element algorithms, the local stability conditions are not accounted for in the formulation of the tangent stiffness matrix. As a result, the loss of the local stability is not adequately related to the onset of the global instability. The phenomenon typically arises with material-type localizations, such as shear bands and plastic hinges. This paper addresses the problem in the context of the planar, finite-strain, rate-independent, materially non-linear beam theory, although the proposed technology is in principle not limited to beam structures. A weak formulation of Reissner's finite-strain beam theory is first presented, where the pseudocurvature of the deformed axis is the only unknown function. We further derive the local stability conditions for the large deformation case, and suggest various possible combinations of the interpolation and numerical integration schemes that trigger the simultaneous loss of the local and global instabilities of a statically determined beam. For practical applications, we advice on a procedure that uses a special numerical integration rule, where interpolation nodes and integration points are equal in number, but not in locations, except for the point of the local instability, where the interpolation node and the integration point coalesce. Provided that the point of instability is an end-point of the beam-a condition often met in engineering practice-the procedure simplifies substantially; one of such algorithms uses the combination of the Lagrangian interpolation and Lobatto's integration. The present paper uses the Galerkin finite element discretization, but a conceptually similar technology could be extended to other discretization methods.

Numerical Integration based on Harmonic Oscillation and Jacobi Iteration for Efficient Simulation of Soft Objects with GPU (GPU를 활용한 고성능 연체 객체 시뮬레이션을 위한 조화진동 모델과 야코비 반복법 기반 수치 적분 기술)

  • Kang, Young-Min
    • Journal of Korea Game Society
    • /
    • v.18 no.5
    • /
    • pp.123-132
    • /
    • 2018
  • Various methods have been proposed to efficiently animate the motion of soft objects in realtime. In order to maintain the topology between the elements of the objects, it is required to employ constraint forces, which limit the size of the time steps for the numerical integration and reduce the efficiency. To tackle this, an implicit method with larger steps was proposed. However, the method is, in essence, a linear system with a large matrix, of which solution requires heavy computations. Several approximate methods have been proposed, but the approximation is obtained with an increased damping and the loss of accuracy. In this paper, new integration method based on harmonic oscillation with better stability was proposed, and it was further stabilized with the hybridization with approximate implicit method. GPU parallelism can be easily implemented for the method, and large-scale soft objects can be simulated in realtime.

Analysis of the Axisymmetric Hydro-Mechanical Deep Drawing Process by Using the Finite Element Method (유한 요소법을 이용한 축대칭 하이드로 미케니칼 디프 드로잉 공정의 해석)

  • 양동열;김한경;이항수;김경웅
    • Transactions of the Korean Society of Mechanical Engineers
    • /
    • v.16 no.5
    • /
    • pp.873-882
    • /
    • 1992
  • The study is concerned with the rigid-plastic element analysis for axisymmetric hydromechanical deep drawing in which the fluid flow influences the metal deformation. Due to the fluid pressure acting on the sheet material hydromechanical deep drawing is distinguished from the conventional deep drawing processes. In considering the pressure effect, the governing equation for fluid pressure is solved and the result is reflected on the global stiffness matrix. The solution procedure consists of two stages ; i.e., initial bulging of the sheet surface before the initiation of steady fluid flow in the flange and fluid-lubricated stage. The problem is decoupled between fluid analysis and analysis of solid deformation by deformation by iterative feedback of mutual computed results. The corresponding experiments are carried out for axisymmetric hydro-mechanical deep drawing of annealled aluminium sheet as well as for deep drawing. It has been shown from the experiments that the limit drawing ratio for hydro-mechanical deep drawing is improved as compared with deep drawing. The computed results are in good agreement with the experiment for variation of punch head and chamber pressure with respect to the punch travel and for distribution of thicknees strain. It is thus shown that the present method of analysis can be effectively applied to the analysis of axisymmetric hydro-mechanical deep drawing processes.

-An Analysis of Pre-Stressed Concrete Farn Sild by the Finite Element Method- (유한요소법에 의한 PC 농업용 사이로의 해석에 관한 연구 -제2보 탄성지반에 놓인 경우-)

  • 조진구;조현영;박병기
    • Magazine of the Korean Society of Agricultural Engineers
    • /
    • v.24 no.3
    • /
    • pp.73-83
    • /
    • 1982
  • study aims to derive a rational method for the analysis of the farm silo supported on an elastic foundation in which it is assumed that the reaction pressure of the soil at a point is proportional to the deflection at that point. In order to investigate the effects of an elastic foundation on the behaviour of the structures on it, the analysis of the farm silo resting on an elastic foundation was compared with the solution that the ground support may be assumed uniform (which was obtained from part I of this paper). To calculate the deformation of an elastic foundation, Boussinesq's solution which allows an interaction of the various parts of ground was adopted. In this case, the foundation was treated as a superparametric element additionally. In the evaluation of an element stiffness matrix, Gauss quadrature' was used. In above numerical integration, 3-point rule for the farm silo wall and the footing was introduced and 2-point rule for the evaluation of a reaction between the footing and the elastic foundation was adopted. The stresses of a farm silo on an elastic foundation were smaller than those which the distribution of contact pressure between the footing and the soil is assumed uniformly. Since the differences of stresses were remarkable in PS structures than RC structures, it is desirable that designers take into account the effect of an elastic foundation for the case of PS structures. It can be noted that while the effect of an elastic foundation was more conspicuously observed in near of the ground, the value of stresses at far from the soil was little affected by an supported soil.

  • PDF

Improved Stability Design of Plane Frame Members (평면프레임 구조의 개선된 좌굴설계)

  • Kim, Moon Young;Song, Ju Young;Kyung, Yong Soo
    • Journal of Korean Society of Steel Construction
    • /
    • v.18 no.2
    • /
    • pp.225-237
    • /
    • 2006
  • Based on the study conducted by Kim et al. (205a, b), an improved stability design method for evaluating the effective buckling lengths of beam-column members is proposed herein, using system elastic/inelastic buckling analysis and second-order elastic analysis. For this purpose, the stress-strain relationship of a column is inversely formulated from the reference load-carrying capacity proposed in design codes, so as to derive the tangent modulus of a column as a function of the slenderness ratio. The tangent stiffness matrix of a beam-column element is formulated using the so-called "stability functions," and elastic/inelastic buckling analysis Effective buckling lengths are then evaluated by extending the basic concept of a single simply-supported column to the individual members as one component of a whole frame structure. Through numerical examples of several structural systems and loading conditions, the possibilities of enhancement in stability design for frame structures are addressed by comparing their numerical results obtained when the present design method is used with those obtained when conventional stability design methods are used.