ANSYS循环搜索命令流文档格式.docx

ANSYS循环搜索命令流文档格式.docx

《ANSYS循环搜索命令流文档格式.docx》由会员分享，可在线阅读，更多相关《ANSYS循环搜索命令流文档格式.docx（16页珍藏版）》请在冰豆网上搜索。

*dim,xlx,,36

xlx

（1）=74.25,72.43,70.61,68.79,66.97,65.15,62.905,60.66,58.415,56.17

xlx（11）=53.925,51.71,49.495,47.28,45.065,42.85,40.665,38.48,36.295,34.11

xlx（21）=31.925,29.77,27.615,25.46,23.305,21.15,19.025,16.9,14.775,12.65

xlx（31）=10.525,8.42,6.315,4.21,2.105,0

*do,i,1,36,1

angle=atan（0.5115*sinh（1.475*xlx（i）/74.25））

h=（（2.5**3）/（（1-（1-0.5225）*xlx（i）/74.25）*cos（angle）））**（1/3）

z1=-xlx（i）-（h/2）*sin（angle）

y1=-（32.474/1.30）*（cosh（1.475*（-xlx（i））/74.25）-1）+（h/2）*cos（angle）

z2=-xlx（i）+（h/2）*sin（angle）

y2=-（32.474/1.30）*（cosh（1.475*（-xlx（i））/74.25）-1）-（h/2）*cos（angle）

k,i,0,y1,z1

k,36+i,0,y2,z2

*do,i,1,31,5

a,i,i+1,i+2,i+3,i+4,i+5,i+41,i+40,i+39,i+38,i+37,i+36

arsym,z,1,7,1,100,0,0

allsel

nummrg,all

k,1001,0,0,100

k,1002,23.9,0,100

l,1001,1002

vdrag,1,2,3,4,5,6,151

vdrag,8,9,10,11,12,13,151

vdrag,7,14,,,,,151

ldel,151,,,1

/view,1,-1,1,-1

vplo

腹拱墩

*dim,zgd,,7

zgd

（1）=65.154,53.929,42.854,31.929,21.154,10.529,0

*dim,kd,,7

kd

（1）=1.9,1.75,1.6,1.45,1.3,1.15,1

*do,i,1,7,1

z=-zgd（i）

k,2000+i,0.1,2.25,z-kd（i）/2

k,2010+i,0.1,2.25,z+kd（i）/2

k,2020+i,23.9,2.25,z-kd（i）/2

k,2030+i,23.9,2.25,z+kd（i）/2

a,2000+i,2010+i,2030+i,2020+i

k,4001,0,-23.123,0

k,4002,0,0,0

l,4001,4002

vdrag,172,,,,,,481

ldel,481,,,1

k,4001,0,-15.296,0

vdrag,174,,,,,,481

k,4001,0,-9.365,0

vdrag,186,,,,,,481

k,4001,0,-5.395,0

vdrag,187,,,,,,481

k,4001,0,-2.735,0

vdrag,188,,,,,,481

k,4001,0,-1.009,0

vdrag,189,,,,,,481

k,4001,0,-0.7,0

vdrag,190,,,,,,481

gplot

挖空1

*dim,zgdx,,4

zgdx

（1）=65.154,53.929,42.854,31.929

*dim,kdx,,4

kdx

（1）=1.9,1.75,1.6,1.45

*dim,dkg,,4

dkg

（1）=19.323,11.396,5.465,1.195

*do,i,1,4,1

z=-zgdx（i）

k,2040+i,3.26,-0.15,z-kdx（i）/2

k,2050+i,3.26,-0.15,z+kdx（i）/2

k,2060+i,6.86,-0.15,z-kdx（i）/2

k,2070+i,6.86,-0.15,z+kdx（i）/2

k,2080+i,3.26,-1.05-dkg（i）,z-kdx（i）/2

k,2090+i,3.26,-1.05-dkg（i）,z+kdx（i）/2

k,2100+i,6.86,-1.05-dkg（i）,z-kdx（i）/2

k,2110+i,6.86,-1.05-dkg（i）,z+kdx（i）/2

a,2040+i,2050+i,2090+i,2080+i

a,2060+i,2070+i,2110+i,2100+i

a,2080+i,2090+i,2110+i,2100+i

a,2040+i,2050+i,2070+i,2060+i

ASEL,,,,226,229,1

VSBA,15,all

ASEL,,,,230,233,1

VSBA,16,all

ASEL,,,,234,237,1

