当前位置：文档库 › (完整版)用Opensees进行IDA分析(桥墩模型命令流)

(完整版)用Opensees进行IDA分析(桥墩模型命令流)

wipe

# Opensees dandun

#Units:kN, m, sec

# -----------------

# Start of model generation

# -----------------

# Create ModeBulider (with two-dimensions and 3 DOF/node) model basic -ndm 2 -ndf 3

# -----------------

# tag X Y

node 1 0.0 0.0

node 2 0.0 0.0

node 3 0.0 2.0

node 4 0.0 4.0

node 5 0.0 6.0

node 6 0.0 8.0

node 7 0.0 10.0

node 8 0.0 12.0

node 9 0.0 14.0

node 10 0.0 16.0

node 11 0.0 18.0

node 12 0.0 20.0

# -----------------

# Fix supports at base of columns

# tag DX DY RZ

fix 1 1 1 1

# ----------------

# Concrete tag fc ec0 fcu ecu # Core concrete (confined)

uniaxialMaterial Concrete01 1 -25600.0 -0.00219 -17780.0 -0.01 #Cover concrete (unconfined)

uniaxialMaterial Concrete01 2 -23400.0 -0.002 -0.0 -0.006

# STEEL

# Reinforcing steel

set fy 400000.0; #Yield stress

set E 200000000.0;# Young's modulus

# tag fy E0 b

uniaxialMaterial Steel02 3 $fy $E 0.01 18.5 0.925 0.15 uniaxialMaterial Elastic 11 29043600

uniaxialMaterial Elastic 12 12326600

uniaxialMaterial Elastic 13 587247596

#Define cross-section for nonlinear columns

# ---------------------

# set some parameters

set colWidth 8.18

set colDepth 4.28

set cover 0.05

set As 0.00049

# some variables derived from the parameters

set y1 [expr $colDepth/2.0]

set z1 [expr $colWidth/2.0]

section Fiber 1 {

# Create the concrete core fibers

patch rect 1 20 30 [expr $cover-$y1] [expr $cover-$z1] [expr $y1-$cover] [expr $z1-$cover]

# Create the concrete cover fibers (top, bottom, left, right)

patch rect 2 20 5 [expr -$y1] [expr $z1-$cover] $y1 $z1

patch rect 2 20 5 [expr -$y1] [expr -$z1] $y1 [expr $cover-$z1]

patch rect 2 5 10 [expr -$y1] [expr $cover-$z1] [expr $cover-$y1] [expr $z1-$cover]

patch rect 2 5 10 [expr $y1-$cover] [expr $cover-$z1] $y1 [expr $z1-$cover]

# Create the reinforcingfibers (left, middle, right)

layer straight 3 175 $As [expr $y1-$cover] [expr $z1-$cover] [expr $y1-$cover] [expr $cover-$z1] layer straight 3 175 $As [expr $cover-$y1] [expr $z1-$cover] [expr $cover-$y1] [expr $cover-$z1] layer straight 3 115 $As [expr $y1-$cover] [expr $z1-$cover] [expr $cover-$y1] [expr $z1-$cover] layer straight 3 115 $As [expr $y1-$cover] [expr $cover-$z1] [expr $cover-$y1] [expr $cover-$z1] }

# Define column elements

# ----------------------

# Geometry of column elements

# tag

geomTransf Linear 1

# Number of integration points along length of elements

set np 5

# Create the columns using Beam-column elements

# tag ndI ndJ secID transfTag

element nonlinearBeamColumn 2 2 3 $np 1 1

element nonlinearBeamColumn 3 3 4 $np 1 1

element nonlinearBeamColumn 4 4 5 $np 1 1

element nonlinearBeamColumn 5 5 6 $np 1 1

element nonlinearBeamColumn 6 6 7 $np 1 1

element nonlinearBeamColumn 7 7 8 $np 1 1

element nonlinearBeamColumn 8 8 9 $np 1 1

element nonlinearBeamColumn 9 9 10 $np 1 1

element nonlinearBeamColumn 10 10 11 $np 1 1

element nonlinearBeamColumn 11 11 12 $np 1 1

equalDOF 1 2 1 2

element zeroLength 111 1 2 -mat 13 -dir 3

set m [expr 3355.7]

set n [expr 175.1]

# tag MX MY RZ

mass 2 $n $n 1e-8

mass 3 $n $n 1e-8

mass 4 $n $n 1e-8

mass 5 $n $n 1e-8

mass 6 $n $n 1e-8

mass 7 $n $n 1e-8

mass 8 $n $n 1e-8

mass 9 $n $n 1e-8

mass 10 $n $n 1e-8

mass 11 $n $n 1e-8

mass 12 $m $m 1e-8

# Set a parameter for the axial load

set P 33557.0; # of axial capacity of columns

set Q 1715.5; # of axial capacity of columns

# Create a Plain load pattern with a Linear TimeSeries pattern Plain 1 "Constant" {

# Create nodal loads at nodes 9

# nd FX FY MZ

load 12 0.0 [expr -$P] 0.0

load 2 0.0 [expr -$Q] 0.0

load 3 0.0 [expr -$Q] 0.0

load 4 0.0 [expr -$Q] 0.0

load 5 0.0 [expr -$Q] 0.0

load 6 0.0 [expr -$Q] 0.0

load 7 0.0 [expr -$Q] 0.0

load 8 0.0 [expr -$Q] 0.0

load 9 0.0 [expr -$Q] 0.0

load 10 0.0 [expr -$Q] 0.0

load 11 0.0 [expr -$Q] 0.0

}

system SparseGeneral -piv

# Create the constraint handler

constraints Transformation

# Create the time integration scheme

# Create the DOF numberer

numberer RCM

test NormDispIncr 1.0e-8 30 5

# Create the solution algorithm

algorithm Newton

# create the transient analysis

