怎么把四边形单元刚度矩阵阵对应到总刚度矩阵

来源：蜘蛛抓取(WebSpider) 时间：2015-12-12 09:48 标签： 单元刚度矩阵
后使用快捷导航没有帐号？
只需一步，快速开始
请教ANSYS中beam3的单元刚度矩阵和质量矩阵是如何定义的？
一根简支梁，长6m，宽0.25m，高0.2m，密度2500，E＝3.2e10Pa。取平面二维单元beam3，全梁分为15个单元分析，计算结构的固有频率和振型。命令流如下： &BR& ＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝ &BR& /filname,Mode analysis &BR& /title,Mode analysis of a Simply Beam.&&/prep7&&et,1,beam3&&mp,ex,1,3.2e10 &BR& mp,dens,1,2500 &BR& mp,nuxy,1,0.20 &BR& r,1,0.05,1.66e-4,0.2&&k, &BR& k,,6 &BR& l,1,2,&&lesize,1,,,15 &BR& latt,1,1,1,,,,1 &BR& lmesh,all &BR& save&&!*************************************** &BR& /OUTPUT,cp,out,, ! 将输出信息送到cp.out文件&&&BR& /debug,-1,,,1 ! 指定输出单元矩阵&&&BR& !***************************************&&/solu &BR& antype,modal &BR& modopt,lanb,10 &BR& mxpand,10 &BR& nsel,s,node,,1 &BR& d,1,ux,0 &BR& d,1,uy,0 &BR& nsel,s,node,,2 &BR& d,2,uy,0 &BR& allsel,all &BR& solve &BR& fini&&/post1 &BR& set,list &BR& ＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝ &BR& 现在的问题是，根据ansys计算的单元刚度矩阵和质量矩阵如下图片1（未消除约束前，各单元的K和M应该是一样的，此处取第8号单元）； &BR&& &&BR& 而很多参考书介绍的二维梁单元（每节点三自由度）的单元刚度和质量矩阵见图片2（质量矩阵为一致质量矩阵，而集中质量矩阵为对角线相应元素为非零的对角阵，在此没有列出）。 &BR&& &&BR&& &&BR& 通过比较本例的ansys计算结果和用图片2的计算公式得出的结果（图片3）可以发现： &BR& （1）二者的单元刚度矩阵刚度矩阵基本一致； &BR& （2）ansys的质量矩阵采用的是一致质量矩阵，但二者的质量矩阵数值结果却存在较大差别。我想这种差别不应该是计算精度误差导致的，所以想请高手帮看看是什么原因？那么在ansys中的质量矩阵究竟是取什么样的形式？ansys帮助文件中能查到各单元的具体矩阵显式表达吗？ &BR&
(100.16 KB, 下载次数: 13)
10:13 上传
图片2 &BR&
(29.15 KB, 下载次数: 30)
21:57 上传
图片3 &BR&
(58.65 KB, 下载次数: 14)
16:05 上传
工作时间：8:00-24:00
Copyright & 2015 -星级专家为您解答难题！
800-820-5114
您尚未登录，请先
遇到问题怎么办？赶快来提问！研发平台专家为您解答！
问题解决中
问题提出时间： 14:09:43
请教刚度矩阵和柔度矩阵的问题
梁的单元刚度矩阵是12×12的，但用mathematica算他的逆矩阵确算不了，我想应该是刚度矩阵是奇异的，k的行列式为0。那么如果我要用梁单元的柔度矩阵该怎么弄？
行业分类：其他专业技术
回答时间： 13:01:26
自由式单刚当然不可逆
只有约束单元才能从刚度矩阵求逆得到柔度矩阵
可以通过振型和频率来求柔度矩阵
共0条评论&
回答时间： 16:59:26
自由式单刚当然不可逆
只有约束单元才能从刚度矩阵求逆得到柔度矩阵
可以通过振型和频率来求柔度矩阵
共0条评论&
回答时间： 19:11:33
自由式单刚当然不可逆
只有约束单元才能从刚度矩阵求逆得到柔度矩阵
可以通过振型和频率来求柔度矩阵
共0条评论&
这个问题已提出一个月以上，回答功能关闭，如需求助可重新您好，请问一下，对公式变分、推导单元刚度、总刚度矩阵，然后迭代，编写的有限元程序可用matlab计算吗？_百度知道
您好，请问一下，对公式变分、推导单元刚度、总刚度矩阵，然后迭代，编写的有限元程序可用matlab计算吗？
您没限元程序呢我知道编写请您帮忙解答谢谢
我有更好的答案
知想要MATLAB编写源限元求解程序或想用MATLAB发限元程序发没问题源简单教程序网想工程应用找商用软件吧
我是需要编写有限元计算程序。因为理论推导出来了几个公式，要把这几个公式应用到有限元计算中去。这几个公式需要变分，并采用分部积分法推导出8节点等参元有限元法的单元刚度矩阵和总刚度矩阵，然后采用Newton-Raphoson迭代法进行求解。我就是不理解这个计算是通过matlab计算还是编写命令流程序使用ansys计算？
其他类似问题
为您推荐：
迭代的相关知识
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁单元刚度矩阵与整体刚度矩阵
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。在ansys中如何得到单元的刚度矩阵？
在ansys中如何得到单元的刚度矩阵？
08-12-12 &
资料来源：最近在计算一个模型管状在地下分析中要施加地层给予的支撑我想简化为弹性支撑可是发现在ansys 中找不到这样的支撑情况 -用弹簧可是问题很多因为是一个管体整个管壁都有支撑简化为弹簧出现的问题很多了比如用多少个合适隔多长合适线性支撑和非线性支撑能等同吗我最近看校验手册发现V135 这个例子可以借鉴可是还是有许多麻烦哪位大侠如果有这方面经验请指点指点 -如何梁单元同时施加两个方向分布荷载竖向侧向 -请问如何打印图形结果 -ansys5.6.2 计算一个压力容器,划分网格9 万余个.运行环境win2000,cpu:1G,内存640M,硬盘40G,模型建立,网格划分均无问题.但是运行到占用硬盘大约9.2G 的时候,就出错.提示内存无法read.请问该如何解决?是内存的问题?还是ansys 本身问题?还是系统问题?硬件问题?模型建立实在无法简化. -我在求解是出现如下的提示 1. real constant 1 referenced by at least element types 1 and 2. 2. 2.coefficient ratio exceeds 1.0e8. 我全部忽略上述提示得出了结果结果能可信吗我有错在哪里呢模型我都检查了数遍了还是和初始文件名有关我是在一个初始名下面建立若干相关的模型 由于我的单元只是板梁的结合所以我重新检查了一下删除了一个影响不大的扁梁很扁尺寸不大系统可能认为是板了运行就通过了上述的提示可能是因为系统对某些单元有默认的尺寸限制 -slab on girder highway bridge 桥板为SOLID 单元由于girder 上部flange 嵌贴在板的下表面所以建模时我将桥板定义为SOLID 单元桥板下表面AREA 定义为SHELL63 单元但SOLID 和SHELL 单元的自由度不同这样建模肯定有不对的地方运算结果与实验结果也相差了40% 有什么好的解决办法 -Q ansys 有没有评价一个网格划分好坏的参数或命令请告知 -Q (1)在gui 下运行宏时,重新定义原来已有的数组,ansys 就会有个提示,想不要这个提示,怎么办?(/uis,msgpop,3 这个行不通,我试过了) (2)ansys 中数组的定义为a*b*c a&10040 b&255 c&7 而我实际上要定义更大的有什么好办法,难道非要写到文件里再读进来? -Q 我现在做一个板壳的屈曲分析想要得到该结构在荷载作用下的极限承载力利用非线性屈曲分析请问极限承载力的荷载比例因子在什么命令下可以查到 -Q 如果材料一定为各向异性我用的材料为钢板该如何选用单元类型 ANSYS使用问答精华Q为问A为回复-------------------------------------------------------------------------------Q:模态分析得到的结果是不是某个方向上的各阶频率啊我要得到各个方向的一阶频率能做到吗A模态分析得到的结果是你所选择的自由度内的振型如果需要得到某一个方向内的振型用RUDECED法选择你所希望的主自由度如UX将所求的频率值设置为一即可Q不过大型复杂结构的振型一般都不是一个方向的应该是一个方向为主其它方向为次如果把所要求的方向设为主自由度会不会跟真是实际结果有出入Q对于复杂的大型实体其主自由度无法给出此时就无法应用reduce法可不可以改约束还是就认为一次得到前3阶频率就是x,y,z方向的一阶频率啊A正是因为复杂结构的主自由度不好确定所以我还是倾向于用SUBSPACE和BLOCK LANCZOS法约束不要随便改它应该描述结构真实的状态不能说前3阶频率即为XYZ方向的第一阶频率你可以好好想想所用方程中[K]和[M]的意义
请登录后再发表评论!

我要回帖

更多关于 单元刚度矩阵 的文章

 

随机推荐