《结构分析中的有限元法》2015-有限元习题-参考答案.docx

本科有限元习题参考答案2015年3月10日作业1、简述力学课程中介绍的各种力学模型的简化条件、基本假设和适用范围(包括有拉压杆模型、弯曲梁模型、平面应力和平面应变模型、轴对称模型、板模型、壳模型等)力学模型简化条件基本假设适用范围拉压杆线弹性，无塑性变形拉压平面假设拉伸或压缩的轴向二力杆弯曲梁弯曲为主要变形的曲杆弯曲平面假设单向受力假设弯曲为主要变形的曲杆平面应力平面应变微体所受所有外力均在同一平面；构件内一点外侧变形均在同一平面理想弹性体只有平面应力分量在仅为x,y函数；只有平面应变存在仅为x,y函数；轴对称具有对称性对称的理想弹性体弹性体几何形状、约束情况、及所受外力都对称于某一轴板板厚远小于其他两个方向的尺寸弹性体板件壳曲面薄板弹,性体壳体2、给出弹性力学问题中平衡方程、几何方程、物理方程的表达式及其意义。(1)平衡方程：二。ayzx-zx÷=+丝物理意义：应力分量与体力分量之间的关系。（2）几何方程：uvw,J=京，z=oxyOZuvvwwuy=T"+7yz+yxOzoyoxz物理意义：应变分量与位移分量之间的关系。（3）物理方程：=v-(v+j匕=4b、i(q+q)J4上一(%+%)E111y=xyvz=yzzx=Tzx物理意义：应变分量与应力分量之间的关系。3、简述最小势能原理的主要内容和主要公式。根据虚功原理得到：-uF）d-uTd=O,由=U（）贝IJmP=义（g-uF）d-urTd=O其中，口产Lg£%i%）dC-rT"即为系统的总势能，它是弹性体变形势能和外力势能之和。上面变分为零式表明：在所有区域内满足几何关系，在边界上满足给定位移条件的可能位移中，真实位移使系统的总势能取驻值（可证明此驻值为最小值）。4、“有限元法”都有哪些名称（包括中文和英文）？有限元法，也叫有限单元法、有限元素法、有限元分析；FEM（FiniteElementMethod）、FEA（FiniteElementAnalysis）5、简述有限元法的发展和现状。近几十年，伴随着计算机科学和技术的快速发展，有限元法作为工程分析的有效方法在理论、方法的研究、计算机程序的开发以及应用领域的开拓者方面均取得了根本性的发展。（1）单元的类型和形式为了扩大有限元法的应用领域，新的单元类型和形式不断涌现（等参元,梁板壳,复合材料）（2）有限元法的理论基础和离散格式将Hellinger-Reissner>Hu一WaShiZU（多场变量变分原理）应用于有限元分析，发展了混合模型、杂交型的有限元表达格式，应研究了各自的收敛条件；将加权余量法用于建立有限元的表达格式；进一步研究发展有限元解的后验误差估计和应力磨平方法。（3）有限元方程的解法（大型复杂工程结构问题一一静态，特征值，瞬态等）（4）有限元法的计算机软件（专用软件，通用软件）6、结构力学分析问题中的三种“非线性”都包含哪些，并解释其含义。非线性可以是由材料性质、变形状态和边界接触条件引起的，分别称为材料、几何、边界非线性。材料非线性就是材料的本构关系不是线性的。几何非线性时结构在载荷作用过程中产生大的位移和转动，如板壳结构的大挠度。边界非线性是指高挠度部件或由多个部件组成的结构组合件，渐进位移将会增大部件自身或是部件之间产生接触的可靠性，以此特征的特定类型几何非线性为边界条件或者接触非线性。7、简述有限元法的未来。以有限元法为代表的计算力学提出一系列新的课题：（1）为了真实地模拟新材料和新结构的行为，需要发展新的材料本构模型和单元型式。例如对于特种合金、复合材料、陶瓷材料、机敏材料、智能材料、生物材料以及纳米材料等。（2）为了分析和模拟各种类型和形式的结构在复杂载荷工况和环境作用下的全寿命过程的响应，需要发展新的数值分析方案。如：多重非线性（材料、几何、边界等）相耦合的分析方法；多场（结构、流体、热、电、化学）耦合作用的分析方法；跨时间/空间多尺度；非确定性（随机/模糊）的分析方法；自适应的分析方法。（3）有限元软件和CAD/CAM等软件系统共同集成完整的虚拟产品开发（VPD）系统。这个系统强烈影响着未来工程系统的设计、制造、和运行，主要体现在：它能提供对所设计的工程系统从加工制造到运行，直至失效和破坏的全寿命过程的更深入认识，从而能更好地识别它的属性和特征。它能够鉴定和评估所设计对象的性能和质量，并允许以最低的费用在设计过程中就对所设计的对象进行修改和优化。它能显著地缩短工程对象设计和投产的周期，降低生产成本，提高市场竞争力。2015年3月17日作业1、简述有限元法的基本思想，并结合简单结构来说明。（1）有限元法，也叫有限单元法，它的基本思想是将一个结构或连续体的求解域离散为若干个子域（单元），并通过它们边界上的结点相互联结成为组合体。（2）有限元法用每一个单元内所假设的近似函数来分片地表示全求解域内待求的未知场变量。而每个单元内的近似函数由未知函数或其导数在单元各个结点上的数值和与其对应的插值函数来表示。由于在联结相邻单元的结点上，场函数应具有相同的数值，因而将它们用作数值求解的基本未知量。这样一来，求解原来待求场函数的无穷自由度问题转换为求解场函数结点值的有限自由度问题。（3）有限元法是通过和原问题数学模型（基本方程、边界条件）等效的变分原理或加权余量法，建立求解基本未知量（场函数的结点值）的代数方程组或微分方程组。此方程组称为有限元求解方程，并表示成规范的矩阵形式。接着用数值方法求解此方程，从而得到问题的解答。2、简述结构离散（或有限元建模）的内容和要求。有限元建模的内容：1）网格划分一即把结构按一定规则分割成有限单元2）边界处理一即把作用于结构边界上约束和载荷处理为结点约束和结点载荷有限元建模的要求：1）离散结构必须与原始结构保形一单元的几何特性2）一个单元内的物理特性必须相同一单元的物理特性3、简述结点力和结点载荷的差别。结点力：单元与单元间通过结点的相互作用力。结点载荷：作用于结点上的外载荷。4、列表给出有限元几类基本单元的图形、结点数、结点自由度数和单元总自由度数（包括杆单元、梁单元、平面三角形单元、平面四边形单元、轴对称问题三角形单元、四边形壳单元、四面体单元）。单元类型单元图形结点数结点自由度杆单元AA225、写出适用于插值下表的基函数，并给出拉格朗日插值多项式。X1/31/41/(%)2-17甘N蚪,J、V4工一11o2459.z1x.,26416基函数：A.()=×-=18x4X,Ik(_)=16xXH31_1L22"43334371a(1)=2x2-x÷-OO根据拉格朗日插值多项式：Pe)=却仆)喈医学力将=2带入：p9(x)=-38x2632015年3月24日作业1、为了保证有限单元法解答的收敛性，位移函数应满足哪些条件？完备协调元、非协调元和完备元分别是什么意思？为了保证有限单元法解答的收敛性，位移函数应满足：1)位移函数必须包括单元的刚性位移(即常量项)；2)位移函数必须包括常量应变(即线性项)；3)位移函数在单元内部必须连续(连续性条件)；4)位移函数应使得相邻单元间的位移协调(协调性条件)注：上述四个条件称为有限元解收敛于真实解的充分条件；前三个条件称为必要条件。满足四个条件的位移函数构成的单元称为完备协调元；满足前三个条件的单元称为非协调元；满足前两个条件的单元称为完备元。2、如下图所示4结点平面应力单元，结点1结点4对应的结点坐标分别为(0,0),(0,1),(2,0),(2,1),结点1结点4对应的结点位移分别为(ul,vl),(u2,v2),(u3,v3),(u4,v4),试基于拉格朗日插值基函数构造如下单元的位移函数。2(0,1),4(2,1)1(0,0)3(2,0)解：根据拉格朗日插值基函数：u(x,y)=I1(x,y)ul÷2(x,y)u2+Z3(,y)u2t+4(x,y)u4I(x9y)=(X2)(y=1(-2)(y-1),(0-2)(0-1)2l2(x,y)=-(x-2)yl4(,y)=-y所以：u(x9y)=-(-2)(y-l)w1-(x-2)yu2-x(y-l)w3+；xyu4同求：v(x,y)O3、说明有限元方法解误差的主要来源？答：影响有限元解的误差：1)离散误差。边界上以直线代曲线导致离散化模型与实际物体的差异；2)位移函数误差。一般情况下单元位移函数不可能与实际单元的位移场一致；3)计算机计算误差。计算机字长的限制、相差悬殊的数值加减运算。4、说明用有限单元法解题的主要步骤。答：研究问题的力学建模；结构离散；单元分析；整体分析与求解；结果分析及后处理。5、推导基于变分原理的总势能泛函极值条件。解：有积分形式确立的标量泛函有其中是未知函数，R和七是特定的算子，。是求解域，是。的边界。"称为未知函数的泛函，随函数的变化而变化。连续介质问题的解使泛函对于微小的变化均取驻值，即泛函的“变分”等于零b=0,此为变分法。将虚功原理用于弹性变形时，总功W要包括外力功(T)和内力功(U)两部分,即：W=T-U;内力功(-U)前面有一负号，是由于弹性体在变形过程中，内力是克服变形而产生的，所有内力的方向总是与变形的方向相反,所以内力功取负值。根据虚功原理，总功等于零得：T-U=O,即外力虚功T=内力虚功U弹性力学中的虚功原理可表达为：在外力作用下处于平衡状态的弹性体，如果发生了虚位移，那么所有的外力在虚位移上的虚功(外力功)等于整个弹性体内应力在虚应变上的虚功(内力功)。根据虚功原理得到JC7b-uF)d-rurTdV=0犯二«J。(g,-uf)J-LuTTdr=o其中的n即为总势能泛函。由上面变分为零式表明：在所有区域内满足几何关系，在边界上满足给定位移条件的可能位移中，真实位移使系统的总势能取驻值(可证明此驻值为最小值)。此即总势能泛函的极值条件。2015年3月31日作业1、给出利用最小势能原理建立单元有限元静力平衡方程的一般推导过程。解：由单元位移函数：U=Nlle式中：N为插值函数(或称形函数)，得到单元内的应变和应力分别为=Be=D=DBue=Sue其中却为弹性矩阵，它完全取决于弹性常数E和将位移、应力和应变代入势能泛函有=t,r(oBDBd)wc-ue(JQVrFd)-"(NTd)根据最小势能原理，势能泛函取驻值的必要条件：=(JQBtDBd)we一(JCNFd)-(fNrTd)=O式中尸ObdC=Ke称为单元刚度矩阵，LN称为体积力等效结点力，LN7Td=6称为面力等效结点力，则上式可写成单元平衡方程：KcUe=fe此式即单元有限元静力平衡方程。2、利用最小势能原理推导结构(假设包含N个单元)有限元静力平衡方程。解：利用弹性体所有单元的总势能最小，将所有单元进行组装后即可形成整个结构的有限元静力平衡方程。为此先需要建立单元结点位移向量和总体结点位移向量之间的对应关系，即%=iu式中：叫表示第i个单元的结点位移向量；是总体结点位移向量；矩阵7；.是第，个单元的结点位移提取矩阵。对于杆单元可以表示为.O1OO/=%,%,%+JT=I八八1nlUU1U利用提取矩阵对单元矩阵进行合同变换与累加后，可得到总体的刚度矩阵K.质量矩阵Af和载荷向量产，即K=Tk,TiM=TmiTiF=Tifi=lil1=1根据最小势能变分原理即可得到m=F,此式即为总体结构的有限元静力平衡方程。2015年4月7日作业1、推导拉（压）杆单元形状函数，利用最小势能原理推导拉（压）杆单元在局部坐标系中的单元刚度方程。解：（1）如图1表示某一杆单元，现约定附属于该单元的局部坐标系（单元坐标系）为。'x'y'z',i点为原点，X轴沿着杆轴线，其正方向为由i指向j,其余各轴按右手螺旋规则确定。设%,vf.,吗,勺，vy,%为杆元结点位移分量q,Fyi,Fzi,Fxj9Fyj9%杆单元结点力分量，一律规定和坐标轴正向一致时为正。设杆的长度为/,弹性模量为E,横截面积为A。图1杆单元结点位移结点力分量对于钱接杆单元，在小变形假设的前提下，与杆垂直方向的位移并不使杆产生应变和应力。于是，对每一个结点只需考虑一个结点位移和结点力，即如图2所示Wq)勺(4)"71X图2二力杆单元二力杆单元的位移函数为(x)=q+x0式中为，牝是两个待定常数，可由3/两结点的位移唯一确定。当X=O,u(0)=ui；X=/,(/)=%时，将其代入位移函数有：ui=ai,uj=al+a2lt从而可以得到将为，出的值代入位移函数得Uj-Uiu(x)=ui+X或写成W(X)M/ IJ ujuiUi>=Nue则有u=Nue=Njui+NjujNi=I-TNj=I式中N,Nj称为i点、j点的形状函数，N成为形函数矩阵。（2）根据位移函数与应变的定义%=半得ax£=£X=处吟=(ddxl)册IlM或写成£=及/,其中5=1-11称为应变矩阵。根据应力与应变的关系=戊，将应变代入有。=E=oar=其中S=OB称为应力矩阵，对于拉（压）杆有S=彳-1lo根据杆单元在局部坐标系单元刚度矩阵的表达式Ke=fB'O&ZV代入矩阵JB=l-1与S=OB=kl1得到拉（压）杆单元在局部坐标系中的单元刚度方程为2、给出将局部坐标系中的拉（压）杆单元刚度矩阵通过坐标变换得到总体坐标系中的单元刚度矩阵推导过程。解：设o3为总体坐标系，为局部坐标系，如图1所示。图1平面杆单元总体坐标系位移在局部轴方向分量规定由总体坐标系平面X轴到局部坐标系X,轴的夹角。逆时针为正。杆单元总体坐标系下的结点位移分量用表示，局部坐标下的位移分量用表示。则平面杆单元结点i在总体坐标系和局部坐标系下的位移分量关系有ui=uicos+匕sinavi=-uisina+vfcosa同理对于/结点有uj=Ujcosa+V；sinavj=-ujsina+vycosa用矩阵表示为,fwJcosasinezulW/=5.r=ulU=i,j)vljL-smacosajvl对两结点杆单元，当用总体坐标系位移。表示局部坐标系中位移（优丫时有转换关系uiUj/,eV；xV.Uj_zjUj.vJ式中T矩阵称为坐标变换矩阵。对于图1中的杆单元就可以很容易将局部坐标系的刚度矩阵转换为总体坐标系的刚度矩阵为Ke=EA(cosaSina一SinaCoSa0000EA1000-100000CoSa-Sina00SinacosaCOSaSina00-SinaCoSa0000COSaSina00-SinaCOSacos2aCOSaSina:osasinasin2a-cos-cos2a-COSaSinac(2cosasina-sinacos:os2asi5S2aczsin-10rnaa010000-cosasin-sin2aCOSaSinCsin2aar3、计算下图所示三杆平面桁架结构的总体坐标系下各单元的单元刚度矩阵。杆单元材料弹性模量E=ZOXloAIl（N/!«八2）,截面积A=LoXloN-4）（111八2）。注：物理量单位均采用国际单位。解：总体坐标系下的单元刚度矩阵为cos2aCOSaSina-cos2a-CoSaSina.t,EACOSaSinasin2a一COSaSina-sin2aKe=-cos2a一COSaSinacos2aCOSaSina一COSaSina-sin24?CoSaSinasin2a对杆有：=0Sinal=Ocosal=1对杆有：2=90sina2=1cos2=0_2.0×10"×1.0×104000100O-1x/m71000-10O-0"2-0.30000-10100-1010000对杆有：Cr3=-arctan34sina3=354cosa3=16121612'16121612'25252525252525251291291291292.0×10,×1.0×1025254x17252530.512164T1V12162525252591291225252525-102.0x10UXLoXl(T40 00.4-1 00 0-101 0 -1 00 00 0 0 0= 5×1071 0-10 100 00 0 0 04、推导平面梁弯曲的挠曲线控制方程。解：如下图为梁单元，结点为，和J因弯曲而引起的轴向位移为(x,y)=-45,其中y是所讨论点离中面的距离,应变步一者。又根据平衡关系有巴二半，其中/是切面惯性矩,=£/dAo根据胡克定律，=Eg有yMM=-唔又根据剪力与弯矩的关系Q=也有dLvd2v-EI-2.d-MddxdrO,现在讨论的梁E,1为常数,故有E/白=小若梁上无分布剪力，则有E/9则可以判断是X的三次函数。设则梁的转角v(x)= b() +bx + b2x2 ÷ b3x3dvo(x)= = b.+ 2b)x+3bx dx将i, /结点位移代入上式有匕=4，4=瓦vy=+ W2+ %=b+2bj + 3bj2用矩阵表示为Vi1004010Vj卜=0/I2%012100 3 3/2b b2从而得1000"%4b2=0_3_1_203F07匕2.F1F2F1F.q1000"v(x)=l,x,x2,0_21_203F0_1"7%、4Vj'=Njivj2_2_13为,FFF7.单元的形状函数矩阵M,为NV=3x22x32x2X33x22x3WE+T,x+F,TT,"+F=%H°jH1.其中rr13x22x3Hoi=I-J+J'设g=则上式可写成/=1-3+2昇Hii=l(-22-y)aH.j=32-2Hy=(3-2)以上四个函数即是二结点梁单元的形状函数.由于每结点有两个位移参数，则每结点有两个形状函数。"o'是指i结点零阶导数(即i点位移)对应的形函数,Hli指i结点一阶导数(即，结点转角)对应的形函数，位移函数用内插多项式表示为U(X)=HoM+Hhq+H0jvj+Hjj=Nv<2015年4月14日作业1、推导平面直梁单元形状函数，利用最小势能原理推导平面直梁单元在局部坐标系中的单元刚度方程。解：先将轴力与剪力、弯矩分开考虑。对于平面直梁单元的两个节点1,2各有一个挠度和转角，故可设其挠度为：vv(x)=bii+blx+b2x2+b3x3=V1(x)w1+N2(x)l+N3(x)w2+TV4(x)=Nuw其中，4(i=0,l,2,3)为待定常数，Nj(X)(i=l,2,3,4)为单元形函数,N.,.=N1(X)N2(X)TV3(X)N4(x)t;w=4b2Ar利用Hermit插值或书中解方程的方法，可求得：7V1(x)= 1-3 +2("，N2(x)= x-2LN3(X)对于轴力引起的位移w*),设:(x)=a0+alx=NljXX1-7,7吗U2W24,进而，位移函数改写为Yl-OOONTN2位移函数求得后，可得到应变和应力的表达式：若忽略剪切影响，分是拉力引起的应变，4是弯曲引起的应变-OO-nf62x(4,6xnLU23)UI2)对应的应力为=<=0,下面求梁单元的总势能。梁应变能为：U=-7JV=-(we)rB,EBuedV2°2°=-()7BEB6Vmc=-(wc)7KV2。2梁上的结点、力F'=MQ,Mj,Nj,Qj,MjTV=-(we)rFf-v(x)q(x)dx总势能=U+V=g(1fF'=g()TKV-(wc)rFc-(wt)rQ式中。=(Nz(X)公是分布载荷的等效结由最小势能原理，e=5ue(Kewe由施IC的任意性，知KV=F°°<9,/、/,、4=Bu(62x(26x)”，"lF"Tj77JJ叱2=Ez=EBut，并有分布力q(x)作用。故外力势为:点栽荷,其中丽舄-Fe-Qe)=Oe+Qe其中刚度矩阵e=EjiBrBdV将B带入上式，得到EAOOEAOOLLO12E/6EIO12E/6EIL3TL3O6EI4F/O6E12E1HTLL2LEAOOEAOOLLO12E/6EIO12EI6EIE13Z?L2O6EI2EIO6EI4E/Lj2L2、计算下表所示直梁单元承受的4种典型载荷分布形式的等效结点载荷(要求计算过程!)表1几种典型载荷分布形式解：等效节点载荷的基本原理是功的互等原理，即对任意虚位移，等效节点力做的功与实际力做的功相等。对于分布载荷夕(X),可推得等效节点载荷的计算公式(书中P28)为:(1)flx)qxdx其中N7(x)为单元形函数向量。对于梁单元:Nl(X)= 1 3 nN2(x)= x-2LXN3(X)M(X)=T-L)设等效节点载荷为方=以MQ2M27(1)当梁发生虚位移为形函数TV1(X)= 1-3+2()时，由形函数定义知，等效节点载荷中只有Ql做功，由功的互等原理，有:故Q=PQ-TV1(O)=P.2V1(a)(1一a)2(L+2a)_出?(3a+力)Lj7同理2JV3(O)=pv3();M1JV2(O)=P-TV2(«);M2JV4(O)=PJV4();求得：CPZ(、Pa2(a+3b)Q2=PN,(a)=p；D八2陷二PM3)=%;M2=PN4(G=-等；对任意虚位移，等效节点力做的功与实际力做的功相等。而对于梁单元，所有容许的位移都可表示为单元形函数的线性组合。所以，当每个形函数作为虚位移，等效节点力与实际力做功相等时，对任意容许位移，等效节点力与实际力做功也必相等。满足功的互等。(2)将夕(X)=4带入公式(1),求得等效节点载荷为:储qL212qL2ql)12M一20处30这20卫20(4)将q(X)= <%+ g;(-<x0)L 2带入公式(1),一千+Z(0<x 勺求得等效节点载荷为:2(V+ 2 -l)QMlQ2此.2x+ 4)x2l1/23xY+L /、.3-2-l)(x-L-L)dx-qL45ql96qL4 5ql96（3）将（7（x）="带入公式（1）,求得等效节点载荷为:3、计算如下所示受一个集中力和一个集中力矩的直梁单元的等效结点载荷。解法一：仿照2中(1)的方法，设等效节点载荷列向量为f；=以MQ2m2o当梁发生虚位移为N(X)=I+：荷中只有0做功，由功的互等原理，彳)+m普Ie同样地，Q=PM(4)+M簧LM=P小"L%=PM()+M鬻Z解法二：利用狄拉克函数5。)的性质，入第2题公式(1)中，得：°x-2lf-l+f-M=P%*-)1"画l)l)-Lf-Y+f-YLUJS=PN()+MN'S)进而可得出与解法一相同的结果。2015年4月21日作业2时，由形函数定义知，等效节点载：P(L-a)2(L+2a)m6b(b-L)L3L3Pa2(3L-2a)Mb(Lb)L31e,Pah2MGb-L)(b-L)"L21U'Pa2bb(3b-2L)L2,L2'可设夕(X)=P5(x-4)+f5'(x-6),带,-F<'、x-2f-l+f-ldx+Mx-bY：公3一2L)l)-<r÷<)/2Z1、请给出由单元刚度矩阵组装形成总体刚度矩阵的基本原理和过程。(补充公式推导)O* 0Oo12(3.4)(3,92)图3,123杆平IE桁姿馅梅*jMl ku 5此心此J 6Of若3夕原UM¼诙鼻粽产构，*，必度知少R的W码布扬Fi f K；Jq I："位向同曳。集发结构检体刚度用汰除入端/ 的阿过矩阵K .然后料单元r的定位向耻分别写住;元刖收臣府, - '柑部,阵,设每个杆单元的单元刚度矩阵分别如F；345怀tn AJt Au Kq3旗铲M,求出单元*住初飞保率IJb.K的上方利右恻，戊左奥：VOL 8 -0 o* o *O 1 ,' o J岛G）e-W>r,U：(w,4”.仆,展把中足杆端位移tt<w1）篇翻就扁需除东*方*.它成为电元的定位向后.利用他元定位向it1的句行或拉列就J或元定位向域的一个分心相对应.这个分K即为黑中加小的'利列在结构总体刚度短阵K中的行码或刖码卜必.按明由单元定付向附中分出的行科相刚码,就加够:将弟元刚度知阵的元ICiE确地势加到结构刚度拉阵K中去.如图3.I所小.：图3.11直接刚度法集成示虐图考感如图3.12所示的平面桁架结构，图中.1（1,2）、2（3,4）、3（5.6）中括号前面的数字代表倍体结点编号括号中的数字即为位移向股的,对于学年杆单元有MF关系式"（2：1,2，）。孝寺七人屉由前面知.f*乐）杆单元刚度矩阵工2个4阶知竹噌也10-2,4G/为）基本原理：结构总体刚度矩阵的元素是由单元刚度矩阵的元素组成的，只要确定了单元刚度矩阵各元素在结构总体刚度矩阵中的位置，就可以由单元刚度矩阵直接集成结构总体刚度矩阵。基本思想是：把单元杆端位移分量（局部码）所对应的结构结点位移向量的序号（整体位移码）组成一向量，它成为单元的定位向量。利用单元定位向量可以完全确定单元刚度矩阵的每个元素在结构（原始）刚度矩阵中的行码和列码。过程：先求出单元在整体坐标系中的刚度矩阵，然后将单元的定位向量分别写在单元刚度矩阵的上方和右侧（或左侧）。这样，的每一行或每一列就与单元定位向量的一个分量相对应，这个分量即为中相应的行和列在结构刚度矩阵中的行码或列码。于是，按照由单元定位向量中分量的行码和列码，就能够将单元刚度矩阵的元素正确地叠加到结构刚度矩阵中去。2、计算下图所示三杆平面桁架结构的总体刚度矩阵。杆单元材料弹性模量E=20×10,1Vm2,截面积A=LOxlO-4W2o注：物理量单位均采用国际单位。解：先计算各单元相对总体坐标系的刚度矩阵。各单元局部坐标系的选取与X到X的角度如下表:单元编号局部坐标系原点cosSina201022013-arctan(4)0.8-0.6由书中P20公式(3.31)得到各单元的刚度矩阵为:10-1O-0000EA0000EA010-12-F-1010000000000-1010.64-0.48-0.640.48-_EA-0.480.360.48-0.36(/=0.Itw)3-1=V-0.640.480.64-0.480.48-0.36-0.48035组装得到总刚:0.378-0.096-0.250-0.1280.096-0.0960.072000.0960.072K=2×107×-0.2500.250000000.3330-0.333-0.1280.096000.128-0.0960.0960.0720-0.333-0.0960.4053、从半带宽最小的角度出发，如下图所示平面桁架或刚架的结点编号哪一种最好，其“半带宽B”是多少？解：对于桁架计算半带宽：(a)单元节点号之差最大的为3(如节点1和4),节点自由度数为2半带宽：8=(3+l)x2=8(b)单元节点号之差最大的为9(如节点1和10),节点自由度数为2半带宽：B=(9+l)×2=20(c)单元节点号之差最大的为5(如节点3和8),节点自由度数为2半带宽：B=(5+l)×2=12对于桁架计算半带宽：(d)单元节点号之差最大的为3(如节点1和4),节点自由度数为3半带宽：B=(3+l)×3=12(e)单元节点号之差最大的为9(如节点1和10),节点自由度数为3半带宽：B=(9+l)×3=30(f)单元节点号之差最大的为5(如节点3和8),节点自由度数为3半带宽：8=(5+l)x3=18由上可知，(a)编号最好2015年4月28日作业1、计算下图所示三杆平面桁架结构的结点位移、各杆轴力、轴向正应力及支座约束反力。1结点处受X和y方向的载荷大小为Px=Py=IOooo.O(N)。杆单元材料弹性模量E=2.0×10,TVm2,截面积A=LoXlOTm?。解：先计算各单元相对总体坐标系的刚度矩阵。各单元局部坐标系的选取与X到的角度如下表:单元编号局部坐标系原点OtCOSCtsin20102/2013-arctan(4)0.8-0.6由书中P20公式(3.31)得到各单元的刚度矩阵为:10-1O-0000"EA0000K一坐010-1K=Tr-10103/000000000-1010.64-0.48-0.640.48-r,EA-0.480.360.48-0.36K3.=(Z=0.1/72)3,5/-0.640.480.64-0.480.48-0.36-0.48035将单元刚度矩阵组成总刚，各节点X、y方向上对应的位移依次设为H1,M2,W3,W4,W5,M6,引入边界条件的=%=%=5=0，用删行删列法处理，删去第2至第5行与第2至第5歹U，总刚度方程变为：EA0.3780.096lwl1_RrlT.O960.405jmJ0J(0.1407)=mmIk-0.0334)下面求各杆轴力、轴向正应力。对于号杆，局部坐标系中=0,w11=w1=0.1407/ZZW,带入局部单元刚度矩阵，轴向应力Cru.,N,Nug*-/T2.0x10"×0.1407p,=E=E=E(U,L+W1=E-=Pa=70.35MP«1x-xlx4/400轴力为：F.=iA=70.35×100V=7.035皿同理可求解号杆：轴向应力：6=E生二包二22.27MPo23/轴力：F2=2A=2.22