integrator LoadControl 1

# -----------------

# End of analysis generation

# -------------------

analysis Static

initialize

analyze 1

# End of static analysis

# -------------------

# --建立文件

# 第1条波---15

set dataDir Dynamic-Output-15;# name of output folder

file mkdir $dataDir; # create output folder

# ------

# 定义峰值加速度（g）

foreach pga {0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5} {

puts [format "The PGA is %3.2f g" $pga]

recorder Node -file [format "$dataDir/zhuanjiao%4.3f.out" $pga] -time -node 2 -dof 3 disp recorder Node -file [format "$dataDir/new%4.3f.out" $pga] -time -node 12 -dof 1 disp

recorder plot [format "$dataDir/new%4.3f.out" $pga] Node_Disp 650 0 710 390 -columns 1 2 recorder Element -file [format "$dataDir/Curvature-%4.3f.out" $pga] -time -ele 2 section 1 deformation

recorder Element -file [format "$dataDir/Force-%4.3f.out" $pga] -time -ele 2 section 1 force wipeAnalysis

setTime 0.0

set g [expr 6.289*9.81*$pga]

set tabasFN "Path -filePath 15.txt -dt 0.01 -factor $g"

#Define the excitation using the 960 ground motion records

# tag dir accel series args

pattern UniformExcitation 2 1 -accel $tabasFN

# add some mass proportional damping

# define DAMPING

#---------------------------------------------------------------------

set xDamp 0.05;# 2% damping ratio

set lambda [eigen 1]; # eigenvalue mode 1

set omega [expr pow($lambda,0.5)];

set alphaM 0.;# M-prop. damping; D = alphaM*M

set betaKcurr 0.; # K-proportional damping; +beatKcurr*KCurrent set betaKcomm [expr 2.*$xDamp/($omega)]; # K-prop. damping parameter; +betaKcomm*KlastCommitt

set betaKinit 0.; # initial-stiffness proportional damping +beatKinit*Kini

# define damping

rayleigh $alphaM $betaKcurr $betaKinit $betaKcomm; # RAYLEIGH damping

#---------------------------------------------------------------------

# Create the convergence test

test EnergyIncr 1.0e-8 30 5

# Create the solution algorithm

algorithm Newton

# Create the system of equation storage and solver

system SparseGeneral -piv

# Create the constraint handler

constraints Transformation

# Create the time integration scheme

integrator Newmark 0.5 0.25

# Create the DOF numberer

numberer RCM

# create the transient analysis

analysis VariableTransient

recorder Element -file [format "$dataDir/ele1Concrete-%4.3f.out" $pga] -time -ele 2 section 1 fiber -2.05 -4.0 1 stressStrain

recorder Element -file [format "$dataDir/ele1sec1StressStraingangjin%4.3f.out" $pga] -time -ele 2 section 1 fiber -2.05 -4.0 3 stressStrain

recorder plot [format "$dataDir/ele1sec1StressStraingangjin%4.3f.out" $pga] strain-stress 0 300 700 400 -columns 3 2

recorder Element -file [format "$dataDir/ele1-%4.3f.out" $pga] -time -ele 1 localForce

recorder Element -file [format "$dataDir/ele2-%4.3f.out" $pga] -time -ele 2 localForce

#Perform the analysis

# numSteps dt

analyze 4000 0.005 0.0000000001 0.005 30

# ------

set PI 3.1415926

set lambda [eigen 1]

set omega [expr pow($lambda,0.5)]

set Tperiod [expr 2*$PI/$omega]

puts "T1: $Tperiod"

remove loadPattern 2

remove recorders;

reset;

}

注：同文目录下需要编号为15.txt的地震记录文件，时间间隔为0.01s，PGA为0.159g，如下：

-0.0124

-0.00268

0.0082

0.0198

0.0377 0.0378 0.0309 0.0201 0.00885 0.00074 -0.00136 0.0027 0.0105 0.0189

0.025 0.0274 0.0267

0.026 0.0282 0.0345 0.0436 0.0532 0.0595 0.0594 0.0521 0.0405 0.0278 0.0164 0.00667 -0.00233 -0.0117 -0.0221 -0.0331 -0.0444 -0.0552 -0.0628 -0.064 -0.0578 -0.0464 -0.0334 -0.0213 -0.0109 -0.00136 0.00811

0.0169

0.0234

0.0268

0.0265

0.0251

0.023

0.0186

0.0106 0.000209 -0.00818 -0.0109 -0.00759 -0.000752 0.00608 0.00956 0.00827

0.0033 -0.00281 -0.00686 -0.00659 -0.00208 0.00454

0.0101

0.0121

0.0107 0.00815 0.00744 0.00959

0.0131

0.0152

0.0145

0.0113 0.00779 0.00629 0.00706 0.00692 0.00174 -0.00993 -0.0262 -0.0439 -0.0605

-0.0741

-0.0814

-0.0803

-0.0724

-0.061

-0.0398 -0.0313 -0.023 -0.014 -0.004 0.00632 0.0141 0.0152 0.00802 -0.00407 -0.0126 -0.0104 0.00346

0.024 0.0453 0.0637

0.077 0.0804 0.0699 0.0468 0.0167 -0.0125 -0.0303 -0.029 -0.0101 0.0178 0.0438 0.0612 0.0698 0.0742 0.0769 0.0757 0.0676 0.0525

0.033 0.0124 -0.00684 -0.0214 -0.0272 -0.0233 -0.0125 0.000534

0.0125

0.0312 0.0397 0.0458

0.046 0.0391 0.0271 0.0134 0.00139 -0.00477 -0.00219 0.00819 0.0226 0.0361 0.0423 0.0376 0.0236 0.00522 -0.0128 -0.0283 -0.0418 -0.0549 -0.0678 -0.0777 -0.0823 -0.0818 -0.0784 -0.0748 -0.073 -0.0736 -0.0756 -0.0774 -0.0765 -0.071 -0.0611 -0.0485 -0.0357 -0.0238 -0.0131 -0.0026 0.00822

