ansys提取刚度阵质量阵的快速方法
ansys分析混凝土的若干问题
1. 讨论两种Ansys求极限荷载的方法 （1）力加载 可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。 （2）位移加载 给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。 希望众高手讨论一下 （1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？ （2）位移法求极限荷载的具体步骤？ 2. 需要注意的问题 1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题； 2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中； 3. 六面体的SOLID 65 单元一般比四面体的单元计算要稳定且收敛性好，因此，只要条件允许，应该尽量使用六面体单元； 4. 正确选择收敛标准，一般位移控制加载最好用位移的无穷范数控制收敛，而用力控制加载时可以用残余力的二范数控制收敛。在裂缝刚刚出现和接近破坏的阶段，可以适当放松收敛标准，保证计算的连续性； 3. 关于下降段的问题 1）在实际混凝土中都有下降段，但是在计算的时候要特别小心下降段的问题。 2）下降段很容易导致计算不收敛，有时为了计算的收敛要避免设置下降段，采用rush模型。 3）利用最大压应变准则来判断混凝土是否破坏。 4. Solid65单元中的破坏准则 1）采用Willam&Warnke五参数破坏准则 2）需要参数： 单轴抗拉强度，单轴，双轴抗压强度，围压压力，在围压作用下双轴，单轴抗压强度
ansys质量矩阵刚度矩阵提取
ansys质量矩阵刚度矩阵提取 看了这么久了都没人回，查了一些质料终于找到答案了，，下面提供三种方法：方便与其他程序进行接口编程1. Which matrix you would like? element stiffness matrix or full stiffness matrix? element stiffness is within file.emat. full stiffness matrix is within file.full A simple way to dump the matrix is as follow: ------------------- /aux2 fileaux2,file,emat form,long dump,all ------------------- 2. 可以使用/DEBUG命令来得到。详细步骤参见下面的宏文件 finish /clear PI=3.1415926 w1=3 w2=10 w3=6 w4=1.2 r=.8 t=0.08 /PREP7 !* ET,1,SHELL63 R,1,t ET,2,MASS21 R,2,500,500,500,2000,2000,2000,
!* UIMP,1,EX, , ,2e11 UIMP,1,NUXY, , ,0.3, UIMP,1,DAMP, , ,0.2, UIMP,1,DENS, , ,7800, BLC4,0,0,w2,w1 ESIZE,1.5,0, AMESH,all NSEL,S,LOC,X,0.0 D,all, , , , , ,ALL, , , , , allsel,all SFA,all,1,PRES,12 FINISH /OUTPUT,cp,out,, ! 将输出信息送到cp.out文件 /debug,-1,,,1 ! 指定输出单元矩阵 /SOLU SOLVE finish /OUTPUT, TERM ! 将输出信息送到output windows中 ! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵 3. 其原理很简单，即使用ansys的超单元即可解决问题。定义超单元，然后列出超单元的刚度矩阵即可。 面是一个小例题，自可明白。 /prep7 k,1 k,2,3000 l,1,2 et,1,beam3 mp,ex,1,2e5 mp,prxy,1,0.3 mp,dens,1,2e3 r,1,5000,2e7,200 lesize,all,,,10 lmesh,all
钢筋混凝土梁的ansys分析
摘要 本文介绍ANSYS 模拟钢筋混凝土梁的过程，讨论了有限元模型的建立以及在 ANSYS 中的实现，给出了用分离式配筋方法对混凝土梁的分析的一般过程。并给出了详细的命令流过程。并在此基础上对混凝土梁进行了分析，讨论了在力的作用下混凝土梁的塑形变形和裂缝的发展过程。 关键词 Ansys 混凝土梁 分离式配筋 The analysis of mechanics of a reinforced concrete based on ANSYS Abstract This paper introduces ANSYS simulation of the reinforced concrete beam process, discusses the establishment of the finite element model and the realization, and gives the ANSYS reinforcement method with separate the analysis of concrete beams of the general process. And gives the detailed command flow process. Based on the analysis of concrete beams, and discussed the concrete beam under the action of forces of the body deformation and fracture process. Keywords Ansys concrete beams reinforced separated 1 引言 由于钢筋混凝上材料性质复杂，使其表现出明显的非线性行为[1]。长期以来采用线弹性理论的设计方法来研究钢筋混凝上结构的应力或内力，显然不太合理，尽管有此理论是基于人量试验数据上的经验公式，还是不能准确反映混凝上的力学性能，特别是受力复杂的重要结构，必须采用三维钢筋混凝上非线性有限元方法才能很好地掌握其力学性能。利用ANSYS 对钢筋混凝上结构弹塑性的仿真分析，可以对结构自开始受荷载直到破坏的全过程进行分析，获得不同阶段的受力性能。本文将以混凝土梁的弹塑性分析为例，介绍在Ansys 中分析材料非线性问题的具体实现方法。 2 问题介绍 如图所示的钢筋混凝土梁[2]，横截面尺寸为200400b h mm mm ?=?,梁的跨度为3.0L m =，支座宽度为250mm 采用C20混凝土，梁内受拉纵筋3φ20，架立筋采用2φ12， 箍筋采用φ6@150，钢筋保护层厚度为25mm 。如图一。 图一 对于梁中所采用的所有钢筋，弹性模量为5 2.110MPa ?，抗拉强度设计值210MPa ， 密度33 7.810/kg m ?，泊松比为0.3。
ANSYS中混凝土的本构关系
一、关于模型 钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、分布式和组合式模型。考虑钢筋和混凝土之间的粘结和滑移，则采用引入粘结单元的分离式模型；假定混凝土和钢筋粘结很好，不考虑二者之间的滑移，则三种模型都可以；分离式和分布式模型适用于二维和三维结构分析，后者对杆系结构分析比较适用。裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型，后者目前尚处研究之中，主要应用的是前两种。离散裂缝模型和分布裂缝模型各有特点，可根据不同的分析目的选择使用。随着计算速度和网格自动划分的快速实现，离散裂缝模型又有被推广使用的趋势。 就ANSYS而言，她可以考虑分离式模型(solid65+link8，认为混凝土和钢筋粘结很好，如要考虑粘结和滑移，则可引入弹簧单元进行模拟，比较困难！)，也可采用分布式模型(带筋的solid65)。而其裂缝的处理方式则为分布裂缝模型。 二、关于本构关系 混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类，其中研究最多的是非线性弹性和弹塑性本构关系，其中不乏实用者。混凝土破坏准则从单参数到五参数模型达数十个模型，或借用古典强度理论或基于试验结果等，各个破坏准则的表达方式和繁简程度各异，适用范围和计算精度差别也比较大，给使用带来了一定的困难。 就ANSYS而言，其问题比较复杂些。 1 ANSYS混凝土的破坏准则与屈服准则是如何定义的？ 采用tb,concr,matnum则定义了W-W破坏准则(failure criterion)，而非屈服准则(yield criterion)。W-W破坏准则是用于检查混凝土开裂和压碎用的，而混凝土的塑性可以另外考虑（当然是在开裂和压碎之前）。理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的，例如在高静水压力下会发生相当的塑性变形，表现为屈服，但没有破坏。而工程上又常将二者等同，其原因是工程结构不容许有很大的塑性变形，且混凝土等材料的屈服点不够明确，但破坏点非常明确。 定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系，即W—W破坏准则、混凝土开裂和压碎前均为线性的应力应变关系，而开裂和压碎后采用其给出的本构关系。但屈服准则尚可另外定义(随材料的应力应变关系，如tb,MKIN,则定义的屈服准则是Von Mises,流动法则、硬化法则也就确定了)。 2 定义tb,concr后可否定义其它的应力应变关系 当然是可以的，并且只有在定义tb,concr后，有些问题才好解决。例如可以定义tb,miso,输入混凝土的应力应变关系曲线(多折线实现)，这样也就将屈服准则、流动法则、硬化法则等确定了。 这里可能存在一点疑问，即ANSYS中的应力应变关系是拉压相等的，而混凝土材料显然不是这样的。是的，因为混凝土受拉段非常短，认为拉压相同影响很小，且由于定义的tb,concr 中确定了开裂强度，所以尽管定义的是一条大曲线，但应用于受拉部分的很小。 三、具体的系数及公式 1 定义tb,concr时候的两个系数如何确定？ 一般的参考书中，其值建议先取为0.3~0.5(江见鲸)，原话是“在没有更仔细的数据时，不妨先取0.3~0.5进行计算”，足见此0.3~0.5值的可用程度。根据我的经验和理由，建议此值取大些，即开裂的剪力传递系数取0.5,（定要>0.2）闭合的剪力传递系数取1.0。支持此说法的还有 现行铁路桥规的抗剪计算理论，以及原公路桥规的容许应力法的抗计剪计算。
提取刚度矩阵
============为什么折腾这个文档======== 我有一个计算线性动力学方程组的瞬态、谐响应和静力学的python程序，现希望开发一个将ANSYS组集好的总体矩阵导入该PYTHON程序中的接口。 该问题可分解为： [STEP1] [ANSYS]->[包含矩阵信息的文件] [STEP2] [包含矩阵信息的文件]->[python通用数据对象] [STEP3] [python通用数据对象]->[程序特定数据对象]->[进行计算] 因此检索了一些帖子，基本上完成了这项工作，本文是对[STEP1]和[STEP2]的整理，并且利用[STEP3]对结果进行了验证 ============主要内容================== 1，了解从ANSYS中提取总体矩阵和载荷向量的方法； 2，了解提取出来的矩阵是怎样表示的； 3，说明在Python中，如何读取这样的矩阵； 4，构造一个简单的算例，说明整个【建模】-【提取】-【读取】过程及其正确性； =========站内检索综述==================== 检索词：提取矩阵 得到21个结果，代表性的帖子有下面这9个： 编号[1] 标题：ansys中怎样提取质量，刚度，阻尼矩阵？ 地址：https://www.wendangku.net/doc/7f12896761.html,/forum-vi ... fromuid-159019.html 要点：pengweicai给出了一段网上最常见的提取代码,该程序以fortran 写成，可以利用.full文件以及一些列约定将ANSYS中的总体矩阵读入FORTRAN中。 编号[2] 标题：如何得知HBMAT命令提取的质量、刚度矩阵对应的自由度？ 地址：https://www.wendangku.net/doc/7f12896761.html,/forum-vi ... fromuid-159019.html 要点：提出了使用HBMAT命令提取稀疏矩阵时常见的问题：我们如何知道提取出来的信息是怎么储存的呢？ 编号[3] 标题：[分享]ANSYS中整体、单元刚度和质量矩阵的提取 地址：https://www.wendangku.net/doc/7f12896761.html,/forum-vi ... fromuid-159019.html 要点：在该帖子的7楼，其实已经给出了帖子[2]中问题的解答，即HBMAT 中提取出来的矩阵是Harwell-Boeing格式的，并且给出了该格式的细节，可惜是英文的，没引起多少关注。 编号[4] 标题：帮我看看提取的刚度与质量矩阵 地址：https://www.wendangku.net/doc/7f12896761.html,/forum-vi ... fromuid-159019.html 要点：这个帖子所示的矩阵并非是使用HBMAT命令提出出来的，而应该是SELIST命令列举出来的未压缩的矩阵，后续楼层的回帖给了大家一个提示，即有可能提取出来的矩阵是引入了边界条件的（即删除了被约束的行和列的）。 编号[5] 标题：提取刚度矩阵的问题 地址：https://www.wendangku.net/doc/7f12896761.html,/forum-vi ... fromuid-159019.html 要点：本帖作者的工作是基于单元刚度矩阵的，因此ANSYS中提取的单元刚度矩阵是否处于总体坐标系就成为问题。该问题并非本文内容，但仍值得关注。 编号[6] 标题：提取刚度矩阵丢失节点的问题 地址：https://www.wendangku.net/doc/7f12896761.html,/forum-vi ... fromuid-159019.html 要点：帖子[5]作者的又一帖，在这里帖子[5]的问题得到了欧阳中华老师的回答。 编号[7] 标题：提取刚度矩阵的ANSYS操作过程 地址：https://www.wendangku.net/doc/7f12896761.html,/forum-vi ... fromuid-159019.html 要点：实际上这就是使用HBMAT从ANSYS中提取总体矩阵的全过程！只是还有一些细节待确定。 编号[8] 标题：提取整体刚度矩阵、质量矩阵及阻尼矩阵的简单方法 地址：https://www.wendangku.net/doc/7f12896761.html,/forum-vi ... fromuid-159019.html 要点：给出了利用“不减缩的”子结构方法来得到总体矩阵的方法（这也是网络上常见的代码之一） 编号[9] 标题：质量矩阵、刚度矩阵如何提取？ 地址：https://www.wendangku.net/doc/7f12896761.html,/forum-vi ... fromuid-159019.html 要点：16443在5楼的回帖中给出了提取刚度矩阵的三种方法 =======站外检索略述======================== 百度检索：提取矩阵 比较好的帖子有： 编号[10] 来源:百度文库 标题:怎样从ansys中提取单元刚度矩阵与质量矩阵 地址:https://www.wendangku.net/doc/7f12896761.html,/view/3cf5e567f5335a8102d220d9.html 要点:这应该就是16443在帖子[9]中回复的内容了，全面的总结了在帖子[3,4,5,9]中涉及的问题。 编号[11] 来源:中华钢结构标题:ansys刚度矩阵Harwell-Boeing格式的具体含义讨论 地址:https://www.wendangku.net/doc/7f12896761.html,/forum/viewthread.php?tid=184007 要点:如题，后续楼层给出了一些将矩阵读入ANSYS的APDL（好不容易读出来，又读进去干嘛呢……） 编号[12] 来源:simwe 标题:关于ANSYS（质量、刚度、阻尼）矩阵Harwell-boeing格式数据的说明 地址:https://www.wendangku.net/doc/7f12896761.html,/archiver/tid-924778.html 要点:比[11]更透彻的HB格式说明! ============================================================= =======1.从ANSYS中提取总体矩阵的方法================================= ============================================================= 1,用/DEBUG命令 2,子结构法
Ansys中节点力提取
Ansys中节点力提取几个问题的说明 对于ansys中节点力提取的命令，一般有如下命令可以用， *GET,Par, NODE, N, RF,FX(FY/FZ/MX/MY/MZ) 这组命令是我们最开始用的，用来提取节点反力，但是有个缺陷，节点反力只在有约束位置才能提取，如果在结构中任何一个节点处提取此节点所受合力，界面操作有两种方法。 Main Menu>General Postproc>Nodal Calcs>Total Force Sum Main Menu>General Postproc>Nodal Calcs>Sum @ Each Node 但是执行上面两个操作有个前提，需要选出对应的单元和节点，下面举例说明： 如下图：800臂架结构 由于要对连接架+塔帽进行单独的详细分析，需要提取旋转架与塔帽连接处铰点对塔帽的作用力。而且为了在详细模型中施加载荷的时候方便，提取结果的坐标系需要是X向沿着主臂的局部坐标系，见示图1。 运用Main Menu>General Postproc>Nodal Calcs>Total Force Sum 或者 Main Menu>General Postproc>Nodal Calcs>Sum @ Each Node这两个操作可以实现。 下面就这两个操作的的结果进行比较。就如上模型，研究塔帽和旋转架连接左侧铰耳处受力提取，见示图2。 将塔帽上与此铰点相连接的两个单元选出来，选择此节点，见示图3。 读取结果文件，设置结果坐标系为要求的局部坐标系（文件中为局部坐标系11）。
执行Main Menu>General Postproc>Nodal Calcs>Total Force Sum，选取Active Rsys，结果界面如下： 执行Main Menu>General Postproc>Nodal Calcs>Sum @ Each Node，结果界面如下： 从结果界面上可以看出，了；两个结果界面中显示的结果来看，“结果显示1”与“结果显示3”的结果数据相同，而“结果显示2”与其他两个在在力的显示数据上相同，而矩的显示数据上并不相同。之所以出现这种差别，原因在于在求矩的作用数据时，矩的作用中心不相同。“结果显示2”中所显示的数据其矩的作用中心为提取载荷点（1600078），而“结果数据1”和“结果数据3”中矩的作用中心为结果坐标系的原点。 要想使上述3个结果数据显示值相同，只需要将结果显示中矩的作用中心设置到提取点上就可以了。 通过Main Menu>General Postproc>Nodal Calcs>Summation Pt>At Node将矩的作用点设置到1600078，再执行命令Main Menu>General Postproc>Nodal Calcs>Total Force Sum，结果
ANSYS 钢筋混凝土建模
ANSYS 钢筋混凝土建模 一、简介 钢筋混凝土有限元建模的方法与结果评价（前后处理），是对钢筋混凝土结构进行数值模拟的重要步骤，能否把握模型的可行性、合理性，如何从计算结果中寻找规律，是有限元理论应用于实际工程的关键一环。Blackeage以自己做过的一组钢筋混凝土暗支撑剪力墙的数值模拟为例，从若干方面提出一些经验与建议。希望大家一起讨论、批评指正(wang.jian@https://www.wendangku.net/doc/7f12896761.html,)。 程序：ANSYS 单元：SOLID65、BEAM188 建模方式：分离 暗支撑剪力墙结构由北京工业大学曹万林所提出，简言之就是一种在普通钢筋配筋情况下，加配斜向钢筋的剪力墙结构。 二、单元选择 以前经常采用的钢筋混凝土建模方法是通过SOLID65模拟混凝土，通过SOLID65的实常数指定钢筋配筋率，后来发现这种整体式的模型并不理想，而且将钢筋周围的SOLID65单元选择出来，再换算一个等效的配筋率，工作量也并不小。最关键的是采用整体式模型之后，得不出什么有意义的结论，弄一个荷载-位移曲线出来又和实验值差距比较大。只有计算的开裂荷载与实验还算是比较接近，但这个手算也算得出来的东西费劲去装模作样的建个模型又有什么意义？ 所以，这次我尝试采用分离式的模型，钢筋与混凝土单元分别建模，采用节点共享的方式。建模时发现，只要充分、灵活地运用APDL的技巧，处理好钢筋与混凝土单元节点的位置，效率还是很高的。 暗支撑剪力墙数值模型 看过很多的资料，分离式模型是用LINK8与SOLID65的组合方式，这样做到是非常直观，因为LINK8是spar类型的单元，每个节点有3个自由度，这与SOLID65单元单节点自由度数量是一致的。但是问题也就由此产生，当周围的混凝土开裂或是压碎时，SOLID65将不能对LINK8的节点提供足够地约束（如
ANSYS中混凝土的计算问题
ANSYS中混凝土的计算问题【精华】 最近做了点计算分析，结合各论坛关于这方面的讨论，就一些问题探讨如下，不当之处敬请指正。 一、关于模型 钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、分布式和组合式模型。考虑钢筋和混凝土之间的粘结和滑移，则采用引入粘结单元的分离式模型；假定混凝土和钢筋粘结很好，不考虑二者之间的滑移，则三种模型都可以；分离式和分布式模型适用于二维和三维结构分析，后者对杆系结构分析比较适用。裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型，后者目前尚处研究之中，主要应用的是前两种。离散裂缝模型和分布裂缝模型各有特点，可根据不同的分析目的选择使用。随着计算速度和网格自动划分的快速实现，离散裂缝模型又有被推广使用的趋势。 就ANSYS而言，她可以考虑分离式模型(solid65+link8，认为混凝土和钢筋粘结很好，如要考虑粘结和滑移，则可引入弹簧单元进行模拟，比较困难！)，也可采用分布式模型(带筋的solid65)。而其裂缝的处理方式则为分布裂缝模型。 二、关于本构关系 混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类，其中研究最多的是非线性弹性和弹塑性本构关系，其中不乏实用者。混凝土破坏准则从单参数到五参数模型达数十个模型，或借用古典强度理论或基于试验结果等，各个破坏准则的表达方式和繁简程度各异，适用范围和计算精度差别也比较大，给使用带来了一定的困难。 就ANSYS而言，其问题比较复杂些。 1 ANSYS混凝土的破坏准则与屈服准则是如何定义的？ 采用tb,concr,matnum则定义了W-W破坏准则(failure criterion)，而非屈服准则(yield criterion)。W-W破坏准则是用于检查混凝土开裂和压碎用的，而混凝土的塑性可以另外考虑（当然是在开裂和压碎之前）。理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的，例如在高静水压力下会发生相当的塑性变形，表现为屈服，但没有破坏。而工程上又常将二者等同，其原因是工程结构不容许有很大的塑性变形，且混凝土等材料的屈服点不够明确，但破坏点非常明确。 定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系，即W—W破坏准则、混凝土开裂和压碎前均为线性的应力应变关系，而开裂和压碎后采用其给出的本构关系。但屈服准则尚可另外定义(随材料的应力应变关系，如tb,MKIN,则定义的屈服准则是V on Mises,流动法则、硬化法则也就确定了)。
用ANSYS对钢筋混凝土梁进行计算模拟
一、用钢筋混凝土简支梁的数值模拟为实例，对ANSYS的使用方法进行说明 钢筋混凝土简支梁，尺寸为长2000mm，宽150mm，高300mm。混凝土采用C30，钢筋全部采用HRB335，跨中集中荷载P作用于一刚性垫板上，垫板尺寸为长150mm，宽100mm。 建立分离式有限元模型，混凝土采用SOLID65单元，钢筋采用LINK8单元，不考虑钢筋和混凝土之间的粘结滑移。创建分离式模型时，将几何实体以钢筋位置切开，划分网格时将实体的边线定义为钢筋即可。加载点以均布荷载近似代替钢垫板，支座处则采用线约束和点约束相结合。单元尺寸以50mm左右为宜。 二、命令流 !钢筋混凝土简支梁数值分析 !分离式模型 FINISH /CLEAR /PREP7 !1.定义单元与材料属性 ET,1,SOLID65,,,,,,,1 ET,2,LINK8 MP,EX,1,13585 !混凝土材料的初始弹模以及泊松比 MP,PRXY,1,0.2 FC=14.3 !混凝土单轴抗压强度和单轴抗拉强度 FT=1.43 TB,CONCR,1 TBDA TA,,0.5,0.95,FT,-1 !定义混凝土材料及相关参数，关闭压碎 TB,MISO,1,,11 !定义混凝土应力应变曲线，用MISO模型 TBPT,,0.0002,FC*0.19 TBPT,,0.0004,FC*0.36 TBPT,,0.0006,FC*0.51 TBPT,,0.0008,FC*0.64 TBPT,,0.0010,FC*0.75 TBPT,,0.0012,FC*0.84 TBPT,,0.0014,FC*0.91 TBPT,,0.0016,FC*0.96 TBPT,,0.0018,FC*0.99 TBPT,,0.002,FC TBPT,,0.0033,FC*0.85 MP,EX,2,2.0E5 !钢筋材料的初始弹模以及泊松比 MP,PRXY,2,0.3 TB,BISO,2 TBDA TA,,300,0 !钢筋的应力应变关系，用BISO模型
提取单元刚度矩阵
单元刚度矩阵的提取 刚度矩阵在有限元求解过程中扮演者非常重要的角色，以最小位能原理求解过程为例最终越是转换为含有结构刚度矩阵的能量泛函的取值问题。有限元过程中涉及到三类刚度：单元刚度矩阵，组合结构刚度矩阵和最终求解刚度矩阵。 其中单元刚度矩阵：仅与单元的自身自由度有关，同一编号的单元矩阵的维数是固定。组合结构刚度：矩阵根据求解的初始变量个数决定刚度矩阵的维数，属于单元组装后的初始刚度，维数和整个单元初始变量个数相等。最终求解刚度矩阵：代入边界条件简化后的刚度。以《Finite Element Analysis-Theory and Application With ANSYS》中的梁单元例子为例，解释刚度提取过程： 此模型的单元刚度矩阵：（学则beam3梁单元后，该单元包含两节点，每个节点具有三个自由度，因此对应单元刚度矩阵为6*6的方阵）
组合结构刚度矩阵：（该结构含有三个节点，每个节点具有三个原始自由度，因此组合结构刚度矩阵具有9*9阶的形式） 最终求解刚度矩阵：（由于边界条件的存在，该结构中，1,3点的自由度不存在，求解参数中有六个参数已知，因此对最终求解刚度矩阵为三阶方阵） 通过最终的刚度矩阵组成的方程，求解出2节点的位移解，再以这些原始解得出应力，应变，支反力的其他的解。 ansys实现过程： 提取思路如下：通过/debug提取单元刚度矩阵，通过filname.full文件提取后两者的矩阵 ansys实现过程如下： finish /clear /filname,k,1 /prep7 N,1 N,2,120 N,3,120,-108 et,1,beam3 mp,ex,1,3.0e7 mp,prxy,1,0.3 R,1,7.65,204,10 E,1,2 E,2,3 /debug,-1,,,1,,,,,
混凝土的ANSYS分析
【原创】钢筋混凝土分离式建模方法（含ANSYS命令流） 钢筋混凝土, 分离式, 建模, ANSYS, 命令钢筋混凝土, 分离式, 建模, ANSYS, 命令 一、简介 钢筋混凝土有限元建模的方法与结果评价（前后处理），是对钢筋混凝土结构进行数值模拟的重要步骤，能否把握模型的可行性、合理性，如何从计算结果中寻找规律，是有限元理论应用于实际工程的关键一环。Blackeage以自己做过的一组钢筋混凝土暗支撑剪力墙的数值模拟为例，从若干方面提出一些经验与建议。希望大家一起讨论、批评指正(wang.jian@https://www.wendangku.net/doc/7f12896761.html,)。 程序：ANSYS 单元：SOLID65、BEAM188 建模方式：分离 暗支撑剪力墙结构由北京工业大学曹万林所提出，简言之就是一种在普通钢筋配筋情况下，加配斜向钢筋的剪力墙结构。 二、单元选择 以前经常采用的钢筋混凝土建模方法是通过SOLID65模拟混凝土，通过SOLID65的实常数指定钢筋配筋率，后来发现这种整体式的模型并不理想，而且将钢筋周围的SOLID65单元选择出来，再换算一个等效的配筋率，工作量也并不小。最关键的是采用整体式模型之后，得不出什么有意义的结论，弄一个荷载-位移曲线出来又和实验值差距比较大。只有计算的开裂荷载与实验还算是比较接近，但这个 手算也算得出来的东西费劲去装模作样的建个模型又有什么意义？ 所以，这次我尝试采用分离式的模型，钢筋与混凝土单元分别建模，采用节点共享的方式。建模时发现，只要充分、灵活地运用APDL 的技巧，处理好钢筋与混凝土单元节点的位置，效率还是很高的。 [center] 暗支撑剪力墙数值模型[/center] 看过很多的资料，分离式模型是用LINK8与SOLID65的组合方式，这样做到是非常直观，因为LINK8是spar类型的单元，每个节点有3个自由度，这与SOLID65单元单节点自由度数量是一致的。但是问题也就由此产生，当周围的混凝土开裂或是压碎时，SOLID65将不能对LINK8的节点提供足够地约束（如下图箭头方向），从而导致总刚矩阵小主元地出现影响计算精度，或者干脆形成瞬变体系导致计 算提前发散。 [center] LINK8+SOLID65的问题[/center] 如果采用梁单元模拟暗钢筋，就算包裹钢筋的混凝土破坏了，钢筋单元本身仍可对连接点提供一定的侧向刚度（其实钢筋本身就是有一定抗弯刚度的），保证计算进行下去。ANSYS中的梁单元比较多，建议选取beam188单元。beam188支持弹塑性分析、自定义截面。可以用内力计算结果按截面插值得出应力结果，这样，SOLID65+beam188不仅解决了SOLID65+beam188的小主元问题，而且可以方便地控 制钢筋单元的划分密度，也扩充了钢筋单元输出信息。 三、单元组合方式 将剪力墙中所有钢筋单元（包括暗柱、梁的纵、箍筋、暗支撑钢筋、暗支撑箍筋、暗分布筋）单独建模，为了能够与混凝土单元节点共享，将混凝土单元细化，单元高度设为暗柱箍筋间距与墙片分布筋间距的最大公约数。 钢筋与混凝土单元节点共享。不考虑粘接-滑移影响。其实由于混凝土单元已经细化过了，钢筋周围的混凝土由于钢筋作用而开裂之
基于matlab的有限元法分析平面应力应变问题刘刚
姓名：刘刚学号：15 平面应力应变分析有限元法 Abstruct:本文通过对平面应力/应变问题的简要理论阐述，使读者对要分析的问题有大致的印象，然后结合两个实例，通过MATLAB软件的计算，将有限元分析平面应力/应变问题的过程形象的展示给读者，让人一目了然，快速了解有限元解决这类问题的方法和步骤！ 一.基本理论 有限元法的基本思路和基本原则以结构力学中的位移法为基础，把复杂的结构或连续体看成有限个单元的组合，各单元彼此在节点出连接而组成整体。把连续体分成有限个单元和节点，称为离散化。先对单元进行特性分析，然后根据节点处的平衡和协调条件建立方程，综合后做整体分析。这样一分一合，先离散再综合的过程，就是把复杂结构或连续体的计算问题转化简单单元分析与综合问题。因此，一般的有限揭发包括三个主要步骤：离散化单元分析整体分析。 二.用到的函数 1. LinearTriangleElementStiffness(E,NU,t,xi,yi,xj,yj,xm,ym,p) (K k I f) (k u) (k u A) (E NU t) 三.实例 例1.考虑如图所示的受均布载荷作用的薄平板结构。将平板离散化成两个线性三角元，假定E=200GPa，v=，t=0.025m,w=3000kN/m. 1.离散化 2.写出单元刚度矩阵
通过matlab 的LinearTriangleElementStiffness 函数，得到两个单元刚度矩阵1k 和2k ，每个矩阵都是6 6的。 >> E=210e6 E = >> k1=LinearTriangleElementStiffness(E,NU,t,0,0,,,0,,1) k1 = +006 * Columns 1 through 5 0 0 0 0 0 0 0 0 Column 6 >> NU= NU = >> t= t = >> k2=LinearTriangleElementStiffness(E,NU,t,0,0,,0,,,1)
ANSYS提取单元内力的方法
ANSYS提取单元或节点内力的方法 方法1：节点荷载（List Results→Nodal Loads） 方法2：节点合力计算（Nodal Cals→Sum @ Each Node） 方法3：单元解中的节点解（List Results→Element Solution→Structural Forces & Moments）方法4：支座反力（List Results→Reaction Solu） 方法5：单元表（List Results→Elem Table Data） 上述各方法提取的结果关系如下： （1）方法1和方法2提取的结果完全相同，但结果为0的项在方法1的结果列表中不显示，而方法2的结果列表则会全部显示。 （2）方法3提取的结果是每个单元各节点在该单元中的内力，针对同一节点，将其在各个单元中的内力求和，其累加结果与方法1和2得到的结果一致。 （3）方法4提取的结果只显示有施加位移约束的节点反力，其数值大小与方法1和2得到的结果相差一个正负号，即节点内力和节点反力刚好是一对作用力与反 作用力。 （4）方法5提取的结果是单元的内力，如果单元的形函数为线性（如BEAM188单元设置“KEYOPT(3)=0”），则ANSYS会取单元中点作为积分点并将其数值代 替单元内的线性变化，因此其输出结果的绝对值等于方法3中对应单元的各节 点相应内力绝对值的平均值；如果单元的形函数为非线性（如BEAM188单元 设置“KEYOPT(3)=2”），则单元各节点的内力不同，其结果与方法3得到的结 果一致。 （5）方法1~4提取的结果都是默认基于整体坐标系的，而方法5提取的结果是基于单元坐标系的，因此提取结果的方向和正负号需特别注意。有限元中力的方向 和结构力学中的方向是有区别的，不论是什么结果坐标系，力的正方向取为对 应结果坐标的正方向，弯矩则是对应坐标轴的顺时针为正。
ANSYS岩土计算例子
ANSYS土工结构计算案例 ANSYS－CHINA广州办事处 2019年6月24日
目录 计算题目及计算要求说明 (1) 题目一 (4) 一、计算说明 (4) 二、计算所用ANSYS邓肯－张的E－B模型说明 (5) 三、计算有限元模型及计算结果 (6) 题目二 (7) 一、用三维有限元模型计算 (7) 二、用三维有限元模型计算 (8) 题目三 (10) 一、计算说明 (10) 二、计算有限元模型及边界条件 (10) 三、强夯地基固结计算 (10) 题目四 (17) 一、计算说明 (17) 二、计算几何模型和有限元模型 (17) 三、计算结果 (18) 1、计算边界条件 (18) 2、计算结果 (19) 3、结论 (20)
计算题目及计算要求说明 题目一：高桩码头桩基与岸坡相互作用的线性有限元和非线性有线元分析 题目二：大圆筒结构、波浪与地基的相互作用分析（大圆筒作为重力式码头结构，波浪为水平动荷载，门吊为竖向动荷载，地基为三层以上地基包括自抛碎石垫层、粘土层、粉细砂层和岩层，粉细砂层可能在波浪动荷载作用下液化造成圆筒倾覆） 题目三：（冲击荷载下）强夯地基固结有限元分析（提供固结方程或固结方程处理方案，孔隙水压力消散计算方案、沉降计算方案及其他一些处理技巧） 题目四：在降雨情况下土工格栅加筋土挡墙边坡上公路稳定分析（由上至下为公路面层，垫层，挡墙，挡墙面板采用预制混凝土块0.6?0.6?0.6m3，混凝土后方为钩挂式土工格栅，边坡比较陡，边坡有一定排水特性）。 具体处理方案包括： 1、提供计算输入界面 2、计算模型或采用本构情况 3、前处理方案及网格划分技巧 4、特殊材料或模型嵌入技术 5、计算技巧及解决方案 6、后处理提供内容
ansys动力学分析全套讲解
第一章模态分析 §模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性（固有频率和振型），即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS产品家族中的模态分析是一个线性分析。任何非线性特性，如塑性和接触（间隙）单元，即使定义了也将被忽略。ANSYS提供了七种模态提取方法，它们分别是子空间法、分块Lanczos法、PowerDynamics法、缩减法、非对称法、阻尼法和QR阻尼法。阻尼法和QR阻尼法允许在结构中存在阻尼。后面将详细介绍模态提取方法。 §模态分析中用到的命令 模态分析使用所有其它分析类型相同的命令来建模和进行分析。同样，无论进行何种类型的分析，均可从用户图形界面（GUI）上选择等效于命令的菜单选项来建模和求解问题。 后面的“模态分析实例（命令流或批处理方式）”将给出进行该实例模态分析时要输入的命令（手工或以批处理方式运行ANSYS时）。而“模态分析实例（GUI方式）” 则给出了以从ANSYS GUI中选择菜单选项方式进行同一实例分析的步骤。（要想了解如何使用命令和GUI选项建模，请参阅<
ANSYS混凝土问题分析实例
ANSYS混凝土问题分析学习资料 1．关于模型 钢筋混凝土有限元模型根据钢筋的处理方式分为三种：分离式、整体式和组合式模型◆分离式模型： 把混凝土和钢筋作为不同的单元来处理，即混凝土和钢筋各自被划分为足够小的单元，两者的刚度矩阵是是分开来求解的，考虑到钢筋是一种细长的材料，通常可以忽略起横向抗剪强度，因此可以将钢筋作为线单元处理。钢筋和混凝土之间可以插入粘结单元来模拟钢筋与混凝土之间的粘结和滑移。一般钢筋混凝土是存在裂缝的，而开裂必然导致钢筋和混凝土变形的不协调，也就是说要发生粘结的失效与滑移，所以此种模型的应用最为广泛。 ◆整体式模型： 将钢筋分布与整个单元中，假定混凝土和钢筋粘结很好，并把单元视为连续均匀材料，与分离式模型不同的是，它求出的是综合了混凝土与钢筋单元的整体刚度矩阵；与组合式不同之点在于它不是先分别求出混凝土与钢筋对单元刚度的贡献然后再组合，而是一次求得综合的刚度矩阵。 ◆组合式模型 组合式模型分为两种：一种是分层组合式，在横截面上分成许多混凝土层和若干钢筋层，并对截面的应变作出某些假设，这种组合方式在钢筋混凝土板、壳结构中应用较广；另一种组合方法是采用带钢筋膜的等参单元。 当不考虑混凝土和钢筋二者之间的滑移，三种模型都可以。分离式和整体式模型使用于二维和三维结构分析。 就ANSYS而言，可以考虑分离式模型：混凝土(SOLID65)+钢筋(LINK单元或PIPE单元)，认为混凝土和钢筋粘结很好。如要考虑粘结和滑移，则可引入弹簧单元进行模拟，如果比较困难也可以采用整体式模型(带筋的SOLID65)。 2．本构关系及破坏准则 ◆本构关系 混凝土本构关系的模型对钢筋混凝土结构的非线性分析有重大影响。混凝土的本构就是表示在各种外荷载作用下的混凝土应力应变的响应关系。在建立混凝土本构关系时一般都是基于现有的连续介质力学的本构理论，在结合混凝土的力学特性，确定甚至调整本构关系中各种所需的材料参数。通常，混凝土的本构关系可以分为线性弹性、非线性弹性、弹塑性及其他力学理论等四类。其中研究最多的是非线性弹性和弹塑性本构关系，其他的不怎么用。线性弹性理论认为应力应变加载、卸载时呈线性关系，服从虎克定律，应力应变关系是相互对应的关系。在实际结构设计中线性弹性仍然是应用很广泛的本构模型。 非线性弹性理论认为应力应变不成正比，但是有一一对应的关系。卸载后没有残余应变，应力状态完全由应变状态决定，而与加载历史无关。非线性弹性本构关系分为全量型（如Ottosen模型）和增量型（如Darwin-Pecknold）两类。 弹塑性本构关系则把屈服面和破坏面分开处理。根据混凝土单轴受压的试验研究结果，混凝土在应力未达到其强度极限以前，应力应变的非线性关系受塑性变形的影响，这可以用屈服面理论来解释。而在曲线的下降阶段，混凝土的非线性关系则主要受混凝土内部微断裂的影响，表现微损伤断裂的关系，可用破坏准则来评判。一般在经典的强度理论中，有Tresca、VonMises和Druck-Prager等屈服准则，此外还有Zienkiewicz-Pande、W.F.Chen、Nilsson 屈服条件，破坏准则有Mohr。 混凝土破坏准则从单参数到五参数多大数十个模型，或借用古典强度理论或基于试验结果等。各个破坏准则的表达式和繁简程度各异，适用范围和计算精度也差别较大，给使用带来了一定的困难。