15 February 2023 1 9K Report

The shear modulus of intact rock is obtained from G=E/2(1+v).

I want to numerically(FLAC3D or 3DEC) calculate the shear modulus of intact rock as the ratio of shear stress to shear strain. Has anyone done this before who can help me?

I have written the FISH related to shear strain and I have also written the FISH related to shear stress as below, but the output is not correct.

def yzstress

yzstr=0

count =0

bh=block_head ;index to list of blocks

loop while bh # 0

bz=b_zone(bh) ;index to block zone list

loop while bz # 0

zy1=z_y(bz) ;y-coordinate of zone centroid

if zy1>0.1 then

yzstr=yzstr+z_syz(bz) ;yz-stress

count=count+1

end_if

bz=z_next(bz) ;index to next zone in block list

end_loop

bh=b_next(bh) ;index of next block in main list

end_loop

yzstress= yzstr/count

end

This code is related to 3dec.

The FISH related to shear strain is correct, but the FISH related to shear stress may be wrong because the amount of stress applied at the boundary to the model is not applied to the output of the FISH. For example, a stress of 2 MPa is applied to the model at the boundary, but the FISH shows no more than 1.3 MPa shear stress.

More Aref Jaberi's questions See All
Similar questions and discussions