0.0196

0.0299

0.034

0.0124 -0.00687 -0.0218 -0.0274 -0.0247 -0.0182 -0.0124 -0.00987 -0.0116 -0.0177 -0.027 -0.0378 -0.0483 -0.0582 -0.0672 -0.0724 -0.0688 -0.0548 -0.0331 -0.00849

0.015 0.0362

0.056 0.0759

0.097

0.119

0.137

0.145

0.141

0.128

0.113 0.0983 0.0875 0.0817 0.0814 0.0858 0.0923 0.0985

0.1

0.0938

0.0797

0.061

0.0417

0.00794 -0.00748 -0.0235 -0.0406 -0.0584 -0.0764 -0.0935 -0.106

-0.11 -0.104 -0.0925 -0.0788 -0.0664 -0.0563 -0.0482 -0.0441 -0.0465 -0.0553 -0.068 -0.0814 -0.0915 -0.0939 -0.0873 -0.0743 -0.059 -0.0449 -0.0328 -0.0219 -0.0107 0.00151 0.0148 0.0284 0.0419 0.0548 0.0673 0.0787 0.0869 0.0906

0.0905

0.0887

0.0865

0.0828

0.0746

0.0421

0.0218

0.0021 -0.0163 -0.0335

-0.048 -0.0566 -0.0575 -0.0524 -0.0447 -0.0375

-0.032 -0.0275 -0.0231

-0.018 -0.0119 -0.0056 -0.000207 0.00316 0.00446

0.0051

0.0077

0.014

0.0236

0.0348

0.0447

0.0491

0.0452

0.0344

0.0208

0.0105 0.00822

0.0139

0.0234

0.0302

0.0297

0.0221

0.011 0.000109 -0.00847

-0.0148

-0.0192

-0.0214

-0.0209 -0.0198 -0.0197 -0.021 -0.0219

-0.02 -0.0144 -0.00597 0.00335 0.0104 0.0116 0.00573 -0.00525 -0.017 -0.0233 -0.0201 -0.00854 0.00724 0.0228

0.036 0.0472

0.058 0.0691 0.0764 0.0728 0.0553 0.0271 -0.00512 -0.0359 -0.0613 -0.0772 -0.0822 -0.0791 -0.0727 -0.0671 -0.0635 -0.0613 -0.0589 -0.0546 -0.0467

-0.035 -0.0208 -0.00583

0.0221

0.0353

0.0485

0.0623

0.0761

0.0857

0.083

0.0639

0.0319 -0.00216 -0.0243 -0.0275 -0.0151 0.00342

0.0197

0.0294

0.0302

0.0226 0.00974 -0.00459 -0.0179 -0.0298 -0.0411

-0.053

-0.066 -0.0788 -0.0873 -0.0872 -0.0782 -0.0637 -0.0482 -0.0344 -0.0227 -0.0119 -0.000533

0.012

0.0252

0.0369

0.0442

0.0462

0.0447

0.0508

0.0767

0.0804

0.0696

0.0471

0.0198 -0.00578 -0.0257

-0.036 -0.0351 -0.0257 -0.0128 -0.000396 0.00988

0.0188

0.0277

0.0369

0.0419

0.0378

0.0236 0.00283 -0.0197 -0.0398 -0.0524 -0.0532 -0.0422 -0.0241 -0.00461

0.0129

0.0279

0.0413

0.0513

0.0545

0.0501

0.043

0.0411

0.0491

0.0656

0.0832

0.0904

0.0796

0.0533

0.0232

0.00148

-0.0075 -0.00916 -0.0192 -0.0384 -0.0623 -0.0858

-0.106

-0.124

-0.14

-0.154

-0.159

-0.15

-0.129

-0.102

-0.074 -0.0493 -0.0272 -0.00574

0.0168

0.0394

0.0558

0.0597

0.0508

0.0342

0.0163 0.000908 -0.0111 -0.0193 -0.0222 -0.0195 -0.0128 -0.00493 0.0000243 -0.00136 -0.00938 -0.0214

CAD常用命令汇总及详解

CAD中有哪些命令？我们可以把它们分为几类。一类是绘图类，二类是编辑类，三类是设置类，四类是其它类，包括标注、视图等。我们依次分析。第一类，绘图类。常用的命令有： Line 直线 Xline 构造线 mline 双线 pline 多义线 rectang 矩形 arc 圆弧 circle 圆 hatch 填充 boundary 边界 block 定义块 insert 插入快第二类，编辑类。常用的命令有： Matchprop 特性匹配 Hatchedit 填充图案编辑 Pedit 多义线编辑 Erase 擦除 Copy 拷贝 Mirror 镜像 Offset 平移 Array 阵列 Move 移动 Rotate 旋转 Scale 缩放 Stretch 拉伸 Lengthen 拉长 Trim 裁减 Extend 延伸 Break 打断 Fillet 倒圆角 Explode 炸裂 Align 对齐 Properties 属性

