exa_hydfracture10.1.5 Hydraulically induced fracture in a well bore 液压诱导井眼破裂
units s, kN and KPa and m
Solid Section, elset=lower_shale, material=Mat_lower_shale
**
*Material, name=Mat_lower_shale下层页岩
*Elastic
0.625e+07, 0.19
*DruckerPrager
30., 1.0, 30.
*DruckerPrager Hardening
32.00e+03 , 0.00
33.000e+03, 0.005
33.500e+03, 0.08
*Permeability, specific=9.8
7.6720e-09, 0.17
7.6720e-09, 0.19
8.12200e-09, 0.24
9.0500e-09, 0.30
1.8800e-08, 0.36
**
*Solid Section, elset=upper_shale, material=Mat_upper_shale
**
*Material, name=Mat_upper_shale上层页岩
*Elastic
0.4137e+07, 0.15
*DruckerPrager
28., 1.0, 28.
*DruckerPrager Hardening
26.207e+03 , 0.00
26.500e+03, 0.005
27.00e+03, 0.08
*Permeability, specific=9.8
9.6720e-09, 0.17
9.6720e-09, 0.19
1.0200e-08, 0.24
2.0500e-08, 0.28
4.8800e-08, 0.32
**
*Solid Section, elset=target_formation, material=Mat_target_formation
**
*Material, name=Mat_target_formation目标地层
*Elastic
1.2e+07, 0.22
*DruckerPrager
36., 0.95, 36.
*DruckerPrager Hardening
37.930e+3, 0.
38.20e+3, 0.005
38.34e+3, 0.02
38.40e+3, 0.08
*Permeability, specific=9.8？？？？？other》pore fluid井眼流动》permeability表示渗透性。。。。。
0.560e-08, 0.2
1.000e-08, 0.22
4.836e-08, 0.25
4.880e-08, 0.29
4.980e-08, 0.31
5.280e-08, 0.33
5.800e-08, 0.35
6.525e-08, 0.40
8.100e-08, 0.48
*Membrane Section, elset=casing, material=casing
0.01,
*Material, name=casing
*Elastic
2.e+06, 0.3
**
*Initial Conditions, type=Ratio,user
*Initial Conditions, type=Pore pressure, user
*Initial Conditions, type=Stress, user
**
*elset,elset=initopen,generate
5183,5186
*INITIAL CONDITIONS,TYPE=INITIAL GAP
initopen
*nset,nset=iniflm,generate
1000076,1000078
*nset,nset=inifle 1000075,1000079
*nset,nset=inifl iniflm,inifle
*Nset, nset=target_zone
exa_hydfracture10.1.5 Hydraulically induced fracture in a well bore
*Elset, elset=__PickedSurf260_S4, generate
5383, 5440, 3
*Surface, type=ELEMENT, name=_PickedSurf260
__PickedSurf260_S5, S5
**
*COHESIVE SECTION,ELSET=upper_inter,MATERIAL=upper_inter_mat,RESPONSE=TRACTION SEPARATION,
THICKNESS=SPECIFIED,CONTROLS=VISCO_upper
*SECTION CONTROLS,NAME=VISCO_upper,VISCOSITY=0.01
*MATERIAL,NAME=upper_inter_mat
*ELASTIC,TYPE=TRACTION
*DAMAGE INITIATION,CRITERION=QUADS
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=BK,POWER=
*GAP FLOW
*FLUID LEAKOFF, USER
**
*COHESIVE SECTION,ELSET=lower_inter,MATERIAL=lower_inter_mat,RESPONSE=TRACTION SEPARATION,
THICKNESS=SPECIFIED,CONTROLS=VISCO_lower
*SECTION CONTROLS,NAME=VISCO_lower,VISCOSITY=0.01
*MATERIAL,NAME=lower_inter_mat
*ELASTIC,TYPE=TRACTION
*DAMAGE INITIATION,CRITERION=QUADS
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=BK,POWER=
*GAP FLOW
*FLUID LEAKOFF, USER
**
*COHESIVE SECTION,ELSET=middle_inter,MATERIAL=middle_inter_mat,RESPONSE=TRACTION SEPARATION,
THICKNESS=SPECIFIED,CONTROLS=VISCO_middle
*SECTION CONTROLS,NAME=VISCO_middle,VISCOSITY=0.01
*MATERIAL,NAME=middle_inter_mat
*ELASTIC,TYPE=TRACTION
*DAMAGE INITIATION,CRITERION=QUADS
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=BK,POWER=
*GAP FLOW
*FLUID LEAKOFF, USER
**
*Solid Section, elset=lower_shale, material=Mat_lower_shale
**
*Material, name=Mat_lower_shale
*Elastic
0.625e+07, 0.19
*DruckerPrager
30., 1.0, 30.
*DruckerPrager Hardening
32.00e+03 , 0.00
33.000e+03, 0.005
33.500e+03, 0.08
*Permeability, specific=9.8
7.6720e-09, 0.17
7.6720e-09, 0.19
8.12200e-09, 0.24
9.0500e-09, 0.30
1.8800e-08, 0.36
**
*Solid Section, elset=upper_shale, material=Mat_upper_shale
**
*Material, name=Mat_upper_shale
*Elastic
0.4137e+07, 0.15
*DruckerPrager
28., 1.0, 28.
*DruckerPrager Hardening
26.207e+03 , 0.00
26.500e+03, 0.005
27.00e+03, 0.08
*Permeability, specific=9.8
9.6720e-09, 0.17
9.6720e-09, 0.19
1.0200e-08, 0.24
2.0500e-08, 0.28
4.8800e-08, 0.32
**
*Solid Section, elset=target_formation, material=Mat_target_formation **
*Material, name=Mat_target_formation
*Elastic
1.2e+07, 0.22
*DruckerPrager
36., 0.95, 36.
*DruckerPrager Hardening
37.930e+3, 0.
38.20e+3, 0.005
38.34e+3, 0.02
38.40e+3, 0.08
*Permeability, specific=9.8
0.560e-08, 0.2
1.000e-08, 0.22
4.836e-08, 0.25
4.880e-08, 0.29
4.980e-08, 0.31
5.280e-08, 0.33
5.800e-08, 0.35
6.525e-08, 0.40
8.100e-08, 0.48
*Membrane Section, elset=casing, material=casing
0.01,
*Material, name=casing
*Elastic
2.e+06, 0.3
**
*Initial Conditions, type=Ratio,user
*Initial Conditions, type=Pore pressure, user
*Initial Conditions, type=Stress, user
**
*elset,elset=initopen,generate
5183,5186
*INITIAL CONDITIONS,TYPE=INITIAL GAP
initopen
*nset,nset=iniflm,generate
1000076,1000078
*nset,nset=inifle
1000075,1000079
*nset,nset=inifl
iniflm,inifle
*Nset, nset=target_zone
5, 6, 7, 8, 9, 10, 25, 28, 65, 66, 67, 68, 69, 70, 71, 72
73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88
89, 90, 91, 92, 93, 94, 95, 96, 173, 174, 175, 176, 177, 178, 179, 363
364, 365, 366, 367, 368, 369, 457, 458, 459, 460, 461, 462, 463, 464, 465, 466
467, 468, 469, 470, 471, 472, 473, 474, 616, 617, 618, 619, 620, 621, 622, 623
624, 625, 626, 627, 628, 629, 630, 631, 632, 633, 634, 635, 636, 637, 638, 639
640, 641, 642, 643, 644, 645, 646, 647, 648, 649, 650, 651, 652, 653, 654, 655
656, 657, 658, 659, 660, 661, 662, 663, 664, 665, 666, 667, 668, 669, 670, 671
672, 673, 674, 675, 676, 677, 678, 2264, 2265, 2266, 2267, 2268, 2269, 2270, 2271, 2272
2273, 2274, 2275, 2276, 2277, 2278, 2279, 2280, 2281, 2282, 2283, 2284, 2285, 2286, 2287, 2288
2289, 2290, 2291, 2292, 2293, 2294, 2295, 2296, 2297, 2298, 2299, 2300, 2301, 2302, 2303, 2304
2305, 2306, 2307, 2308, 2309, 2310, 2311, 2312, 2313, 2314, 2315, 2316, 2317, 2318, 2319, 2320
2321, 2322, 2323, 2324, 2325, 2326
iniflm,
inifle,
*elset,elset=one
5184
*elset,elset=QA
3113,3105,3097,3089
5185,5186,5187,5188
*nset,nset=QA
1001015,1001016,1001017,1001018
**
** BOUNDARY CONDITIONS
**
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
_PickedSet4, XSYMM
** Name: BC-2 Type: Displacement/Rotation
*Boundary
_PickedSet5, 1, 1
_PickedSet5, 2, 2
** Name: BC-3 Type: Displacement/Rotation
*Boundary
_PickedSet6, 3, 3
**
*Amplitude, name=volumerate
0.0,0.0, 200.0,-1.0
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1, nlgeom=YES, unsymm=YES
*Geostatic
**
** LOADS
**
** Name: Load-1 Type: Pressure
*Dsload
_PickedSurf260, P, 42000.
** Name: Load-2 Type: Body force
*Dload
_PickedSet8, BZ, -20.
** Name: Load-3 Type: Pressure
*Dsload
well_bore, PNU, 1.
**
*Boundary,user
**_PickedSet5, 8, 8, 1
TOP, 8, 8, 1
BOT, 8, 8, 1
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
** MAKE CHANGES HERE
** use the variable=PRESELECT version of the *output keyword line
** *Output, field, variable=PRESELECT
*Output, field
*node output
POR,
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: Pump Stage
**
*Step, name="Pump Stage", nlgeom=YES, inc=200, unsymm=YES *Soils, consolidation, end=PERIOD, utol=500000.
** MAKE CHANGES HERE
** Use the 1200s data line in place of the 50s dataline
**1., 50., 0.1, 30.,
1., 1200., 0.01, 30.,
**
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary,op=new
_PickedSet4, XSYMM
** Name: BC-2 Type: Displacement/Rotation
*Boundary,op=new
_PickedSet5, 1, 1
_PickedSet5, 2, 2
** Name: BC-3 Type: Displacement/Rotation
*Boundary,op=new
_PickedSet6, 3, 3
**
*Boundary,user,op=new
TOP, 8, 8, 1
BOT, 8, 8, 1
**Controls, PARAM=FIELD, FIELD=PORE
**, , 1.0e-3
**
*cflow,amplitude=volumerate
iniflm,,
inifle,,
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
** MAKE CHANGES HERE
** Remove the nset=QA, and elset=QA parameter settings **Node Output,nset=QA
*Node Output
POR,U
**element output,elset=QA
*element output
S,LE,PE
GFVR
SDEG
DMICRT
PFOPEN
LEAKVRT
LEAKVRB
ALEAKVRT
ALEAKVRB
FLVEL
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*node output,nset=inifl
por
*element output,elset=one
LEAKVRT
LEAKVRB
ALEAKVRT
ALEAKVRB
SDEG
*node print,nset=inifl
por
*End Step
** ----------------------------------------------------------------
**
** STEP: Hold
**
*Step, name=Hold, nlgeom=YES, unsymm=YES
hold
** MAKE CHANGES HERE
** Use the 7200s data line in place of the 100s dataline
**Soils, consolidation, end=SS, utol=1.E6
**100., 100., 0.005, 1500., 10.0
*Soils, consolidation, end=SS, utol=1.0e+6
100., 7200., 0.005, 1500., 10.0
**
** BOUNDARY CONDITIONS
**
** Name: BC-4 Type: Velocity/Angular velocity
*Boundary, type=VELOCITY
_PickedSet273_1, 1
_PickedSet273_2, 2
_PickedSet273_3, 3, 3
**
*cflow
inifl,,0.0
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=1
**
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
** MAKE CHANGES HERE
** Remove the nset=QA, and elset=QA parameter settings **Node Output,nset=QA
*Node Output
POR,U
**element output,elset=QA
*element output
S,LE,PE
GFVR
SDEG
DMICRT
PFOPEN
LEAKVRT
LEAKVRB
ALEAKVRT
ALEAKVRB
FLVEL
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*node output,nset=inifl
por
*element output,elset=one
LEAKVRB
ALEAKVRT
ALEAKVRB
SDEG
*node print,nset=inifl
por
*End Step
** ----------------------------------------------------------------
**
** STEP: Draw down
**
*Step, name="Draw down", nlgeom=YES, unsymm=YES drawdown
** MAKE CHANGES HERE
** Use the 864000s data line in place of the 200s dataline **Soils, consolidation, end=SS, utol=1.E6
**100., 200., 100., 200., 0.05
*Soils, consolidation, end=SS, utol=1.E6
100., 864000., 100., 43200., 0.05
**
** BOUNDARY CONDITIONS
**
** Name: BC-4 Type: Velocity/Angular velocity
*Boundary, type=VELOCITY,op=new
_PickedSet273_1, 1
_PickedSet273_2, 2
_PickedSet273_3, 3, 3
**
*Boundary,op=new
target_zone,8,8,20000.
**
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre *Boundary,op=new
_PickedSet4, XSYMM
** Name: BC-2 Type: Displacement/Rotation
*Boundary,op=new
_PickedSet5, 1, 1
_PickedSet5, 2, 2
** Name: BC-3 Type: Displacement/Rotation
*Boundary,op=new
_PickedSet6, 3, 3
**
*Boundary,user,op=new
BOT, 8, 8, 1
**
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
** MAKE CHANGES HERE
** Remove the nset=QA, and elset=QA parameter settings **Node Output,nset=QA
*Node Output
POR,U
**element output,elset=QA
*element output
S,LE,PE
GFVR
SDEG
DMICRT
PFOPEN
LEAKVRT
LEAKVRB
ALEAKVRT
ALEAKVRB
FLVEL
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*node output,nset=target_zone
por,rvf,rvt
*element output,elset=one
LEAKVRT
LEAKVRB
ALEAKVRT
ALEAKVRB
SDEG
*node print,nset=target_zone
*End Step
Planar model of the oil wellbore1.1.22 Erosion of material (sand production) in an oil wellbore
出沙油井井眼材料腐蚀
*Elset, elset=__PickedSurf11_S2, internal, instance=Formation-1, generate
200, 288, 8
*Elset, elset=__PickedSurf11_S1, internal, instance=Formation-1, generate
289, 296, 1
*Surface, type=ELEMENT, name=_PickedSurf11, internal
__PickedSurf11_S2, S2
__PickedSurf11_S1, S1
*Elset, elset=__PickedSurf12_S2, internal, instance=Formation-1, generate
296, 352, 8
*Elset, elset=__PickedSurf12_S1, internal, instance=Formation-1, generate
353, 364, 1
*Surface, type=ELEMENT, name=_PickedSurf12, internal
__PickedSurf12_S2, S2
__PickedSurf12_S1, S1
*Elset, elset=Formation_Face_P2, instance=Formation-1
898, 927, 956, 985, 1014, 1043, 1072, 1101, 1130, 1159, 1188, 1217, 1246, 1275, 1304, 1333
1362, 1391, 1420
*Elset, elset=Formation_Face_P4, instance=Formation-1
1421, 1450, 1479, 1508, 1537, 1566, 1595, 1624, 1653, 1682, 1711
*End Assembly
**
** MATERIALS
**
*Material, name=Rock
*DruckerPrager
40., 1.,35.
*DruckerPrager Hardening
0.500, 0.
0.72, 0.005
1.0, 0.05
1.05, 0.1
1.2, 0.2
*Elastic
** units: MPa and meter
1000., 0.22
*Permeability, specific=0.0076
2.34833e-08, 0.04
1.17341e-07, 0.07
3.36324e-07, 0.1
7.45578e-07, 0.13
1.42563e-06, 0.16
2.4774e-06, 0.19
4.02852e-06, 0.22
6.24154e-06, 0.25
9.3249e-06, 0.28
1.35478e-05, 0.31
1.92602e-05, 0.34
2.69213e-05, 0.37
3.71384e-05, 0.4
5.0723e-05, 0.43
6.87716e-05, 0.46
9.27846e-05, 0.49
0.000124845, 0.52
0.000167889, 0.55
0.000226136, 0.58
0.000305778, 0.61
0.000416122, 0.64
0.000571566, 0.67
0.000795141, 0.7
0.0011252, 0.73
0.00162879, 0.76
0.00243043, 0.79
0.0037799, 0.82
0.0062315, 0.85
0.0112, 0.88
**
** BOUNDARY CONDITIONS
**
** Name: Dir-1 Type: Symmetry/Antisymmetry/Encastre *Boundary
_PickedSet8, XSYMM
** Name: Dir-2 Type: Symmetry/Antisymmetry/Encastre *Boundary
_PickedSet9, YSYMM
** Name: Initial Pore Pressure Type: Pore pressure
*Boundary
"All Rock", 8, 8
*Initial Conditions, type=ratio
"All Rock", 0.4286
**"All Rock", 0.0001
** initial geostatic stress: sig_V 16.25 MPasig_H 24.7 MPasig_h 19.5 MPa *Initial Conditions, type=stress
"All Rock", -24.7E0, -19.5E0, -16.25E0,
** ----------------------------------------------------------------
**
** STEP: Geostatic
**
*Step, name=Geostatic, nlgeom
Equilibrium w. Initial Stress
*Geostatic
**
** LOADS
**
** Name: Sig_H Type: Pressure
*Dsload
_PickedSurf11, P, 24.7
** Name: sig_h_min Type: Pressure
*Dsload
_PickedSurf12, P, 19.5
**
** OUTPUT REQUESTS
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*Output, history
*El Print, freq=999999
*Node Print, freq=999999
*End Step
** ----------------------------------------------------------------
**
** STEP: Drilling
**
*Step, name=Drilling, nlgeom, amplitude=RAMP
Remove Well Bore Elements
*Soils, utol=1.0
0.001, 0.001, 1e-08, 0.001,
*Model Change, remove
"Well Bore",
**
** Drilling Mud Load. 500 psi underbalanced. Formation pressure 1500 psi. *DLOAD
Formation_face_p2, P2, 6.895
Formation_face_p4, P4, 6.895
**
** OUTPUT REQUESTS
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*Node Output
*Element Output
**
*End Step
** ----------------------------------------------------------------
**
** STEP: Underbalanced for 32 hours
**
*Step, name="Underbalanced for 32 hours", nlgeom, inc=200
500 psi underbalanced for 32 hours
*Soils, consolidation, end=PERIOD, utol=10., creep=none
5.e3, 55.e3, 2e-07, 10.e3,
**
** BOUNDARY CONDITIONS
**
** Name: Dir-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary, op=NEW
_PickedSet8, XSYMM
** Name: Dir-2 Type: Symmetry/Antisymmetry/Encastre
*Boundary, op=NEW
_PickedSet9, YSYMM
** Name: Far Field Pore Type: Pore pressure
*Boundary, op=NEW
_PickedSet10, 8, 8
** Name: Initial Pore Pressure Type: Pore pressure
*Boundary, op=NEW
** Name: underbalance drilling Type: Pore pressure
*Boundary, op=NEW
"Formation Face", 8, 8, -3.4475
*ADAPTIVE MESH, ELSET=amc_region, FREQUENCY=1, MESH SWEEP=5
*ADAPTIVE MESH CONSTRAINT, type=velocity,user
"Formation Face",1,,1.1e-2
**
** OUTPUT REQUESTS
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, freq=9999
*Node Output,nset=QA_TEST_FORMATION_FACE
U,
*Node Output
U, RF, CF, POR, RVT, RVF
*Element Output
S, LE, VOIDR, SAT, FLVEL, PE, PEEQ, PEMAG
**
*restart,write
*End Step
Three-dimensional model of the oil wellbore perforation tunnel1.1.22 Erosion of material (sand production) in an oil wellbore
** MATERIALS
** Units: inch, lbf, psi
**
**
*Material, name=Rock
*DruckerPrager
45., 0.8 ,40.
*DruckerPrager Hardening
900 , 0
1320 , 0.003
1600 , 0.006
1800 , 0.009
1900 , 0.012
2000 , 0.018
2080 , 0.025
2150 , 0.035
2220 , 0.05
**2490 , 0.2
*Elastic
1.32e6, 0.22
*Permeability, specific=0.67
150 ,0.271597
180 ,0.285892
225 ,0.300939
275 ,0.316778
325 ,0.33345
400 ,0.351
440 ,0.35802
490 ,0.36518
570 ,0.372484
1000 ,0.4
**
*Material, name="Low Stifness"
*Elastic
0.001, 0.2
**
**
** Porosity 26%
*Initial Conditions, type=Ratio
"All Rock" , 0.3513514
**
*Initial Conditions, type=Pore
"All Rock" , 5500.00
*Initial Conditions, type=Stress, user
** ---------------------------------------------------------------- **
** STEP: geostatic
**
*Step, name=geostatic, nlgeom
Establish geostatic equilibrium
*Geostatic
*Boundary
Face-1, 1, 2, 0.00
Face-2, 1, 2, 0.00
Face-Radial, 1, 1, 0.00
Face-Bottom, 3, 3, 0.00
"All Rock", 8, 8, 5500.
"Center Line", 1, 1, 0.000
**
*DSLOAD
TOP_Face , P, 13000.00
**
** OUTPUT REQUESTS
**
***Restart, write, overlay, frequency=1
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*ELEMENT OUTPUT
PE, PEEQ, S
*NODE OUTPUT
U,POR
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*El Print, freq=999999
*Node Print, freq=999999
*End Step
*********************************************************** ** STEP 2: Drilling Step
**
**
** Elset _I5 is casing ID
** _I6 is perforation face
** _I7 is well bore face
**
*********************************************************** *Step, name="Drilling", inc=6, nlgeom=yes
*SOILS, utol=300.0
0.3, 1.0, 1e-4,
***
*** Borehole and Perf removal
***
*Model Change, remove
"Well Bore", Perf_vol_Rock, Formation-1._I5
***
**
*Boundary
"Well Bore Face", 1, 1
*************************************************
*** LOADS
****
*** Borehole pressure applied to borehole face
*Dload
ANSYS命令流中文说明（2） 默认分类 2009-10-02 10:28 阅读106 评论0 字号：大大中中小小 KB、KE: 待划分线的定向关键点起始、终止号 SECNUM: 截面类型号 u SECPLOT,SECID,MESHKEY 画梁截面的几何形状及网格划分 SECID:由SECTYPE命令分配的截面编号 MESHKEY：0：不显示网格划分 1：显示网格划分 u /ESHAPE, SCALE 按看似固体化分的形式显示线、面单元 SCALE: 0:简单显示线、面单元 1：使用实常数显示单元形状 u esurf, xnode, tlab, shape 在已存在的选中单元的自由表面覆盖产生单元 xnode: 仅为产生surf151 或surf152单元时使用 tlab: 仅用来生成接触元或目标元 top 产生单元且法线方向与所覆盖的单元相同，仅对梁或壳有效，对实体单元无效 Bottom产生单元且法线方向与所覆盖的单元相反，仅对梁或壳有效，对实体单元无效Reverse 将已产生单元反向 Shape: 空与所覆盖单元形状相同 Tri 产生三角形表面的目标元 注意：选中的单元是由所选节点决定的，而不是选单元，如同将压力加在节点上而不是单元上 u Nummrg,label,toler, Gtoler,action,switch 合并相同位置的item label: 要合并的项目 node: 节点，Elem,单元，kp: 关键点（也合并线，面及点） mat: 材料，type: 单元类型，Real: 实常数 cp：耦合项，CE：约束项，CE: 约束方程，All：所有项 toler: 公差 Gtoler：实体公差 Action: sele 仅选择不合并 空合并 switch: 较低号还是较高号被保留（low, high） 注意：可以先选择一部分项目，再执行合并。如果多次发生合并命令，一定要先合并节点，再合并关键点。合并节点后，实体荷载不能转化到单元，此时可合并关键点解决问题。 u Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线 type: s 从全部线中选一组线 r 从当前选中线中选一组线 a 再选一部线附加给当前选中组 au none u(unselect) inve: 反向选择 item: line 线号 loc 坐标
A N S Y S命令流解释大 全 Document serial number【LGGKGB-LGG98YT-LGGT8CB-LGUT-
一、定义材料号及特性 mp,lab, mat, co, c1,…….c4 lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens） ex: 弹性模量 nuxy: 小泊松比 alpx: 热膨胀系数 reft: 参考温度 reft: 参考温度 prxy: 主泊松比 gxy: 剪切模量 mu: 摩擦系数 dens: 质量密度 mat: 材料编号（缺省为当前材料号） c 材料特性值，或材料之特性，温度曲线中的常数项 c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数二、定义DP材料： 首先要定义EX和泊松比：MP，EX，MAT，…… MP，NUXY，MAT，…… 定义DP材料单元表（这里不考虑温度）：TB，DP，MAT 进入单元表并编辑添加单元表：TBDATA，1，C TBDATA，2，ψ TBDATA，3，…… 如定义：EX=1E8，NUXY=，C=27，ψ=45的命令如下：
MP，EX，1，1E8 MP，NUXY，1， TB，DP，1 TBDATA，1，27 TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg 三、单元生死载荷步 !第一个载荷步 TIME,... !设定时间值（静力分析选项） NLGEOM,ON !打开大位移效果 NROPT,FULL !设定牛顿－拉夫森选项 ESTIF,... !设定非缺省缩减因子（可选） ESEL,... !选择在本载荷步中将不激活的单元 EKILL,... !不激活选择的单元 ESEL,S,LIVE !选择所有活动单元 NSLE,S !选择所有活动结点 NSEL,INVE !选择所有非活动结点（不与活动单 元相连的结点） D,ALL,ALL,0 !约束所有不活动的结点自由度（可 选） NSEL,ALL !选择所有结点 ESEL,ALL !选择所有单元
Finish(退出四大模块，回到BEGIN层) /clear (清空存，开始新的计算) 1．定义参数、数组，并赋值. 2．/prep7(进入前处理) 定义几何图形：关键点、线、面、体 定义几个所关心的节点，以备后处理时调用节点号。 设材料线弹性、非线性特性 设置单元类型及相应KEYOPT 设置实常数 设置网格划分，划分网格 根据需要耦合某些节点自由度 定义单元表 3．/solu 加边界条件 设置求解选项 定义载荷步 求解载荷步 4./post1（通用后处理） 5./post26 （时间历程后处理） 6.PLOTCONTROL菜单命令 7.参数化设计语言 8.理论手册 Finish(退出四大模块，回到BEGIN层) /clear (清空存，开始新的计算) 1.定义参数、数组，并赋值. dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组 par: 数组名 type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符） table imax,jmax, kmax 各维的最大下标号 var1,var2,var3 各维变量名，缺省为row,column,plane(当type为table时) 2./prep7(进入前处理) 2.1 设置单元类型及相应KEYOPT ET, itype, ename, kop1……kop6, inopr 设定当前单元类型 Itype：单元号 Ename：单元名设置实常数 Keyopt, itype, knum, value itype: 已定义的单元类型号 knum: 单元的关键字号
为方便大家的交流和学习,特推出"跟我学命令流"课程本课程分为三部分:前处理,加载求解,后处理 每部分的学习时间:10天,共计30天 每天学习大约10个命令 希望本课程对大家能有所帮助 第一天 目标:熟悉ANSYS基本关键字的含义 k --> Keypoints 关键点 l --> Lines 线 a --> Area 面 v --> Volumes 体 e --> Elements 单元 n --> Nodes 节点 cm --> component 组元 et --> element type 单元类型 mp --> material property 材料属性 r --> real constant 实常数 d --> DOF constraint 约束 f --> Force Load 集中力 sf --> Surface load on nodes 表面载荷 bf --> Body Force on Nodes 体载荷 ic --> Initial Conditions 初始条件 第二天 目标:了解命令流的整体结构,掌握每个模块的标识 !文件说明段 /BATCH /TITILE,test analysis !定义工作标题 /FILENAME,test !定义工作文件名 /PREP7 !进入前处理模块标识 !定义单元,材料属性,实常数段 ET,1,SHELL63 !指定单元类型 ET,2,SOLID45 !指定体单元 MP,EX,1,2E8 !指定弹性模量 MP,PRXY,1,0.3 !输入泊松比 MP,DENS,1,7.8E3 !输入材料密度
R,1,0.001 !指定壳单元实常数-厚度 ...... !建立模型 K,1,0,0,, !定义关键点 K,2,50,0,, K,3,50,10,, K,4,10,10,, K,5,10,50,, K,6,0,50,, A,1,2,3,4,5,6, !由关键点生成面 ...... !划分网格 ESIZE,1,0, AMESH,1 ...... FINISH !前处理结束标识 /SOLU !进入求解模块标识 !施加约束和载荷 DL,5,,ALL SFL,3,PRES,1000 SFL,2,PRES,1000 ...... SOLVE !求解标识 FINISH !求解模块结束标识 /POST1 !进入通用后处理器标识 ...... /POST26 !进入时间历程后处理器 …… /EXIT,SAVE !退出并存盘 以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助 /ANGLE !指定绕轴旋转视图 /DIST !说明对视图进行缩放
Ansys命令流 第一天 目标:熟悉ANSYS基本关键字的含义 k --> Keypoints 关键点 l --> Lines 线 a --> Area 面 v --> Volumes 体 e --> Elements 单元 n --> Nodes 节点 cm --> component 组元 et --> element type 单元类型 mp --> material property 材料属性 r --> real constant 实常数 d --> DOF constraint 约束 f --> Force Load 集中力 sf --> Surface load on nodes 表面载荷 bf --> Body Force on Nodes 体载荷 ic --> Initial Conditions 初始条件 第二天 目标:了解命令流的整体结构,掌握每个模块的标识 !文件说明段 /BATCH /TITILE,test analysis !定义工作标题 /FILENAME,test !定义工作文件名 /PREP7 !进入前处理模块标识 !定义单元,材料属性,实常数段 ET,1,SHELL63 !指定单元类型 ET,2,SOLID45 !指定体单元 MP,EX,1,2E8 !指定弹性模量 MP,PRXY,1,0.3 !输入泊松比 MP,DENS,1,7.8E3 !输入材料密度 R,1,0.001 !指定壳单元实常数-厚度...... !建立模型 K,1,0,0,, !定义关键点
K,2,50,0,, K,3,50,10,, K,4,10,10,, K,5,10,50,, K,6,0,50,, A,1,2,3,4,5,6, !由关键点生成面 ...... !划分网格 ESIZE,1,0, AMESH,1 ...... FINISH !前处理结束标识 /SOLU !进入求解模块标识 !施加约束和载荷 DL,5,,ALL SFL,3,PRES,1000 SFL,2,PRES,1000 ...... SOLVE !求解标识 FINISH !求解模块结束标识 /POST1 !进入通用后处理器标识 ...... /POST26 !进入时间历程后处理器 …… /EXIT,SAVE !退出并存盘 以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助 /ANGLE !指定绕轴旋转视图 /DIST !说明对视图进行缩放 /DEVICE !设置图例的显示,如:风格,字体等 /REPLOT !重新显示当前图例 /RESET !恢复缺省的图形设置 /VIEW !设置观察方向
ANSYS最常用命令流+中文注释 VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes，用于2个solid相减操作，最终目的是要nv1-nv2=?通过后面的参数设置，可以得到很多种情况：sepo项是2个体的边界情况，当缺省的时候，是表示2个体相减后，其边界是公用的，当为sepo的时候，表示相减后，2个体有各自的独立边界。keep1与keep2是询问相减后，保留哪个体？当第一个为keep时，保留nv1,都缺省的时候，操作结果最终只有一个体，比如：vsbv,1,2,sepo,,keep,表示执行1-2的操作，结果是保留体2，体1被删除，还有一个1-2的结果体，现在一共是2个体（即1-2与2），且都各自有自己的边界。如vsbv,1,2,,keep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。同理，将v换成a 及l是对面和线进行减操作！ mp,lab, mat, co, c1,…….c4 定义材料号及特性 lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens） ex: 弹性模量 nuxy: 小泊松比 alpx: 热膨胀系数 reft: 参考温度 reft: 参考温度 prxy: 主泊松比 gxy: 剪切模量 mu: 摩擦系数 dens: 质量密度 mat: 材料编号（缺省为当前材料号） co: 材料特性值，或材料之特性，温度曲线中的常数项 c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数 定义DP材料： 首先要定义EX和泊松比：MP，EX，MA T，…… MP，NUXY，MAT，…… 定义DP材料单元表（这里不考虑温度）：TB，DP，MA T 进入单元表并编辑添加单元表：TBDATA，1，C TBDATA，2，ψ TBDATA，3，…… 如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8 MP，NUXY，1，0.3 TB，DP，1 TBDATA，1，27 TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg VSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、反选（inv）等,其余方式不常用 Item, Comp 是选取的原则以及下面的子项 如volu 就是根据实体编号选择， loc 就是根据坐标选取，它的comp就可以是实体的某方向坐标！ 其余还有材料类型、实常数等 MIN, VMAX, VINC，这个就不必说了吧！ ,例：vsel,s,volu,,14 vsel,a,volu,,17,23,2 上面的命令选中了实体编号为14，17，19，21，23的五个实体 VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体 nv1:初始体号 nv2:最终的体号 ninc:体号之间的间隔 kswp=0:只删除体 kswp=1:删除体及组成关键点,线面 如果nv1=all,则nv2,ninc不起作用 其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令，输入后直接回车，就可以显示刚刚选择了的体、面或节点，很实用的哦！ Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备 Type: S: 选择一组新节点（缺省） R: 在当前组中再选择 A: 再选一组附加于当前组 U: 在当前组中不选一部分 All: 恢复为选中所有 None: 全不选 Inve: 反向选择 Stat: 显示当前选择状态 Item: loc: 坐标 node: 节点号
一、定义材料号及特性 mp,lab, mat, co, c1,…….c4 lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens） ex: 弹性模量 nuxy: 小泊松比 alpx: 热膨胀系数 reft: 参考温度 reft: 参考温度 prxy: 主泊松比 gxy: 剪切模量 mu: 摩擦系数 dens: 质量密度 mat: 材料编号（缺省为当前材料号） c 材料特性值，或材料之特性，温度曲线中的常数项 c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数二、定义DP材料： 首先要定义EX和泊松比：MP，EX，MAT，…… MP，NUXY，MAT，…… 定义DP材料单元表（这里不考虑温度）：TB，DP，MAT 进入单元表并编辑添加单元表：TBDATA，1，C TBDATA，2，ψ TBDATA，3，……
如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下： MP，EX，1，1E8 MP，NUXY，1，0.3 TB，DP，1 TBDATA，1，27 TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg 三、单元生死载荷步 !第一个载荷步 TIME,... !设定时间值（静力分析选项） NLGEOM,ON !打开大位移效果 NROPT,FULL !设定牛顿－拉夫森选项 ESTIF,... !设定非缺省缩减因子（可选） ESEL,... !选择在本载荷步中将不激活的单元 EKILL,... !不激活选择的单元 ESEL,S,LIVE !选择所有活动单元 NSLE,S !选择所有活动结点 NSEL,INVE !选择所有非活动结点（不与活动单 元相连的结点） D,ALL,ALL,0 !约束所有不活动的结点自由度（可 选） NSEL,ALL !选择所有结点
1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9 此命令用已知的一组关键点点（P1~P9）来定义面（Area）， 最少使用三个点才能围成面，同时产生转围绕些面的线。 点要依次序输入，输入的顺序会决定面的法线方向。 如果超过四个点，则这些点必须在同一个平面上。 Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPs 2、*ABBR，Abbr,String－－定义一个缩略语． Abbr:用来表示字符串＂String＂的缩略语，长度不超过８个字符． String：将由＂Abbr＂表示的字符串，长度不超过６０个字符．3、ABBRES，Lab，Fname，Ext－从一个编码文件中读出缩略语． Lab：指定读操作的标题， NEW：用这些读出的缩略语重新取代当前的缩略语（默认）CHANGE：将读出的缩略语添加到当前缩略语阵列，并替代现存同名的缩略语． Ext：如果＂Fname＂是空的，则缺省的扩展命是＂ABBR＂．4、ABBSA V，Lab，Fname，Ext－将当前的缩略语写入一个文本文件里 Lab：指定写操作的标题，若为ALL，表示将所有的缩略语都写入文件（默认） 5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc 将ia,ib,ic变量相加赋给ir变量 ir, ia,ib,ic：变量号 name: 变量的名称 6、Adele,na1,na2,ninc,kswp ！kswp=0时只删除掉面积本身，＝1时低单元点一并删除。 7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 ！面积的建立，沿某组线段路径，拉伸而成。 8、Afillt,na1,na2,rad ！建立圆角面积，在两相交平面间产生曲面，rad为半径。 9、*AFUN，Lab 在参数表达式中，为角度函数指定单位． Lab：指定将要使用的角度单位．有３个选项． RAD：在角度函数的输入与输出中使用弧度单位（默认）DEG：在角度函数的输入与输出中使用度单位． STAT：显示该命令当前的设置（即是度还是弧度）． 10、Agen, itime,na1,na2,ninc,dx,dy,dz,kinc,noelem,imove ！面积复制命令。itime包含本身所复制的次数；na1,na2,ninc为现有的坐标系统下复制到其他位置(dx,dy,dz)；kinc为每次复制时面积号码的增加量。 11、AINV, NA, NV 面与体相交生成一个相交面． NA, NV ：分别为指定面，指定体的编号．其中NA可以为P．说明：面与体相交生成新面．如果相交的区域是线，则生成新线． 指定源实体的单元属性和边界条件不会转换到新生成的实 体上． 12、AL,L1,L2,L3,L4,L5,L6,L7,L8,L9,L10 此命令由已知的一组直线（L1,…L10）围绕成面（Area）， 至少须要3条线才能形成面，线段的号码没有严格的顺序限制，只要它们能完成封闭的面积即可。 同时若使用超过4条线去定义面时，所有的线必须在同一平面上，以右手定则来决定面积的方向。如果L1为负号，则反向。Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>By Lines 13、ALLSEL, LABT, ENTITY 选中所有项目 LABT: ALL: 选所有项目及其低级项目 BELOW: 选指定项目的直接下属及更低级项目 ENTITY: ALL: 所有项目（缺省） VOLU:体高级 AREA:面 LINE :线 KP:关键点 ELEM:单元 NODE:节点低级 14、Amesh, nA1,nA2,ninc 划分面单元网格nA1,nA2,ninc 待划分的面号，nA1如果是All,则对所有选中面划分 15、ANORM, ANUM, NOEFLIP 修改面的正法线方向． ANUM：面的编号，改变面的正法线方向与面的法线方向相同．NOEFLIP：确定是否要改变重定向面上单元的正法线方向，这样可以使他们与面的正法线方向一致 若为0，改变单元的正法线方向； 若为1，不改变已存在单元的正法线方向； 说明：重新改变面的方向使得他们与指定的正法线方向相同. 不能用＂ANORM＂命令改变具体或面载荷的任何单元的正法线方向． 16、数学函数 ABS(X) 求绝对值 ACOS(X) 反余弦 ASIN(X) 反正弦 ATAN(X) 反正切 ATAN2(X,Y) 反正切, ArcTangent of (Y/X) , 可以考虑变量X,Y 的符号 COS(X) 求余弦 COSH(X) 双曲余弦 EXP(X) 指数函数 GDIS(X,Y) 求以X为均值,Y为标准差的高斯分布,在使用蒙地卡罗法研究随机荷载和随机材料参数时,可以用该函数处理计算结果 LOG(X) 自然对数 LOG10(X) 常用对数(以10为基) MOD(X,Y) 求X/Y的余数. 如果Y=0, 函数值为0 NINT(X) 求最近的整数 RAND(X,Y) 取随机数,其中X 是下限, Y是上限 SIGN(X,Y) 取X的绝对值并赋予Y的符号. Y>=0, 函数值为|X|, Y<0, 函数值为-|X|,. SIN(X) 正弦 SINH(X) 双曲正弦 SQRT(X) 平方根 TAN(X) 正切 TANH(X) 双曲正切 17、antype, status, ldstep, substep, action 声明分析类型，即欲进行哪种分析，系统默认为静力学分析。 antype: static or 1 静力分析 buckle or 2 屈曲分析 modal or 3 模态分析 trans or 4 瞬态分析 status: new 重新分析（缺省），以后各项将忽略 rest 再分析，仅对static,full transion 有效 ldstep: 指定从哪个荷载步开始继续分析，缺省为最大的，runn数（指分析点的最后一步） substep: 指定从哪个子步开始继续分析。缺省为本目录中，runn文件中最高的子步数 action, continue: 继续分析指定的ldstep,substep 说明：继续以前的分析（因某种原因中断）有两种类型singleframe restart: 从停止点继续 需要文件：jobname.db 必须在初始求解后马上存盘 jobname.emat 单元矩阵 jobname.esav 或.osav : 如果.esav坏了，将.osav 改为.esav results file: 不必要，但如果有，后继分析的结果也将很好地附加到它后面 注意：如果初始分析生成了.rdb, .ldhi, 或rnnn 文件。必须删除再做后继分析 步骤：（1）进入anasys 以同样工作名 （2）进入求解器，并恢复数据库 （3）antype, rest （4）指定附加的荷载 （5）指定是否使用现有的矩阵（jobname.trl）（缺省重新生成） kuse: 1 用现有矩阵
ANSYS结构分析单元功能与特性 /可以组成一一些命令，一般是一种总体命令（session),三十也有特殊，比如是处理/POST1 ! 是注释说明符号，，与其他软件的说明是一样的，ansys不作为命令读取， * 此符号一般是APDL的标识符，也就是ansys的参数化语言，如*do ,,,*enddo等等 NSEL的意思是node select，即选择节点。s就是select，选择。 DIM是定义数组的意思。array 数组。 MP命令用来定义材料参数。 K是建立关键点命令。K,关键点编号,x坐标,y坐标，z坐标。K, NPT, X, Y, Z是定义关键点，K是命令，NPT是关键点编号，XYZ是坐标。 NUMMRG, keypoint 用这个命令，要保证关键点的位置完全一样，只是关键点号不一样的才行。这个命令对于重复的线面都可以用。这个很简单，压缩关键。 Ngen 复制节点 e,节点号码：这个命令式通过节点来形成单元 NUMCMP,ALL：压缩所有编号，这样你所有的线都会按次序重新编号~你要是需要固定的线固定的标号NSUBST,100,500,50：通过指定子步数来设置载荷步的子步 LNSRCH线性搜索是求解非线性代数方程组的一种技巧，此法会在一段区间内，以一定的步长逐步搜索根，相比常用的牛顿迭代法所要耗费的计算量大得多，但它可以避免在一些情况下牛顿迭代法出现的跳跃现象。LNSRCH激活线性搜索 PRED 激活自由度求解预测 NEQIT指定一个荷载步中的最大子步数 AUTOTS 自动求解控制打开自动时间步长. KBC -指定阶段状或者用跳板装载里面一个负荷步骤。 SPLINE:P1，P2，P3，P4，P5，P6，XV1，YV1，ZV1，XV6，YV6，ZV6（生成分段样条曲线） *DIM，Par，Type，IMAX，JMAX，KMAX，Var1，Var2，Var3（定义载荷数组的名称） 【注】Par: 数组名 Type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省） char 字符串组（每个元素最多8个字符） table IMAX，JMAX，KMAX各维的最大下标号 Var1，Var2，Var3 各维变量名，缺省为row,column,plane(当type为table时) /config是设置ansys配置参数的 命令格式为/CONFIG, Lab, V ALUE Lab为参数名称value为参数值 例如：/config，MXEL，10000的意思是最大单元数为10000 杆单元:LINK1、8、10、11、180 梁单元：BEAM3、4、23、24，44，54，188，189 管单元:PIPE16，17，18，20，59，60 2D实体元:PLANE2，25，42，82，83，145，146，182，183 3D实体元:SOLID45，46，64，65，72，73，92，95，147，148，185，186，187，191 壳单元:SHELL28，41，43，51，61，63，91，93，99，143，150，181，208，209 弹簧单元:COMBIN7，14，37，39，40 质量单元:MASS21 接触单元:CONTAC12，52，TARGE169，170，CONTA171，172，173，174，175，178 矩阵单元:MATRIX27，50
全套完整版ANSYS命令流教学手册
————————————————————————————————作者：————————————————————————————————日期：
目录 1ANSYS结构分析单元功能与特性 0 1.1杆单元: LINK1、8、10、11、180 0 1.2梁单元:BEAM3、4、23、24，44，54，188，189 (1) 1.3管单元:PIPE16，17，18，20，59，60 (1) 1.42D实体单元:PLANE2，25，42，82，83，145，146，182， 183 2 1.53D实体元:SOLID45，46，64，65，72，73，92，95，147，148，185，186，187，191 (4) 1.6壳单元:SHELL28，41，43，51，61，63，91，93，99，143，150，181，208，209 (5) 1.7弹簧单元:COMBIN7，14，37，39，40 (5) 1.8质量单元:MASS21 (5) 1.9接触单元:CONTAC12，52，TARGE169，170，CONTA171，172，173，174，175，178 (5) 1.10 矩阵单元:MATRIX27，50 (6) 1.11 表面效应元:SURF153，154 (6) 1.12 粘弹实体元:VISCO88，89，106，107，108， (7) 1.13 超弹实体元:HYPER56，58，74，84，86，158 (7) 1.14 耦合场单元:SOLID5，PLANE13，FLUID29，30，38，SOLID62，FLUID79，FLUID80，81，SOLID98，FLUID129，INFIN110，111，FLUID116，130 (7) 1.15 界面单元:INTER192，193，194，195 (7) 1.16 显式动力分析单元:LINK160，BEAM161，PLANE162，SHELL163，SOLID164，COMBI16 (7) 1.17 预紧、多点约束、网分单元 (7) 2ANSYS 的基本使用 (9)
ansys——ANSYS命令流（Ⅰ） 1. A，P1，P2，…，P17，P18（以点定义面） 2. AADD，NA1，NA2，…NA8，NA9（面相加） 3. AATT，MAT，REAL，TYPE，ESYS，SECN（指定面的单元属性） 【注】ESYS为坐标系统号、SECN为截面类型号。 4. *ABBR，Abbr，String（定义一个缩略词） 5. ABBRES，Lab，Fname，Ext（从文件中读取缩略词） 6. ABBSAVE，Lab，Fname，Ext（将当前定义的缩略词写入文件） 7. ABS，IR，IA，--，--，Name，--，--，FACTA（取绝对值） 【注】************* 8. ACCAT，NA1，NA2（连接面） 9. ACEL，ACEX，ACEY，ACEZ（定义结构的线性加速度） 10. ACLEAR，NA1，NA2，NINC（清除面单元网格） 11. ADAMS，NMODES，KSTRESS，KSHELL 【注】************* 12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC 【注】************* 13. ADD，IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC（变量加运算） 14. ADELE，NA1，NA2，NINC，KSWP（删除面） 【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。 15. ADRAG，NL1，NL2，…，NL6，NLP1，NLP2，…，NLP6（将既有线沿一定路径拖拉成面） 16. AESIZE，ANUM，SIZE（指定面上划分单元大小）
ANSYS常用命令 Fini(退出四大模块，回到BEGIN层) /cle (清空存，开始新的计算) 1．定义参数、数组，并赋值. 2． /prep7(进入前处理) 定义几何图形：关键点、线、面、体 定义几个所关心的节点，以备后处理时调用节点号。 设材料线弹性、非线性特性 设置单元类型及相应KEYOPT 设置实常数 设置网格划分，划分网格 根据需要耦合某些节点自由度 定义单元表 3．/solu 加边界条件 设置求解选项 定义载荷步 求解载荷步 4./post1（通用后处理） 5./post26 （时间历程后处理） 6.PLOTCONTROL菜单命令 7.参数化设计语言 8.理论手册 Fini(退出四大模块，回到BEGIN层) /cle (清空存，开始新的计算) 1定义参数、数组，并赋值. u dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组 par: 数组名 type： array 数组，如同fortran,下标最小号为1，可以多达三维（缺省） char 字符串组（每个元素最多8个字符） table imax,jmax, kmax 各维的最大下标号 var1,var2,var3 各维变量名，缺省为row,column,plane(当type为table时) 2 /prep7(进入前处理) 2.1 定义几何图形：关键点、线、面、体 u csys,kcn kcn , 0 迪卡尔zuobiaosi 1 柱坐标 2 球 4 工作平面 5 柱坐标系（以Y轴为轴心） n 已定义的局部坐标系 u numstr, label, value
Ansys命令流大全(整理)
1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9 此命令用已知的一组关键点点（P1~P9）来定义面（Area）， 最少使用三个点才能围成面，同时产生转围绕些面的线。 点要依次序输入，输入的顺序会决定面的法线方向。 如果超过四个点，则这些点必须在同一个平面上。 Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPs 2、*ABBR，Abbr,String－－定义一个缩略语． Abbr:用来表示字符串＂String＂的缩略语，长度不超过８个字符． String：将由＂Abbr＂表示的字符串，长度不超过６０个字符．3、ABBRES，Lab，Fname，Ext－从一个编码文件中读出缩略语． Lab：指定读操作的标题， NEW：用这些读出的缩略语重新取代当前的缩略语（默认）CHANGE：将读出的缩略语添加到当前缩略语阵列，并替代现存同名的缩略语． Ext：如果＂Fname＂是空的，则缺省的扩展命是＂ABBR＂．4、ABBSA V，Lab，Fname，Ext－将当前的缩略语写入一个文本文件里 Lab：指定写操作的标题，若为ALL，表示将所有的缩略语都写入文件（默认） 5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc 将ia,ib,ic变量相加赋给ir变量 ir, ia,ib,ic：变量号 name: 变量的名称 6、Adele,na1,na2,ninc,kswp ！kswp=0时只删除掉面积本身，＝1时低单元点一并删除。 7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 ！面积的建立，沿某组线段路径，拉伸而成。 8、Afillt,na1,na2,rad ！建立圆角面积，在两相交平面间产生曲面，rad为半径。 9、*AFUN，Lab 在参数表达式中，为角度函数指定单位． Lab：指定将要使用的角度单位．有３个选项． RAD：在角度函数的输入与输出中使用弧度单位（默认）DEG：在角度函数的输入与输出中使用度单位． STAT：显示该命令当前的设置（即是度还是弧度）． 10、Agen, itime,na1,na2,ninc,dx,dy,dz,kinc,noelem,imove ！面积复制命令。itime包含本身所复制的次数；na1,na2,ninc为现有的坐标系统下复制到其他位置(dx,dy,dz)；kinc为每次复制时面积号码的增加量。 11、AINV, NA, NV 面与体相交生成一个相交面． NA, NV ：分别为指定面，指定体的编号．其中NA可以为P．说明：面与体相交生成新面．如果相交的区域是线，则生成新线． 指定源实体的单元属性和边界条件不会转换到新生成的实 体上． 12、AL,L1,L2,L3,L4,L5,L6,L7,L8,L9,L10 此命令由已知的一组直线（L1,…L10）围绕成面（Area）， 至少须要3条线才能形成面，线段的号码没有严格的顺序限制，只要它们能完成封闭的面积即可。 同时若使用超过4条线去定义面时，所有的线必须在同一平面上，以右手定则来决定面积的方向。如果L1为负号，则反向。Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>By Lines 13、ALLSEL, LABT, ENTITY 选中所有项目 LABT: ALL: 选所有项目及其低级项目 BELOW: 选指定项目的直接下属及更低级项目 ENTITY: ALL: 所有项目（缺省） VOLU:体高级 AREA:面 LINE :线KP:关键点 ELEM:单元 NODE:节点低级 14、Amesh, nA1,nA2,ninc 划分面单元网格 nA1,nA2,ninc 待划分的面号，nA1如果是All,则对所有选中面划分 15、ANORM, ANUM, NOEFLIP 修改面的正法线方向． ANUM：面的编号，改变面的正法线方向与面的法线方向相同．NOEFLIP：确定是否要改变重定向面上单元的正法线方向，这样可以使他们与面的正法线方向一致 若为0，改变单元的正法线方向； 若为1，不改变已存在单元的正法线方向； 说明：重新改变面的方向使得他们与指定的正法线方向相同. 不能用＂ANORM＂命令改变具体或面载荷的任何单元的 正法线方向． 16、数学函数 ABS(X) 求绝对值 ACOS(X) 反余弦 ASIN(X) 反正弦 ATAN(X) 反正切 ATAN2(X,Y) 反正切, ArcTangent of (Y/X) , 可以考虑变量X,Y 的符号 COS(X) 求余弦 COSH(X) 双曲余弦 EXP(X) 指数函数 GDIS(X,Y) 求以X为均值,Y为标准差的高斯分布,在使用蒙地卡罗法研究随机荷载和随机材料参数时,可以用该函数处理计 算结果 LOG(X) 自然对数 LOG10(X) 常用对数(以10为基) MOD(X,Y) 求X/Y的余数. 如果Y=0, 函数值为0 NINT(X) 求最近的整数 RAND(X,Y) 取随机数,其中X 是下限, Y是上限 SIGN(X,Y) 取X的绝对值并赋予Y的符号. Y>=0, 函数值为|X|, Y<0, 函数值为-|X|,. SIN(X) 正弦 SINH(X) 双曲正弦 SQRT(X) 平方根 TAN(X) 正切 TANH(X) 双曲正切 17、antype, status, ldstep, substep, action 声明分析类型，即欲进行哪种分析，系统默认为静力学分析。 antype: static or 1 静力分析 buckle or 2 屈曲分析 modal or 3 模态分析 trans or 4 瞬态分析 status: new 重新分析（缺省），以后各项将忽略 rest 再分析，仅对static,full transion 有效 ldstep: 指定从哪个荷载步开始继续分析，缺省为最大的，runn数（指分析点的最后一步） substep: 指定从哪个子步开始继续分析。缺省为本目录中，runn文件中最高的子步数 action, continue: 继续分析指定的ldstep,substep 说明：继续以前的分析（因某种原因中断）有两种类型singleframe restart: 从停止点继续 需要文件：jobname.db 必须在初始求解后马上存盘 jobname.emat 单元矩阵 jobname.esav 或.osav : 如果.esav坏了，将.osav 改为.esav results file: 不必要，但如果有，后继分析的结果也将很好地附加到它后面 注意：如果初始分析生成了.rdb, .ldhi, 或rnnn 文件。必须删除再做后继分析
有限元分析软件ANSYS命令流中文说明 Command VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes，用于2个solid相减操作，最终目的是要nv1-nv2=?通过后面的参数设置，可以得到很多种情况：sepo项是2个体的边界情况，当缺省的时候，是表示2个体相减后，其边界是公用的，当为sepo的时候，表示相减后，2个体有各自的独立边界。keep1与keep2是询问相减后，保留哪个体？当第一个为keep时，保留nv1,都缺省的时候，操作结果最终只有一个体，比如：vsbv,1,2,sepo,,keep,表示执行1-2的操作，结果是保留体2，体1被删除，还有一个1-2的结果体，现在一共是2个体（即1-2与2），且都各自有自己的边界。如vsbv,1,2,,keep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。同理，将v换成a及l是对面和线进行减操作！ mp,lab, mat, co, c1,…….c4 定义材料号及特性 lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens） ex: 弹性模量 nuxy: 小泊松比 alpx: 热膨胀系数 reft: 参考温度 reft: 参考温度 prxy: 主泊松比 gxy: 剪切模量 mu: 摩擦系数
dens: 质量密度 mat: 材料编号（缺省为当前材料号） co: 材料特性值，或材料之特性，温度曲线中的常数项 c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数 定义DP材料： 首先要定义EX和泊松比：MP，EX，MAT，…… MP，NUXY，MAT，…… 定义DP材料单元表（这里不考虑温度）：TB，DP，MAT 进入单元表并编辑添加单元表：TBDATA，1，C TBDATA，2，ψ TBDATA，3，…… 如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下： MP，EX，1，1E8 MP，NUXY，1，0.3 TB，DP，1 TBDATA，1，27 TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg VSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、反选（inv）等,其余方