VSBA,17,all

ASEL,,,,238,241,1

VSBA,18,all

VDELE,15

VDELE,16

VDELE,17

VDELE,22

挖空2

k,2140+i,10.06,-0.15,z-kdx（i）/2

k,2150+i,10.06,-0.15,z+kdx（i）/2

k,2160+i,13.66,-0.15,z-kdx（i）/2

k,2170+i,13.66,-0.15,z+kdx（i）/2

k,2180+i,10.06,-1.05-dkg（i）,z-kdx（i）/2

k,2190+i,10.06,-1.05-dkg（i）,z+kdx（i）/2

k,2200+i,13.66,-1.05-dkg（i）,z-kdx（i）/2

k,2210+i,13.66,-1.05-dkg（i）,z+kdx（i）/2

a,2140+i,2150+i,2190+i,2180+i

a,2160+i,2170+i,2210+i,2200+i

a,2180+i,2190+i,2210+i,2200+i

a,2140+i,2150+i,2170+i,2160+i

ASEL,s,,,207,209,2

ASEL,a,,,252,253,1

VSBA,23,all

ASEL,s,,,254,257,1

VSBA,24,all

ASEL,,,,258,261,1

VSBA,25,all

ASEL,,,,262,265,1

VSBA,26,all

VDELE,24

挖空3

k,2240+i,16.86,-0.15,z-kdx（i）/2

k,2250+i,16.86,-0.15,z+kdx（i）/2

k,2260+i,20.46,-0.15,z-kdx（i）/2

k,2270+i,20.46,-0.15,z+kdx（i）/2

k,2280+i,16.86,-1.05-dkg（i）,z-kdx（i）/2

k,2290+i,16.86,-1.05-dkg（i）,z+kdx（i）/2

k,2300+i,20.46,-1.05-dkg（i）,z-kdx（i）/2

k,2310+i,20.46,-1.05-dkg（i）,z+kdx（i）/2

a,2240+i,2250+i,2290+i,2280+i

a,2260+i,2270+i,2310+i,2300+i

a,2280+i,2290+i,2310+i,2300+i

a,2240+i,2250+i,2270+i,2260+i

ASEL,s,,,250,251,1

ASEL,a,,,276,277,1

ASEL,s,,,278,281,1

ASEL,,,,282,285,1

ASEL,,,,286,289,1

VDELE,18

VDELE,23

numcmp,all

vsymm,z,15,17,1

vsymm,z,19,21,1

wpoff,0.1,-0.15,0

wprota,,90,,

vsbw,15

vsbw,19

vsbw,20

vsbw,21

vsbw,22

vsbw,25

vsbw,26

vsbw,27

wpcsys,-1

wpoff,0.1,-1.05-19.323,0

vsbw,33

vsbw,29

wpoff,0.1,-1.05-11.396,0

vsbw,34

vsbw,30

wpoff,0.1,-1.05-5.465,0

vsbw,35

vsbw,31

wpoff,0.1,-1.05-1.195,0

vsbw,36

vsbw,32

wpoff,3.26,0,0

wprota,,,90,

vsbw,40

vsbw,48

vsbw,56

vsbw,64

vsbw,68

vsbw,60

vsbw,52

vsbw,28

vsbw,44

wpoff,6.86,0,0

vsbw,69

vsbw,70

vsbw,71

vsbw,72

vsbw,73

vsbw,74

vsbw,75

vsbw,76

vsbw,77

vsbw,78

vsbw,79

vsbw,80

vsbw,81

vsbw,82

vsbw,83

vsbw,84

wpoff,10.06,0,0

vsbw,85

vsbw,86

vsbw,87

vsbw,88

vsbw,89

vsbw,90

vsbw,91

vsbw,92

vsbw,93

vsbw,94

vsbw,95

vsbw,96

vsbw,97

vsbw,98

vsbw,99

vsbw,100

wpoff,13.66,0,0

vsbw,101

vsbw,102

vsbw,103

vsbw,104

vsbw,105

vsbw,106

vsbw,107

vsbw,108

vsbw,109

vsbw,110

vsbw,111

vsbw,112

vsbw,113

vsbw,114

vsbw,115

vsbw,116

wpoff,16.86,0,0

vsbw,117

vsbw,118

vsbw,119

vsbw,120

vsbw,121

vsbw,122

vsbw,123

vsbw,124

vsbw,125

vsbw,126

vsbw,127

vsbw,128

vsbw,129