绘图工具栏: 直线（L）：全称(line) 在屏幕上指定两点可画出一条直线。也可用相对坐标或者在正交模式打开的情况下，直接给实际距离鼠标拖动来控制方向构造线（XL）：全称（xline） H为水平V为垂直O为偏移A为角度B为等分一个角度。多段线（PL）：全称（pline）首先在屏幕上指定一点，然后有相应提示：指定下一个点或[圆弧(A)/半宽(H)/长度(L)/放弃(U)/宽度(W)]。可根据需要来设置。其中“圆弧”指定宽度可画任意角度圆弧；“半宽”指多段线的一半宽度，即如要高线宽为10，则5；“长度”给相应的值，则画出相应长度的多段线；“放弃”指放弃一次操作；“宽度”指多段线的宽度多边形（pol）：全称（polygon）所绘制多边形为正多边形，边数可以自己设 E：根据边绘制多边形也可根据圆的半径利用外切和内接来画正多边形矩形（REC）：全称（rectang）点击矩形工具后出现下列提示：指定第一个角点或[倒角(C)/标高(E)/圆角(F)/厚度(T)/宽度(W)] 其中“倒角”是将90度直角的两条边割去一点。变成一个斜角。“标高”是空间上的意义可以在三视图当中展现出来,标高是相对的；“圆角”：即是将四个直角边倒成半径为X的圆角；“厚度”：空间上的意义，可在Z轴上表现出来“宽度”：平面空间的概念，指矩形四边的宽度。圆弧（ARC或A）：默认为3点画圆弧，成弧方向为逆时针，画优弧半径给负值。绘图菜单中有如下选项：起点、圆心、端点；起点、圆心、角度；起点、圆心、长度；起点、端点、角度；起点、端点、方向；起点、端点、半径；圆心、起点、端点；圆心、起点、角度；圆心、起点、长度；

OpenSEES重点笔记

1、利用零长单元模拟阻尼， uniaxialMaterial Elastic 1 6.8098e6; uniaxialMaterial Viscous 2 3.24e5 1; uniaxialMaterial Parallel 3 3 5; element zeroLength 1 $iNode $jNode -mat 3 -dir 1; 通常有两种方式：（1）truss element and viscous material.（桁架单元和阻尼材料）（2）force-based beam-column element and Maxwell material（基于力的梁柱单元和Maxwell 材料）。－、如何运行OpenSEES 有三种方法可以执行OpenSees/Tcl命令： 1、interactive交互式直接将命令输入Prompt。 2、执行文件输入这种方法是最常用的一种，以source inputfile.tcl方式执行已写好的外部命令文件。 3、Batch模式即以Opensees inputFile.tcl方式在MS－DOS/Unix promt中运行。二、定义单位和常数在编写一个较大的Opensees命令时。最好先定义好单位及常数。在Opensees中，编译器不能自行转换单位。所以一开始就要先定义好。单位定义包括两部分：首先定义基本单位；再定义合成单位。其中基本单位要相互独立。同时，在定义单位时，既可以按国际公制单位，也可以按私制单位。因些在单位定义文件中可能是混合的。我个人建议，还是采用国际公制单位较好。像国外常用英制单位。很不习惯。对于一些常数，如和g等常数要事先定义好。在定义这些单位时所用的命令是“set”。

ANSYS APDL命令流学习参数化建模

第一天目标:熟悉ANSYS基本关键字的含义 k --> Keypoints 关键点 l --> Lines 线 a --> Area 面 v --> V olumes 体 e --> Elements 单元 n --> Nodes 节点 cm --> component 组元 et --> element type 单元类型 mp --> material property 材料属性 r --> real constant 实常数 d --> DOF constraint 约束 f --> Force Load 集中力 sf --> Surface Force 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 !前处理结束标识

ORACLE SQLPLUS 常用命令及解释

Oracle SQLPlus常用命令及解释 1.@ 执行位于指定脚本中的SQLPlus语句。可以从本地文件系统或Web服务器中调用脚本。可以为脚本中的变量传递值。在iSQL*Plus中只能从Web服务器中调用脚本。 2.@@ 执行位于指定脚本中的SQL*Plus语句。这个命令和@（“at”符号）命令功能差不多。在执行嵌套的命令文件时它很有用，因为它会在与调用它的命令文件相同的路径或url中查找指定的命令文件。在iSQL*Plus中只支持url形式。 3./ 执行保存在SQL缓冲区中的最近执行的SQL命令或PL/SQL块。在SQL*Plus命令行中，可在命令提示符或行号提示符使用斜线（/）。也可在iSQL*Plus的输入区中使用斜线（/）。斜线不会列出要执行的命令。 4.ACCEPT 可以修改既有变量，也可定义一个新变量并等待用户输入初始值,读取一行输入并保存到给出的用户变量中。ACCEPT在iSQL*Plus中不可用。 5.APPEND 把指定文本添加到SQL缓冲区中当前行的后面。如果text的最前面包含一个空格可在APPEND和text间输入两个空格。如果text的最后是一个分号，可在命令结尾输入两个分号（SQL*Plus会把单个的分号解释为一个命令结束符）。APPEND 在iSQL*Plus中不可用。 6.ARCHIVE LOG 查看和管理归档信息。启动或停止自动归档联机重做日志，手工（显示地）归档指定的重做日志，或者显示重做日志文件的信息。 7.ATTRIBUTE 为对象类型列的给定属性指定其显示特性，或者列出单个属性或所有属性的当前显示特性。 8.BREAK 分开重复列。指定报表中格式发生更改的位置和要执行的格式化动作（例如，在列值每次发生变化时跳过一行）。只输入BREAK而不包含任何子句可列出当前的BREAK定义。 9.BTITLE 在每个报表页的底部放置一个标题并对其格式化，或者列出当前BTITLE定义。

ansys旋转经典命令流

1 旋转摩擦 (1) 2. 电磁三d命令流实例（论坛看到） (11) 3. 帮助感应加热例子induction heating of a solid cylinder billet (15) 4. 感应加热温度场的数值模拟（论文）inducheat30命令流 (19) 5. 如何施加恒定的角速度？Simwe仿真论坛 (24) 6. 旋转一个已经生成好的物体 (27) 7. 产生这样的磁力线 (28) 8. 旋转摩擦生热简单例子（二维旋转） (32) 8.1. 原版 (32) 8.2. 部分gui操作 (35) 9. VM229 Input Listing (39) 10 轴承---耦合+接触分析 (47) 11. 板的冲压仿真 (52) 1 旋转摩擦 FINISH /FILNAME,Exercise24 !定义隐式热分析文件名 /PREP7 !进入前处理器 ET,1,SOLID5 !选择单元类型 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7800 !定义材料1的密度 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,460 !定义材料1的比热 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,KXX,1,,66.6 !定义材料1的热传导系数 MPTEMP,,,,,,,, MPTEMP,1,0 UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度 MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,206e9 !定义材料1的弹性模量 MPDATA,PRXY,1,,0.3 !定义材料1的泊松比 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,8900 !定义材料2的密度 MPTEMP,,,,,,,, MPTEMP,1,0

