声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1694|回复: 7

新手急求老师指点

[复制链接]
发表于 2006-4-24 10:53 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
<P>下面是我作得一个小小例子但是老是有问题。我想在立方体的侧面施加变载荷压力,用得表面效应单元。但是每次运行只能施加到一个单元上。不知道那个地方的问题,请老师帮忙看看。我在这里表示感谢。<BR>/BATCH  <BR>/COM,ANSYS RELEASE  9.0    UP20041104       08:50:38    04/23/2006<BR>!*  <BR>/NOPR   <BR>/PMETH,OFF,0<BR>KEYW,PR_SET,1   <BR>KEYW,PR_STRUC,1 <BR>KEYW,PR_THERM,0 <BR>KEYW,PR_FLUID,0 <BR>KEYW,PR_ELMAG,0 <BR>KEYW,MAGNOD,0   <BR>KEYW,MAGEDG,0   <BR>KEYW,MAGHFE,0   <BR>KEYW,MAGELC,0   <BR>KEYW,PR_MULTI,0 <BR>KEYW,PR_CFD,0   <BR>/GO <BR>!*  <BR>/COM,   <BR>/COM,Preferences for GUI filtering have been set to display:<BR>/COM,  Structural   <BR>!*  <BR>/PREP7  <BR>!*  <BR>ET,1,SOLID92<BR>!*  <BR>ET,2,SURF154<BR>!*  <BR>!*  <BR>MPTEMP,,,,,,,,  <BR>MPTEMP,1,0  <BR>MPDATA,EX,1,,2e11   <BR>MPDATA,PRXY,1,,0.35 <BR>BLOCK,0,2,0,2,0,2,  <BR>TYPE,   1   <BR>MAT,       1<BR>REAL,   <BR>ESYS,       0   <BR>SECNUM, <BR>TSHAP,LINE  <BR>!*  <BR>MSHKEY,0<BR>MSHAPE,1,3d <BR>CM,_Y,VOLU  <BR>VSEL, , , ,       1 <BR>CM,_Y1,VOLU <BR>CHKMSH,'VOLU'   <BR>CMSEL,S,_Y  <BR>!*  <BR>VMESH,_Y1   <BR>!*  <BR>CMDELE,_Y   <BR>CMDELE,_Y1  <BR>CMDELE,_Y2  <BR>!*  <BR>SAVE</P>
<P>ASEL,S,LOC,X,2<BR>type,2<BR>amesh,all<BR>ESEL,S,TYPE,,2</P>
<P>*GET,ECOUNT,ELEM,,COUNT<BR>*DIM,ENUM,ARRAY,ECOUNT<BR>*GET,ENUM(1),ELEM,,NUM,MIN<BR>*DO,I,2,ECOUNT<BR>   ENUM(I)=ELNEXT(ENUM(I-1))<BR>*ENDDO</P>
<P>C2=5<BR>*DO,I,1,ECOUNT<BR>  ESEL,S,,,ENUM(I)<BR>  NSLE,R<BR>  *GET,N_E,NODE,,COUNT<BR>   <BR>  <BR>  YN=NY(NELEM(ENUM(I),1))  <BR>  PN1=C2*YN<BR>  <BR>  <BR>  YN=NY(NELEM(ENUM(I),2))  <BR>  PN2=C2*YN<BR>   <BR>  <BR>  YN=NY(NELEM(ENUM(I),3))  <BR>  PN3=C2*YN<BR> <BR>  PN4=0<BR>  *IF,N_E,EQ,4,THEN<BR>    <BR>     YN=NY(NELEM(ENUM(I),4))     <BR>     PN4=C2*YN<BR>   *ENDIF<BR>  <BR>   SFE,ENUM(I),1,PRES,,PN1,PN2,PN3,PN4<BR>*ENDDO</P>
<P>/PSF,PRES,NORM,2,1</P>
回复
分享到:

使用道具 举报

 楼主| 发表于 2006-4-27 08:55 | 显示全部楼层

大家帮帮吧真的很急。

帮帮忙吧很急用。真的很感谢
 楼主| 发表于 2006-4-28 13:37 | 显示全部楼层
<P>高手给看看吧,指点一下吧<BR>[em04]</P>
发表于 2006-4-30 16:37 | 显示全部楼层
我很同情你,但没办法,我是菜鸟
发表于 2006-5-1 20:53 | 显示全部楼层
我也是菜鸟,你的心情我也理解
发表于 2006-5-2 07:40 | 显示全部楼层

回复:(wz1105)新手急求老师指点

