请ABAQUS的高手赐教,关于刚度矩阵的问题 输出单元刚度矩阵的方法。输出单元刚度矩阵的方法。1.用命令:*ELEMENT MATRIX OUTPUT只设定Required parameter:ELSET的话,由于结果文件(*.fil)是二进制文件,用文本编辑器打开是一堆我们看不明白的乱码,所以有必要设置一下文件格式。2.用命令:*FILE FORMAT设定文件格式将文件格式设为ASCII,那么结果文件我们就能够看明白了。3.命令:*ELEMENT MATRIX OUTPUT的其他参数设定:首先FREQUENCY要设定,否则会生成每一步的结果文件,这个文件会很大,举例来说,1800个单元,共50步,全部输出的单元刚度矩阵有600M,这是一个很大的数据,所以要把不需要的删除。并且在打开这个600M大的文件时发现,UntraEdit比记事本和写字板的打开速度快得多。假如共50步,我只想要最后一步的单元刚度矩阵,那么可以设定FREQUENCY=50,设定FREQUENCY=0表示不输出,默认为每一步都输出。其次可以设定结果文件名,设定时要使用OUTPUT FILE=USER DEFINED这个选项,将FILE NAME的值设定为想要的文件名,注意不要扩展名。例如:Elset,Elset=All,Generate1,1800,1File Format,AsciiElement Matrix Output,Elset=All,File Name=abc,Frequency=50,Output File=User Defined,。
区域古构造应力场研究 区域构造背景分析和构造形迹分析是重塑古构造应力场的通用方法,e799bee5baa6e997aee7ad94e59b9ee7ad9431333433616236前者可以定性获得工区的区域构造应力状态,后者则得到局部构造应力场分布的详细信息,两者相结合可以得到特定地质历史时期构造应力场的方向。1.古构造应力确定构造活动的强度反映构造应力数值的大小,通过类比可近似估算出古应力的大小。声发射方法是求取古应力大小较准确的方法(丁元辰,2000,2001),在边界条件已知的情况下,数值模拟的方法是确定古应力场最有效的方法,它不仅可以确定盆地内部详细的古应力数值,还可以确定古应力方向。为了进一步弄清济阳坳陷北带馆陶组的构造特征及其形成机理,需扩展区域范围,以济阳坳陷作为对象来进行研究。(1)古应力方向的确定①区域构造和动力学机制反映的应力场方向。深部地幔物质上涌控制着浅层地壳伸展减薄,产生拉张应力场。在地幔物质持续上升的过程中,若上涌强度过大,地壳上隆幅度大于地壳伸展减薄的幅度,又逢区域构造挤压环境,就会遭受短暂的抬升剥蚀;之后若地幔上涌导致的地壳减薄作用和盆地裂陷伸展作用占较大优势,盆地在短暂抬升之后会发生伸展裂陷作用。地幔上涌过程中可能会形成多次。
ANSYS中怎样提取刚度矩阵 用!命令正确建模并且加上约束条件,然后使用下面的命令流:soluantype,7。substructuring分析类型seopt,matname,2。设置文件名称和刚度矩阵类型(刚度,质量,阻尼等)!matname是设置的输出矩阵名称,可以用其他字母替换。数字1代表输出刚度矩阵,2代表输出刚度和质量矩阵,3代表输出刚度、质量和阻尼矩阵。这些都是整体的。nsel,all。选择所有节点m,all,all。定义所有节点自由度为主自由度solve。求解selist,matname,3。列出所需要的整体刚度矩阵
我获得ansys分析 桁架结构的单元刚度矩阵和质量刚度矩阵,用贡献矩阵法组装整体矩阵,如何得到贡献矩阵? 这里就是要把单元坐标系,换到总体坐标系,用变换矩阵做一下相应的变化,就可以通过单元刚度矩阵,得到贡献矩阵。最后再叠加就可以了。
跪求 求解某一桁架或者刚架(随便某一问题的都可以,简单点的也可以)的MATLAB程序,要详细的程序 你提出的问题我之前刚好做过,使用有限元方法来进行桁架结构分析。Matlab编程实现平面杆单元分析首先,明确Matlab程序要实现的5个重要模块分别为:单元刚度矩阵的求解、单元组装、节点位移的求解、单元应力的求解、节点力的求解。下面给出这5个模块的实现。1.单元刚度矩阵求解定义函数Bar2D2Node_Stiffness,该函数计算单元的刚度矩阵,输入弹性模量E,横截面积A,两个节点坐标输出单元刚度矩阵k(4X4)。具体代码如下:function k=Bar2D2Node_Stiffness(E,A,x1,y1,x2,y2)L=sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));x=acos((x2-x1)/L);C=cos(x);S=sin(x);k=E*A/L*[C*C C*S-C*C-C*S;C*S S*S-C*S-S*S;C*C-C*S C*C C*S;C*S-S*S C*S S*S];2.单元组装定义函数Bar2D2Node_Assembly,该函数进行单元刚度矩阵的组装,输入单元刚度矩阵k,单元的节点编号i、j。输出整体刚度矩阵KK,具体代码如下:function z=Bar2D2Node_Assembly(KK,k,i,j)DOF(1)=2*i-1;DOF(2)=2*i;DOF(3)=2*j-1;DOF(4)=2*j;for n1=1:4for n2=1:4KK(DOF(n1),DOF(n2))=KK(DOF(n1),DOF(n2))+k(n1,n2);endendz=KK;3.节点位移的求解定义函数Bar2D2Node_Disp(KK,num,p),该函数输入KK为总体刚度矩阵;num为活动。
你好,请问一下桁架结构带未知参数的刚度矩阵如何在matlab中实现? 你提出的问题我之前刚好做过,使用有限元方法来进行桁架结构分析.Matlab编程实现平面杆单元分析首先,明确Matlab程序要实现的5个重要模块分别为:单元刚度矩阵的求解、单元组装、节点位移的求解、单元应力的求解、节点力的求解.下面给出这5个模块的实现.1.\\x09单元刚度矩阵求解定义函数Bar2D2Node_Stiffness,该函数计算单元的刚度矩阵,输入弹性模量E,横截面积A,两个节点坐标输出单元刚度矩阵k(4X4).具体代码如下:function k=Bar2D2Node_Stiffness(E,A,x1,y1,x2,y2)L=sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));x=acos((x2-x1)/L);C=cos(x);S=sin(x);k=E*A/L*[C*C C*S-C*C-C*S;C*S S*S-C*S-S*S;C*C-C*S C*C C*S;C*S-S*S C*S S*S];2.\\x09单元组装定义函数Bar2D2Node_Assembly,该函数进行单元刚度矩阵的组装,输入单元刚度矩阵k,单元的节点编号i、j.输出整体刚度矩阵KK,具体代码如下:function z=Bar2D2Node_Assembly(KK,k,i,j)DOF(1)=2*i-1;DOF(2)=2*i;DOF(3)=2*j-1;DOF(4)=2*j;for n1=1:4for n2=1:4KK(DOF(n1),DOF(n2))=KK(DOF(n1),DOF(n2))+k(n1,n2);endendz=KK;3.\\x09节点位移的求解定义函数Bar2D2Node_Disp(KK,num,p),该函数输入KK为总体刚度矩阵;num。
ABAQUS怎么输出刚度矩阵? 楼上说得有道理。这是个8*8矩阵,根据你的模型有25个单元来看,这种矩阵应该有25个。为啥是8*8矩阵呢,因为一个单元四个节点,每个结点两个自由度4*2=8.
ABAQUS做模态分析,怎样提取模态振型? 不是要振型云图,而是模态振型矢量。最好是GUI操作,关于命令流我不是很熟悉可以的话再说一下质量和刚度矩阵的提取方法。谢谢大家了