OPENSEES

OPENSEES opensees中的单元问题梁柱单元 1. Nonlinear BeamColumn 基于有限单元柔度法理论。允许刚度沿杆长变化，通过确定单元控制截面各自的截面抗力和截面刚度矩阵,按照Gauss-Lobatto积分方法沿杆长积分计算出整个单元的抗力与切线刚度矩阵。NonlinearBeamColumn单元对于截面软化行为，构件反应由单元积分点数控制，为保证不同积分点数下构件反应的一致性，可以通过修正材料的应力-应变关系来实现，但同时会造成截面层次反应的不一致，因此需要在截面层次进行二次修正。一根构件不需要单元划分，使用1个单元即可，建议单元内使用4个截面积分点，截面上使用6*6的纤维积分点。 [5] 2. Displacement – Based BeamColumn 基于有限单元刚度法理论。允许刚度沿杆长变化，按照Gauss -Legendre积分方法沿杆长积分计算出整个单元的抗力与切线刚度矩阵。 Displacement - BasedBeam- Column单元对于截面软化行为，构件反应由遭受软化行为的单元长度控制，为保证计算结果的精确性，一般需要将构件离散为更多的单元，而截面层次的反应与构件的单元离散数无关，可以较为准确地反应截面的软化行为。建议一根构件划分为5个单元，单元内使用4个截面积分点，截面上使用6*6的纤维积分点。[5] 3. Beam With Hinges 基于有限单元柔度法理论。假定单元的非弹性变形集中在构件的两端，在杆件端部设置2个积分控制截面,并设定恰当的塑性铰长度，按照Gauss - Radau积分方法沿塑性铰长度积分来模拟构件和整体结构的非线性反应特点，而杆件中部的区段仍保持弹性。

ANSYS命令流学习笔记10-利用APDL在WorkBench中进行非线性屈曲分析

!ANSYS命令流学习笔记10-利用APDL在WorkBench中进行非线性屈曲分析 !学习重点： !1、强化非线性屈曲知识首先了解屈曲问题。在理想化情况下,当F < Fcr时, 结构处于稳定平衡状态，若引入一个小的侧向扰动力，然后卸载, 结构将返回到它的初始位置。当F > Fcr时, 结构处于不稳定平衡状态, 任何扰动力将引起坍塌。当F = Fcr时,结构处于中性平衡状态，把这个力定义为临界载荷。在实际结构中, 几何缺陷的存在或力的扰动将决定载荷路径的方向。在实际结构中, 很难达到临界载荷，因为扰动和非线性行为, 低于临界载荷时结构通常变得不稳定。要理解非线性屈曲分析，首先要了解特征值屈曲。特征值屈曲分析预测一个理想线弹性结构的理论屈曲强度，缺陷和非线性行为阻止大多数实际结构达到理想的弹性屈曲强度，特征值屈曲一般产生非保守解, 使用时应谨慎。 !理论解，根据Euler公式。其中μ取决于固定方式。 !有限元方法，已知在特征值屈曲问题：求解，即可得到临界载荷而非线性屈曲问题：其中为结构初始刚度，为有缺陷的结构刚度，为位移矩阵，为载荷矩阵。非线性屈曲分析时考虑结构平衡受扰动（初始缺陷、载荷扰动）的非线性静力分析，该分析时一直加载到结构极限承载状态的全过程分析，分析中可以综合考虑材料塑性、几何非线性、接触、大变形。非线性屈曲比特征值屈曲更精确，因此推荐用于设计或结构的评价。 !2、熟悉WB中非线性屈曲分析流程 (1) 前处理，施加单元载荷，进行预应力静力分析。 (2) 基于预应力静力分析，指定分析类型为特征值屈曲分析，完成特征值屈曲分析。 (3) 在APDL模块将一阶特征屈曲模态位移乘以适当系数，将此变形后的形状当做非线性分析的初始模型。

用Opensees进行IDA分析(桥墩模型命令流)

wipe # Openseesdandun # #Units:kN, m, sec # ----------------- # Start of model generation # ----------------- # CreateModeBulider (with two-dimensions and 3 DOF/node) model basic -ndm 2 -ndf 3 # ----------------- # tag X Y node 1 0.0 0.0 node 2 0.0 0.0 node 3 0.0 2.0 node 4 0.0 4.0 node 5 0.0 6.0 node 6 0.0 8.0 node 7 0.0 10.0 node 8 0.0 12.0 node 9 0.0 14.0 node 10 0.0 16.0 node 11 0.0 18.0 node 12 0.0 20.0 # ----------------- # Fix supports at base of columns # tag DX DY RZ fix 1 1 1 1 # ---------------- # Concrete tag fc ec0 fcuecu # Core concrete (confined) uniaxialMaterial Concrete01 1 -25600.0 -0.00219 -17780.0 -0.01 #Cover concrete (unconfined) uniaxialMaterial Concrete01 2 -23400.0 -0.002 -0.0 -0.006 # STEEL # Reinforcing steel setfy 400000.0; #Yield stress set E 200000000.0;# Young's modulus # tag fy E0 b uniaxialMaterial Steel02 3 $fy $E 0.01 18.5 0.925 0.15 uniaxialMaterial Elastic 11 29043600 uniaxialMaterial Elastic 12 12326600 uniaxialMaterial Elastic 13 587247596 #Define cross-section for nonlinear columns # ---------------------

