资源描述：
《黏弹性人工边界算例》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、!研究问题模型：验证Deeks和刘晶波教授提出的黏弹性人工边界精度（吸能效果）【平直边界】!二维线弹性全空间内部作用瞬态正弦荷载.边界物理元件的弹簧系数：Kt=G/2r,Kn=E/2rFINISH/CLEAR/TITLE,InnerFountainQuestion/PREP7!定义单元类型ET,1,PLANE42!定义波动传输的介质单元KEYOPT,1,3,2!平面应变ET,2,COMBIN14!定义法向2—D纵向弹簧阻尼单元KEYOPT,2,3,2!设定土体被假定为线弹性材料时的材料参数density=2000!密度poisson=0.25!
2、泊松比Cs=500!剪切波速G=density*Cs**2!剪切模量E=2*G*(1+poisson)!弹性模量K=E/(3.*(1.-2*poisson))!体积模量Cp=sqrt((K+4*G/3)/density)!纵波波速!定义土体材料,假定为线性材料mp,ex,1,E!土体弹性模量mp,dens,1,density!土体密度mp,prxy,1,0.25!土体的主泊松比!以下定义一些与边界元件有关的参数alfan=0.5!法向边界参数alfat=0.5!切向边界参数【可能用到】R=100!先【不考虑R值的影响】as=10!以下建立几何模
3、型!定义关键点k,1,100,0,0k,2,100,-100,0k,3,-100,-100,0k,4,-100,0,0!建立面a,1,2,3,4!以下划分网格FLST,5,2,4,ORDE,2FITEM,5,1FITEM,5,3CM,_Y,LINELSEL,,,,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1,,,10,,,,,1FLST,5,2,4,ORDE,2FITEM,5,2FITEM,5,4CM,_Y,LINELSEL,,,,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1,,,20,,,,,1
4、MSHAPE,0,2DMSHKEY,0CM,_Y,AREAASEL,,,,1CM,_Y1,AREACHKMSH,'AREA'CMSEL,S,_YAMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2!底边界上法向和切向的弹簧阻尼单元!检索模型底边界上的所有节点numnode1=100!定义一个参数即以下一维数组的上界*dim,node1,,numnode1!定义一个一维数组nsel,s,loc,y,-100!选择一组节点子集满足y坐标为-100*get,np,node,,count!得到选中的结点数,存入np*get,no
5、de1(1),node,,NUM,MIN!提取选中的节点编号的最小值*DO,k,2,np!此循环是对数组中的每一元素赋上节点编号,赋值方法有意思*GET,node1(k),node,node1(k-1),NXTH*ENDDO!正式施加边界元件*get,npmax,node,,num,maxd!得到已经定义的最大结点编号,存入npmax,对所有节点而言*do,i,1,npx=nx(node1(i))y=ny(node1(i))z=nz(node1(i))npnum=node(x,y,z)!节点编号npmax=npmax+1!新增节点编号n,npm
6、ax,x,y-4,z!定义底边界外法向的结点,以便与土体边界点形成法向单元,先假定单元长度为4*if,abs(x),eq,100,then!底边界上的角节点KbN=alfan*G/(R)*as*0.5!法向弹簧系数CbN=density*Cp*as*0.5!法相阻尼系数r,2,KbN,CbNnreal=2*else!底边界上的其他节点KbN=alfan*G/(R)*asCbN=density*Cp*asr,3,KbN,CbNnreal=3*endiftype,2!设定成编号为2的单元real,nreale,npnum,npmax!创建单元d,n
7、pmax,all!约束新生成的节点npmax=npmax+1!新增节点编号n,npmax,x-4,y,z!定义底边界外切向x向的结点,以便与土体边界节点形成切向单元*if,abs(x),eq,100,thenKbT=alfat*G/(R)*as*0.5CbT=density*Cs*as*0.5r,4,KbT,CbTtreal=4*elseKbT=alfat*G/(R)*asCbT=density*Cs*asr,5,KbT,CbTtreal=5*endiftype,2real,treale,npnum,npmaxd,npmax,all!约束新生成
8、的点*enddo!以下建立左边界上法向和切向的弹簧阻尼单元!先检索模型左边界上的所有节点numnode3=100!定义一个参数即以下一维数组的上界*d
