有限元命令流与部分基础知识.docVIP

一、命令流举例： 有一长为100mm的矩形截面梁，截面为10X1mm，与一规格为 20mmX7mmX10mm的实体连接，约束实体的端面，在梁端施加大小为3N的y方向的压力，梁与实体都为一材料，弹性模量为30Gpa，泊松比为0.3。本例主要讲解梁与实体连接处如何利用耦合及约束方程进行处理。 命令流如下：FINI/CLE/FILNAME,BEAM_AND_SOLID_ELEMENTS_CONNECTION !定义工作文件名/TITLE,COUPLE_AND_CONSTRAINT_EQUATION??? !定义工作名/PREP7 !进入前处理ET,1,SOLID95??????? !定义实体单元类型为SOLID95ET,2,BEAM4??????? !定义梁单元类型为BEAM4MP,EX,1,3E4??????? !定义材料的弹性模量MP,PRXY,1,0.3??????? !定义泊松比R,1????????? !定义实体单元实常数R,2,10.0,10/12.0,1000/12.0,10.0,1.0?? !定义梁单元实常数BLC4,,,20,7,10??????? !创建矩形块为实体模型WPOFFS,0,3.5??????? !将工作平面向Y方向移动3.5WPROTA,0,90???????? !将工作平面绕X轴旋转90度VSBW,ALL???????? !将实体沿工作平面剖开WPOFFS,0,5???????? !将工作平面向Y方向移动5??WPROTA,0,90???????? !将工作平面绕X轴旋转90度VSBW,ALL??????? !将实体沿工作平面剖开WPCSYS,-1??????? !将工作平面设为与总体笛卡儿坐标一致K,100,20,3.5,5?????? !创建关键点K,101,120,3.5,5?????? !创建关键点L,100,101???????? !连接关键点生成梁的线实体LSEL,S,LOC,X,21,130????? !选择梁线LATT,1,2,2??????? !指定梁的单元属性LESIZE,ALL,,,10?????? !指定梁上的单元份数LMESH,ALL??????? !划分梁单元VSEL,ALL???????? !选择所有实体VATT,1,1,1??????? !设置实体的单元属性ESIZE,1???????? !指定实体单元尺寸MSHAPE,0,2D??????? !设置实体单元为2DMSHKEY,1???????? !设置为映射网格划分方法VMESH,ALL??????? !划分实体单元ALLS???????? !全选FINI????????? !退出前处理!------------------------ /SOLU????????? !进入求解器ASEL,S,LOC,X,0???????? !选择实体的端面DA,ALL,ALL????????? !约束实体端面ALLS???????? !全选FK,101,FY,-3.0??????? !在两端施加Y向压力CP,1,UX,1,21??????? !耦合节点1和节点21X方向自由度CP,2,UY,1,21?????? !耦合节点1和节点21Y方向自由度CP,3,UZ,1,21??????? !耦合节点1和节点21Z方向自由度CE,1,0,626,UX,1,2328,UX,-1,1,ROTY,-ABS(NZ(626)-NZ(2328)) !设置约束方程CE,2,0,67,UX,1,4283,UX,-1,1,ROTZ,-ABS(NY(67)-NY(4283))?? !设置约束方程CE,3,0,67,UZ,1,4283,UZ,-1,1,ROTX,-ABS(NY(67)-NY(4283))??? !设置约束方程ALLS???????? !全选SOLVE???????? !保存 FINI???????? !退出求解器!------------------------ /POST1???????? !进入通用后处理PLNSOL, U,Y, 0,1.0????? !显示Y方向位移PLNSOL, S,EQV, 0,1.0???? !显示等效应力ETABLE,ZL1,SMISC,1????? !读取梁单元上I节点X方向的力ETABLE,ZL2,SMISC,7????? !读取梁单元上J节点X方向的力ETABLE,MZ1,SMISC,6????? !读取梁单元上I节点Z方向的力矩ETABLE,MZ2,SMISC,12???? !读取梁单元上J节点Z方向的力矩PLETAB,ZL1?????? !显示梁单元X方向的力PLETAB,MZ1????