个人总结ansys命令流

Q235 属性：弹性模量E=2.1e5 N/mm2 密度=7.85e-6kg/mm3 泊松比=0.3 mp,ex,1,2.1e5 mp,prxy,1,0.3 mp,dens,1,7.85e-6 1,ksymm 镜像点 2,arsym 镜像面 3,kgen 复制点 4.adele删除面 6，kdist,k1,k2 测量两关键点的距离 7,adele,a,,,1 删除area and below 8,创建圆柱面： circle 创建圆然后创建直线然(轴线) 利用拉伸命令创建圆柱面creat__areas__by Lines adrag 线拉伸成面modeling>operate>extrude>lines>>along lines VDRAG 面拉伸成体modeling>operate>extrude>areas>>along lines ！创建空心圆柱体这个命令 CYLIND, RAD1, RAD2, Z1, Z2, THETA1, THETA2 Main Menu>Preprocessor>Modeling>Create>Volumes>Cylinder>By Dimensions Main Menu>Preprocessor>Trefftz Domain>TZ Geometry>Create>Volume>Cylinder>By Dimensions 9,aptn 分割面 10,asbw 用工作平面切割面 11.wpoffs 12.wprota

https://www.wendangku.net/doc/d03014447.html,ng 过圆外一点做圆的切线(0°或180°) 14，nummrg 将重复的点消除 15，asba 面减去面 16,两个圆柱面的相贯线作法：做出两个相穿的圆柱面，利用APTN命令 17，选择面，不选择一部分面 asel,u,loc,z,kz(735) 18.在工作平面上生成一个矩形面 RECTING,X1,X2,Y1,Y2 X1，X2——矩形在工作平面X方向坐标值的变化范围 Y1,Y2——矩形在工作平面Y方向坐标值的变化范围 18,圆阵列建立工作平面与圆柱的横截面平行，在工作平面情况下建立局部坐标系（柱坐标系），然后利用agen命令复制。 19，转换成局部柱坐标系 20,kfill 在两个关键点之间生成一个或多个关键点 21.网格划分 aatt,1,14,1, !aatt,mat,real,type,esys,secn aesize,all,1000 !aesize,anum,size, 单元尺寸 mshape,0,2d !mshape,key,dimension 指定划分单元形状amesh,all k,1,24000,33000,2230 k,2,24000,33000,-2230 k,3,-24000,33000,-2230 k,4,-24000,33000,2230 kfill,2,3,23,5,1,1 kfill,1,4,23,28,1,1 *do,i,5,26 l,i,i+1 *enddo

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/TILE,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,,ALLSFL,3,PRES,1000SFL,2,PRES, 1000......SOLVE!求解标识FINISH!求解模块结束标识/POST1!进入通用后处理器标识....../POST26!进入时间历程后处理器……/EXIT,SAVE!退出并存盘以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助/ANGLE!指定绕轴旋转视图/DIST!说明对视图进行缩放/DEVICE!设置图例的显示,如：风格,字体等/REPLOT!重新显示当前图例/RESET!恢复缺省的图形设置/VIEW!设置观察方向/ZOOM!对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分 1.生成关键点K,关键点编号,X坐标,Y坐标,Z坐标例：

ANSYS APDL命令流建模及模态分析实例相关内容

本文介绍了轮毂的ANSYS APDL命令流建模及模态分析实例相关内容。 ANSYS命令流及注释五个辐条的轮毂 ! !初始化ANSYS环境 ! FINISH /CLEAR !清空内存 /FILNAM,WHEEL5 !文件名 /TITILE,WHEEL5 PARAMETER MODELING !工作名 ! !定义几何尺寸参数 ! R1=180 R2=157 R3=75 R4=75 R5=30 R6=28 R7=20 R8=90 R9=60 S_HOLE=5 TH1=48 TH2=23 TH3=11 TH4=180 TH5=40 TH6=45 TH7=105

TH8=25 TH9=15 TH10=25 TH11=13 /VIEW,1,1,1,1 !改变视图/ANG,1 /PNUM,LINE,1 /PNUM,AREA,1 /PNUM,VOLU,1 /NUMBER,1 ! !关键点 ! /PREP7 k,1,r5,r7,0 k,2,r4-ky(1),ky(1),0 k,3,r4,0,0 k,4,r1,0,0 k,5,kx(4),th5-th9,0 k,6,r1-th8,ky(5),0 k,7,kx(6),th4/2,0 k,8,kx(7)+th11,ky(7)+th10,0 k,9,kx(8),th4-th3,0 k,10,kx(4),ky(9),0 k,11,kx(4),th4,0 k,12,r2,ky(11),0 k,13,kx(12),ky(8),0 k,14,kx(7)-th3,ky(7),0 k,15,kx(14),th5,0 k,16,r3+r6,ky(15),0

Cad常用命令及使用方法

