梁的有限元分析原理.pptVIP

有限元程序设计;§1. 介绍.;§2. 经典梁单元 (Bernoulli-Euler) Beam;平衡方程;最小势能原理;其中;Pj —— 集中荷载;;§3. 铁木辛柯梁理论;其中 ψ (x) 为只考虑梁弯曲理论中的线性单元转角.;内部力;实际上τxz采用以下形式：;其中k为与截面及泊松比μ相关的函数,可从弹性理论推导得到;;;;;;几种方法避免产生剪切闭锁;;;;Timoshenko 梁 （采用精确积分）;采用缩减积分;24;结构离散;坐标系; 杆系结构单元主要有铰接杆单元和梁单元两种类型。它们都只有2个节点i、j。;平面桁架杆单元（2D LINK1） 空间杆单元（3D LINK8） 平面刚架，BEAM3 空间梁单元(BEAM4) ;2-D Elastic Beam;BEAM3 is a uniaxial element with tension, compression, and bending capabilities ;31;3-D Elastic Beam;33;34;BEAM188 3-D Linear Finite Strain Beam;36;BEAM189 3-D Quadratic Finite Strain Beam;38;有限元程序设计方法简介;2、单元分析 ;总刚存贮;;方阵存贮和半带宽存贮地址关系;举例 ;例：计算下图半带宽。;（2） 变带宽存贮（一维压缩存贮）;;; 变带宽存贮：按列存贮方式。从左到右，逐列存放；对每一列，先存主对角线元素，然后由下而上顺序存放，直到顶线下第一个元素为止。为避免混淆，我们把存贮[K]的一维数组称为[A]。 实现变带宽存贮的关键问题是：总刚中元素Kij在一维数组A中的地址是什么？为此，需要知道主元Kii在A中的位置和相应列高hi。 主元位置：采用一个一维数组MAXA存主元在A中位置。 MAXA =[1,2,4,6,10,12,16,18，22]。; 从第j列的主对角线元素起到该列上方第一个非零元素为止，所含元素的个数称为第j列的列高，记为hj ；如果把第j列上方第1个非零元素的行号记为mj，则第j列的列高为 hj = j - mj + 1 其实，hj就是第j行的左带宽，因而必有 UBW= max(hj) j=1,2, …,N;1; 连接数组： 1号单元： LM=[0, 0, 1, 0, 0, 2] 2号单元： LM=[0, 0, 2, 3, 4, 5] 3号单元： LM=[3, 4, 5, 6, 7, 8] 4单元： LM=[0, 0, 1, 6, 7, 8]; a) 如果 ID（i, j ) = 0 则表明j号节点第i个自由度受有约束。;主元在一维数组[A]中的地址; MAXA(j) = h1 +…+ hj-2 - hj-1 + 1 =（h1 +…+ hj-2 + 1）+ hj-1 = MAXA（j - 1）+ hj-1 因为永远有 MAXA(1)= 1， MAXA(2)= 2 故计算主元地址的公式可写为 MAXA(j+1)= MAXA(j)+ hj （2-1） 式中， j = 2，3，…，N； hj——刚度矩阵[K]第j列的列高。;Ki,j在一维数组[A]中的地址 ;Magnetic Resonance Imaging;发生事件;MR成像基本原理;实现人体磁共振成像的条件:; 人体内的H核子可看作是自旋状态下的小星球。 自然状态下， H核进动杂乱无章，磁性相互抵消;; 三、弛豫（Relaxation） 回复“自由”的过程 ?1. 纵向弛豫（T1弛豫）： M0（MZ）的恢复 ，“量变” 高能态1 H → 低能态1 H 自旋—晶格弛豫、热弛豫;T1弛豫时间： MZ恢复到 M0的2/3所需的时间 T1愈小、M0恢复愈快 ;T2弛豫时间： MXY 丧失2/