(2010) 基于有限元模型修正的施工参数识别
基于有限元模型修正的施工参数识别
作者：陈宇， CHEN Yu
作者单位：中铁大桥勘测设计院有限公司,湖北武汉,430050
刊名：
桥梁建设
英文刊名：BRIDGE CONSTRUCTION
年，卷(期)：2010(5)
参考文献(5条)
1.徐君兰大跨度桥梁施工控制 2000
2.Jungwon Huh;Achintya Haldar Stochastic Finite-Element-Based Seismic Risk of Nonlinear Structures [外文期刊] 2001(03)
3.陈宇有限元模型修正技术在桥梁工程中的应用研究 2006
4.薛嘉庆最优化原理与方法 1992
5.向天宇;赵人达;刘海波基于静力测试数据的预应力混凝土连续梁结构损伤识别[期刊论文]-土木工程学报2003(11)
本文链接：https://www.360docs.net/doc/b618372818.html,/Periodical_qljs201005009.aspx
等效电路模型参数在线辨识
第四章 等效电路模型参数在线辨识 通过第三章函数拟合的方法可以确定钒电池等效电路模型中的参数，但是在实际运行过程中模型参数随着工作环境温度、充放电循环次数、SOC 等因素发生变化，根据离线试验数据计算得到的参数值估算电池SOC 可能会造成较大的估计误差。因此，在实际运行时，应对钒电池等效电路模型参数进行在线辨识，做出实时修正，提高基于模型估算SOC 的精度。 4.1 基于遗忘因子的最小二乘算法 参数辨识是根据被测系统的输入输出来，通过一定的算法，获得让模型输出值尽量接近系统实际输出值的模型参数估计值。根据能否实时辨识系统的模型参数，可以将常用的参数辨识方法分为离线和在线两类，离线辨识只能在数据采集完成后进行，不能对系统模型实时地在线调整参数，对于具有非线性特性的电池系统往往不能得到满意的辨识结果；在线辨识方法一般能够根据实时采集到的数据对系统模型进行辨识，在线调整系统模型参数。常用的辨识方法有最小二乘法、极大似然估计法和Kalman 滤波法等。因最小二乘法原理简明、收敛较快、容易理解和掌握、方便编程实现等特点，在进行电池模型参数辨识时采用了效果较好的含遗忘因子的递推最小二乘法。 4.1.1 批处理最小二乘法简介 假设被辨识的系统模型： 12121212()()()1n n n n b z b z b z y z G z u z a z a z a z ------+++==++++L L （4-1） 其相应的差分方程为： 1 1 ()()()n n i i i i y k a y k i b u k i ===--+-∑∑（4-2） 若考虑被辨识系统或观测信息中含有噪声，则被辨识模型式（4-2）可改写为： 1 1 ()()()()n n i i i i z k a y k i b u k i v k ===--+-+∑∑（4-3） 式中， ()z k 为系统输出量的第k 次观测值；()y k 为系统输出量的第k 次真值，()y k i -为系统输出量的第k i -次真值；()u k 为系统的第k 个输入值，()u k i -为 系统的第k i -个输入值；()v k 为均值为0的随机噪声。
有限元复习题答案
1、何为有限元法？其基本思想是什么？ 有限元法是一种基于变分法而发展起来的求解微分方程的数值计算方法，该方法以计算机为手段，采用分片近似，进而逼近整体的研究思想求解物理问题。 基本思想是化整为零集零为整。 2、为什么说有限元法是近似的方法，体现在哪里？ 有两点：用离散单元的组合体来逼近原始结构，体现了几何上的近似；而用近似函数逼近未知变量在单元内的真实解，体现了数学上的近似。 3、单元、节点的概念？ 节点：表达实际结构几何对象之间相互连接方式的概念 单元：网格划分中的每一个小部分称为单元，网格间相互联结点称为节点 4、有限元法分析过程可归纳为几个步骤？ 结构离散化、单元分析、整体分析 5、有限元方法分几种？本课程讲授的是哪一种？ 位移法、力法、混合法本课程讲授位移法 6、弹性力学的基本变量是什么？何为几何方程、物理方程及虚功方程？弹性矩阵的特点？ 弹性力学变量:外力、应力、应变和位移。 描述弹性体应变分量与位移分量之间的方程称为几何方程；物理方程描述应力分量与应变分量之间的关系；弹性体上外力在虚位移发生过程中所做的虚功与储存在弹性体内的需应变能相等。 弹性矩阵由材料的弹性模量和泊松比确定，与坐标位置无关。 7、何为平面应力问题和平面应变问题？ 平面应力问题：在结构上满足a几何条件：研究对象是等厚度薄板。b载荷条件：作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布，而在两板面无外力作用。 平面应变问题：满足a几何条件：长柱体，即长度方向的尺寸远远大于横截面的尺寸，且横截面沿长度方向不变。b载荷条件：作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力两条件的弹性力学问题。 1、何为结构的离散化？离散化的目的？何为有限元模型？ ①离散化:把连续的结构看成由有限个单元组成的集合体。②目的：建立有限元计算模型③通常把由节点,单元及相应的节点载荷和节点约束构成的模型称为有限元模型2、结构离散化时，划分单元数目的多少以及疏密分布，将直接影响到什么？确定单元数量的原则？通常如何设置节点？
北京科技大学有限元试题及答案
一 判断题（20分） （×）1. 节点的位置依赖于形态，而并不依赖于载荷的位置 （√）2. 对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元 （×）3. 不能把梁单元、壳单元和实体单元混合在一起作成模型 （√）4. 四边形的平面单元尽可能作成接近正方形形状的单元 （×）5. 平面应变单元也好，平面应力单元也好，如果以单位厚来作模型化 处理的话会得到一样的答案 （×）6. 用有限元法不可以对运动的物体的结构进行静力分析 （√）7. 一般应力变化大的地方单元尺寸要划的小才好 （×）8. 所谓全约束只要将位移自由度约束住，而不必约束转动自由度 （√）9. 同一载荷作用下的结构，所给材料的弹性模量越大则变形值越小 （√）10一维变带宽存储通常比二维等带宽存储更节省存储量。 二、填空（20分） 1．平面应力问题与薄板弯曲问题的弹性体几何形状都是 薄板 ，但前者受力特点是： 平行于板面且沿厚度均布载荷作用 ，变形发生在板面内； 后者受力特点是： 垂直于板面 的力的作用，板将变成有弯有扭的曲面。 2．平面应力问题与平面应变问题都具有三个独立的应力分量： σx ，σy ，τxy ，三个独立的应变分量：εx ，εy ，γxy ，但对应的弹性体几何形状前者为 薄板 ，后者为 长柱体 。3．位移模式需反映 刚体位移 ，反映 常变形 ，满足 单元边界上位移连续 。 4．单元刚度矩阵的特点有：对称性 ， 奇异性 ，还可按节点分块。 5．轴对称问题单元形状为：三角形或四边形截面的空间环形单元 ，由于轴对称的特性，任意一点变形只发生在子午面上，因此可以作为 二 维问题处理。 6．等参数单元指的是：描述位移和描述坐标采用相同的形函数形式。等参数单元优点是：可以采用高阶次位移模式，能够模拟复杂几何边界，方便单元刚度矩阵和等效节点载荷的积分运算。 7．有限单元法首先求出的解是 节点位移 ，单元应力可由它求得，其计算公式为 {}{} [][]e D B σδ=。（用符号表示即可） 8．一个空间块体单元的节点有 3 个节点位移： u ，v ，w 9．变形体基本变量有位移应变应力 基本方程 平衡方程 物理方程 几何方程 10.实现有限元分析标准化和规范化的载体就是单元
西工大-有限元试题(附答案)
1．针对下图所示的3个三角形元，写出用完整多项式描述的位移模式表达式。 2．如下图所示，求下列情况的带宽: a)4结点四边形元； b)2结点线性杆元。 3．对上题图诸结点制定一种结点编号的方法，使所得带宽更小。图左下角的四边形在两种不同编号方式下，单元的带宽分别是多大 4．下图所示，若单元是2结点线性杆单元，勾画出组装总刚后总刚空间轮廓线。系统的带宽是多大按一右一左重新编号（即6变成3等）后，重复以上运算。
5．设杆件1－2受轴向力作用，截面积为A，长度为L，弹性模量为E，试写出 杆端力F 1，F 2 与杆端位移 2 1 ,u u之间的关系式，并求出杆件的单元刚度矩阵)(] [e k 6．设阶梯形杆件由两个等截面杆件○1与○2所组成，试写出三个结点1、2、3的 结点轴向力F 1，F 2 ，F 3 与结点轴向位移 3 2 1 , ,u u u之间的整体刚度矩阵[K]。 7．在上题的阶梯形杆件中，设结点3为固定端，结点1作用轴向载荷F 1 =P，求各结点的轴向位移和各杆的轴力。
8． 下图所示为平面桁架中的任一单元，y x ,为局部坐标系，x ，y 为总体坐标系，x 轴与x 轴的夹角为θ。 （1） 求在局部坐标系中的单元刚度矩阵 )(][e k （2） 求单元的坐标转换矩阵 [T]； （3） 求在总体坐标系中的单元刚度矩阵 )(][e k 9．如图所示一个直角三角形桁架，已知27/103cm N E ?=，两个直角边长度 cm l 100=，各杆截面面积210cm A =，求整体刚度矩阵[K]。
10．设上题中的桁架的支承情况和载荷情况如下图所示，按有限元素法求出各结点的位移与各杆的力。 11．进行结点编号时，如果把所有固定端处的结点编在最后，那么在引入边界条件时是否会更简便些 12．针对下图所示的3结点三角形单元，同一网格的两种不同的编号方式，单元的带宽分别是多大
西工大-有限元试题(附答案)汇总
1．针对下图所示的 3 个三角形元，写出用完整多项式描述的位移模式表达式 2．如下图所示，求下列情况的带宽 a) 4 结点四边形元； b) 2 结点线性杆元。 3．对上题图诸结点制定一种结点编号的方法，使所得带宽更小。图左下角的四 边形在两种不同编号方式下，单元的带宽分别是多大？ 4．下图所示，若单元是 2结点线性杆单元， 勾画出组装总刚后总刚空间轮廓线。 系统的带宽是多大？按一右一左重新编号(即 6变成 3等)后，重复以上运算
5．设杆件1－2 受轴向力作用，截面积为A，长度为L，弹性模量为E，试写出杆端力F1，F2与杆端位移u1, u2之间的关系式，并求出杆件的单元刚度矩阵[k](e) 6．设阶梯形杆件由两个等截面杆件○ 1 与○2 所组成，试写出三个结点1、2、3 的结点轴向力F1，F2，F3与结点轴向位移u1, u2, u3之间的整体刚度矩阵[K] 。 7．在上题的阶梯形杆件中，设结点3 为固定端，结点1作用轴向载荷 F1=P，求各结点的轴向位移和各杆的轴力。 8．下图所示为平面桁架中的任一单元， x, y 为局部坐标系，x，y 为总
体坐标系， x 轴与x 轴的夹角为。 1) 求在局部坐标系中的单元刚度矩 阵 [k]( e) 2) 求单元的坐标转换矩阵[T] ； 3) 求在总体坐标系中的单元刚度矩 阵[k] (e)
9．如图所示一个直角三角形桁架，已知E 3 107N / cm2，两个直角边长度 2 l 100cm ，各杆截面面积 A 10cm2，求整体刚度矩阵[K] 。 10．设上题中的桁架的支承情况和载荷情况如下图所示，按有限元素法求出各结点的位移与各杆的内力。
浅析电力系统模型参数辨识
浅析电力系统模型参数辨识 （贵哥提供） 一、现状分析 随着我国电力事业的迅猛发展, 超高压输电线路和大容量机组的相继投入, 对电力系统稳定计算、以及其安全性、经济性和电能质量提出了更高的要求。现代控制理论、计算机技术、现代应用数学等新理论、新方法在电力系统的应用，正在促使电力工业这一传统产业迅速走向高科技化。 我国大区域电网的互联使网络结构更复杂，对电力系统安全稳定分析提出了更高的要求，在线、实时、精确的辨识电力系统模型参数变得更加紧迫。由于电力系统模型的基础性、重要性，国外早在上世纪三十年代就开始了这方面的分析研究，[1,2]国内外的电力工作者在模型参数辨识方面做了大量的研究工作。[3]随后IEEE相继公布了有关四大参数的数学模型。1990年全国电网会议上的调查确定了模型参数的地位，促进了模型参数辨识的进一步发展，并提出了研究发电机、励磁、调速系统、负荷等元件的动态特性和理论模型，以及元件在极端运行环境下的动态特性和参数辨识的要求。但传统的测量手段，限制了在线实时辨识方法的实现。 同步相量测量技术的出现和WAMS系统的研究与应用，使实现在线实时的电力系统模型参数辨识成为可能。同步相量是以标准时间信号GPS作为同步的基准，通过对采样数据计算而得的相量。相量测量装置是进行同步相量测量和输出以及动态记录的装置。PMU的核心特征包括基于标准时钟信号的同步相量测量、失去标准时钟信号的授时能力、PMU与主站之间能够实时通信并遵循有关通信协议。 自1988年Virginia Tech研制出首个PMU装置以来，[4]PMU技术取得了长足发展，并在国内外得到了广泛应用。截至2006年底，在我国范围内，已有300多台P MU装置投入运行，并且可预计，在不久的将来PMU装置会遍布电力系统的各个主要电厂和变电站。这为基于PMU的各种应用提供了良好的条件。 二、系统辨识的概念 系统模型是实际系统本质的简化描述。[5]模型可分为物理模型和数学模型两大类。物理模型是根据相似原理构成的一种物理模拟，通过模型试验来研究系统的
有限元试题及答案
有限元试题及答案
一判断题（20分） （×）1. 节点的位置依赖于形态，而并不依赖于载荷的位置 （√）2. 对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元 （×）3. 不能把梁单元、壳单元和实体单元混合在一起作成模型 （√）4. 四边形的平面单元尽可能作成接近正方形形状的单元 （×）5. 平面应变单元也好，平面应力单元也好，如果以单位厚来作模型化 处理的话会得到一样的答案 （×）6. 用有限元法不可以对运动的物体的结构进行静力分析 （√）7. 一般应力变化大的地方单元尺寸要划的小才好 （×）8. 所谓全约束只要将位移自由度约束住，而不必约束转动自由度 （√）9. 同一载荷作用下的结构，所给材料的弹性模量越大则变形值越小（√）10一维变带宽存储通常比二维等带宽存储更节省存储量。 二、填空（20分） 1．平面应力问题与薄板弯曲问题的弹性体几何形状都是薄板，但前者受力特点是：平行于板面且沿厚度均布载荷作用，变形发生在板面内； 后者受力特点是：垂直于板面的力的作用，板将变成有弯有扭的曲面。 2．平面应力问题与平面应变问题都具有三个独立的应力分量：σx，σy，τxy ，三个独立的应变分量：εx，εy，γxy，但对应的弹性体几何形状前者为薄板，后者为长柱体。3．位移模式需反映刚体位移，反映常变形，满足单元边界上位移连续。 4．单元刚度矩阵的特点有：对称性，奇异性，还可按节点分块。
5．轴对称问题单元形状为：三角形或四边形截面的空间环形单元，由于轴对称的特性，任意一点变形只发生在子午面上，因此可以作为二维问题处理。6．等参数单元指的是：描述位移和描述坐标采用相同的形函数形式。等参数单元优点是：可以采用高阶次位移模式，能够模拟复杂几何边界，方便单元刚度矩阵和等效节点载荷的积分运算。 7．有限单元法首先求出的解是节点位移，单元应力可由它求得，其计算公式为。（用符号表示即可） 8．一个空间块体单元的节点有 3 个节点位移： u，v，w 9．变形体基本变量有位移应变应力基本方程平衡方程物理方程几何方程 10.实现有限元分析标准化和规范化的载体就是单元
GARTEUR 有限元模型修正与确认研究
收稿日期:2003207207;修订日期:2004203225 基金项目:教育部博士学科点专项基金(20010227012)资助项目 文章编号:100026893(2004)0420372204 GARTEUR 有限元模型修正与确认研究 费庆国,张令弥,郭勤涛 (南京航空航天大学振动工程研究所,江苏南京 210016) Case Study of FE Model Updating and Validation via an Air craft Model Structur e FEI Qing 2guo,Z HANG Ling 2mi,G UO Qin 2tao (Institu te o f Vi brati on Engi neering,Nanjing University of Aeronau tics and Astro nautics,Nanjing 210016,China)摘 要:待修正参数的选择以及修正后模型的质量评估是有限元模型修正的两个重要问题。以欧洲学术界广泛采用的GA RTEUR 飞机模型为例,利用基于灵敏度分析的模型修正方法,通过仿真算例研究参数选择对模型修正质量的影响,并以试验数据为目标值对有限元模型进行修正与确认。为全面评估模型的修正质量,引入三级标准对修正后有限元模型进行确认。 关键词:固体力学;模型确认;有限元法;模型修正;参数选择中图分类号:O 248121 文献标识码:A Abstr act:Parameter selection and quali ty validation are of g reat i mpo rtance in fini te element model updating.This paper presents so me results which demonstrate the relationship betw een parameter selection and updated model .s quality throu gh si mulation cases.Three q uali ty levels w ith corresponding validation criteria are emplo yed with an emphasis o n updated mod 2el .s predictio n ability.Results of updating based on exper i mental modal test data are sho w n as an application example.A n aircraft test structure,GA RTEUR,which is g enerally utilized in Europe,is employed in bo th the si mulation case and the exper i mental case.Sensi tivity 2based model updating appro ach is applied. Key wor ds:solid mechanics;model validation;finite element method;model updating;parameter selection 在航空工程中,准确的有限元模型对于动态响应预测以及动态设计至关重要。建模过程中的不确定因素,如离散化误差、材料物理参数的不确定性、边界条件的近似等,导致有限元模型必然存在误差。设计规范规定,有限元模型必须通过振动模态试验或者地面共振试验来检验[1]。 近30年来,有限元模型修正技术得到了长足的发展[2~6]。根据修正对象的不同可将修正方法分为矩阵型方法和设计参数型方法。后者物理意义明确,更具工程应用价值。本文采用基于灵敏度分析的设计参数型修正方法。 基于灵敏度分析的设计参数型修正方法主要包括待修正设计参数选择,灵敏度分析,参数修正以及模型确认等环节。 待修正设计参数的选择是模型修正的起始环节。通常,候选参数是有限元模型存在不确定性因素的参数。近20年虽然发展了很多种参数选择或者误差定位的策略与算法,工程应用中仍然难以准确无遗漏地确定误差参数。因此,有必要讨论参数选择对模型修正质量的影响。 模型确认是模型修正的检验环节。在当前的研究及工程应用中,通常只要求修正后模型的计算结果能够复现修正过程中利用的试验数据。事实上,为全面评估模型的质量,模型的复现能力与预测能力应予以同等重视[7]。本研究引入了三级质量标准对修正后的有限元模型进行确认。 本文采用G ARTE UR 飞机模型为研究对象,通过仿真算例来研究参数选择对模型修正质量的影响,并给出了利用振动模态测试结果对G AR 2TEUR 飞机模型的有限元模型进行修正与确认的结果。 1 模型修正方法与模型确认准则 (1)模型修正方法 模型修正可归结为以下的优化问题[8] Min p +R(p )+2 2,R(p )=f E -f A (p )s.t V L [p [V U (1) 其中:p 代表设计参数;f E ,f A 是结构动态特性试验与分析结果;R 代表残差;V L ,V U 是设计参数的下、上限。 令设计参数的初始值为p 0,动态特性f 是设计参数的隐函数,其泰勒展开式为 第25卷 第4期航 空 学 报 Vol 125N o 14 2004年 7月ACT A AERO NA U TICA E T AS TRO NA U TICA SINICA July 2004
板结构有限元分析实例详解
板结构有限元分析实例详解1：带孔平板结构静力分析本节介绍带孔平板结构静力分析问题，同时介绍布尔操作的基本用法。 8.3.1 问题描述与分析 有孔的矩形平板，左侧边缘固定，长400mm，宽200 mm，厚度为10 mm，圆孔在板的正中心，半径为40 mm，左侧全约束，右侧边缘均布应力1MPa，如图8.7所示。求板的变形、位移及应力变化情况。（材料的材料属性为：弹性模量为300000 MPa，剪切模量为0.31。） 图8.7 带孔的矩形平板 由于小孔处边缘不规则，本文采用PLANE82高阶平面单元进行分析。 8.3.2 求解过程 8.3.2.1 定义工作目录及文件名 启动ANSYS Mechanical APDL Product Launcher窗口，如图8.8所示。在License下 拉选框中选择ANSYS Multiphysics产品，在Working Directory输入栏中输入工作目 录：C:\ANSYS12.0 Structural Finite Elements Analysis and Practice\Chapter 8\8-1，在Job Name一栏中输入工作文件名：Chapter8-1。以上参数设置完毕后，单 击Run按钮运行ANSYS。
图8.8 ANSYS设置窗口菜单 可以先在目标文件位置建立工作目录，然后单击Browse按钮选择工作目录；也 可以通过单击Browse按钮选择工作文件名。 8.3.2.2 定义单元类型和材料属性 选择Main Menu>Preferences命令，出现Preferences for GUI Filtering对话框， 如图8.9所示，在Individual discipline(s) to show in the GUI中勾选Structural，过滤掉ANSYS GUI菜单中与结构分析无关的选项，单击OK按钮关闭该对话框。 图8.9 Preferences for GUI Filtering对话框
西工大有限元试题附答案68872
1.针对下图所示的3个三角形元,写出用完整多项式描述的位移模式表达式。 ２.如下图所示,求下列情况的带宽: a)4结点四边形元; b)2结点线性杆元。 ３、对上题图诸结点制定一种结点编号的方法,使所得带宽更小。图左下角的四边形在两种不同编号方式下,单元的带宽分别就是多大? 4、下图所示,若单元就是2结点线性杆单元,勾画出组装总刚后总刚空间轮廓线。系统的带宽就是多大?按一右一左重新编号(即6变成3等)后,重复以上运算。
５. 设杆件1-2受轴向力作用,截面积为Ａ,长度为Ｌ,弹性模量为Ｅ,试写出杆端力Ｆ１,F 2与杆端位移21,u u 之间的关系式,并求出杆件的单元刚度矩阵)(][e k 6、设阶梯形杆件由两个等截面杆件＼o ＼a ｃ(○,1)与错误!所组成,试写出三个结点1、2、3的结点轴向力F 1,F 2,Ｆ3与结点轴向位移321,,u u u 之间的整体刚度矩阵［K]。 7. 在上题的阶梯形杆件中,设结点3为固定端,结点1作用轴向载荷F 1=Ｐ,求各结点的轴向位移与各杆的轴力。 8、 下图所示为平面桁架中的任一单元,y x ,为局部坐标系,ｘ,y 为总体坐标系,x 轴与x 轴的夹角为 。 (1) 求在局部坐标系中的单元刚度矩阵 )(][e k (2) 求单元的坐标转换矩阵 ［T]; (３) 求在总体坐标系中的单元刚度矩阵 )(][e k
９.如图所示一个直角三角形桁架,已知27/103cm N E ?=,两个直角边长度cm l 100=,各杆截面面积210cm A =,求整体刚度矩阵［K ］ 。 10. 设上题中的桁架的支承情况与载荷情况如下图所示,按有限元素法求出各结点的位移与各杆的内力。
基于最小二乘模型的Bayes参数辨识方法
基于最小二乘模型的Bayes 参数辨识方法 王晓侃1，冯冬青2 1 郑州大学电气工程学院，郑州（450001） 2 郑州大学信息控制研究所，郑州（450001） E-mail ：wxkbbg@https://www.360docs.net/doc/b618372818.html, 摘 要：从辨识定义出发,首先介绍了Bayes 基本原理及其两种常用的方法,接着重点介绍了基于最小二乘模型的Bayes 参数辨识,最后以实例用MATLAB 进行仿真,得出理想的辨识结果。 关键词:辨识定义；Bayes 基本原理；Bayes 参数辨识 中国图书分类号：TP273+.1 文献标识码：A 0 概述 系统辨识是建模的一种方法。不同的学科领域，对应着不同的数学模型，从某种意义上讲，不同学科的发展过程就是建立它的数学模型的过程。建立数学模型有两种方法:即解析法和系统辨识。L. A. Zadehll 于1962年曾对”辨识”给出定义[1]:系统辨识是在对输入和输出观测的基础上，在指定的一类系统中，确定一个与被识别的系统等价的系统。一般系统输出y(n)通常用系统过去输出y(n-m)和现在输入u(n)及过去输入u(n-m)的函数描述 y(n)=f(y(n-1),y(n-2)，...，y(n-m y ), u(n),u(n-1),... ,u(n-m u ))=f(x(n),n) x(n)=[y(n-1),y(n-2),...y(n-m y ), u(n),u(n-1),...,u(n-m u )]’ 这里f(,)为未知函数关系，一般情况为泛函数，可以是线性函数或非线性函数，分别对应于线性或非线性系统，通常这个函数未知，但是局部输入输出数据可以测出，系统辨识的任务就是根据这部分信息寻找确定函数或确定系统来逼近这个未知函数。但实际上我们不可能找到一个与实际系统完全等价的模型。从实用的角度来看，系统辨识就是从一组模型中选择一个模型，按照某种准则，使之能最好地拟合由系统的输入输出观测数据体现出的实际系统的动态或静态特性。接下来本文就以最小二乘法为基础的Bayes 辨识方法为例进行分析介绍并加以仿真[4]。 1 Bayes 基本原理 Bayes 辨识方法的基本思想是把所要估计的参数看做随机变量,然后设法通过观测与该参数有关联的其他变量,以此来推断这个参数。 设μ是描述某一动态系统的模型，θ是模型μ的参数，它会反映在该动态系统的输入输出观测值中。如果系统的输出变量z(k)在参数θ及其历史纪录(1) k D ?条件下的概率密度函 数是已知的，记作p(z(k)｜θ，(1) k D ?)，其中(1) k D ?表示（k-1）时刻以前的输入输出数据集 合，那么根据Bayes 的观点参数θ的估计问题可以看成是把参数θ当作具有某种先验概率密 度p （θ，(1) k D ?）的随机变量，如果输入u(k)是确定的变量，则利用Bayes 公式，把参数θ 的后验概率密度函数表示成[2] p (θ，k D ）= p (θ｜z (k ),u(k ), (1) k D ?)=p (θ｜z (k ),(1) k D ?) = (k-1) (k-1) p(z(k)/,D )p(/D ) (k-1)(k-1)p(z(k)/,D )p(/D )d θθθθθ∞∫?∞ (1) 在式(1)中,参数θ的先验概率密度函数p(θ｜(1) k D ?)及数据的条件概率密度函数p(z(k)｜θ，
结构有限元及其应用软件
中国海洋大学本科生课程大纲 课程属性：公共基础/通识教育/学科基础/专业知识/工作技能，课程性质：必修、选修 一、课程介绍 1.课程描述（中英文）： 本课程是一门重要的结构计算分析课程，通过多媒体教学和上机练习，系统学习结构有限元FEM的基本原理和方法，熟悉掌握通用有限元应用软件ANSYS进行结构静力和动力分析的方法和步骤，并初步掌握使用ANSYS进行海工典型结构强度计算的方法。 Structural finite element method and its application software is an important course of structural calculation and analysis. Through multimedia teaching and computer practice, the basic principles and methods of Finite Element Method (FEM) are learned systematically. The general finite element application software ANSYS for the methods and procedures of structural static and dynamic analysis are mastered.At the same time, the strength calculation method of typical ocean engineering structures using ANSYS is preliminarily mastered. 2.设计思路： 有限元方法是一种现代设计方法，应用于结构设计中，是一种具有重要经济意义和巨大潜力的先进结构设计技术。因此选择该课程作为结构设计方面的一门必修课程，主要介绍结构有限元的基本原理和方法，还选择了通用的有限元软件ANSYS进行示例分析。包括要求掌握有限元法的基本思想和基本原理、平面刚架结构的有限元法、弹
最新有限元法基础试题
有限元法基础试题（A ） 一、填空题（5×2分） 1.1单元刚度矩阵e T k B DBd Ω = Ω? 中，矩阵B 为__________，矩阵D 为___________。 1.2边界条件通常有两类。通常发生在位置完全固定不能转动的情况为_______边界，具体指定有限的非零值位移的情况，如支撑的下沉，称为_______边界。 1.3内部微元体上外力总虚功： ()(),,,,e x x xy y bx xy x y y by d W F u F v dxdy δστδτσδ??=+++++??+(),,,,x x y y xy y x u v u u dxdy σδσδτδδ??+++??的表达式中，第一项为____________________的虚功，第二项为____________________的虚功。 1.4弹簧单元的位移函数1N +2N =_________。 1.5 ij k 数学表达式：令j d =_____，k d =_____，k j ≠，则力i ij F k =。 二、判断题（5×2分） 2.1位移函数的假设合理与否将直接影响到有限元分析的计算精度、效率和可靠性。（ ） 2.2变形体虚功原理适用于一切结构（一维杆系、二维板、三位块体）、适用于任何力学行为的材料（线性和非线性），是变形体力学的普遍原理。 （ ） 2.3变形体虚功原理要求力系平衡，要求虚位移协调，是在“平衡、协调”前提下功的恒等关系。 （ ） 2.4常应变三角单元中变形矩阵是x 或y 的函数。 （ ） 2.5 对称单元中变形矩阵是x 或y 的函数。 （ ） 三、简答题（26分） 3.1列举有限元法的优点。（8分） 3.2写出有限单元法的分析过程。（8分） 3.3列出3种普通的有限元单元类型。（6分） 3.4简要阐述变形体虚位移原理。（4分） 四、计算题（54分） 4.1对于下图所示的弹簧组合，单元①的弹簧常数为10000N/m ，单元②的弹簧常数为20000N/m ，单元③的弹簧常数为10000N/m ，确定各节点位移、反力以及单元②的单元力。（10分） 4.2对于如图所示的杆组装，弹性模量E 为10GPa ，杆单元长L 均为2m ，横截面面积A 均为2×10-4m 2，弹簧常数为2000kN/m ，所受荷载如图。采用直接刚度法确定节点位移、作用力和单元②的应力。（10分）
机械结构有限元分析
机械结构有限元分析 有限元分析软件ANSYS在机械设计中的应用 摘要：在机械设计中运用ANSYS软件进行有限元分析是今后机械设计发展的必然趋势，将有限元方法引入到机械设计课程教学中，让学生参与如何用有限元法来求解一些典型零件的应力，并将有限元结果与教材上的理论结果进行对照。这种新的教学方法可以大大提高学生的学习兴趣，增强学生对专业知识的理解和掌握，同时还可以培养学生的动手能力。在机械设计课程教学中具有很强的实用价值。 关键词：机械设计有限元 Ansys 前言：机械设计课程是一门专业基础课，其中很多教学内容都涉及到如何求取零件的应力问题，比如齿轮、v带、螺栓等零件。在传统的教学过程中，都是根据零件的具体受力情况按材料力学中相应的计算公式来求解。比如，在求解齿轮的接触应力时，是把齿轮啮合转化为两圆柱体的接触，再用公式求解。这些公式本身就比较复杂，还要引入各种修正参数，因此我们在学习这些内容时普遍反映公式难记，学习起来枯燥乏味，而且很吃力。 近年来有限元法在结构分析中应用越来越广泛，因此如果能将这种方法运用到机械设计课程中，求解一些典型零件的应力应变，并将分析结果和教材上的理论结果进行对比，那么无论是对于提高学生学习的热情和积极性，增强对重点、难点知识的理解程度，还是加强学生的计算机水平都是一件非常有益的事情。 由于直齿圆柱齿轮的接触强度计算是机械设计课程中的一个重要内容，齿轮强度的计算也是课程中工作量最繁琐的部分。下面就以渐开线直齿圆柱齿轮的齿根弯曲疲劳强度的计算为例，探讨在机械设计课程中用ANSYS软件进行计算机辅助教学的步骤和方法，简述如何将有限元方法应用到这门课程的教学中。 1．传统的直齿圆柱齿轮齿根弯曲疲劳强度的计算 传统方法把轮齿看作宽度为b的矩形截面的悬臂梁。因此齿根处为危险剖面，它可用30。切线法确定。如图l所示。 作与轮齿对称中心线成30。角并与齿根过渡曲线相切的切线，通过两切点作平行与齿轮轴线 的剖面，即齿根危险剖面。理论上载荷应由同时啮合的多对齿分担，但为简化计算，通常假设全部载荷作用于齿顶来进行分析，另用重合度系数E对齿根弯曲应力予以修正。 由材料力学弯曲应力计算方法求得齿根最大弯曲应力为：
Bouc-Wen 滞回模型的参数辨识
上海交通大学 硕士学位论文 Bouc-Wen滞回模型的参数辨识及其在电梯振动建模中的应用 姓名：周传勇 申请学位级别：硕士 专业：机械设计及理论 指导教师：李鸿光 20080201
Bouc-Wen滞回模型的参数辨识 及其在电梯振动建模中的应用 摘 要 电梯导靴是连接轿箱系统与导轨的装置，它能起到导向和隔振减振的作用。同时，在电梯的运行过程中它又将导轨由于制造或安装所造成的表面不平顺度传递给轿箱系统，从而引起轿箱系统的水平振动。国内外学者在电梯水平振动的建模和分析中，往往把导靴视为线性弹簧-阻尼元件来建模而忽略了非线性因素。事实上导靴与导轨之间存在非线性的迟滞摩擦力，本文通过实验的方法，采用Bouc-Wen 滞回模型来建立导靴-导轨非线性摩擦力模型。 Bouc-Wen滞回模型因其微分形式的非线性表达式而使得其参数辨识存在较大的困难，本文利用模型中部分参数的不敏感性，通过数学变换将非线性参数辨识问题转化为线性参数辨识问题，从而使得问题大大简化，参数辨识的效果也能满足要求。 基于以上导靴-导轨间摩擦力模型，本文进而建立了轿箱-导轨耦合水平振动动力学模型，该模型将轿箱系统等效为2自由度的平面运动刚体，将导靴等效为质量-弹簧-阻尼单元，同时考虑了导靴-导轨间的非线性摩擦力，以及导靴靴衬与导轨间接触的不连续性等。 在建立了轿箱-导轨耦合水平振动动力学模型后，利用Matlab/Simulink，建立了相应的仿真模型，开展了几种典型导轨不
平顺度激励（弯曲、失调和台阶）下的仿真分析。研究结果表明，这些分析对于电梯结构优化设计和动力学建模与分析有理论指导意义。 关键词：迟滞，参数辨识，非线性，动力学建模，系统仿真
结构有限元分析的形状处理方法_杜平安
结构有限元分析的形状处理方法 杜平安 摘要 介绍结构形状处理的各种方法,包括类型简化、细节简化、形式变换、局部结构和利用对称性等。 关键词 形状处理 有限元分析 建模 Abstract The processing method is intro-duced in the paper ,including ty pe simplifica tion 、details simplifica tio n 、fo rm tra nsfo rmatio n 、local structure a nd symm etry utiliza tion . Key words Shape processing Finite element analysis Modelling 收稿日期:1999-08-18 1 结构类型简化 根据结构形状、载荷和约束条件的特点,结构类型可分为空间问题、平面问题、轴对称问题、板壳问题和杆件问题等。其中平面问题和轴对称问题的几何模型是一平面图形,在平面上划分网格比在空间内划分要容易得多,单元数量也少得多。因此将空间问题作适当近似,使其按平面问题来处理,则可使分析过程大为简化。在图1a 中,计算轮毂与轴过盈配合的接触压力时,由于辐孔尺寸较小且远离接触面,因此可以不考虑辐孔而将轮毂简化为轴对称结构。同样,在计算图1b 中螺栓与螺母螺纹面上的接触压力时,由于螺旋升角较小,也可以不考虑升角的影响,而将螺栓与螺母简化为轴对称结构 。 图1 结构类型简化结构 2 结构细节简化 细节是结构中相对尺寸很小的局部,如倒圆、倒角、退刀槽和加工凸台等。根据网格划分特点,一条直线或曲线至少要划分一个单元边;一个平面或曲面至少要划分一个单元面;一个圆至少要用三个单元边离散,因此几何模型中的细节将限制细节处及其附近的网格大小,从而影响整个结构的网格分布和增加网格数量。图2是有无细节时自动划分出的网格,从中可以看出细节对网格划分的影响 。 图2 细节对网格划分的影响 因此,建立几何模型时应尽量忽略一些不必要的细节。在静力分析中,高应力区域中的细节会引起应力集中,细节大小和形状对应力影响很大,这些细节不能忽略。而处于结构低应力区的细节一般可以忽略。在动力计算中,由于结构固有频率和模态振型主要取决于结构的质量分布和刚度,因此细节一般可以忽略。在热分析中,细节不会在结构中引起局部高温,这时也可以考虑较少的细节。 3 结构形式变换 有些结构尽管形状不是很复杂,但划分网格却很困难。如果对结构形式作适当变换,则可使网格划分变得容易,划分出的单元更少。例如图3a 所示的带肋板,划分网格时需要用板单元和梁单元组合,且两类单元为偏心连接,自动分网难以满足这种要求。如果将带肋板变换为平板(图3b 所示),则在平板上划分网格要容易得多。 由于带肋板用于焊接而成支撑箱式立柱,其特 性要求主要是刚度,因此可按等刚度条件作为变换 · 26·《机械与电子》2000(1)
有限元考试试题及答案
一、 简答题（共40分，每题10分） 1. 论述单元划分应遵循的原则。 2. 说明形函数应满足的条件。 3. 说明四边形等参数单元中“等参数”的含义，即为什么要引入等参数单元。 4. 阐述边界元法的主要优缺点。 二、 计算题（共60分，每题20分） 1. 一杆件如图3所示，杆件上方固定后，在下方受垂直向下的集中力作用，已 知：杆件材料的杨氏模量2 721/100.3in lbf E E ?==，截面积2125.5in A =， 2275.3in A =，长度in L L 1221==，集中力lbf P 100=，用有限元方法求解B 点 和C 点位移。备注：（1）1 lbf （磅力，libra force ） = N 。（2）杨氏模量、弹性 模量、Young 氏弹性模量具有相同含义（10分） 2. 如图2 所示，有一正方形薄板，沿对角承受压力作用，厚度t=1m ，载荷 F=20KN/m ，设泊松比μ=0，材料的弹性模量为E ，试求它的应力分布。（15分） 学院 专业 学号 姓名 y 图1
图2 3. 图示结点三角形单元的124边作用有均布侧压力q，单元厚度为t，求单元的等效结点荷载。 图3
一、简答题 1. 答： 1）合理安排单元网格的疏密分布 2）为突出重要部位的单元二次划分 3）划分单元的个数 4）单元形状的合理性 5）不同材料界面处及荷载突变点、支承点的单元划分 6）曲线边界的处理，应尽可能减小几何误差 7）充分利用结构及载荷的对称性，以减少计算量 2. 答： 形函数应满足的三个条件： a.必须能反映单元的刚体位移，就是位移模式应反映与本单元形变无关的由 其它单元形变所引起的位移。 b.能反映单元的常量应变，所谓常量应变，就是与坐标位置无关，单元内所 有点都具有相同的应变。当单元尺寸取小时，则单元中各点的应变趋于相 等，也就是单元的形变趋于均匀，因而常量应变就成为应变的主要部分。 c.尽可能反映位移连续性；尽可能反映单元之间位移的连续性，即相邻单元 位移协调。 3. 答： 含义：所谓的等参数单元，就是在确定单元形状的插值函数和确定单元位移场的插值函数中采用了完全相同的形函数。 意义：构造出一些曲边地高精度单元，以便在给定地精度下，用数目较少地单元，解决工程实际地具体问题。 4. 答： 有限单元法是基于变分原理的里兹（Ritz）法的另一种形式，从而使里兹法分析的所有理论基础都适用子有限单元法，确认了有限单元法是处理连续介质问题的一种普遍方法．利用变分原理建立有限元方程和经典里兹法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的，面且事先不要求满足任何边界条件，因此它可以用来处理很复杂的连续介质问题。有限单元法中所利用的主要是伽辽金（Galerkin）法。它可以用于已经知道问题的微分方程和
结构有限元分析
中国海洋大学本科生课程大纲 一、课程介绍 1.课程描述： 本课程是船舶与海洋工程专业重要的结构计算分析课程，通过多媒体教学和上机练习，系统学习结构有限元FEM的基本原理和方法，熟悉掌握通用的有限元软件ANSYS进行结构静力和动力分析的方法和步骤。 2.设计思路： 有限元方法是一种现代设计方法，应用于结构设计中，是一种具有重要经济意义和巨大潜力的先进结构设计技术。因此选择该课程作为结构设计方面的一门必修课程，主要介绍结构有限元的基本原理和方法，还选择了通用的有限元软件ANSYS进行示例分析。包括要求掌握有限元法的基本思想和基本原理、平面刚架结构的有限元法、弹性力学平面问题以及结构动力学问题的有限元基本理论，并通过通用的有限元软件ANSYS了解解决相关问题的过程，同时掌握ANSYS进行结构静力和模态分析的基本步骤和方法，了解ANSYS进行结构瞬态动力分析的基本步骤。 3.课程与其他课程的关系 先修课程：结构力学、弹性力学。本课程与结构力学和弹性力学相关，在掌握了结构分析的基本概念和方法之后才能很好地学习结构有限元分析。在后续课程中结构 - 3 -
有限元分析为学生在海洋平台设计课程设计及毕业设计中提供了结构分析的方法和软件工具。 二、课程目标 本课程的目标是学习掌握现代结构分析方法FEM，初步掌握通用的有限元软件ANSYS，为船舶与海洋工程结构设计、强度校核提供计算结果，为海洋结构动力响应分析提供建模基础。 三、学习要求 结构有限元分析是一门理论和实践性都很强的课程，在机房上课，人手一台计算机，强调实际ANSYS操作能力的培养。要达到以上学习任务，学生必须： （1）按时上课，上课认真听讲，积极参与结构分析典型案例分析。本课程将包含较多的课堂有限元ANSYS作业练习和课后结构有限元计算作业。 （2）保质保量地按时完成课堂ANSYS作业练习和课后结构有限元计算作业，每位学生一个账号通过网络提交课堂ANSYS作业，只有在各项作业中认真练习才能够不断提高ANSYS的操作水平和结构分析的技能。 四、教学内容 - 3 -