Cad常用命令及使用方法一、绘图命令直线：L 用法：输入命令L/回车/鼠标指定第一点/输入数值（也就是指定第二点）/回车（这时直线就画出来了）/回车（结束命令）射线：RAY 用法：输入命令RAY/回车/鼠标指定射线起点/指定通过点/回车（结束命令）构造线：XL 用法：输入命令XL/回车/鼠标指定构造线起点/指定通过点/回车（结束命令）多段线：PL 用法1：同直线命令用法2：输入命令PL/回车/指定起点/输入W（绘制带有宽度的线）/回车/指定线起点宽度/回车/指定线结束点宽度/回车/输入数值（线的长度值）/回车（结束命令）正多边形：POL 用法：输入命令POL/回车/指定边数/回车/鼠标指定正多边形的中心点/输入选项（C外切于圆；I内接于圆）/回车/输入半径/回车（结束命令）矩形：REC 用法1：输入命令REC/回车/鼠标指定第一角点/指定第二角点用法2：输入命令REC/回车/输入C（绘制带有倒角的矩形）/回车/输入第一倒角值/回车/输入第二倒角值/回车/鼠标指定第一角点/指定第二角点用法3：输入命令REC/回车/输入F（绘制带有圆角的矩形）/回车/输入圆角半径/回车/指定第一角点/指定第二角点圆弧：A 用法：输入命令A/回车/指定圆弧起点/指定圆弧中点/指定圆弧结束点（绘制圆弧的方法有11种，可参考绘图菜单---圆弧选项）圆：C 用法：输入命令C/回车/鼠标指定圆心/输入半径值/回车（命令结束）（绘制圆的方法有6种，可参考绘图菜单---圆选项）样条曲线：SPL 用法：输入命令SPL/回车/鼠标指定要绘制的范围即可/需要三下回车结束命令椭圆：EL

opensees总结

1、定义梁柱单元局部坐标轴的命令流为： geomTransf Linear $transfTag $vecxzX $vecxzY $vecxzZ 其中，$transfTag 代表局部坐标轴矢量的编号，$vecxzX $vecxzY $vecxzZ 表示局部坐标轴的方向矢量值。 2、OPENSEES 的刚性隔板假定命令流格式为： rigidDiaphragm $perpDirn $masterNodeTag $slaveNodeTag1 $slaveNodeTag2 ... 其中，$perpDirn 表示刚性隔板的方法，如实例中楼板的刚性隔板的平移方向为U1（X 方向）与U2（Y 方向），即1-2 平面，该值应为3。$masterNodeTag 为主结点，$slaveNodeTag1 为从结点。主结点一般为刚性隔板刚心。实例中：rigidDiaphragm 3 35 2，表示刚性隔板平动方向为1-2 平面，刚心主节点为35 点，2号结点为从结点。 3、弹性梁柱单元的命令流： element elasticBeamColumn $eleTag $iNode $jNode $A $E $G $J $Iy $Iz $transfTag 需要提供截面的截面积A、截面Y 轴惯性矩Iy，截面Z 轴惯性矩Iz，截面扭转矩，截面材料的弹性模量E 及剪切模量G。其中：$transfTag 与$eleTag 是一致的，表示一个单元有自已特定的坐标轴向量，为了编程的方便。陈：例题三 4、非线性材料模型的定义 (1)uniaxialMaterial Steel01 1 335 200000 0.00001 表示，钢筋的屈服强度为335MPa，弹性模量为200000MPa，硬化系数为0.00001，即屈服平台基本上为水平段。将混凝土材料本构C40 改为非线性混凝土本构【Concrete01】，命令流如下： (2)uniaxialMaterial Concrete01 2 -26.8 -0.002 -10 -0.0033材料参数意见参考图所示。注意：混凝土本构Concrete01 是最简单的混凝土本构，注意数值是负数，即表示受压段。该本构没有受拉段，即受拉强度为0，表示结构一分析即进入弹塑性。 5、采用纤维单元，需要定义纤维截面，纤维截面的定义如下面代码所示： section Fiber 1 { fiber -1.125E+002 -2.700E+002 4.500E+003 2 ……… fiber 1.150E+002 -2.650E+002 4.900E+002 1 } 以上命令流表示，纤维截面编号为1，{}内部为子命令流，表示每一个纤维的信息，每一个纤维的定义格式如下： fiber $Y $Z $Area $Mat 命令中，$Y 表示每个纤维的截面Y 坐标（截面中心为原点0）；$Z 表示每个纤维的截面Z 坐标（截面中心为原点0）；$Area 表示每个纤维的贡献面积；$ Mat 表示每个纤维使用的非线性材料本构的编号。注意：纤维的坐标与材料切线模量可以组装成截面的刚度，而纤维的坐标与材料的应力可以组装成截面的内力（抗力），那么每个纤维的应变可以通过截面的变形与坐标求出。采用纤维截面的单元，即为基于平截面假定。截面变形求解应变是基于平截面假定的。 6、采用的单元为非线性梁柱单元，即基于柔度法的纤维单元（Nonlinear BeamColumn Element or Force Beam Column Element），需要输入命令流如下： element nonlinearBeamColumn $eleTag $iNode $jNode $numIntgrPts $secTag $transfTag 其中，$eleTag 为单元编号；$iNode 为开始结点；$jNode 为结束结点；$numIntgrPts 为积分点数量；$secTag 为纤维截面编号，$transfTag 为局部坐标轴编号。积分点数量，也就纤维单元的计算截面数量，纤维单元的刚度与抗力是由截面刚度与抗力沿杆件长度积分所得，显然，不能将全部截面积分，只能采用

ANSYS的基本使用

2ANSYS的基本使用；2.1ANSYS环境简介；ANSYS有两种模式：一种是交互模式（Inter；运行该程序一般采用Interactive进入，这；进入系统后会有6个窗口，提供使用者与软件之间的交；各窗口的功能如下：；1.应用命令菜单（UtilityMenu）：包含；设定（WorkPlane）、参数化设计（Para；及辅助说明（Help）等；2.主菜单（M 2 ANSYS 的基本使用 2.1 ANSYS环境简介 ANSYS有两种模式：一种是交互模式（Interactive Mode），另一个是非交互模式（Batch Mode）。交互模式是初学者和大多数使用者所采用，包括建模、保存文件、打印图形及结果分析等，一般无特别原因皆用交互模式。但若分析的问题要很长时间，如一、两天等，可把分析问题的命令做成文件，利用它的非交互模式进行分析。运行该程序一般采用 Interactive 进入，这样可以定义工作名称，并且存放到指定的工作目录中。若使用 Run Interactive Now 进入还需使用命令定义工作文件名或使用默认的文件名，使用该方式进入一般是为恢复上一次中断的分析。所以在开始分析一个问题时，建议使用 Interactive 进入交互模式。进入系统后会有6个窗口，提供使用者与软件之间的交流，凭借这6个窗口可以非常容易的输入命令、检查模型的的建立、观察分析结果及图形输出与打印。整个窗口系统称为GUI(Graphical User Interface).如图2-1所示。各窗口的功能如下： 1. 应用命令菜单（Utility Menu）：包含各种应用命令，如文件控制（Fi le）、对象选择（Select）、资料列式（List）、图形显示（Pplot）、图形控制（PlotCtrls）、工作界面

