Ansys桥梁计算.pdf

上传人:赵** 文档编号:60871843 上传时间:2022-11-18 格式:PDF 页数:19 大小:966.34KB
返回 下载 相关 举报
Ansys桥梁计算.pdf_第1页
第1页 / 共19页
Ansys桥梁计算.pdf_第2页
第2页 / 共19页
点击查看更多>>
资源描述

《Ansys桥梁计算.pdf》由会员分享,可在线阅读,更多相关《Ansys桥梁计算.pdf(19页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、-桥梁计算桥梁计算(常用的计算方法常用的计算方法)在ns单元库中,有近 200 种单元类型,在本章中将讨论一些在桥梁工程中常用到的单元,包括一些单元的输人参数,如单元名称、节点、自由度、实常数、材料特性、表面荷载、体荷载、专用特性、关键选项 KEYOP等。*关于单元选择问题这是一个大问题,方方面面很多,主要是掌握有限元的理论知识。首先当然是由问题类型选择不同单元,二维还是三维,梁,板壳,体,细梁,粗梁,薄壳,厚壳,膜等等,再定义你的材料:各向同性或各向异性,混凝土的各项参数,粘弹性等等。接下来是单元的划分与网格、精度与求解时间的要求等选择,要对各种单元的专有特性有个大概了解。使用s,还要了解

2、Ansys 的一个特点是笼统与通用,因此很多东西被掩盖到背后去了。比如单元类型,在Solid 里面看到十几种选择,Solid45,Slil,Sold5 等,看来区别只是节点数目上。但是实际上每种类型里还有eyopt 分成多种类型,比如最常用的线性单元oid5,其 Keot(1):includeorxludeetradispacement shae,就分为非协调元和协调元,Keopt:)(fullinegtion。ducedingratio其实又是两种不同的单元,这样不同组合一下这个 S5 实际上是包含了 6 种不同单元,各有各的不同特点和用处。因此使用 Any要注意各单元的 Keyop选项。不

3、同的选项会产生不同的结果。举例来说:对线性元例如 Sold5,要想把弯曲问题计算得比较精确,必须要采用非协调模式。采用完全积分会产生剪切锁死,减缩积分又会产生零能模式(ZM),非协调的线性元可以达到很高的精度,并且计算量比高阶刷、很多,在变形较大时,用 EhancdStrain 比非协调位移模式(EnhacedDilacmnt)更好(Soi85)。但是这些非协调元都要求网格比较规则才行,网格不规则的话,精度会大大下降,所以如何划分网格也是一门实践性很强的学问。采用高阶单元是提高精度的好办法,拿不定主意时采用高阶元是个比较保险的选择,但是高阶单元在某些情况下也会出现剪切锁死,并且很难发现,因此用

4、减缩积分的高阶元通常是最保险的选择,但是在大位移时,网格扭曲较大,减缩积分就不适用。不同结构形式的桥梁具有不同的力学行为,必须针对性地创建其模型,选择维数最低的单元去获得预期的效果(尽量做到能选择点而不选择线,能选择线而不选择平面,能选择平面而不选择壳,能选择壳而不选择三维实体)。下面的几节介绍一下桥梁工程计算中经常会用到的单元。*桥梁仿真单元类型一、建议选用的单元类型在桥梁用 Anys 建立模型时,可参照以下建议用的单元进行桥梁模型的建立。.1 梁(配筋)单元:桥墩、箱梁、纵横梁。-2板壳(配筋)单元:桥面系统。实体(配筋)单元:桥墩系统、基础结构。4拉杆单元:拱桥的系杆、吊杆。5.拉索单元

5、:斜拉桥的索、悬索桥的钢丝绳。.预紧单元:索力控制、螺栓铆钉连接。7.连接单元:支座、地基。二、常见桥梁连接部位在桥梁建立模型时要准确模拟边界条件,因此要准确分析连接部位的固有特性。(一)常见桥梁连接部位.固定支座、铰支、可滑移支座等空间支座系统。.带减振和隔振措施的减振支座系统。3地基主体之间桩-基系统。4刚构之间的螺栓连接、铆接等。5梁管之间的球接和铰接等。(二)连接部分解决方法Asy在解决桥梁不同的连接部位时可选用如下的方法:.Comin7、Cmbi4、Lnkll、ontact52、obne3弹簧(阻尼、间隙元):可用来模拟支座、绳索、拉杆等桥梁部件。.预紧单元可解决螺栓、铆钉连接的问题

6、。,.二力杆拉杆、索可解决拉索问题。.耦合与约束方程可解决梁与塔横梁的边界约束关系。5.接触单元如 Contact2 可模拟滑动支座、销接等部件的真实情况。(三)常见桥梁接触问题桥梁各个部分之间可能存在如下三种接触方式。1.滑移连接:点点接触。2.绑定连接:点面接触。3转动连接:面面接触。用接触单元可模拟如:滑移支座接触、挡块与其他部件的接触、振动时不同构件的碰撞等问题,这里不再一一赘述。三、桥梁基础的处理方式为了真实的模拟桥梁的实际情况,需要真实模拟桥梁的基础受力、变形及约束情况,建议建立模型时采用如下方案。1.基础平台与桩基:用实体模型、预应力配筋。2基础与岩土系统:有限区域实体模型、预应

7、力配筋。*桥梁常见模型处理一、桥梁中常用的模型可以用相应的单元.刚构桥、拱桥:梁与杆单元组合模型。2.钢管混凝土:复合截面梁模型。3.连续梁:梁模型。4.斜拉桥/悬索桥:梁、板壳、索或杆单元组合模型。5,立交桥:实体墩、板壳桥面和加强梁混合模型。6局部详细计算:实体(考虑配筋)或板模型,以便考虑模型细节特征,如结构尺寸构造倒角、厚薄或粗细过度、凹凸部分以及配筋等。二、桥梁建模要综合运用各种合适的单元对桥梁进行总体分析应该遵循如下原则:1.支座系统采用弹簧阻尼系统;2连接部位采用耦合与约束方程;3桥墩系统采用截面梁、配筋梁;4.桥面系统采用截面梁、配筋梁、板壳、梁板组合。对桥梁进行局部分析应该遵

8、循如下原则:支座系统采用实体模型:粘弹(粘塑、超弹、塑性)大变形(位移);-连接部位采用接触模型:实体、板壳、梁或组合模型;3.桥墩系统采用实体模型:配筋与混凝土破坏;.桥面系统采用实体或板壳:配筋与混凝土破坏,组合梁之间的耦合与约束方程。三、选用合适的分析方法、在对桥梁进行建模计算时对不同的计算目的要采取不同的计算步骤。(一)静态计算1根据分析类型承载特点建立合理梁、板、实体、拉杆(模拟索)模型;2.材料与几何非线性效应;3连接部位与支座的正确处理。(二)动力分析1尽量采用梁、板壳或二者组合模型;附属结构简化为质点,建立与总体结构耦合关系;3.连接部位与支座自由度协调合理;4应当考虑大变形、

9、初应力以及预张力的动力影响;.必须正确考虑阻尼效应;.材料与几何非线性效应。施加预应力的方式一、预应力的模拟方式Any里加预应力有几种方式:1.直接在单元中加,(Lkl0 等单元可以通过ea实常数来加)。2.用加力,然后在分析时打开 Prstlss,加。用温度变化模拟。在常用的软件系统中,预应力混凝土分析根据作用不妨分为两类:分离式和整体式。所谓分离式就是将混凝土和力筋的作用分别考虑(脱离体),以荷载的形式取代预应力钢筋的作用,典型.的如等效荷载法;而整体式则是将二者的作用一起考虑,典型的如 Ansys 中用 Ln单元模拟力筋的方法。(一)线性或非线性的考虑对于预应力混凝土结构,只要是开裂前阶

10、段的应力分析,完全可以将混凝土视为弹性材料,当然钢筋也是弹性材料,这主要在使用荷载阶段的应力 分析。假如要进行开裂和极限分析,则必须考虑二者的非线性特性。(二)分离式方法(等效荷载法)的特点主要优点是建模简单,不必考虑力筋的位置而可直接建模,当然网格划分也简单,对结构在预应力作用下的整体效应可比较快捷地掌握。其缺点是比较明显的:不便模拟细部,例如力筋所在位置对结构的影响显然是不同的,假如一定要模拟,则荷载必须施加在力筋的位置上,故其建模的方便性就消失了;等效荷载法没有考虑力筋对混凝土的作用分布和方向,力筋对混凝土作用显然在各处是不同的,而等效荷载法则没有计及此点;对张拉过程无法模拟;,在其他外

11、荷载作用下的共同作用不便考虑,否则要加入力筋(其建模则同整体式),不能确定力筋在外荷载作用下的应力增量;可以模拟应力损失的影响。但是对于只关注预应力混凝土结构的基本性能时,可以考虑采用等效荷载法。(三)整体式方法的特点将混凝土和力筋划分为不同的单元逛考虑,而模拟预应力可以采用降温方法和初应变方法。降温方法比较简单,同时可以模拟力筋的损失,单元和实常数几种即可;而采用初应变又要模拟力筋各处不同的应力时,每个 单元的实常数各不相等,工作量较大。所以比较而言,采用整体式时考虑降温方法为宜。主要缺点是建模不便,尤其是当力筋较多且曲线布筋时(可以采用DL 解决);其优点也比较明显:力筋的具体位置一定,对

12、结构的影响可全面的考虑;-力筋对混凝土的作用近似的得到考虑(在结点处);可以模拟张拉不同的力筋,以优化张拉顺序;不管何种荷载,都是力筋和混凝土共同承担的,可以得到力筋在任何荷载下的应力;可以模拟应力损失的影响。但在后张法中有几个问题是应该考虑的(当然可以不予理睬):力筋的滑动问题。在张拉过程中,力筋与混凝土之间没有粘结,存在接触和滑动,而张拉完毕后,一般又都建立了粘结。这个问题可以这样考虑,因为分析总是张拉完毕(哪怕是某一束),这时显然没有滑动问题了,即可以按有粘结处理;而在荷载作用下有了粘结,自然可以按有粘结处理。在张拉完毕后力筋的应力是已知的,在分析时输入降温也是按张拉应力反算的,计算后力

13、筋的应力显然不等于张拉应力。这里有弹性压缩的问题,即降温应该计人混凝土弹性压缩损失,你可以考虑增大一定的比例,然后降温计算,二者相符或差别合适时认可。综上,类似计算分析,建议采用整体式之降温模拟方法。二、建立预应力的模型在模拟预应力钢筋时,传统的方式是把预应力钢筋作为体积的边界,把混凝土体积分割开来,Glu后划分混凝土单元,边界就作为 Lik 或 Ba单元了。普通钢筋可以用 Solid65 的分布钢筋模拟,其实常数数据很简单。Sold单元模拟精度很高,只要各材料参数取的正确就可以十分好地 模拟混凝土构件加载到变形开裂破坏的全过程,与实验数据十分接近。但是要注意混凝土的取值,弹性模量抗压强度不能

14、取规范 值,要用实验公式。在模拟钢筋混凝土的时候,一般的建模做法是:先建立体,然后使用面去切割体,把体切割成几个部分,在各个部分之间有共用的边界线,把这些边界线定义成 Lik 单元,这样就在体单元内建立于 Link 单元。因为 Lnk 单元就在old 体单元的边界线上,这样在网格划分的时候两种单元会产生共同的节点,也就可以共同工作了。土和砼的模拟*土弹簧的模拟Likl相当于用杆件提供约束,只不过此类单元可以设为只能受拉或只能受压。如果混凝土与节点底板分离,单元将不起作用,否则 Hk0 单元要承受拉力。如果用 Linkl0 的话,你可以把不与结构相联接的一端所有自由度都约束,与结构连接二端是否约

15、束则看结构的具体情况,与Link0 无关。Lik0 单元只能提供轴向位移的约束,不可抗剪。刚度可以这样算,xAL,不与结构连接的节点可取沿约束方向上的任意位置。两节点确定后,A是定的。你只要保证 K 的值没错,愿意怎样组合都行。至于特定问题,可以定义两种 Linl单元,第一种只能受拉,第二种只能受压。当然 ElemnTyeNumer(单元类型指代号码)是不同的。如果你担心自己两种容易搞错,可以只定义一种 Lik0 单元,把第二种的节点取约束的反方向即可。*混凝土的模拟-Ansy的 Siid65 单元是专为混凝土、岩石等抗压能力远大于抗拉能力的非均匀材料开发的单元。它可以模拟混凝土中的加强钢筋,

16、以及材料的拉裂和压溃现象。钢筋混凝土有限元有三种基本模型:分离式、组合式和整体式。Sd5单元采用Q6 是整体式有限元模型,即将钢筋弥散于整个单元中,将加筋混凝土视为连续均匀材料,求出的是一个统一的刚度矩阵。Solid6单元中的钢筋采用实常数的方法进行添加,钢筋的尺寸由混凝土的体积比确定。可以得到三种钢筋,这三种钢筋可以具有不同的材料,不同的方向。.从抗剪的角度出发,箍筋在截面的位置可以是任意的,因此这种方法对于钢筋混凝土中均匀分布的箍筋的设置比较适合。但与纵筋的实际情况却有一定的距离,下面这两种方法可以更好地模拟纵筋的受力情况:)1(将纵筋密集的区域设置为不同的体,使用带筋的 Sol65 单元

17、,而无纵筋区则设置为无筋5 单元,这样就可以将钢筋区域缩小,接近真实的工程情况。()采用杆单元来模拟纵筋,即采用分离式的有限元模型。为了建模方便,可将实体分为几个部分,使其交线为纵筋位置,这样就可以对交线划分籽单元。此时,还可以对杆施加预应力来模拟预应力混凝土。工况组合荷载工况组合是荷载工况之间的运算,典型情况为当前在数据库中的荷载工况和在另外一结果文件中的荷载工况间的运算,运算结果将改写数据库中的结果数据部分,可以显示及列出荷载工况组合。一、典型的荷载工况步骤组合包括以下几步:1.用 kcde命令定义荷载工况,kcdef,cno,Lep,Sbtep,Km.用_cs命令将荷载工况一读人数据库

18、Lese,cno3用Lcpr 命令执行所需的运算 Lopr,Oper,Lcasel,Ope,Lcas2下面举例说明一下工况的具体应用。假设结果文件包括针对几个荷载步的结果,若想比较荷载步 5 和荷载步 7,并将最大值存人内存,做法如下:l_cdf,1,5 将荷载工况 1 指向荷载步kcef,2,7 将荷载工况 2 指向荷载步se,1 将荷载工况读人内存I_cver,Max,与荷载工况 2 比较数据库并将最大值存人内存。Lwite,1写当前荷载工况到文件 Jobname.L12Lcae,3 将荷载工况 3 读人内存Icope,Add,1在 JobnameL2 文件中将数据库追加到荷载工况中二、存

19、储组合后的荷载工况缺省情况下,荷载工况组合的结果存在内存中,并覆盖数据库中的结果部分。要保存这些结果,作为以后浏览或以后的荷载工况组合,用下列方法之一:1.将数据写到荷载工况文件中。-2将数据追加到结果文件中用 kwri,Lcno,Fnme,Ext,Di命令把当前内存中的荷载工况写到荷载工况文件中。用 kcwrite 命令把当前内存中的荷载工况写到荷载工况文件中。文件名为bnamLnn。这里皿为分配的荷载工况号。在后续的荷载工况组合中皿指的是存人荷载工况文件的荷载工况。例如有这样的问题:工况是自重,工况 2 是集中荷载,工况 3 二 15 倍的工况 1+11 倍的工况 2,这个工况如何定义?可

20、通过荷载步定义。在 NO3STEP 定义倍的工况 1;OSTEP 定义 1.倍的工况 2;NOSTP 定义 NO+O4+,其中 11 和 15 的系数可在荷载步中的选项中定义。也可对每种荷载分别定义为一种loadcae,然后在通用后处理器potl中进行组合。风荷载的确定首先要确定场地的风特性、平均风速、谱特性等,将紊流风速分解为定常平均风分量和相应的紊动分量 U+u()。平均风速可由地表条件根据对数律得出。“(c)就是要模拟的紊流风速序列。有了其谱特性,可以生成大气边界层紊流的人工风速序列,这是做结构风振计算的重要步骤。目前的随机序列拟合法有基于 FORIER 分析的波叠加法和时间序列理论的M

21、A 模型法,都是生成具有目标风速谱的高斯平稳序列。有实测得来的风速序列最好,但通常是用模拟的风速序列作为输入。模拟出的人工紊流序列要尽可能地符合实际大气边界层紊流特性,才能保证计算结果的合理性。目前为止都是以紊流目标谱或相关作为拟合目标,对多点则是相关矩阵或互谱矩阵。作随机响应分析 PS是很好的办法,基于大量实测资料而统计出的谱本身就是作随机响应分析的很好输入,比单纯的一条序列更有代表性。但由于随机计算难以考虑非线性,故生成具体的随机序列作为输入,用有限元法计算动力响应是更普遍的方法。常用于结构设计的谱公式有:avnt 谱,Kmal 谱,Kara谱等。谱拟合:不仅在风工程中有应用,在人工地震波

22、和人工海浪波生成中也有应用。相 应方法与风序列的生成大同小异。地震波的输入对于地震波的输入,可以把荷载记录做成文件,利用 Apl 的读取功能读人数据库中。下面的例子是自己编的一个小文件,修改一下可以更简洁。地震波时程记录分成了个文件,每个文件是一列,分别记录、y、z方向的加速度。这样就可以把加速度记录读取到 Asy数据库中作为数组。也可以把加速度记录做成一个文件,这样程序就简单多了。下面是计算部分语句:olu!进入求解模块Antye,ns!求解类型为瞬态mStrt=001 1 开始时间TmEnd=15001结束时间TmIncr=1 1 时间步长D,TM,trt,M nd,M ncr!循环ie,

23、0.05!指定时间-alpha,!指定质量矩际系数iBad,!指定刚度矩际系数Ael,Acx(),Ay(),Aez()!指定力 p 速度大小Sole!求解Edd!结束循环exit!退出求解模块初应力荷载在作桥梁设计时,为了验证结构的可靠度,往往需要对结构施加初应力荷载。在进行结构分析时,Ansys 中可以使用输入文件来把初应力指定为一种荷载。初应力荷载只许用于静态和完全瞬态分析中(分析可以是线性或非线性的)。初应力只能在分析的第一个荷载步中施加,用stle 命令来指定、列表和删除初应力。该命令只能用于/Soltion 处理器中。fi命令的oc 变元用于指定这些初应力的位置。初应力可指定在单元的

24、中心或单元积分点处,kc 的缺省值对应于单元中心,kc:对应于单元积分点。对于网格中的每个个别单元,也可以通过c=2 来指定不同的初应力位置。在这种情况下,每个单元的初应力位置将用个别单元的局部位置标志记录在初应力文件中。如二指定网格中的每一个单兀的初应力状态都是相同的。对于这种情况,对所有单元只需指定一个应力张量。只有单元类型lane、Pan42、Sol45、lane2、Slid2、Soid9、SheHl1、l2、Plne83、oll85、oid86、Slidl87,al8 和 Beaml89 支持初应力输入功能。Beam 单元和 Shdl 单元的初应力必须在所有域段的积分点处指定。要使用s

25、tie 命令,初应力必须列在一个外部 ASCII 文件中,初应力文件中的注释用“!”标记在注释行的第一个字符处指定。每个单元记录的第一行应该由字符串LS开头,后面跟单元号和任意的局部位置标志,这些项必须用逗号隔开。如果si命令的变元的值为、1 或,则局部位置标志将被忽略。如果oc 二 2,则必须对每一个单元指定局部位置标志。局部位置标志必须是下面的一个值:对于单元中心(缺省)为 0,对于积分点为。任何的其他值都会产生错误并使得sfile 命令被忽略。每个单元记录的第一行后面紧跟的一行指定单元的每一个应力点的单元应力记录。当oc:0 时,只需指定每个单元中心处的一个应力记录,当如1 时,每个单元

26、的应力记录的数目等于单元积分点的数目。nsys 要求每个应力记录中有 6 个应力张量分量,当 Loc:3 时,初应力文件中的第一个单元的应力记录将被用于指定所有单元的相同的初应力。如果对一个单元定义了单元坐标系(Eys),则初应力必须在这个坐标系中指定。如何实现铰接Ans可采用两种方法来实现铰接:1在同一位置用 2 个 Node,然后 C。2.只用一个 Ne,然后根据需要用一些可以释放某些自由度的单元,如 Bea4(eanl4 可以定义成 PINI,IN-FIX.)t,22,beam4,,,,,1keyopt,X,11!eamelemeti-fi(端铰结 J 端固结)-et,23,ea44,1

27、keop,3,8,1!beamelmentfix-p(端铰结,I 端固结)keopt,24,7,11!beamelemetn-(1 端铰结,J 端固结)下面举例说明铰接的应用。如上图所示为一中间铰接的结构,两端固结,均布荷载,作出弯矩图。命令流如下:/pr7!进入前处理e,1,am44!定义混凝土单元类型e,ea44!定义钢筋单元类型e,3,bam44!定义混凝土弹性模量keyt,,8,11!第一种单元的 J 节点 X、Y 方向转动放松keyopt,2,11!第一种单元的 I 节点 X、Y 方向转动放松mp,d,1,260!单元一的质量密度mp,ex,1,37!单元一的弹性模量r,3,4,5,

28、1,,0!单元一的实常数,1,!定义关键点k,2,10k,,10,8k,4,5,512,1,!通过关键点连接成直线1,2,3,liv,1,10,!将直线等分为 1份ldiv,2,8,!将直线 2 等分为 10 份lsel,11!选择直线 11LATF,,1,1,,4!赋予直线相关的特性sel,s,2!选择直线 2LATF,,1,2,,!赋予直线相关的特性s,s,,,0,1!选择直线 3 到 10lsel,a,!添加直线 1lsd,a,,12,18,!添加直线 1到 1T,,1,3,,4!赋予直线相关的特性allsel!全选les,dl!划分所有直线,1,!约束节点d,104,dlese,,,1

29、0,30,1!选择单元 10 到0sfeam,all,1,prs,e2!施加均布荷载lseesel,s,,35,5,1!选择单元5 到 50sfbeam,ll,1,re,2e!施加均布荷载lsellu!进入求解模块 slve!求解POSl!进入后处理tble,m,smi,5!建立弯矩单元表pls,,!绘制弯矩图,如下图 4 所示AUTOCAD 模型输入可以充分利用 AtoA强大的绘图功能,在 AtoCA中建立模型后,再输入nys 中进行计算。AutoCD 建立的模型可以通过以下两种方法传人 Ansys。1.对于三维实体(3d oject)AuD:Fie-Export 保存类型选IS(x.st)

30、-输入文件名-选实体(选d obct)ys:FlIpt-at.输入即可优点:用 Sa文件转换方便,而且一般不会有转换问题.缺点:只能转换三维实体或面域2.用 Is 格式文件交换AutoCADl2 自带输出s 格式文件工具,其他可通过 Alr 软件将 Df格式的模型转换为 Iges 的格式文件然后再转入 Anss:FleIort-iges输入即可优点:各种实体类型都能转换查询函数的使用-在sys 操作过程或条件语句中,常常需要知道有关模型的许多参数值,如选择集中的单元数、节点数、最大节点号等。此时,一般可通过Get命令来获得这些参数。现在,对于此类问题,我们有了一个更为方便的选择,那就是查询函数

31、一 Inquiry Fncto。Inqiuncton 类似于 Any的*Get 命令,它访问 Asys 数据库并返回要查询的数值,方便后续使用。nsys 每执行一次查询函数,便查询一次数据库,并用查询值替代该查询函数。假如你想获得当前所选择的单元数,并把它作为Do 循环的上界。传统的方法是使用并 Get 命令来获得所选择的单元数并把它赋给一个变量,则此变量可以作为*Do 循环的上界来确定循环的次数。Gt,Elmx,elm,,outD,I,1,Elmx*Enddo现在可以使用查询函数来完成这件事,把查询函数直接放在 xDO 循环内,它就可以提供所选择的单元数。*,I,lmior(0,13)*En

32、ddo这里的 Elmir 并不是一个数组,而是一个查询函数,它返回的是现在所选择的单元数。括弧内的数是用来确定查询函数的返回值。第一个数是用来标识你所想查询的特定实体(如单元、节点、线、面号等),括弧内的第二个数是用来确定查询函数返回值的类型(如选择状态、实体数量等)。同本例一样,通常查询函数有两个变量,但也有一些查询函数只有一个变量,而有的却有三个变量。查询函数的种类和数量很多,下面是一些常用、方便而快速快捷的查询函数。1.Area-Arnr(Aad,y)Arid-查询的面,对于 Key:2,1,4 可取为 0;Key,标识关于reidr 的返回信息,=选择状态 12=,定义的数目=13,选

33、择的数目=14,定义的最大数,1材料号2=,单元类型3=,实常数=4,节点数,=单元数Arinr(Araid,Key)的返回值对于 Key:10=,eaid 未定义1=,reaid 未被选择=1,reaid 被选择2.Keypoints-inr(id,Key)Kpid-查询的关键点,对于 K:12,13,14 为 0Ke标识关于 Kd 的返回信息=1,选择状态=12,定义的数目=13,选择的数目=14,定义的最大数目1=,数料号,=单元类型=,实常数4=,节点数,如果已分网-7,单元数,如果已分网Kpinq(pid,Key)的返回值对于 Ke=1=,未选择=0,未定义1=,选择3.Lne-Ls

34、inq(Lsd,e)Lsid-查询的线段,对于ey=12,13,4 为 0Ky,标识关于 Lsid 的返回信息=,选择状态=,长度=1,定义的数目=3,选择的数目=14,定义的最大数1=,材料号2=,单元类型,3实常数4=,节点数=6,单元数Nde-Ndiqr(Nod,Key)ode-节点号,对于 Ke:2,1,14 为ey标识关于oe 的返回信息1=,选择状态=12,定义的数目,3=选择的数目1=,定义的最大数,超单元标记3=,主自由度4,激活的自由度,5附着的实体模型Nnqr(No,Ke)的返回值对于 y:11=,未选择,=未定义,1选择5。omsVln(Vni,Key)nm一查询的体,对

35、于 Key12,3,1为 0y标识关于 Vi 的返回信息1=,选择状态,21定义的数目13=,选择的数目14=,定义的最大数目1=,数料号2=,单元类型=3,实常数=,节点数=6,单元数=8,单元形状9=,中节点单元,=单元坐标系lir(Vni,Ke)的返回值.对于 Key:11=,未选择=0,未定义=1,选择用 ANSYS 作桥梁计算十三(其他文件网格划分)原创导人的图形联系在一起,但是没有粘合在一起,如果这样进行网格划分,各个图形就不会通过单元节点联系在一起,因此载荷也没有办法进行传递,必须运用布尔操作进行处理,其处理方法有以下几种。.直接进行 Gue 布尔操作,如果图形形状简单,布尔操作

36、就可能成功,这样可以直接进行下一步的处理。2如果进行布尔操作失败,可根据提示进行修正。有两种提示错误:一是说这种情况不适合用 Glue 操作,建议用 Overl命令,那就用 Overlap 命令进行处理,Overla完成后,再进行 Glu处理。如果提示拓扑退化,参见二说明。-二是说拓扑退化,建议放宽布尔操作的精度,这时应该放宽精度,放宽精度的时候应注意,放宽的精度要至少比图形中的最小尺寸精度小一个数 量级,否则可能导致图形划分出现大的误差,甚至导致错误。放宽精度后,进行ue 操作,一般都会成功(注意放宽精度进行布尔操作完成后,必须把精度调回来)。3如果仍然不成功,可以采取 Add 操作,但是不

37、能把所有的图形都进行 Ad 操作,如果模型的材料不同,将来就无法给各个模型分配不同的材料属性。这时要采取的步骤是把材料相同的有连接面或线的模型也就是体积 Add 在一起,使他们成为一个独立的实体,把所有的模型都进行了各自相关的dd 操作后,对新生成的几个体积进行 Gl操作,按照 1、2 的提示进行,就会成功。一种比较没有依据的方法,就是把体积分开,一个一个的进行布尔操作,不行的话,记着换一换顺序,有时就可以顺利的通过了。载荷步选项载荷步选项是各选项的总称,这些选项包括控制载荷如何在求解过程 中使用的选项以及其他诸如输出控制、阻尼特性设置和响应频谱数据等选项。载荷步选项随载荷步的不同而不同,有以

38、下几种类型的载荷步选项:1通用选项2动力学选项输出控制 4.ot-Saart 选项5.谱选项其中,通用选项包括:瞬态或静态分析中载荷步结束的时间,子步数或时间步大小,载荷阶跃或递增,以及热应变计算的参考温度。以下是对每个选项的简要说明。(一)时间选项Tie 命令用于指定在瞬态或静态分析中载荷步结束的时间。在瞬态或其他与速率相关的分析中,ime 命令指定实际的、按先后顺序的时间,且要求指定一时间值。在与速率无关的分析中,时间作为一跟踪参数。时间不能设置为 0,如果瞬态分析中,可指定一个非常小的值,如 Time,le-6(二)子步数和时间步大小对于非线性或瞬态分析,要指定一个载荷步中需要的子步数。

39、指定子步数的方法如下:Dehi,Dtm,Dtm,Dmx,Car和 Nubst,Nsb,Nsbm,bmn,aryNsubst 命令指定子步数,Dhm 命令指定时间步的大小。在缺省情况下,nys 程序在每个载荷步中使用一个子步。(三)自动时间步长Autos,Key 命令激活时间步自动阶跃,在时间步自动阶跃时,根据结构或组件对施加的载荷的响应,程序计算每个子步结束时最优的时间步。(四)阶跃或递增载荷c,Ke在一个载荷步中指定多个子步时,需要指明载荷是逐渐递增还是阶跃形式。Kc 命令用于此目的:Kbc,0 指明载荷是逐渐递增;K,指明载荷是阶跃载荷。关于阶跃载荷和逐渐递增载荷的说明:.如果指定阶跃载荷

40、,程序按相同的方式处理所有载荷(约束、集中载 荷、表面载荷、体载荷和惯性载荷)。根据情况,可分步施加、分步改变或分步移去这些载荷。2.如果指定递增载荷,那么在第一个载荷步施加的所有载荷,除了薄膜系数外都是从 0 开始逐渐变化的。在随后的载荷步中,所有载荷的变化都是从先前的值开始逐渐变化,-但要注意在完全谐(typ,HarmWitHopt,ll)分析中,表面载荷和体载荷的逐渐变化与在第一个载荷步中的变化相同,且不是从先前的值开始逐渐变化。3对于表格式边界条件,载荷不会逐渐变化,而是在当前时间点计算。如果载荷是在一个载荷步内用表格格式给定的而在下一个载荷步内变为非 表格格式,则该载荷被作为新引入的

41、载荷从 0 开始渐变。在随后的载荷步中被删除的所有载荷,除了体载荷和惯性载荷外,都是分步移去的。5在相同的载荷步中,不应删除或重新指定载荷,在这种情况下,逐渐变化不会按用户期望的方式进行。用 ANSYS 作桥梁计算十五(动力学分析选项)原创一、用于动态和瞬态分析的命令以下命令主要是用于动态和其他瞬态分析的选项,包括:1.mint,ey,La激活或撤消时间积分效应2.arfq,Freq,Fpeqe 在谐波响应分析中指定载荷的频率范围3 Alpa,Vale 指定考虑阻尼时的质量矩阵系数4Betad,Valu指定考虑阻尼时的刚度矩阵系数5.mprat,ato 指定结构动态分析的阻尼系数二、非线性选项

42、指定每个子步最大平衡叠代的次数(缺省二 25)Nt,Nuber指定收敛公差 Cnto,Lab,Value,Tle,Nor,MOnref为终止分析提供选项 Ncnv,Kso,Dli,Itlim,Et,Clim三、输出控制输出控制用于控制分析输出的数量和特性,有两个基本输出控制:.Ootres,Item,Feq,Cne控制nsys 写人数据库和结果文件的内容以及写入的频率。Oup,Item,Fr,nae控制打印的内容以及写入的频率。下例说明 Outrs 和 Opr 命令的使用:Oe,All,5!写入所有数据,每到第 5 子步写人数据utp,Nso,as!仅打印最后子步的节点解求解多载荷步定义和求解

43、多载荷步有三种方法:1.多步求解法2.载荷步文件法3矩阵参数法所有载荷和载荷步选项一起构成一个载荷步,程序可用其进行计算求解,如果有多个载荷步,可将每个载荷步存人一个文件,调用该载荷步文件,并在后面的求解中读人。Lswrie命令写载荷步文件(每个载荷步一个文件,以 Jbnames01,Jbnams0,Jbnames3 等标识),使用Lswrie,anum 命令将荷载步数据存人文件。所有载荷步文件写入后,可以使用命令在文件中顺序读取数据,并求得每个载荷步的解。下列所示的命令组定义多个载荷步:So!载荷步 1D,Nsubst,.bc,tres,Otpr,srit!写载荷步文件:Jobame01!载

44、荷步 2D,F,-Nsubst,.Kbc,Out,Op,swri!写载荷步文件:bname.s02载荷步数据用 Ansys 命令写入文件,ite 命令不捕捉实常数或材料特性的变化。Lswrie 命令自动地将实体模型载荷转换到有限元模型,因此所有载荷按有限元载荷命令的形式被写入文件,特别地,表面载荷总是按Se 或 Sbeam 命令的形式被写入文件,而不论载荷是如何施加的。要修改载荷步文件序号为 n 的数据,执行命令,Lsread,Lnum 在文件中读取数据,作所需的修改,然后执行 Lwrite,n 命令mdele 命令允许你从 Ansys 程序中删除载荷步文件。使用多步求解法是最直接的,它包括在

45、每个载荷步定义好后执行 Solve命令。主要缺点是,在交互使用时必须等到每一步求解结束后才能定义下一载荷步,比较费时。典型的多步求解法命令流如下所示:/lu!载荷步 1D,.SF,.Solve!求解载荷步1 载荷步D,,,S,,ole!求解载荷步 2tm_en50m_inc*DO,t,t_start,_end,tmicTime,tm,1,fy,force(m)!在第 100 号节点上加载随时间变化的力Slve*Endd重新启动一个分析有时在初始运算完成后也许要重新启动分析过程,例如想将更多的载荷步加入到分析中,也许在线性分析中要加入别的载荷条件,或在瞬态分析中加入另外的时间历程加载曲线,或者在

46、非线性分析收敛失败时需要恢复。一、重启动条件要重启动一个分析,模型必须满足下列条件:1.分析类型必须是静态(稳态),谐波或瞬态(只是u方法)分析,其他分析不能重启动。2.在初始运算中,必须已完成至少一次叠代。3.初始运算应该不是由于KiHed作业、系统中断或崩溃而停止的。一般重启动要求作业初始运算产生的某些文件存在,并且要求在 Solve命令之前对输入做一些修改,初始运算产生的下列文件必须可用:Jobnae 曲,onam ea单元矩阵文件(如果已经建立),obname.s或 Jobname.osv 保存的单元数据(esav)或保存的旧的单元数据(.osa),只有当 Job-ameesav 文件

47、丢失、不完整或者由于求解发散、位移超限或主元为负引起Jobnmesa文件的其他方面的破坏时,才需要 Jbnm.osav 文件。如果先前的运算附带产生了.rdb,lh,man 文件,必须在进行一般启动之前删除它们。常用操作分析见表 4-1。二、一般重启动的步骤一般重启动的步骤如下:进入 Anss 程序,用命令给定与第一次运行时相同的文件-名。2.用Su 命令进人求解处理器,用/Resme 命令恢复数据库文件。3执行 Antpe,Rst 命令指出这是一个重启动分析。4按需要规定修正的载荷或增加的载荷。调整先前坡道载荷的起始 点值,新加的坡道载荷从零开始变化,新施加的体载荷从初始值开始。重严加的和删

48、掉的载荷可视为新施加的,而没有修改。待删除的表面载荷和体载荷,必须减小到零或初始值,以保持obnmesav 文件和 Jobnae.o文件的数据库一致。.指定是否要重新使用三角化矩阵(Jonametd),可使用命令 Kuse,通过执行 Kuse 命令可以使程序重新形成单元矩阵,这样对调试分析和处理错误情况是有利的。6发出 Slve 命令初始化重启动步骤。7对增加的载荷步(如果有的话)重复步骤 46,或使用载荷步文件法产生和求解多载荷步(ssolv)。8.按需要进行后处理,然后退出 Ansys。重新启动输入实例:,,,,/ResumSouAt,,Rest!指定新载荷、新载荷步选项等!对非线性分析,

49、采用恰当的操作SveSolveSe/xt三、边界条件重建要为重启动重建正确的边界条件,首先要运行“虚拟”载荷步,过程如下:.将bnae.o文件改为 Jbname.asav 文件。2.进入 Ass 程序,用/命令给定与第一次运行时相同的文件名。3.用Sol命令进入求解处理器,用Rsme 命令恢复数据库文件。4执行 Anype,est 命令指出这是一个重启动分析。5从上一次已成功求解过的子步开始重新规定边界条件。6执行lv命令。7.按需要施加最终载荷及载荷步选项。如载荷步为前面(在虚拟前)加载荷步的延续,需调整子步的数量(或时间步步长)。时间步长编号可能会发生变化,与初始意图不同。如果你需要保留时

50、间步长编号,可在步骤 6中使用一个小的时间增量。8继续前面所述的一般重启动的操作步骤。正确的进行后处理建立有限元模型并求解后,要分析结果的正确性。如果结果不正确可能是由于出现了奇异解。奇异解出现于可能产生不定解或非惟一解的分析中,求解方程的主元为负或零会产生这样的奇异解。有些情况下,尽管遇到主元为负或零,仍可继续进行分析,出现这种情况时,可用Pvhe命令指定是否要停止分析。下述条件会引起求解过程出现奇异性:约束不够。2.模型中有非线性单元(如间隙单元、滑动单元、铰链单元、线缆单元 h:-等),结构的一部分可能已经塌陷或分散了。3材料特性值为负,如密度。4连接点无约束,单元排列可能会引起奇异性。

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 教育专区 > 高考资料

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号© 2020-2023 www.taowenge.com 淘文阁