用普通单元来模拟复合材料

zhaoqite

通过对不同体(面)中的单元赋以不同的单元或不同的单元常数,来模拟一些各向同性,只是强度刚度不同的复合材料.

   (下边是命令) 还望指正!


*set,alphaex,120000                                                
*set,alphaprxy,0.3
*set,alphayield,600
*set,alphatangmode,10000
*set,betaex,80000
*set,betaprxy,0.3
*set,betayield,900
*set,betatangmode,8000
*dim,rectanglearea,array,2,1,1   
rectanglearea(1)=7.8,7.794
*set,angle1,-30
*set,angle2,30
*set,aspectratio,1.8
*set,radio,0.5
*set,rightmove,1.3
*set,up,1.5588
*set,deform,-1.5588

/prep7
ET,1,VISCO108                                                         
KEYOPT,1,3,1
MPTEMP,1,0  
MPDATA,EX,1,,alphaex
MPDATA,PRXY,1,,alphaprxy  
TB,BISO,1,1,2,
TBTEMP,0
TBDATA,,alphayield,alphatangmode,,,,
MPTEMP,1,0   
MPDATA,EX,2,,betaex  
MPDATA,PRXY,2,,betaprxy  
TB,BISO,2,1,2,  
TBTEMP,0
TBDATA,,betayield,betatangmode,,,,   
wpro,-angle1,,   
cswpla,11,1,aspectratio,1
/prep7  
k,1,-radio
k,2,radio
l,1,2   
csys,4  
lsymm,y,1,,,,0,0
LGLUE,all  
AL,all
WPCSYS,-1,0
AGEN,2,1, , ,rightmove, , , ,0
AGEN,2,2, , ,rightmove, , , ,0
AGEN,2,3, , ,rightmove, , , ,0
AGEN,2,4, , ,rightmove, , , ,0
AGEN,2,5, , ,rightmove, , , ,0
AGEN,2,6, , ,rightmove, , , ,0
AGEN,2,7, , ,rightmove, , , ,0
wpro,-angle2,,   
cswpla,11,1,aspectratio,1
/prep7  
k,17,-radio
k,18,radio
l,17,18  
csys,4  
lsymm,y,17,,,,0,0
LGLUE,17,18
AL,17,19

WPCSYS,-1,0
AGEN,0,9, , ,-0.65,up, ,,,1
AGEN,2,9, , ,rightmove, , , ,0
AGEN,2,10, , ,rightmove, , , ,0
AGEN,2,11, , ,rightmove, , , ,0
AGEN,2,12, , ,rightmove, , , ,0
AGEN,2,13, , ,rightmove, , , ,0
AGEN,2,14, , ,rightmove, , , ,0
AGEN,2,15, , ,rightmove, , , ,0
AGEN,2,1,16 , ,,2*up , , ,0
AGEN,2,1,16 , ,,4*up , , ,0
BLC4,0,0,rectanglearea(1),rectanglearea(2)

asbl,1,all
asbl,2,all
asbl,3,all
asbl,4,all
asbl,5,all
asbl,6,all
asbl,7,all
asbl,17,all
asbl,33,all
asbl,42,all
asbl,43,all
asbl,44,all
asbl,45,all
asbl,46,all
asbl,47,all
asbl,39,all
asbl,23,all
adele,1,5,,1
adele,7,9,,1
adele,40,46,,1
adele,50,,,1
adele,57,,,1
adele,25,,,1
adele,17,,,1
adele,33,,,1
adele,48,,,1
adele,66,,,1
adele,32,,,1
adele,67,,,1
adele,24,,,1
adele,16,,,1
aovlap,all
ASEL,all
AATT,1
ASEL,s,,,1
AATT,2
ASEL,all
AESIZE,ALL,0.1
ASEL,all
AMESH,all
FINISH
/SOL
lsel,,line,,3
lsel,a,line,,129
lsel,a,line,,15
lsel,a,line,,126
lsel,a,line,,5
lsel,a,line,,104
DL,all, ,SYMM
lsel,,line,,103
lsel,a,line,,6
lsel,a,line,,107
lsel,a,line,,7
lsel,a,line,,110
lsel,a,line,,8
lsel,a,line,,113
lsel,a,line,,9
lsel,a,line,,116
lsel,a,line,,10
lsel,a,line,,119
lsel,a,line,,13
lsel,a,line,,11
DL,all, ,UY
ksel,,kp,,97
DK,all, , , ,0,ALL
lsel,,line,,4
lsel,a,line,,132
lsel,a,line,,16
lsel,a,line,,135
lsel,a,line,,17
lsel,a,line,,138
lsel,a,line,,19
lsel,a,line,,141
lsel,a,line,,31
lsel,a,line,,144
lsel,a,line,,32
lsel,a,line,,147
lsel,a,line,,1
DL,all, ,UY,deform
allsel,all
nlgeom,on   
TIME,1  
AUTOTS,-1   
NSUBST,20,20,20,1   
KBC,0   
TSRES,ERASE
/STATUS,SOLU

[ 本帖最后由 zhaoqite 于 2008-11-7 08:34 编辑 ]

zhoulon

你这个问题连个模型都没有上传，只是自己说了几句话，别人如何来指正啊。

zhaoqite

已经把命令写出来了

sh_lin30

来模拟一些各向同性,只是强度刚度不同的复合材料

lz的意思不太明白
既然是复合材料，怎么可能是各向同性

可以细致阐明一下，你的想法？

zhaoqite

一些弥散强化材料,不就是各向同性的吗?