CISCO 常用命令解释

视图模式介绍：普通视图 router> 特权视图 router# /在普通模式下输入enable 全局视图 router(config)# /在特权模式下输入config t 接口视图 router(config-if)# /在全局模式下输入int 接口名称例如int s0或int e0 路由协议视图 router（config-route）# /在全局模式下输入router 动态路由协议名称 1、基本配置： router>enable /进入特权模式 router#conf t /进入全局配置模式 router(config)# hostname xxx /设置设备名称就好像给我们的计算机起个名字 router(config)#enable password /设置特权口令 router(config)#no ip domain lookup /不允许路由器缺省使用DNS解析命令 router(config)# Service password-encrypt /对所有在路由器上输入的口令进行暗文加密router(config)#line vty 0 4 /进入设置telnet服务模式 router(config-line)#password xxx /设置telnet的密码 router(config-line)#login /使能可以登陆 router(config)#line con 0 /进入控制口的服务模式 router(config-line)#password xxx /要设置console的密码 router(config-line)#login /使能可以登陆 2、接口配置： router(config)#int s0 /进入接口配置模式 serial 0 端口配置（如果是模块化的路由器前面加上槽位编号，例如serial0/0 代表这个路由器的0槽位上的第一个接口） router(config-if)#ip add xxx.xxx.xxx.xxx xxx.xxx.xxx.xxx /添加ip 地址和掩码router(config-if)#enca hdlc/ppp 捆绑链路协议 hdlc 或者 ppp 思科缺省串口封装的链路层协议是HDLC所以在show run配置的时候接口上的配置没有，如果要封装为别的链路层协议例如PPP/FR/X25就是看到接口下的enca ppp或者enca fr router(config)#int loopback /建立环回口(逻辑接口)模拟不同的本机网段 router(config-if)#ip add xxx.xxx.xxx.xxx xxx.xxx.xxx.xxx /添加ip 地址和掩码给环回口在物理接口上配置了ip地址后用no shut启用这个物理接口反之可以用shutdown管理性的关闭接口 3、路由配置： (1)静态路由 router(config)#ip route xxx.xxx.xxx.xxx xxx.xxx.xxx.xxx 下一条或自己的接口router(config)#ip route 0.0.0.0 0.0.0.0 s 0 添加缺省路由 (2)动态路由 rip协议 router(config)#router rip /启动rip协议 router(config-router)#network xxx.xxx.xxx.xxx /宣告自己的网段 router(config-router)#version 2 转换为rip 2版本 router(config-router)#no auto-summary /关闭自动汇总功能，rip V2才有作用 router(config-router)# passive-int 接口名 /启动本路由器的那个接口为被动接口

opensees总结

o p e n s e e s总结 -标准化文件发布号：（9556-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

ANSYS常用命令解释

1，/PREP7！加载前处理模块 2，/CLEAR,NOSTART！清除已有的数据,不读入启动文件的设置（不加载初始化文件）初始化文件是用于记录用户和系统选项设置的文本文件/CLEAR, START！清除系统中的所有数据,读入启动文件的设置/FILENAME, EX10.5！定义工程文件名称/TITLE, EX10.5 SOLID MODEL OF AN AXIAL BEARING！指定标题4，F,2,FY,-1000！在2号节点上施加沿着-Y方向大小为1000N的集中力6，FINISH ！退出模块命令 7，/POST1！加载后处理模块 8，PLDISP,2！显示结构变形图,参数“2”表示用虚线绘制出原来结构的轮廓9，ETABLE,STRS,LS,1！用轴向应力SAXL的编号”LS,1”定义单元表STRS ETABLE, MFORX,SMISC,1！以杆单元的轴力为内容,建立单元表MFORX ETABLE, SAXL, LS, 1！以杆单元的轴向应力为内容,建立单元表SAXL ETABLE, EPELAXL, LEPEL, 1！以杆单元的轴向应变为内容,建立单元表EPELAXL ETABLE,STRS_ST,LS,1！以杆件的轴向应力“LS,1”为内容定义单元表STRS_ST ETABLE, STRS_CO, LS,1！以杆件的轴向应力“LS,1”定义单元表STRS_CO ETABLE,STRSX,S,X ！定义X方向的应力为单元表STRSX ETABLE,STRSY,S,Y ！定义Y方向的应力为单元表STRSY *GET,STRSS_ST,ELEM,STEEL_E, ETAB, STRS_ST！从单元表STRS_ST中提取STEEL_E单元的应力结果，存入变量STRSS_ST； *GET, STRSS_CO,ELEM,COPPER_E,ETAB,STRS_CO”从单元表STRS_CO中提取COPPER_E单元的应力结果，存入变量STRSS_CO 10 FINISH ！退出以前的模块 11, /CLEAR, START！清除系统中的所有数据,读入启动文件的设置 12/UNITS, SI！申明采用国际单位制 14/NUMBER, 2！只显示编号,不使用彩色

ANSYS-结构稳态(静力)分析之经典实例-命令流格式

ANSYS 结构稳态(静力)分析之经典实例-命令流格式.txt两人之间的感情就像织毛衣，建立的时候一针一线，小心而漫长，拆除的时候只要轻轻一拉。。。。/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT