耦合场分析中如何加初始应力？

TNC

<P>源程序如下：</P>
<P>/prep7, Silicon beam deflection from an applied voltage</P>
<P>et,1,121                    ! PLANE121 element for air region<BR>et,2,121                    ! Temporary element for beam region<BR>emunit,epzro,8.854e-6       ! Free-space permittivity, μMKSV units<BR>mp,perx,2,1                 ! Relative permittivity for air</P>
<P>bl=150                      ! Beam length (um)<BR>bh=2.0                      ! Beam height <BR>glc=bl/2                    ! Center location of ground<BR>gl=90                       ! Ground conductor length<BR>gh=1.5                      ! Ground conductor thickness<BR><BR>gap=4.5                     ! Air gap<BR>vltg=120                    ! Applied voltage</P>
<P>rectng,0,bl,gap,bh+gap      ! Create model<BR>rectng,glc-gl/2,glc+gl/2,-gh,0<BR>rectng,-10,170,-20,30<BR>aovlap,all</P>
<P>asel,s,area,,1              ! Area for beam<BR>aatt,1,,1<BR>           <BR>asel,s,area,,4              ! Area for air elements<BR>cm,air,area                 ! Group air area into component<BR>aatt,2,,2</P>
<P>allsel,all<BR>smrtsiz,2<BR>amesh,1                     ! Mesh beam<BR>mshape,1<BR>amesh,4                     ! Mesh air with triangles</P>
<P>asel,s,area,,1<BR>lsla,s<BR>dl,all,,volt,vltg           ! Apply voltage to beam</P>
<P>asel,s,area,,2<BR>lsla,s<BR>dl,all,,volt,0              ! Ground conductor (not meshed)</P>
<P>allsel,all<BR>et,1,0                      ! Set structure to null element type<BR>physics,write,ELECTROS      ! Write electrostatic physics file<BR>physics,clear               ! Clear Physics</P>
<P>et,1,82,,,2                 ! Define beam elementy type<BR>et,2,0                      ! Set air to null element type</P>
<P>mp,ex,1,170e3               ! Set Modulus  uN/(um)**2<BR>mp,nuxy,1,0.34<BR>mp,alpx,1,2.6e-6</P>
<P>dl,4,,ux,0                  ! Apply beam constraints<BR>dl,4,,uy,0<BR></P>
<P>allsel,all<BR>finish</P>
<P>physics,write,STRUCTURE     ! Write structural physics file</P>
<P>ESSOLV,'ELECTROS','STRUCTURE',2,0,'air',,,,10  ! Solve coupled-field problem</P>
<P>finish</P>
<P>physics,read,ELECTROS       ! Read electrostatic physics file<BR>/post1<BR>set,last                    ! Retrieve results<BR>esel,s,mat,,2               ! Select air elements<BR>etable,fx,fmag,x            ! Retrieve electrostatic forces<BR>etable,fy,fmag,y<BR>ssum                        ! Sum forces<BR>finish</P>

TNC

<P>如果要研究初始应力对结构响应的影响，那么初始应力该怎么加？</P>
<P>尝试过在</P>
<P>physics,write,STRUCTURE     ! Write structural physics file</P>
<P>的前面加了如下四句</P>
<P>nsel,s,loc,y,gap,bh+gap<BR>nsel,r,loc,x,0,bl<BR>istress,100,0,0,0,0,0<BR>allsel,all</P>
<P>加了以后出现问题，程序运行到istress命令就提示错误，并且强制退出了。</P>

seagar

因为这是平面分析，并没有指定厚度啊。