matlab刚度矩阵求逆,ANSYS求结构整体刚度矩阵逆矩阵APDL命令流
虽然采用matlab等相关软件对于大型矩阵的求解非常方便，但是这里涉及到ANSYS与matlab接口的处理问题，比较繁琐，本文编制一个直接在ANSYS中提取结构整体刚度矩阵和对整体刚度求逆的APDL程序，能在ANSYS中直接操作刚度矩阵，比较实用。
利用ANSYS求结构刚度矩阵逆矩阵，模型采用单层球面网壳！
相关的APDL命令流如下所示。
finish
/clear
/filname,modal
/prep7
et,1,beam188
mp,ex,1,2.1e11
mp,prxy,1,0.28
mp,dens,1,7850
*afun,deg
!用户界面设计，输入几何参数
multipro,'start',4
*cset,1,3,f,'Enter the f',6.8
*cset,4,6,span,'Enter the span',30
*cset,7,9,Kn,'Enter the radial number',24
*cset,10,12,Nx,'Enter the node circle number',3
*cset,61,62,'please enter the geometry parameters'
multipro,'end'
!计算节点坐标位置，并定义节点
csys,2
r=(f*f+span*span/4)/(2*f)
Dalfa=atn(span/2/sqrt(r*r-(span/2)*(span/2)))/Nx
n,1,r,0,90
*do,i,1,Nx
*do,j,1,Kn
x=r
y=(j-1)*360/Kn
z=90-i*Dalfa
n,1+Kn*(i-1)+j,x,y,z,
*enddo
*enddo
!定义单元连接
ro=0.219/2
ri=(0.219-0.010*2)/2
sectype,1,beam,ctube,,,
secdata,ri,ro,8
type,1
mat,1
secnum,1
!定义环向杆连接
*do,i,1,Nx
*do,j,1,Kn-1,1
e,1+Kn*(i-1)+j,1+Kn*(i-1)+j+1
*enddo
e,1+Kn*i,1+i*Kn-Kn+1
*enddo
!定义径向杆连接
*do,j,1,Kn,1
e,1,1+j,
*enddo
*do,i,1,Nx-1
*do,j,1,Kn
e,1+Kn*(i-1)+j,1+Kn*i+j
*enddo
*enddo
!定义斜杆连接
*do,i,1,Nx-1
*do,j,1,Kn-1,1
e,1+Kn*(i-1)+j,1+Kn*i+j+1,1
*enddo
e,1+i*Kn,1+i*Kn+1
*enddo
!定义边界约束节点和节点荷载
*do,i,1,1+Kn*Nx
*if,i,le,1+Kn*(Nx-1),then
f,i,fz,-1000
*else
d,i,all
*endif
*enddo
finish
/solu
solve
/aux2
file,modal,full
hbmat,modal,txt,,ascII,stiff,yes !转换hb.file文件hb.txt
finish
*DIM,CONTLINE,,5 !定义一维数组
*VREAD,CONTLINE(1),MODAL,TXT,,,5,,,1 !跳过第1行后读入5个数据
(5F14.0)
PTRCRD=CONTLINE(2) !保存列指针总行数
INDCRD=CONTLINE(3) !保存行索引总行数
VALCRD=CONTLINE(4) !保存矩阵元素总行数
RHSCRD=CONTLINE(5) !保存右边项总行数
*VREAD,CONTLINE(1),MODAL,TXT,,,4,,,2
matlab刚度矩阵求逆,ANSYS求结构整体刚度矩阵逆矩阵APDL命令流相关推荐
- 桁架结构的APDL命令流
ANSYS参数化设计语言(APDL)是一种用来完成有限元常规分析操作或通过参数化变量方式建立分析模型的脚本语言,它用智能化分析的手段,为用户提供了自动完成有限元分析过程的功能,即程序的输入可根 ...
- ANSYS APDL 命令流输入方法_51CAE_新浪博客
APDL的全称是ANSYS Parametric Design Language(ANSYS参数化设计语言),是ANSYS的一种高级应用,具有很多图形界面无法实现的功能,广泛应用于优化设计.焊接模拟等 ...
- ANSYS apdl命令流笔记13-------路径的定义与显示
路径的定义与结果数据显示 一.路径的定义 1.PATH 2.PPATH 二.结果数据在路径上显示 1.PDEF 2.PLPATH 三.举个例子 总结 一.路径的定义 1.PATH 使用功能:定义一个路 ...
- ANSYS apdl命令流笔记8---载荷步与载荷子步的运用
载荷步与载荷子步的运用以及输出控制 前言 一.有关命令 1.TIME命令 2.NSUBST命令 3.KBC命令 4.OUTRES命令 二.需要获得结果的几种情况 1.获取不同加载情况下的结果 2.获取 ...
- ANSYS apdl命令流笔记7---对称与反对称、几何与有限元模型属性的赋予、单元表指定的意义
三个基础知识点 一.对称边界条件与反对称边界条件 1.对称边界条件 2.反对称边界条件 二.几何模型与有限元模型属性指定命令 2.1.几何模型 2.1.1"KATT"命令 2.1. ...
- ANSYS apdl命令流笔记4----非线性材料的定义
定义材料的非线性或特殊材料 前言 一.和温度相关的材料 1.MP命令 2.MPTEMP命令 3.MPDATA 命令 4.举例说明 二.非线性材料及相关命令 1.TB命令 2.TBTEMP命令 3.TB ...
- ANSYS apdl命令流耦合场分析案例4--------三维感应加热
三维电磁感应加热 前言 完整感应加热命令 总结 前言 感应加热的激励源为365000HZ的交流电,线圈电流密度为2.04e8A/m^2,线圈和管子的几何模型如下图所示,另外再在外面在加上一个空气域. ...
- ANSYS apdl命令流瞬态动力学分析案例5--------凸轮从动件运动分析
凸轮从动件分析 问题描述 一.命令分解 1.选择单元和材料属性 2.建立模型和划分网格 3.施加载荷 4.后处理 5.完整命令 总结 问题描述 图1所示为一对心直动尖底从动件盘形凸轮机构,从动件位移s ...
- ANSYS apdl命令流笔记6--------生成线的10种方法
线的生成 1.LSTR命令 2.L命令 3.LAREA命令 4.LTAN与L2TAN命令 4.1 LTAN命令 4.2 L2TAN命令 5.LANG命令与L2ANG命令 5.1 LANG命令 5.2 ...
- 【ANSYS Workbench仿真】如果后处理添加命令流显示渗透压力？
如果后处理添加命令流显示渗透压力? 命令流 /show,png set,2,1 !2表示选择step 2,1或last代表该step的第一个或最后一个substep plesol,cont,fprs ...
最新文章
- iOS 标签自动布局
- PowerDesigner导入MySQL数据库表
- PYTORCH：DenseNet做新冠肺炎CT照片是否确诊分类
- 92行JavaScript代码实现的类似Excel外观的UI
- SAP CRM WebClient UI配置的加载逻辑
- matlab怎么实现OFDM仿真,OFDMmatlab实现仿真源代码
- 疑似谷歌Pixel 4真机谍照曝光：边框宽到没朋友
- 1000道Python题库系列分享25（40道Pandas客观题）
- 四川大专计算机专业排行,四川排名前三中等职业技术学院专业排行
- 点击按钮弹出Outlook（新建邮件状态）
- 总结---JavaScript数组
- FLTK--轻量级C++跨平台GUI库
- HTML的字体10种酷炫效果
- 18年12月英语六级第二套听力单词
- 关于普中科技程序烧录软件PZ-ISP与驱动安装程序SETUP.exe之间血与泪的教训
- 【HTML+CSS】移动端布局之流式布局
- c语言感叹号和逗号标红,感叹号不够用？试试感叹逗号
- 2014年总结——回眸过去
- python编写的软件可以申请专利吗_Python爬虫 | 爬取同一公司用不同名字申请专利的那些Assignees...
- k叉树的性质_二叉树的性质总结.doc