vsbw,130

vsbw,131

vsbw,132

wpoff,20.46,0,0

vsbw,133

vsbw,134

vsbw,135

vsbw,136

vsbw,137

vsbw,138

vsbw,139

vsbw,140

vsbw,141

vsbw,142

vsbw,143

vsbw,144

vsbw,145

vsbw,146

vsbw,147

vsbw,148

/eshape,1

vplot

腹拱

wpoff,0.1,-0.519,70.8

cyl4,,,5.455,27,6.055,153,23.7

vgen,2,164,,,,,-11.3

vgen,2,165,,,,,-11.15

vgen,2,166,,,,,-11

vgen,2,167,,,,,-10.85

vgen,2,168,,,,,-10.7

vgen,2,169,,,,,-10.55

vsymm,z,164,170,1

/eshape

墙座

v,468,466,439,440,212,226,107,106

v,516,514,487,488,222,236,112,111

v,548,546,519,520,232,246,117,116

v,580,578,551,552,242,256,122,121

v,476,474,471,472,252,266,127,126

v,484,482,479,480,262,288,132,131

v,311,312,309,310,285,271,35,135

v,307,308,305,306,268,199,31,32

v,303,304,301,302,195,189,26,27

v,299,300,297,298,185,179,21,22

v,295,296,293,294,175,99,16,17

v,291,292,289,290,95,89,11,12

v,279,280,267,278,85,78,6,7

wpoff,0,2.25,0

wprota,,-90,,

vsbw,all

vdelet,191,218,1,1

赋属性

et,1,plane42

et,2,solid65

1#材料（主拱）

mp,ex,2,1.5e10

mp,prxy,2,0.255

mp,dens,2,2780

tb,mkin,2

tbtemp,,strain

tbdata,,10e-4,15e-4,17.5e-4,20e-4,30e-4

tbtemp,,

tbdata,,15e6,18.75e6,19.68e6,20e6,17e6

real,1

tb,concr,2

tbdata,,0.20,0.8,1e6,20e6

asel,s,,,166

aplo

lesize,73,,,5

lesize,426,,,5

lesize,425,,,8

lesize,427,,,8

mshkey,1

amesh,166

vsel,s,,,1,14,1

type,2

mat,2

vsweep,all

aclear,all

2#材料（腹拱与主墩）

MPTEMP,,,,,,,, 

MPTEMP,1,0 

MPDATA,EX,3,,1.45e10

MPDATA,PRXY,3,,0.255 

MPDATA,DENS,3,,2750

TB,CONC,3,1,9, 

TBTEMP,0

tbdata,,0.20,0.8,1e6,20e6

TBDATA,,,,,,, 

vsel,,,,164,176,1

vatt,3,,2,

mshape,0,3d

vmesh,164,176,1

vsel,,,,15,163,1

vmesh,15,163,1

vsel,,,,177,190,1

vmesh,177,190,1

耦合1

asel,s,,,416 

选择第一个墩上部的所有节点，并编组

asel,a,,,625,881,64

asel,a,,,815

asel,a,,,463

asel,a,,,633,889,64

asel,a,,,823

asel,a,,,491

asel,a,,,641,897,64

asel,a,,,831

asel,a,,,519

asel,a,,,649,905,64

asel,a,,,839

asel,a,,,320

asel,a,,,326

asel,a,,,190

asel,a,,,189

asel,a,,,188

asel,a,,,381

asel,a,,,657,913,64

asel,a,,,846

asel,a,,,210

asel,a,,,665,921,64

asel,a,,,854

asel,a,,,205

asel,a,,,673,929,64

asel,a,,,862

asel,a,,,200

asel,a,,,681,934,64

asel,a,,,935

nsla

选择墩上部节点

取出墩顶节点总数

定义墩顶节点编号数组,节点坐标数组

选择墩顶附近的节点（去除墩顶节点）,找到墩顶的最近节点

nsel,u,cp,,all

/plopts,date,0

/triad,ltop

eplo

/solu

asel,s,,,26

asel,a,,,88

asel,a,,,1004

asel,a,,,1018

nsla,,1

d,all,all

antype,static

nropt,full,,1

autots,on

pred,on

neqit,50

kbc,0

cnvtol,f,100000,0.001,2

cnvtol,u,10,0.001,0

time,1

nsubst,50

acel,,9.81

outpr,all,all

outres,all,all

solve

可以采用这种办法，将每个腹拱上需要耦合的节点定义为一个选