这是一个加表面效应单元的例子,你看一下有没有帮助。 <BR>FINISH  <BR>/CLEAR <BR>!!!gauss加载 <BR>/title, the simulate of laser cladding process <BR>*SET,P,2000,       !the power of the laser <BR>*SET,V,0.008       !scanning speed   <BR>!*SET, RADIUS,0.003       !radius of the laser beam <BR>*SET, R,0.003 <BR>*set,k,0.6          !the coefficient of the absorb <BR>*SET,L,0.05          !the length of the block <BR>/prep7 <BR>!熔覆层高5mm=0.005m <BR>!基体高10mm=0.01m,长x=50mm=0.05,宽=20mm=0.02m <BR>BLOCK,0,0.05,0, -0.0105,-0.01,0, <BR>Vsel,all <BR>vplot <BR>! <BR>/USER, 1 <BR>/VIEW, 1, -0.198026943995E-01,  0.382677277552    ,  0.923669829831  <BR>/ANG,  1, -0.107453019382    <BR>/REPLO <BR>! <BR>ET,1,SOLID70 <BR>ET,2,Surf152  <BR>DOF,TEMP <BR>KEYOPT,2,3,0 <BR>KEYOPT,2,4,1   <BR>KEYOPT,2,5,0 <BR>KEYOPT,2,6,0 <BR>KEYOPT,2,7,0 <BR>KEYOPT,2,8,4      !根据固体表面温度计算对流系数 <BR>KEYOPT,2,9,0 <BR>!用温度表方式定义材料性能 <BR>Mp,kxx,1,82.9 <BR>Mp,dens,1,8900   <BR>Mp,c,1,471                <BR>mptemp,1,20,300,600,750,900,1200,1500,1800,2000         <BR>mpdata,hf,1,1,6,50,120,180,200,250,378,700,850 <BR>! <BR>! <BR>Mshkey,0        !自由网格化 <BR>Mshape,1,3d <BR>type,1 <BR>esize,0.002,          <BR>Vmesh,all <BR>eplot <BR>!  <BR>! <BR>asel,s,loc,y,0          !对流边界条件要施加在5个面的周围 <BR>asel,a,loc,z,-0.01 <BR>asel,a,loc,x,0 <BR>asel,a,loc,x,0.05 <BR>type,2         !定义属性 <BR>amesh,all <BR>esel,s,type,,2         !选择表面效应单元 <BR>sf,all,conv,-1,25         !确定对流边界条件,表格加载 <BR>!sf,all,conv,200,25         !确定对流边界条件 <BR>!!底面也施加轻微的对流 <BR>allsel,all <BR>asel,s,loc,y, -0.0105 <BR>type,2         !定义属性 <BR>amesh,all <BR>esel,s,type,,2 <BR>sf,all,conv,3,25 <BR>allsel,all <BR>! <BR>nsel,all <BR>/psf,hflux,,2,0        !以箭头方式显示载荷 <BR>! <BR>finish <BR>/SOLU  <BR>!!!!! <BR>antype,transient,new <BR>TRNOPT,FULL  <BR>Nropt,,,off            !关闭自适应下降 <BR>LUMPM,0 <BR>Solcontrol,on <BR>Autots,on <BR>Kbc,1 <BR>Outres,all,last, <BR>Lnsrch,on <BR>neqit,50          !每一子步中方程的迭代次数限值  <BR>!load step1,initial conditions 25 <BR>time,0.001 <BR>deltim,0.001,0.001,0.001             <BR>allsel,all <BR>timint,off <BR>tunif,25   !均匀的节点初始温度25,只对分析的第一个子步有效   <BR>solve <BR>EPLOT <BR>! <BR>! load step2………,apply moving heat flux <BR>allsel <BR>nsel,all <BR>t=0 <BR>dt=0.5 <BR>local,11,1, 0.001+0.005,0,0,0,-90,0        <BR>CSYS,11 <BR>esel,s,type,,1 <BR>nsel,r,loc,x, 0,0.003 <BR>nsel,r,loc,y,0,180 <BR>nsel,r,loc,z,0 <BR>sf,,hflux,7000000 <BR>!sf,,hflux,%flux2%     <BR>csys,0  <BR>t=t+dt <BR>time,t   !每一载荷步结束时的时间 <BR>deltim,0.1,0.05,0.1,on    <BR>kbc,1        !设定为阶越载荷 <BR>!neqit,25       !每个子步的默认次数 <BR>autots,on        !打开自动时间步长 <BR>timint,on,therm    !时间积分效果 <BR>TINTP,,,,1     <BR>Lnsrch,on <BR>!outres,all,all <BR>outres,all,last    <BR>!outres,all,10 <BR>allsel <BR>solve <BR>save
发表于 2006-5-2 07:41 | 显示全部楼层

回复:(wz1105)新手急求老师指点

在帮助文件里,search "VM38 Input Listing",里面有两个例子,是surf153.surf154的,你也可以参考一下
发表于 2006-5-2 09:39 | 显示全部楼层
哈哈,我也是面瓜,有好多问题想请教呢~~~~~~~~~~~~~~~~~~~~~
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 05:46 , Processed in 0.052962 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表