梁实体接触.docVIP

梁实体接触.txt36母爱是一缕阳光，让你的心灵即便在寒冷的冬天也能感受到温暖如春；母爱是一泓清泉，让你的情感即使蒙上岁月的风尘仍然清澈澄净。梁实体接触.txt我退化了，到现在我还不会游泳，要知道在我出生之前，我绝对是游的最快的那个/PREP7 ET,1,SOLID95 !定义实体单元类型为SOLID95 ET,2,BEAM4 !定义梁单元类型为BEAM4 MP,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.5 WPROTA,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 MSHKEY,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!显示梁单元Z方向力矩 /prep7 et,1,solid95 !实体单元 mp,ex,1,3.5e4 mp,nuxy,1,0.2 mp,dens,1,2600 et,2,beam44 !梁单元 mp,ex,2,3.2e4 mp,dens,2,2500 mp,prxy,2,0.1667 et,3,175 !contact175 接触单元，准备做梁-实体的MPC KEYOPT,3,2,2 !使用MPC多点约束算法 KEYOPT,3,12,5 !设置接触行为是；绑定粘结 et,4,170 !Target170 r,1 !空实常数，用于赋给实体 r,2 !空实常数，用于识别接触对 block,-3,0,0,2,-4,4 !创建实体 type,1, mat,1 real,1 esize,0.5 mshape,0 mshkey,1 vmesh,all *get,nodemax,node,0,num,max !创建梁 n,nodemax+1,0.1,1,0 ngen,6,1,nodemax+1,,,1 type,2 mat,2 sectype,1,beam,csolid,pier secdata,0.3 secnum,1 *do,i,1,5 e,nodemax+i,nodemax+i+1 *enddo