module MConcrete !混凝土模板

发布日期:[08-12-04 15:24:26] 浏览人次:[]

d0*SINA*COSA,
1 COSA**2-SINA**2/);
return
end subroutine Con_Get_N
subroutine Con_Get_De(C) !得到弹性本构矩阵
type(typ_Concrete) :: C
C%G0=C%E0/(2.d0*(1.d0+C%MU0))
C%De(1,:)=(/1.d0,C%MU0,0.d0/)
C%De(2,:)=(/C%MU0,1.d0,0.d0/)
C%De(3,:)=(/0.d0,0.d0,(1.d0-C%MU0)/2.d0/)
C%De=C%De*C%E0/(1.d0-2.d0*C%MU0**2)
return
end subroutine Con_Get_De
SUBROUTINE MAXMIN (STRESS,P,AG) !得到主应力(应变方向)
implicit real*8 (A-H,O-Z)
real*8 STRESS(3),P(3) !标量, 主方向
real*8 T(3,3) !转换矩阵
PI=atan(1.0d0)*4.0d0 !得到PI
CC = (STRESS(1)+STRESS(2)) * 0.5
BB = (STRESS(1)-STRESS(2)) * 0.5
CR = SQRT(BB**2 + STRESS(3)**2)
AG=PI/4.d0
IF(BB.NE.0.0d0) Then
AG = 0.5d0* ATAN2(-STRESS(3),BB)
end if
SINA=SIN(AG); COSA=COS(AG)
T(1,:)=(/COSA**2,SINA**2,SINA*COSA/);
T(2,:)=(/SINA**2,COSA**2,-SINA*COSA/);
T(3,:)=(/-2d0*COSA*SINA,2.0d0*SINA*COSA,COSA**2-SINA**2/);
P=matmul(transpose(T),STRESS)
if(P(1) CR=P1; P1=P2; P2=CR;
AG=PI/2+AG;
end if
if(P(1)==0.0.and.P(2)==0.0) then
AG=0;
end if
RETURN
end subroutine MAXMIN
function matinv(A) result (B)
real(8) ,intent (in)::A(:,:)
!real(8) , allocatable::B(:,:)
real(8) , pointer::B(:,:)
integer(4):: N,I,J,K
real(8)::D,T
real(8), allocatable::IS(:),JS(:)
N=size(A,dim=2)
allocate(B(N,N))
allocate(IS(N));allocate(JS(N))
B=A
do K=1,N
D=0.0D0
do I=K,N
do J=K,N
if(abs(B(I,J))>D) then
D=abs(B(I,J))
IS(K)=I
JS(K)=J
end if
end do
end do
do J=1,N
T=B(K,J)
B(K,J)=B(int(IS(K)),J)
B(int(IS(K)),J)=T
end do
do I=1,N
T=B(I,K)
B(I,K)=B(I,int(JS(K)))
B(I,JS(K))=T
end do
B(K,K)=1/B(K,K)
do J=1,N
if(J.NE.K) then
B(K,J)=B(K,J)*B(K,K)
end if
end do
do I=1,N
if(I.NE.K) then
do J=1,N
if(J.NE.K) then
B(I,J)=B(I,J)-B(I,K)*B(K,J)
end if
end do
end if
end do
do I=1,N
if(I.NE.K) then
B(I,K)=-B(I,K)*B(K,K)
end if
end do
end do
do K=N,1,-1
do J=1,N
T=B(K,J)
B(K,J)=B(int(JS(K)),J)
B(int(JS(K)),J)=T
end do
do I=1,N
T=B(I,K)
B(I,K)=B(I,int(IS(K)))
B(I,int(IS(K)))=T
end do
end do
return
end function matinv
end module Mconcrete

|<< << < 1 2 3 > >> >>|
www.mapeng.net 马棚网
www.mapeng.net
文章作者:未知 | 文章来源:网络 | 责任编辑:admin | 发送至邮箱: | 加入收藏:
本文关键字:module  MConcrete  混凝土  模板
本文所属专题:ANSYS技术 
相关资讯
热点资讯
推荐资讯

关于我们 | 站点导航 | 使用帮助 | 友情链接 | 广告服务 | 免责声明 | 新手上路
设为首页 | 加入收藏 | 在线留言 | 马棚网QQ群:{92562572}{102901272}{333259257} | 交流QQ: 客户服务 客户服务 客户服务